怎么样用MATLAB绘制三维曲线图

ygkonline

我的毕业设计有数据处理！在一个圆里分割出很多方块，每一个方块都有一个值，现在需要在每一个对应的方块把值显示出来！！
数据如下：
                 

X轴坐标	-20	-15	-10	-5	0	5	10	15	20
Y轴坐标
20			0.335	0.357	0.364	0.356	0.335
15		0.357	0.403	0.429	0.438	0.430	0.403	0.358
10	0.350	0.435	0.491	0.520	0.529	0.520	0.490	0.436	0.350
5	0.435	0.545	0.604	0.630	0.637	0.631	0.604	0.544	0.434
0	0.580	0.709	0.750	0.757	0.757	0.757	0.751	0.708	0.581
-5	0.878	0.958	0.925	0.889	0.875	0.888	0.924	0.957	0.877
-10	1.562	1.284	1.092	0.995	0.966	0.994	1.091	1.283	1.563
-15		1.471	1.166	1.037	1.000	1.037	1.166	1.472
-20			1.091	0.994	0.966	0.995	1.092

请高手将其绘制成三维曲面，谢谢！附上代码！！

老师给的参考代码看不懂（我从来没有学过MATLAB）

[ 本帖最后由 eight 于 2007-6-12 16:11 编辑 ]

wy558558558

不知道你要显示的是图象还是三维曲面？如果要显示图象就用image这个命令．如果其他可以用mesh,surf ,plot3.如果是云图可以用contour．

[ 本帖最后由 eight 于 2007-6-12 19:06 编辑 ]

ygkonline

要求是要三维曲面，能不将上述的代码和截图发过来吗？谢谢老兄！急！

上述的数据是对应于一个以0.0265m为半径的圆，将其网格化为0.005m（也就是上述数据x，y轴坐标的5个单位）其是以方块中心为圆心的，如下图：

[ 本帖最后由 ChaChing 于 2009-9-6 14:38 编辑 ]

w89986581

不明白你的意思,圆形曲面?先解解急,参考一下.
x = -20:5:20;
y = 20:-5:-20;
[x,y] = meshgrid(x,y);
z = [NaN NaN 0.335 0.357 0.364 0.356 0.335 NaN NaN;
    NaN 0.357 0.403  0.429  0.438  0.430  0.403  0.358  NaN;
    0.350  0.435  0.491  0.520  0.529  0.520  0.490  0.436  0.350  ;
      0.435  0.545  0.604  0.630  0.637  0.631  0.604  0.544  0.434  ;
      0.580  0.709  0.750  0.757  0.757  0.757  0.751  0.708  0.581  ;
      0.878  0.958  0.925  0.889  0.875  0.888  0.924  0.957  0.877  ;
      1.562  1.284  1.092  0.995  0.966  0.994  1.091  1.283  1.563  ;
      NaN 1.471  1.166  1.037  1.000  1.037  1.166  1.472 NaN;
      NaN NaN 1.091  0.994  0.966  0.995  1.092 NaN NaN];

surf(x*0.001,y*0.001,z); shading interp

ygkonline

兄弟，我试过了，首先谢谢你哦
还好
就是它不是一个圆曲面
我给的网格化视图中有，靠近圆周的方块没有构成一个完整的方块，故我们没有去取值，但是要求画图时按照一个圆来画，不知道用什么方法比较好，你现在给出的是不靠近圆周的数值图，圆周没有取值的也要算上，最终结果就可以显示一个圆曲面，你再看看行吗？
老师给了段参考程序
我看不懂
你帮忙参考一下
或许对你的理解有好处：


[X1,Y1] = meshgrid(-0.0265:.0001:0.0265);
x=data(:,1); y=data(:,2); z=data(:,3);
dx=reshape(x,1,210); dy=reshape(y,1,210); dz=reshape(z,1,210);
[dxx]=meshgrid(1:210); [dyy]=meshgrid(1:210); [dzz]=meshgrid(1:210);
for i=1:210;dxx(i,:)=dx(1:210);end;
for i=1:210;dyy(i,:)=dy(1:210);end;
for i=1:210;dzz(i,:)=dz(1:210);end;
Z1=griddata(dxx,dyy,dzz,X1,Y1)
meshc(X1,Y1,Z1)

