mesh画三维曲线图出现问题

银杏树之舞

X=0:0.01:1;Z=590:1:(590-500*x);
[x,z]=meshgrid(X,Z);
s=-192912544908704749/17918603558912+40*(z-20/1043*(40000*x+40600).*x-1/26075*(60600-60000*x).*(500-500*x)-12690/1043)./(-1227626434375/447-1/2*(-1600/1043*(40000*x+40600).*x-16/5215*(60600-60000*x).*(500-500*x)-1015200/1043).*(348100-(590-500*x).^2)-20000*(-20/1043*(40000*x+40600).*x-1/26075*(60600-60000*x).*(500-500*x)-12690/1043).^2.*x-1/2*(-1600/1043*(40000*x+40600).*x-16/5215*(60600-60000*x).*(500-500*x)-1015200/1043).*((590-500*x).^2-8100)-40*(-20/1043*(40000*x+40600).*x-1/26075*(60600-60000*x).*(500-500*x)-12690/1043).^2.*(500-500*x)+2092500/149*(40000*x+40600).*x+4185/149*(60600-60000*x).*(500-500*x)-6075*(-20/1043*(40000*x+40600).*x-1/26075*(60600-60000*x).*(500-500*x)-12690/1043).^2).*(-24114068113588093625/4479650889728*x.*(755-250*x)+(-24114068113588093625/4479650889728+24114068113588093625/4479650889728*x).*(415-250*x)+62937360000/149);
mesh(x,z,s)
hold on
plot3(x(1,1),z(1,1),s(1,1),'o','markersize',8)
??? Error using ==> set
View [azimuth elevation] must be a two element vector.

Error in ==> mesh at 128
      set(cax,'View',y);

liuyalong008

1. x=0:0.01:1;
   
2. z=590-500.*x;
   
3. [x,z]=meshgrid(x,z);
   
4. s=-192912544908704749/17918603558912+40*(z-20/1043*(40000*x+40600).*x-1/26075*(60600-60000*x).*(500-500*x)-12690/1043)./(-1227626434375/447-1/2*(-1600/1043*(40000*x+40600).*x-16/5215*(60600-60000*x).*(500-500*x)-1015200/1043).*(348100-(590-500*x).^2)-20000*(-20/1043*(40000*x+40600).*x-1/26075*(60600-60000*x).*(500-500*x)-12690/1043).^2.*x-1/2*(-1600/1043*(40000*x+40600).*x-16/5215*(60600-60000*x).*(500-500*x)-1015200/1043).*((590-500*x).^2-8100)-40*(-20/1043*(40000*x+40600).*x-1/26075*(60600-60000*x).*(500-500*x)-12690/1043).^2.*(500-500*x)+2092500/149*(40000*x+40600).*x+4185/149*(60600-60000*x).*(500-500*x)-6075*(-20/1043*(40000*x+40600).*x-1/26075*(60600-60000*x).*(500-500*x)-12690/1043).^2).*(-24114068113588093625/4479650889728*x.*(755-250*x)+(-24114068113588093625/4479650889728+24114068113588093625/4479650889728*x).*(415-250*x)+62937360000/149);
   
5. mesh(x,z,s)
   
6. hold on
   
7. plot3(x(1,1),z(1,1),s(1,1),'o','markersize',8)

复制代码

liuyalong008

Z=590:1:(590-500*x);  是嘛意思？？？

银杏树之舞

应该是Z=(590-500*x):1:590，不过错误是一样的

银杏树之舞

也就是说Z的取值范围随x而变化，即Z=[(590-500*x),590]

银杏树之舞

好像是我的z定义的有问题，我再想想

银杏树之舞

liuyalong008 发表于 2012-2-18 12:38
Z=590:1590-500*x);  是嘛意思？？？

可能是我写的有问题，x=0:0.1:1;590-500*x=<z<=590;,麻烦大虾再给算一下。:)

银杏树之舞

我算了一下，z=590-500*x与z=90：1:590的结果是一样的，如何才能体现出z=590-500*x  ？？

