Matlab图片功能图像.doc

Matlab画函数图像 HYPERLINK /crazyzcs@126/blog/ \l m=0t=1c=fks_087069080086083075084081083095085084088068080084086070084 代码算法 2010-06-09 13:32:54 阅读234 评论2 字号：大中小 HYPERLINK /blog/static/2052877920086250294907/ /blog/static/2052877920086250294907/ 前几天做家教,教Matlab画函数图像,原想学了四年应该没什么问题,可真去教起来,觉得四年都白学了,画几张图费了我两天时间,虽说最后做出来了,可是依然一知半解的. 下面将这几天所做的图像及程序小记一下（望大虾指教） 一、螺旋线 1.静态螺旋线 a=0:0.1:20*pi;h=plot3(a.*cos(a),a.*sin(a),2.*a,b,linewidth,2);axis([-50,50,-50,50,0,150]);grid onset(h,erasemode,none,markersize,22);xlabel(x轴);ylabel(y轴);zlabel(z轴);title(静态螺旋线); ? 2.动态螺旋线 t=0:0.1:10*pi; i=1; h=plot3(sin(t(i)),cos(t(i)),t(i),*,erasemode,none); grid on axis([-2 2 -2 2 0 35]) for i=2:length(t) ??? set(h,xdata,sin(t(i)),ydata,cos(t(i)),zdata,t(i)); ??? drawnow ??? pause(0.01) end title(动态螺旋线); (图略) ? 3.圆柱螺旋线 t=0:0.1:10*pi; x=r.*cos(t); y=r.*sin(t); z=t; plot3(x,y,z,h,linewidth,2); grid on axis(square) xlabel(x轴);ylabel(y轴);zlabel(z轴); title(圆柱螺旋线) ? 二、旋转抛物面 b=0:0.2:2*pi; [X,Y]=meshgrid(-6:0.1:6); Z=(X.^2+Y.^2)./4; meshc(X,Y,Z); axis(square) xlabel(x轴);ylabel(y轴);zlabel(z轴); title(旋转抛物面) 或直接用：ezsurfc((X.^2+Y.^2)./4)??? 三、椭圆柱面 load clown ezsurf((2*cos(u)),4*sin(u),v,[0,2*pi,0,2*pi]) view(-105,40)????????????? %视角处理 shading interp????????????%灯光处理 colormap(map)????????? %颜色处理 grid on??????????????????????? %添加网格线 axis equal????????????????? %使x,y轴比例一致 xlabel(x轴);ylabel(y轴);zlabel(z轴);?????????????????? %添加坐标轴说明 title(椭圆柱面)????????? %添加标题 四、椭圆抛物面 b=0:0.2:2*pi; [X,Y]=meshgrid(-6:0.1:6); Z=X.^2./9+Y.^2./4; meshc(X,Y,Z); axis(square) xlabel(x轴);ylabel(y轴);zlabel(z轴); title(椭圆抛物面) 或直接用：ezsurfc(X.^2./9+Y.^2./4) 五、双叶双曲面 ezsurf(8*tan(u)*cos(v),8.*tan(u)*sin(v),2.*sec(u),[-pi./2,3*pi./2,0,2*pi]) axis equal grid on axis square xlabel(x轴);ylabel(y轴);zlabel(z轴); title(双叶双曲面) 六、双曲柱面 load clown ezsurf(2*sec(u),2*tan(u),v,[-pi/2,pi/2,-3*pi,3*pi]) hold on????????????? %在原来的图上继续作图 ezsurf(2*sec(u),2*tan(u),v,[pi/2,3*pi/2,-3*pi,3*pi]) colormap(map) shading interp view(-15,30) axis equal grid on axis equal