[ 本帖最后由 ChaChing 于 2009-11-13 22:48 编辑 ]

w89986581

clear all
close all

z1 = interp1([-10:5:10],[0.335 0.357 0.364 0.356 0.335],[-30:5:30],'cubic','extrap');
z2 = interp1([-15:5:15],[0.357 0.403  0.429  0.438  0.430  0.403  0.358],[-30:5:30],'cubic','extrap');
z3 = interp1([-20:5:20],[0.350  0.435  0.491  0.520  0.529  0.520  0.490  0.436  0.350],[-30:5:30],'cubic','extrap');
z4 = interp1([-20:5:20],[0.435  0.545  0.604  0.630  0.637  0.631  0.604  0.544  0.434],[-30:5:30],'cubic','extrap');
z5 = interp1([-20:5:20],[0.580  0.709  0.750  0.757  0.757  0.757  0.751  0.708  0.581],[-30:5:30],'cubic','extrap');
z6 = interp1([-20:5:20],[0.878  0.958  0.925  0.889  0.875  0.888  0.924  0.957  0.877],[-30:5:30],'cubic','extrap');
z7 = interp1([-20:5:20],[1.562  1.284  1.092  0.995  0.966  0.994  1.091  1.283  1.563],[-30:5:30],'cubic','extrap');
z8 = interp1([-15:5:15],[1.471  1.166  1.037  1.000  1.037  1.166  1.472],[-30:5:30],'cubic','extrap');
z9 = interp1([-10:5:10],[1.091  0.994  0.966  0.995  1.092],[-30:5:30],'cubic','extrap');
zz = [z1;z2;z3;z4;z5;z6;z7;z8;z9];
for ii = 1:13,
    zz1(ii,:) = interp1([20:-5:-20],zz(:,ii),[30:-5:-30],'cubic','extrap');
end
zz1 = zz1.';
x = -30:5:30;
y = 30:-5:-30;
[x,y] = meshgrid(x,y);

x1 = -27:0.1:27;
y1 = 27:-0.1:-27;
[x1,y1] = meshgrid(x1,y1);
z11 = interp2(x,y,zz1,x1,y1,'cubic');
index = find(x1.^2+y1.^2 > 26.5^2);
z11(index) = NaN;
surf(x1*0.001,y1*0.001,z11); shading interp;

=============================
由于LZ表述的相当含糊不清，甚至错误的情况下，
能猜测并想到该方法,赞一个．
By xjzuo
=============================

[ 本帖最后由 xjzuo 于 2007-6-14 15:24 编辑 ]

ygkonline

楼上的兄弟你太强了！！谢谢你啊！请斑竹给予奖励！！:@)

[ 本帖最后由 xinyuxf 于 2007-6-16 08:58 编辑 ]

w89986581

客气客气,对于外延插值讲究很多的.我没有找到对应的二维外延插值命令,只有一维一维的插啦.
希望以后不要答辩前抓瞎,呵呵.

pengweicai

不错，长见识了。

特别是 z11(index) = NaN;的应用，我以前都是，在别的软件里面把多余的那块去掉的，很麻烦。。。。

w89986581

嘿嘿,我也是前两天在这个版面知道的.

ygkonline

我以前很少接触MATLAB，就这次毕业设计，本来要求用EXCEL做数据处理就可以了，谁知道还要用MATLAB，现在算是长见识，学到知识了！！呵呵
谢谢斑竹，谢谢帮助我的w89986581 ，也谢谢关注振动的各位朋友！谢谢大家！:@)

mpzhou520

好碟子
收获不少

appleseed05

z11(index) = NaN;确实不错，今天新学到一招

wondersoul

原帖由 w89986581 于 2007-6-14 11:41 发表

                               
登录/注册后可看大图

clear all
close all

z1 = interp1([-10:5:10],[0.335 0.357 0.364 0.356 0.335],[-30:5:30],'cubic','extrap');
z2 = interp1([-15:5:15],[0.357 0.403  0.429  0.438  0.430  0.403  0.358],[-30:5:30], ...

我问个比较菜的问题：w89986581的新建BMP图像是程序直接产生的吗？还是生成后手工旋转了下？

qiuxiaoju

非常感谢！
太强了！