1、Matlab畫函數(shù)圖像代碼算法 2010-06-09 13:32:54 閱讀234 評(píng)論2 字號(hào)：大中小 前幾天做家教,教Matlab畫函數(shù)圖像,原想學(xué)了四年應(yīng)該沒什么問題,可真去教起來,覺得四年都白學(xué)了,畫幾張圖費(fèi)了我兩天時(shí)間,雖說最后做出來了,可是依然一知半解的. 下面將這幾天所做的圖像及程序小記一下（望大蝦指教）一、螺旋線1.靜態(tài)螺旋線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,'er

2、asemode','none','markersize',22);xlabel('x軸');ylabel('y軸');zlabel('z軸');title('靜態(tài)螺旋線'); 2.動(dòng)態(tài)螺旋線t=0:0.1:10*pi;1 / 18i=1;h=plot3(sin(t(i),cos(t(i),t(i),'*','erasemode','none');grid onaxis(-2 2 -2 2 0 35)for i=2:length(t)&

3、#160;   set(h,'xdata',sin(t(i),'ydata',cos(t(i),'zdata',t(i);    drawnow    pause(0.01)endtitle('動(dòng)態(tài)螺旋線');(圖略) 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 onaxis('

4、square')xlabel('x軸');ylabel('y軸');zlabel('z軸');title('圓柱螺旋線') 二、旋轉(zhuǎn)拋物面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('旋轉(zhuǎn)拋物面')或直接用：ezsurfc('(

5、X.2+Y.2)./4')    三、橢圓柱面load clownezsurf('(2*cos(u)','4*sin(u)','v',0,2*pi,0,2*pi)view(-105,40)              %視角處理shading interp           

6、; %燈光處理colormap(map)          %顏色處理grid on                        %添加網(wǎng)格線axis equal         

7、;         %使x,y軸比例一致xlabel('x軸');ylabel('y軸');zlabel('z軸');                   %添加坐標(biāo)軸說明title('橢圓柱面')      

8、60;   %添加標(biāo)題四、橢圓拋物面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

9、(v)','2.*sec(u)',-pi./2,3*pi./2,0,2*pi)axis equalgrid onaxis squarexlabel('x軸');ylabel('y軸');zlabel('z軸');title('雙葉雙曲面')六、雙曲柱面load clownezsurf('2*sec(u)','2*tan(u)','v',-pi/2,pi/2,-3*pi,3*pi)hold on     

10、0;        %在原來的圖上繼續(xù)作圖ezsurf('2*sec(u)','2*tan(u)','v',pi/2,3*pi/2,-3*pi,3*pi)colormap(map)shading interpview(-15,30)axis equalgrid onaxis equalxlabel('x軸');ylabel('y軸');zlabel('z軸');title('雙曲柱面')七、雙曲拋物面（馬鞍面）X,

11、Y=meshgrid(-7:0.1:7);Z=X.2./8-Y.2./6;meshc(X,Y,Z);view(85,20)axis('square')xlabel('x軸');ylabel('y軸');zlabel('z軸');title('雙曲拋物面')或直接用：ezsurfc('X.2./8-Y.2./6') 八、拋物柱面X,Y=meshgrid(-7:0.1:7);Z=Y.2./8;h=mesh(Z);rotate(h,1 0 1,180)   

12、0;            %旋轉(zhuǎn)處理%axis(-8,8,-8,8,-2,6);axis('square')xlabel('x軸');ylabel('y軸');zlabel('z軸');title('拋物柱面')或直接用：ezsurfc('Y.2./8') 九、環(huán)面ezmesh('(5+2*cos(u)*cos(v)','(5+2*cos(u)*sin(v

13、)','2*sin(u)',0,2*pi,0,2*pi)axis equalgrid onxlabel('x軸');ylabel('y軸');zlabel('z軸');title('環(huán)面')十、橢球ezsurfc('(5*cos(u)*sin(v)','(3*sin(u)*sin(v)','4*cos(v)',0,2*pi,0,2*pi)axis equalgrid onxlabel('x軸');ylabel('y軸');zla

14、bel('z軸');title('橢球')十一、單葉雙曲面ezsurf('4*sec(u)*cos(v)','2.*sec(u)*sin(v)','3.*tan(u)',-pi./2,pi./2,0,2*pi)axis equalgrid onxlabel('x軸');ylabel('y軸');zlabel('z軸');title('單葉雙曲面')十二、旋轉(zhuǎn)單葉雙曲面load clownezsurf('8*sec(u)*cos(v)',

15、'8.*sec(u)*sin(v)','2.*tan(u)',-pi./2,pi./2,0,2*pi)colormap(map)view(-175,30)%alpha(.2)                   %透明處理axis equalgrid onaxis squarexlabel('x軸');ylabel('y軸');zlabel('z

16、軸');title('旋轉(zhuǎn)單葉雙曲面')十三、圓柱面subplot(1,2,1)ezsurf('(2*cos(u)','2*sin(u)','v',0,2*pi,0,2*pi)grid onshading interpaxis equalxlabel('x軸');ylabel('y軸');zlabel('z軸');title('圓柱面')subplot(1,2,2)cylinder(30)shading interpaxis squaretitle('

17、調(diào)用cylinder函數(shù)所得圓柱面')下面給出用colormap()改變圖像顏色的例子:(用了燈光效果"shading interp") colormap();     %hot/cool/copper/gray/hsv/spring/summer/winter.                     

18、60;              colormap(hsv)                                  &

19、#160;            colormap(hot)                                     colormap(gray)           &#