|
R数据分析:方法与案例详解中P49练习
求解二元函数g(x1,x2)=(x1^2+x2-11)^2+(x1+x2^2-7)^2
在画图部分,我敲进电脑的程序和书上写的一模一样,如下:
x1<-x2<-seq(-10,10,length=100)
fr2=function(x){
x1=x[1]
x2=x[2]
(x1^2+x2-11)^2+(x1+x2^2-7)^2
}
grr=function(x){
x1=x[1]
x2=x[2]
c(2*(x1^2+x2-11)*2*x1+2*(x1+x2^2-7),2*(x1^2+x2-11)+2*(x1+x2^2-7)*2*x2)
}
z=outer(x1,x2,fr2)
image(x1,x2,z)
contour(x1,x2,z,add=T)
persp(x1,x2,z)
persp(x1,x2,z,box=T,border=T,theta = 45,phi=35)
当程序运行到“z=outer(x1,x2,fr2)”时就会出现问题,显示如下:
Error in FUN(X, Y, ...) : unused argument (Y)
后面无法进行了……
我确实和书上打的一样,然后也使用了rm函数将之前可能存在内存中的变量都去掉了,还是不行。感谢大师援手!
|
|