Giáo trình cơ sở matlab v5 2 1 phần 2 bài tập ứng dụng i

Phần bi tập ví dụ v lời giải Bi tập ứng dụng phần B i Xây dựng hm bậc y = ax + b với tham số a,b đợc đa vo từ bn phím Truy xuất kết lên mn hình đồ hoạ Bi giải: %A.1 Vẽ theo phơng trình hm bậc % y = ax + b clg a=0;b=0;c=0;d=0;e=0; disp('Khong gian hai chieu') disp('Ve thi ham bac nhat y = ax + b'); a=input('Vao he so bac nhat ; a = '); b=input('Vao he so tu : 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 clc B i Xây dựng hm bậc hai y = ax^2 + bx + c với tham số a, b, c đợc đa vo từ bn phím Truy xuất kết lên mn hình đồ hoạ Phần - Cơ sở 120 Phần bi tập ví dụ v lời giải Bi giải: %B.1 Vẽ theo phơng trình hm bậc % y = ax^2 + bx + c disp('Next : Ham so bac hai') pause clg disp('Ve 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 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 clc B i Xây dựng hm bậc hai y = 1/( ax + b ) với tham số a, b đợc đa vo từ bn phím Truy xuất kết lên mn hình đồ hoạ % Vẽ theo phơng trình hm Ham so % y = 1/( ax + b ) disp('Next : Ham so y=1/(ax+b)') pause clg disp('Ve thi ham y =1/(ax + b)'); a=input('Vao he so bac nhat ; a = '); b=input('Vao he so tu : b = '); x=-5:0.1:5; y=1./(a*x+b); Phần - Cơ sở 121 Phần bi 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 clc B i Xây dựng hm r = a* phi với tham số a đợc đa vo từ bn phím Truy xuất kết lên mn hình đồ hoạ với hệ toạ độ dùng l hệ toạ độ cực Bi giải: % Ví dụ hệ toạ độ cực disp('Next : He toa 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 Phần - Cơ sở 122 Phần bi tập ví dụ v lời giải B i Xây dựng hm r = a*cos ( phi ) + b với tham số a,b đợc đa vo từ bn phím Truy xuất kết lên mn hình đồ hoạ với hệ toạ độ dùng l hệ toạ độ cực Bi 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 end B i Xây dựng hm Astroit với tham số a đợc đa vo từ bn phím Truy xuất kết lên mn hình đồ hoạ với hệ toạ độ dùng l hệ toạ độ cực Bi 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 Phần - Cơ sở 123 Phần bi tập ví dụ v lời giải for m=1:8 hold on r1=r*m; polar(tt,r1) hold off end B i Xây dựng phơng trình đờng Lemniscat Becnulli với tham số a đợc đa vo từ bn phím Truy xuất kết lên mn hình đồ hoạ với hệ toạ độ dùng l hệ toạ độ cực Bi 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 end B i Dùng hm bucky để xây dựng hình giả 3chiều Truy xuất kết lên mn hình đồ hoạ Bi giải: %Không gian 3D Phần - Cơ sở 124 Phần bi tập ví dụ v lời giải disp('Khong gian ba chieu ') pause clg %e.1 Vẽ hình 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 off B i Vẽ hm đồ thị không gian chiều Dùng plot3() Bi 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); Bi 10 Vẽ số bề mặt ví dụ không gian chiều với tham số tuỳ chọn Mặt parabolloit, mặt trụ Bi giải %e.3 Vẽ mặt không gian 3D disp('Next: Ve mat khong gian ba chieu') disp('Ve Parabolloit') pause clg Phần - Cơ sở 125 Phần bi 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') pause Bi 11 Xây dựng menu môi trờng Matlab v thực số thao tác xây dựng hm đồ hoạ đơn giản Bao gồm: Vẽ hình cầu, phơng trình đờng sin(x)^2, sin(x^2)*exp(-x), sin(1/x)^2/x v tich phân xác định hm Bi giải function Thuctap(action); %Thuctap Chuong trinh ve mot 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 : Tich phan xac dinh if nargin1 s1=s1+1; s=input('Neu tiep tuc thi s[...]... W=(g -1) :0 .1: (g-0 .1) t=t +1; Phần 1 - Cơ sở 13 6 Phần bi tập ví dụ v l i gi i for i =2: h if i= =2 for l =1: (m+h -2) if w(l +i -1) ==w(l) if w(l +i) ==w(l +1) M(l ,i) =0; else M(l ,i) =(w(l +i) -W)*M(l +1, i -1) /(w(l +i) -w(l +1) ); end else if w(l +i) ==w(l +1) M(l ,i) =(W-w(l))*M(l ,i -1) /(w(l +i -1) -w(l)); else M(l ,i) =(W-w(l))*M(l ,i -1) /(w(l +i -1) -w(l)) +(w(l +i) W)*M(l +1, i -1) /(w(l +i) -w(l +1) ); end end end end if i= =3 for l =1: (m+h-3) if w(l +i -1) ==w(l)... v=input('vao khoang can ve cua u (1 ,2 n)='); g=input('vao khoang can ve cua w (1 ,2 n)='); for i =1: (n+k -1) if u (i) < u (i +1) if u (i) ==v -1 N (i ,1) =1; else N (i ,1) =0; end else N (i ,1) =0; end end for i =1: (n+k -1) t=N (i ,1) ; end for i =1: (m+h -1) if w (i) < w (i +1) if w (i) ==g -1 M (i ,1) =1; else M (i ,1) =0; end else M (i ,1) =0; end end for i =1: (m+h -1) t1=M (i ,1) ; end X1=[];Y1=[];Z1=[]; for U=(v -1) :0 .1: (v-0 .1) t=0; for W=(g -1) :0 .1: (g-0 .1) ... x(3 ,1) =1; x(3 ,2) =1; x(3,3) =1; x(3,4) =1; x(4 ,1) =3;x(4 ,2) =3;x(4,3)=3;x(4,4)=3; y (1, 1)=0;y (1 ,2) =3;y (1, 3)=3;y (1, 4)=0; y (2, 1) =3;y (2, 2)=5;y (2, 3)=5;y (2, 4)=3; Phần 1 - Cơ sở 13 5 Phần bi tập ví dụ v l i gi i y(3 ,1) =3;y(3 ,2) =5;y(3,3)=5;y(3,4)=3; y(4 ,1) =5;y(4 ,2) =5;y(4,3)=5;y(4,4)=5; z (1, 1)=5;z (1 ,2) =3;z (1, 3)=-3;z (1, 4)=-5; z (2, 1) =5;z (2, 2)=3;z (2, 3)=-3;z (2, 4)=-5; z(3 ,1) =5;z(3 ,2) =3;z(3,3)=-3;z(3,4)=-5; z(4 ,1) =5;z(4 ,2) =3;z(4,3)=-3;z(4,4)=-5;... if w(l +i -1) ==w(l) if w(l +i) ==w(l +1) M(l ,i) =0; else M(l ,i) =(w(l +i) -W)*M(l +1, i -1) /(w(l +i) -w(l +1) ); end else if w(l +i) ==w(l +1) M(l ,i) =(W-w(l))*M(l ,i -1) /(w(l +i -1) -w(l)); else M(l ,i) =(W-w(l))*M(l ,i -1) /(w(l +i -1) -w(l)) +(w(l +i) W)*M(l +1, i -1) /(w(l +i) -w(l +1) ); end end end end end for j =2: k if j= =2 for l =1: (n+k -2) Phần 1 - Cơ sở 13 7 Phần bi tập ví dụ v l i gi i if u(l+j -1) ==u(l) if u(l+j)==u(l +1) N(l,j)=0; else... l i gi i k=input('k='); for i =1: (n+4) if i n u (i) =n-k +1; else u (i) =i- k; end end x=input('Nhap vao n toa do Px='); y=input('Nhap vao n toa do Py='); z=input('Nhap vao n toa do Pz='); m=input('vao khoang can ve (1 ,2 n)='); for i =1: (n+3) if u (i) < u (i +1) if u (i) ==m -1 N (i ,1) =1; else N (i ,1) =0; end else N (i ,1) =0; end end for i =1: (n+3) t=N (i ,1) ; end t=0; for U=(m -1) :0. 12 5 :m %U=0. 12 5 *i; ... q=input('Neu tiep tuc thi danh y nguoc lai la n ='); for i =1: (n+k) if i n u (i) =n-k +1; else u (i) =i- k; end end for i =1: (m+h) if i m w (i) =m-h +1; else w (i) =i- h; end end %x=input('Nhap vao n.m toa do Px='); %y=input('Nhap vao n.m toa do Py='); %z=input('Nhap vao n.m toa do Pz='); x (1, 1)=-3;x (1 ,2) =-3;x (1, 3)=-3;x (1, 4)=-3; x (2, 1) = -1; x (2, 2)= -1; x (2, 3)= -1; x (2, 4)= -1; ... 'FaceColor', 'interp') colormap(cool (10 0)) lighting phong; set(gca, 'CameraViewAngle', 7); set(gcf, 'color', [1 1 1] ); Bi 13 Xây dựng v vẽ hình đờng B-Spline trong không gian 2D v 3D từ các i m kiểm soát đợc vo từ bn phím hay các file dữ liệu Trên cơ sở đờng cong phát triển thnh mặt B-spline Bi gi i s =2; s1=0; while s >1 s1=s1 +1; s=input('Neu tiep tuc thi s ... ve (1 ,2 n)='); for i =1: (n+3) if u (i) < u (i +1) if u (i) ==m -1 N (i ,1) =1; else N (i ,1) =0; end else N (i ,1) =0; end end for i =1: (n+3) t=N (i ,1) ; end t=0; for U=(m -1) :0. 12 5 :m %U=0. 12 5 *i; t=t +1; for j =2: 4 if... w (i +1) if w (i) ==g -1 M (i ,1) =1; else M (i ,1) =0; end else M (i ,1) =0; end end for i =1: (m+h -1) t1=M (i ,1) ; end X1=[];Y1=[];Z1=[]; for U=(v -1) :0 .1: (v-0 .1) t=0; for W=(g -1) :0 .1: (g-0 .1) t=t +1; Phần - Cơ sở. .. sở 13 6 Phần bi tập ví dụ v l i gi i for i =2: h if i= =2 for l =1: (m+h -2) if w(l +i -1) ==w(l) if w(l +i) ==w(l +1) M(l ,i) =0; else M(l ,i) =(w(l +i) -W)*M(l +1, i -1) /(w(l +i) -w(l +1) ); end else if w(l +i) ==w(l +1)