Giao trinh matlab v5.2 P13

Các câu lệnh có thể được kết hợp với nhau theo trật tự nào đó. Điều này cho phép thiết lập được các chương trình thực hiện được nhiều chức năng. Xa hơn, ngoài các câu lệnh để thay đổi và điều chỉnh dữ liệu, còn có những kiểu câu lệnh dùng để điều khiển dòng xử lý của máy tính như là phân nhánh, định nghĩa bởi nhiều trường hợp, vòng lặp, hay kết hợp các chức năng. Đây là các thành tố không thể thiếu của một ngôn ngữ lập trình....
Nội dung trích xuất từ tài liệu:
Giao trinh matlab v5.2 P13PhÇn bμi tËp vÝ dô vμ lêi gi¶i hold on plot(x,y,m-) plot(y,zeros(x),c-) plot(zeros(x),x,c-) text(-1,-1.5,O) text(-0.05,max(y),^) text(max(x),0,>) title(Ham y=1/(ax+b)) hold off clcBμ i 4 X©y dùng hμm r = a* phi víi c¸c tham sè a ®−îc ®−a vμo tõ bμn phÝm. Truy xuÊtkÕt qu¶ lªn mμn h×nh ®å ho¹ víi hÖ to¹ ®é dïng lμ hÖ to¹ ®é cùcBμi gi¶i: % VÝ dô vÒ hÖ to¹ ®é cùc disp(Next : He toa do cuc) pause clg % D.1 VÏ ®−êng xoan èc % r = a* phi disp(Ve duong xoan oc : r = a*tt) pause clg a=input(Vao he so a = ); tt=0:0.1:8*pi; r=a*tt; axis(equal,off) polar(tt,r) title(Duong xoan oc) disp(Ve nhieu lan) pause axis(equal,off) for m=1:8 hold on r1=r*m; polar(tt,r1) hold off end 122PhÇn 1 - C¬ sëPhÇn bμi tËp vÝ dô vμ lêi gi¶iBμ i 5 X©y dùng hμm r = a*cos ( phi ) + b víi c¸c tham sè a,b ®−îc ®−a vμo tõ bμn phÝm. Truy xuÊt kÕt qu¶ lªn mμn h×nh ®å ho¹ víi hÖ to¹ ®é dïng lμ hÖ to¹ ®é cùcBμi gi¶i: %D.2 §−êng èc sªn r = a*cos ( phi ) + b disp(Next :duong oc sen r=a*cos(tt)+b) pause clg a=input(Vao he so a = ); b=input(Vao he so b = ); tt=0:0.1:8*pi; r=a*cos(tt)+b; axis(equal,off) polar(tt,r) title(Duong oc sen) disp(Ve nhieu lan) pause for m=1:8 hold on r1=r*m; polar(tt,r1) hold off endBμ i 6 X©y dùng hμm Astroit víi c¸c tham sè a ®−îc ®−a vμo tõ bμn phÝm. Truy xuÊt kÕt qu¶ lªn mμn h×nh ®å ho¹ víi hÖ to¹ ®é dïng lμ hÖ to¹ ®é cùcBμi gi¶i: %D.3 §−êng astroit disp(Next :duong Astroit ) pause clg a=input(Vao he so a = ); tt=0:0.1:8*pi; r=a*sqrt(abs(1-sin(3*tt)/4)); polar(tt,r) title(Duong Astroit) disp(Ve nhieu lan) pause 123PhÇn 1 - C¬ sëPhÇn bμi tËp vÝ dô vμ lêi gi¶i for m=1:8 hold on r1=r*m; polar(tt,r1) hold off endBμ i 7 X©y dùng ph−¬ng tr×nh ®−êng Lemniscat Becnulli víi c¸c tham sè a ®−îc ®−a vμotõ bμn phÝm. Truy xuÊt kÕt qu¶ lªn mμn h×nh ®å ho¹ víi hÖ to¹ ®é dïng lμ hÖ to¹ ®écùcBμi gi¶i: % D.4 §−êng Lemniscat Becnulli disp(Next :duong Lemniscat Becnulli) pause clg a=input(Vao he so a = ); tt=0:0.1:8*pi; r=a*sqrt(abs(2*cos(2*tt))); axis(equal,off) polar(tt,r) title(Duong xoan oc) disp(Ve nhieu lan) pause for m=1:8 hold on r1=r*m; polar(tt,r1) hold off endBμ i 8 Dïng hμm bucky ®Ó x©y dùng h×nh gi¶ 3chiÒu. Truy xuÊt kÕt qu¶ lªn mμn h×nh ®åho¹Bμi gi¶i: %Kh«ng gian 3D 124PhÇn 1 - C¬ sëPhÇn bμi tËp vÝ dô vμ lêi gi¶i disp(Khong gian ba chieu ) pause clg %e.1 VÏ h×nh qu¶ bãng disp(Ve qua bong da) [B,V]=bucky; H=sparse(60,60); k=31:60; H(k,k)=B(k,k); x=V(:,1); y=V(:,2); gplot(H,V,m-) axis(equal,off); hold on gplot(B-H,V,c-) hold offBμ i 9 VÏ hμm ®å thÞ trong kh«ng gian 3 chiÒu. Dïng plot3()Bμi gi¶i %e.2 VÏ ®−êng cã h×nh ¶nh kh«ng gian disp(Ve duong co hinh anh khong gian) pause clg t=0:pi/50:8*pi; plot3(sin(t),cos(t),t);Bμi 10 VÏ mét sè bÒ mÆt vÝ dô trong kh«ng gian 3 chiÒu víi c¸c tham sè tuú chän. MÆtparabolloit, mÆt trô.Bμi gi¶i %e.3 VÏ mÆt kh«ng gian 3D disp(Next: Ve mat khong gian ba chieu) disp(Ve Parabolloit) pause clg 125PhÇn 1 - C¬ sëPhÇn bμi tËp vÝ dô vμ lêi gi¶i t=-5:0.1:5; [x,y]=meshdom(t,t); z=x.^2+y.^2; mesh(z) title(Paraboloit) disp(Next: Mat tru sinh boi y=x^2) pause clg z=sqrt(x.^4+y.^2); mesh(z) title(Mat tru) pauseBμi 11 X©y dùng menu trong m«i tr−êng Matlab vμ thùc hiÖn mét sè c¸c thao t¸c x©ydùng c¸c hμm ®å ho¹ ®¬n gi¶n. Bao gåm: VÏ mét h×nh cÇu, ph−¬ng tr×nh ®−êngsin(x)^2, sin(x^2)*exp(-x), sin(1/x)^2/x vμ tich ph©n x¸c ®Þnh cña hμm bÊt kú.Bμi gi¶ifunction Thuctap(action);%Thuctap Chuong trinh nay ve mot do hoa bao gom chuc % nang ve mot so ham va tich phan%Nhung viec can lam:% Nut1 :Sphere (Hinh cau)% Nut2,3,4 : Phuong trinh cac ham co ban% Nut 5 : Tich phan xac dinhif narginPhÇn bμi tËp vÝ dô vμ lêi gi¶i% Khoang cach giua nut va nhan cua lenh tiep theo kc=0.03;% Khung nen cho cac phim chuc nang :The CONSOLE frame frmBorder=0.01; yPos=0.01; frmPos=[xPos+0.02 yPos-frmBorder btnLen+4*frmBorder0.9+11*frmBorder]; uicontrol( ... Style,frame, ... Units,normalized, ... Position,frmPos, ... Bac ...