银杏树之舞

这样不行吗？

银杏树之舞

银杏树之舞 发表于 2012-2-20 18:11
可能是我写的有问题，x=0:0.1:1;590-500*x=

如何定义z才能得到的我想要图，或其他方法，如何才能实现呢？

liuyalong008

银杏树之舞 发表于 2012-2-21 00:36
我算了一下，z=590-500*x与z=90：1:590的结果是一样的，如何才能体现出z=590-500*x  ？？ ...

表述不清，你到底是要怎么做？？？

银杏树之舞

直接说吧，主要是想做两个图：
X=0:0.01:1; 90<Z<590-500.*X;
s(1,3)=-82*(-428467000*x+858624940)./(-21000*x+47075)-82*(z-(500*(40000*x+40600).*x+(124230-123000*x).*(500-500*x)+317250)./(-21000*x+47075))./(-8235217625/3-14*(590-500*x).^3+40*(500*(40000*x+40600).*x+(124230-123000*x).*(500-500*x)+317250)./(-21000*x+47075).*(348100-(590-500*x).^2)-20000*(500*(40000*x+40600).*x+(124230-123000*x).*(500-500*x)+317250).^2./(-21000*x+47075).^2.*x+82*(500*(40000*x+40600).*x+(124230-123000*x).*(500-500*x)+317250)./(-21000*x+47075).*((590-500*x).^2-8100)-82*(500*(40000*x+40600).*x+(124230-123000*x).*(500-500*x)+317250).^2./(-21000*x+47075).^2.*(500-500*x)+732375*(500*(40000*x+40600).*x+(124230-123000*x).*(500-500*x)+317250)./(-21000*x+47075)-6075*(500*(40000*x+40600).*x+(124230-123000*x).*(500-500*x)+317250).^2./(-21000*x+47075).^2).*(500*(40*(-428467000*x+858624940)./(-21000*x+47075)-11521576189034495/17179869184).*x.*(755-250*x)+(82*(-428467000*x+858624940)./(-21000*x+47075)-1527578).*(500-500*x).*(415-250*x)+732375/2*(-428467000*x+858624940)./(-21000*x+47075)-5618578050)+ 1548898
surf(x,z,s(1,3))
shading interp

X=0:0.01:1; 590-500.*X<Z<590;
s(1,4)=-40*(-428467000*x+858624940)./(-21000*x+47075)-40*(z-(500*(40000*x+40600).*x+(124230-123000*x).*(500-500*x)+317250)./(-21000*x+47075))./(-8235217625/3-14*(590-500*x).^3+40*(500*(40000*x+40600).*x+(124230-123000*x).*(500-500*x)+317250)./(-21000*x+47075).*(348100-(590-500*x).^2)-20000*(500*(40000*x+40600).*x+(124230-123000*x).*(500-500*x)+317250).^2./(-21000*x+47075).^2.*x+82*(500*(40000*x+40600).*x+(124230-123000*x).*(500-500*x)+317250)./(-21000*x+47075).*((590-500*x).^2-8100)-82*(500*(40000*x+40600).*x+(124230-123000*x).*(500-500*x)+317250).^2./(-21000*x+47075).^2.*(500-500*x)+732375*(500*(40000*x+40600).*x+(124230-123000*x).*(500-500*x)+317250)./(-21000*x+47075)-6075*(500*(40000*x+40600).*x+(124230-123000*x).*(500-500*x)+317250).^2./(-21000*x+47075).^2).*(500*(40*(-428467000*x+858624940)./(-21000*x+47075)-11521576189034495/17179869184).*x.*(755-250*x)+(82*(-428467000*x+858624940)./(-21000*x+47075)-1527578).*(500-500*x).*(415-250*x)+732375/2*(-428467000*x+858624940)./(-21000*x+47075)-5618578050)+680004；
surf(x,z,s(1,4))
shading interp

银杏树之舞

少了一步[x,z]=meshgrid(X,Z);