Giáo trình cơ sở Matlab v5.2-1 - Phần 2 Bài tập ứng dụng I

Tài liệu tham khảo Giáo trình cơ sở Matlab v5.2-1 - Phần 2 Bài tập ứng dụng I
Nội dung trích xuất từ tài liệu:
Giáo trình cơ sở Matlab v5.2-1 - Phần 2 Bài tập ứng dụng IPhÇn bμi tËp vÝ dô vμ lêi gi¶i Bμi tËp øng dông phÇn 1Bμ i 1 X©y dùng hμm bËc nhÊt y = ax + 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¹Bμi gi¶i: %A.1 VÏ theo ph−¬ng tr×nh hμm bËc nhÊt % y = ax + b clg a=0;b=0;c=0;d=0;e=0; disp(Khong gian hai chieu) disp(Ve do thi ham bac nhat y = ax + b); a=input(Vao he so bac nhat ; a = ); b=input(Vao he so tu do : b = ); x=-5:0.1:5; y=a*x+b; 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 bac nhat) hold off clcBμ i 2 X©y dùng hμm bËc hai y = ax^2 + bx + c víi c¸c tham sè a, b, c ®−îc ®−a vμo tõbμn phÝm. Truy xuÊt kÕt qu¶ lªn mμn h×nh ®å ho¹ 120PhÇn 1 - C¬ sëPhÇn bμi tËp vÝ dô vμ lêi gi¶iBμi gi¶i: %B.1 VÏ theo ph−¬ng tr×nh hμm bËc 2 % y = ax^2 + bx + c disp(Next : Ham so bac hai) pause clg disp(Ve do thi ham bac hai y = ax^2 + bx +c); a=input(Vao he so bac hai ; a = ); b=input(Vao he so bac nhat : b = ); c=input(Vao he so tu do c = ); x=-3:0.1:3; y=a*(x.^2)+b*x+c; 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 bac hai) hold off clcBμ i 3 X©y dùng hμm bËc hai y = 1/( ax + b ) víi c¸c tham sè a, b ®−îc ®−a vμo tõ bμnphÝm. Truy xuÊt kÕt qu¶ lªn mμn h×nh ®å ho¹ % VÏ theo ph−¬ng tr×nh hμm Ham so % y = 1/( ax + b ) disp(Next : Ham so y=1/(ax+b)) pause clg disp(Ve do thi ham y =1/(ax + b)); a=input(Vao he so bac nhat ; a = ); b=input(Vao he so tu do : b = ); x=-5:0.1:5; y=1./(a*x+b); 121PhÇn 1 - C¬ sëPhÇ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 ...