求助 怎么用MATLAB画出下面函数的曲面图形

sonia001

z = (p1+p2*Ln(x)+p3*(Ln(x))^2+p4*(Ln(x))^3+p5*y+p6*y^2)/(1+p7*Ln(x)+p8*(Ln(x))^2+p9*y)
其中：X=(1:50:3031),y=(30:5:90);
  各参数值如下：
p1=26.248358078106;
p2=-24.1298122868307;
p3=6.91270559265974;
p4=-0.411821386110327;
p5=0.129226276697087;
p6=0.000615572200627997;
p7=-0.442965192301227;
p8=0.0598765826900065;
p9 =0.000802780902706521;

谢谢各位了 ，不胜感激
希望能够给出程序代码，以便能够学习 :@)

ch_j1985

原帖由 sonia001 于 2008-8-2 10:57 发表

                               
登录/注册后可看大图

z = (p1+p2*Ln(x)+p3*(Ln(x))^2+p4*(Ln(x))^3+p5*y+p6*y^2)/(1+p7*Ln(x)+p8*(Ln(x))^2+p9*y)
其中：X=(1:50:3031),y=(30:5:90);
  各参数值如下：
p1=26.248358078106;
p2=-24.1298122868307;
p3=6.9127055926 ...

请问LZ：Ln(x)是不是代表对x取自然对数？
LZ可以到绘图技巧看一下
是否可以参考此帖：http://www.chinavib.com/forum/viewthread.php?tid=6362&extra=page%3D23%26amp%3Bfilter%3Dtype%26amp%3Btypeid%3D23

[ 本帖最后由 ch_j1985 于 2008-8-2 11:34 编辑 ]

sonia001

是求自然对数
谢谢啦
待我瞧瞧内容再讨论  ：）

sonia001

嗯 看过了 自己尝试了一下
但是:@(
比较复杂
最多只能作出三维曲线，但是却得不到曲面图
可否再指点下 谢谢 :@)

ch_j1985

原帖由 sonia001 于 2008-8-2 15:50 发表

                               
登录/注册后可看大图

嗯 看过了 自己尝试了一下
但是:@(
比较复杂
最多只能作出三维曲线，但是却得不到曲面图
可否再指点下 谢谢 :@)

将X=(1:50:3031),y=(30:5:90); 改为X=(1:1:3031),y=(30:1:90); 试试
亦即修改X、Y的步长

guiguihong

我吧x,y的取值范围改变了一下，用如下的代码绘制，得出图形，不晓得是不是你所需要的。我是初学的，
x=(30:1:60);y=(30:1:60); %保证x,y的范围长度相等
[X,Y]=meshgrid(x,y);%产生格点矩阵
p1=26.248358078106;
p2=-24.1298122868307;
p3=6.91270559265974;
p4=-0.411821386110327;
p5=0.129226276697087;
p6=0.000615572200627997;
p7=-0.442965192301227;
p8=0.0598765826900065;
p9 =0.000802780902706521;
Z = (p1+p2*log(X)+p3*(log(X))^2+p4*(log(X))^3+p5*Y+p6*Y^2)/(1+p7*log(X)+p8*(log(X))^2+p9*Y);
surf(X,Y,Z)