Bài tập môn học matlab

Trng H Kinh t k thut cụng nghiờp Phc Lc Chng I : Matlab c bn Chng II : Lý thuyt iu t ng Chng III:V th Phn IV Gii tớch mch in bng Simulink SV:Nguyn Ngc Tnh in1A Trng H Kinh t k thut cụng nghiờp Chng I : Matlab c bn Bai 1: Lp chng trỡnh dang function gii phng trỡnh bc >> %vi du >> bt1 Nhap he so a=3 Nhap he so b=6 Nhap he so c=-9 Cac Nghiem so x1 = x2 = -3 Bi 2: Tớnh tớch phõn ca I bit I= * x^ * cos( x) * dx Gii: Chng trỡnh vit M-file nh sau SV:Nguyn Ngc Tnh in1A Trng H Kinh t k thut cụng nghiờp Chng trỡnh thc hin comman Window: Bi [1 4]H y lập chơng trình cho toán nội suy Lagrange chiều khoảng cách H y tìm kết khoảng đó: S liu; x F 0.67 0.620 0.567 0.512 Gii: Chng trỡnh vit M-file nh sau SV:Nguyn Ngc Tnh in1A Trng H Kinh t k thut cụng nghiờp f= 0.6700 0.6200 0.5670 0.5120 n=4 nhap gia tri cua xA[1 4]=3 diem A can tim thuoc day ngoai suy bang= 1.0000 5.0000 0.6700 2.0000 6.0000 0.6200 3.0000 7.0000 0.5670 4.0000 8.0000 0.5120 Kq= 4.9368 >> Bi 4: Cú n cụng vic phi lm liờn tip ,cụng vic phi lm xong thi gian t1 cụng vic phi hon thnh trc thi gian c1.nu cụng vic I oc hon thnh ti thi im d1 thỡ nú s b pht vi h s lng pht nh sau:pi(di-ci)nu hon thnh ti di sau ci cng nhiu thỡ pht cng nng.Ngcj li nu hon thnh trc (di-ci)>bt5 nhap ma tran so hang la m=2 so cot la n=3 nhap phan tu nhap phan tu nhap phan tu nhap phan tu nhap phan tu nhap phan tu sai = dai=1 SV:Nguyn Ngc Tnh in1A Trng H Kinh t k thut cụng nghiờp rong= -1 sai=0 dai= -1 >> Bài 6: Cho mt qu cu nng cú lng l m=2kg c buc vo mt si dõy di l=1.5m mt ngi nm vo u dõy v quay trũn qu cu theo phng thng ng Vn tc di ca qu cu v trớ cao nht l v=6m/s Lc kộo T trờn si dõy thay i theo gúc nghiờng ca dõy vi phng nm ngang nh theo cụng thc sau: T=(mv2/l)+mg(2-3sin ) Cho bit lc kộo T>145N thỡ si dõy b t Hóy xỏc nh v trớ dõy b t Giải : Kt qu common windown Function y= tinh luc (theta,T) M=2 ;l=1,5;Tmax=145;v=6;g=9,8; Theta=0:10:360; Rad=pi*theta/180; T=m*(v^2/l+g*(2-3*sin(rad))); Bang=[theta;T>Tmax] Plot(theta,T); Grid on >> Tính lực Bang = Columns through 12 10 20 30 40 50 60 0 0 0 70 80 90 100 110 Columns 13 through 24 120 130 140 150 160 170 180 190 200 210 220 230 0 0 0 0 0 0 Columns 25 through 36 240 250 26 270 280 290 300 310 320 330 340 350 0 1 0 0 0 Columns 37 360 T kt qu trờn ta thy ti cỏc v trớ tng ng vi cỏc gúc 260,270.280, thỡ biu thc logic T>T max cú giỏ tr bng tc ti v trớ ú cú th b t Bài 7: SV:Nguyn Ngc Tnh in1A v Trng H Kinh t k thut cụng nghiờp Gii bi mch in di dng script file: Cho mch in nh hỡnh v L1 R L C CHOKE HASH-IRON U Gi i: Chng trỡnh vit M-file nh sau SV:Nguyn Ngc Tnh in1A Trng H Kinh t k thut cụng nghiờp Kt qu commonwindown : SV:Nguyn Ngc Tnh in1A Trng H Kinh t k thut cụng nghiờp Bi Lp chng trỡnh kim tra im s ca mt hc sinh xp loi xut sc hay khỏ gii, khỏ trung bỡnh Gii Chng trỡnh vit M-file nh sau Kt qu Commanwindows: Bi 9: Cho mt dóy s quỏ trỡnh nhp s dng li s a vo l s Hóy xp xp SV:Nguyn Ngc Tnh 10 in1A Trng H Kinh t k thut cụng nghiờp set(gca,'Zscale','log'); view([125,7]); grid on title('ve 3D voi truc Z logarit'); xlabel('x'); zlabel('z log'); Chng trỡnh thc hin comman Window >> Bi 32 bi 32 V hm Z(x,y)=x.e-x2-y2 Chng trinh vit M-file nh sau function Z=nice(X,Y); dx=0.5;dy=1.0; x=[-2:dx:2]; y=[-2:dy:2]; SV:Nguyn Ngc Tnh 30 in1A Trng H Kinh t k thut cụng nghiờp [X,Y]=meshgrid(x,y); Z=X.*exp(-X.^2-Y.^2) mesh(X,Y,Z); view([10,35]) Chng trỡnh thc hin comman Window >> bai32 Z= Columns through -0.0007 -0.0135 -0.0366 -0.0135 -0.0007 -0.0029 -0.0582 -0.1581 -0.0582 -0.0029 -0.0067 -0.1353 -0.3679 -0.1353 -0.0067 -0.0071 -0.1433 -0.3894 -0.1433 -0.0071 0 0 0.0071 0.1433 0.3894 0.1433 0.0071 -0.0071 -0.1433 -0.3894 -0.1433 -0.0071 0 0 0.0071 0.1433 0.3894 0.1433 0.0071 Columns through 0.0067 0.1353 0.3679 0.1353 0.0067 0.0029 0.0582 0.1581 0.0582 0.0029 0.0007 0.0135 0.0366 0.0135 0.0007 ans = Columns through -0.0007 -0.0135 -0.0366 -0.0135 -0.0007 -0.0029 -0.0582 -0.1581 -0.0582 -0.0029 -0.0067 -0.1353 -0.3679 -0.1353 -0.0067 Columns through 0.0067 0.1353 0.3679 0.1353 0.0067 0.0029 0.0582 0.1581 0.0582 0.0029 0.0007 0.0135 0.0366 0.0135 0.0007 SV:Nguyn Ngc Tnh 31 in1A Trng H Kinh t k thut cụng nghiờp bi 33: Dựng cõu lnh Interp1 tớnh ni suy cho hm sau: Y=sin(x), x thuc khong [0:10] Hóy v hm ú bng lnh plot Chng trỡnh M-file nh- sau x=0:10; y=sin(x); xi=0:0.25:10; yi=interp1(x,y,xi); plot(x,y,'o',xi,yi) Chng trỡnh common window >> bai33 SV:Nguyn Ngc Tnh 32 in1A Trng H Kinh t k thut cụng nghiờp bi 34: Cho mch iờn phi tuyn sau: R U0 U1 U2 I(U ) Vi R1= 500 ohm, I(U) = 2*10-4U2A, U0=20V Hóy tớnh I1, U1, U2=U Chng trỡnh vit M-file nh sau R=500; SV:Nguyn Ngc Tnh 33 in1A Trng H Kinh t k thut cụng nghiờp Uo=20; deltaU=1; U=0:deltaU:Uo I1=(Uo/R)*ones(1,length(U))-(1/R)*U I2=2*0E-04*U.^2 plot(U,I1,U,I2) grid Chng trỡnh thc hin comman Window Bi 34: U= Columns through 10 15 16 17 18 Columns 11 through 20 10 11 12 13 14 19 Column 21 20 I1 = Columns through 0.0400 0.0380 0.0360 0.0340 0.0320 0.0300 0.0220 0.0200 0.0180 0.0100 0.0080 0.0060 Columns through 12 0.0280 0.0260 0.0240 Columns 13 through 18 0.0160 0.0140 0.0120 Columns 19 through 21 0.0040 0.0020 I2 = SV:Nguyn Ngc Tnh 34 in1A Trng H Kinh t k thut cụng nghiờp Columns through 10 0 0 0 0 0 0 0 Columns 11 through 20 0 0 Column 21 SV:Nguyn Ngc Tnh 35 in1A Trng H Kinh t k thut cụng nghiờp Phn IV Gii tớch mch in bng Simulink Bi:35: Minh chc nng ca cỏc Product (phộp nhõn), Gain (khuch i tớn hiu u vo), Math Function (trong ny cú mt lng khỏ ln cỏc hm toỏn ó c chun b sn ti ụ Function) Hóy mụ phng phng trỡnh sau: F(t)=80exp(-1/80t)sin(0.25+pi/3) Tham s ca Sine wave Amplitude:1; Frequency (rad/sec):0.25, Phase:pi/3, Sample Time:0 Bi lm SV:Nguyn Ngc Tnh 36 in1A Trng H Kinh t k thut cụng nghiờp bi 36: Hóy hin th c lp ng thi trờn ta XY hai tớn hiu súng hỡnh sin rng ca Bi lm: SV:Nguyn Ngc Tnh 37 in1A Trng H Kinh t k thut cụng nghiờp SV:Nguyn Ngc Tnh 38 in1A Trng H Kinh t k thut cụng nghiờp bi 37: Lm VD s v th mụ phng ca h thng iu khin nh hỡnh v Bi lm SV:Nguyn Ngc Tnh 39 in1A Trng H Kinh t k thut cụng nghiờp bi 41: To súng Bi lm bi 39: To súng hỡnh sin SV:Nguyn Ngc Tnh 40 in1A Trng H Kinh t k thut cụng nghiờp bi 40: Bi lm bi 41: SV:Nguyn Ngc Tnh 41 in1A Trng H Kinh t k thut cụng nghiờp Bi lm SV:Nguyn Ngc Tnh 42 in1A Trng H Kinh t k thut cụng nghiờp Bi 43:Cho s nh hỡnh v Bi lm bi 44: Tuyn tớnh húa h thng: Bi lm SV:Nguyn Ngc Tnh 43 in1A Trng H Kinh t k thut cụng nghiờp SV:Nguyn Ngc Tnh 44 in1A [...]... Chương trình thực hiện trong comman Window >> Bài 28 Bài 29 Vẽ đồ thị hàm z=sin(x)+cos(y) trong 3-D ch-¬ng tr×nh ®-îc viÕt trong M-file nh- sau a=linspace(0,10,100); b=linspace(0,6,100); [x,y]=meshgrid(a,b); z=sin(x)+cos(y); plot3(x,y,z) Chương trình thực hiện trong comman Window SV:Nguyễn Ngọc Tỉnh 27 Điện1A Trường ĐH Kinh tế kỹ thuật công nghiêp >> bai29 bài 30: Vẽ đồ thị hàm z=sin(x)*cos(y) trong... Kinh tế kỹ thuật công nghiêp Kết quả trong Commanwindows: *Bài 20: Cho ma trận sau: A=[3 2/3;2 3/2] h y tÝnh gi¸ trÞ riªng,vÐct¬ riªng ch−¬ng tr×nh ®−îc thùc hiÖn trªn cöa sæ comman Window >> A=sym([3 2/3;2/3 2]); >> [u,v]=eig(A) u= [ 1, 2] [ -2, 1] SV:Nguyễn Ngọc Tỉnh 16 Điện1A Trường ĐH Kinh tế kỹ thuật công nghiêp v= [ 5/3, 0] [ 0, 10/3] *Bài2 1 Giải phương trình sau: Chương trình thực hiện trong... 0.206 s^2 + 0.11 s + 0.625 SV:Nguyễn Ngọc Tỉnh 21 Điện1A Trường ĐH Kinh tế kỹ thuật công nghiêp bài 25: Cho mạch điện như sau; Với R=1.4 ohm, L=2H, C=0.32F, uc(0)=0.5V, iL(0)=0 Tại i=0 mạch điện được đóng vào điện áp u(t)=1(t) Hãy tính các đáp ứng quá độ uc(0), iL(0) và uRL(t), và vẽ các đồ thị đó Sử dụng Matlab để tính quá trình quá độ trong mạch điện tuyến tính ch-¬ng tr×nh ®-îc viÕt trong M-file... u=ones(1,length(t)); [x,y]=lsim(A,B,C1,D,u,t,x0); plot(t,x,t,y); grid SV:Nguyễn Ngọc Tỉnh 22 Điện1A Trường ĐH Kinh tế kỹ thuật công nghiêp Chương trình thực hiện trong comman Window >> bai25 C= 0.3200 bài 26: Sử dụng Matlab để tính và vẽ quá trình quá độ trong mạch điện tuyến tính Cho mạch điện như sau: R U SV:Nguyễn Ngọc Tỉnh U(C) I(R I(L) C 23 L Điện1A Trường ĐH Kinh tế kỹ thuật công nghiêp Với R=10ohm ,L=1h,C=10-3F,... Chương III:Vẽ đồ thị bài 27: Vẽ hàm sin(x) và cos(x) và tan (x),cotan(x) trên cùng một đồ thị Chương trình trong M-file như sau x=-pi:1:pi; y1=sin(x); y2=cos(x); y3=tan(x); plot(x,y1,x,y2,x,y3); title('do thi cac ham luong giac sinx,cosx,tanx'); set(gca,'xtick',-pi:pi/2:pi) set(gca,'Xticklabel',{'-pi','-pi/2','0','pi/2','pi'}) grid on Chương trình thực hiện trong comman Window >> bai27 Bài 28: Vẽ đồ thị... SV:Nguyễn Ngọc Tỉnh 11 Điện1A Trường ĐH Kinh tế kỹ thuật công nghiêp Bài 10: Lập chương trình giải phương trình bậc 2: ax2 + bx +c = 0 khi có nghiệm kép, vô nghiệm, hai nghiệm phân biệt Giải: Chương trình viết trong M-file như sau SV:Nguyễn Ngọc Tỉnh 12 Điện1A Trường ĐH Kinh tế kỹ thuật công nghiêp Kết quả chương trình trong Commanwindows: Bài 11: SV:Nguyễn Ngọc Tỉnh 13 Điện1A Trường ĐH Kinh tế kỹ thuật... Trường ĐH Kinh tế kỹ thuật công nghiêp set(gca,'Zscale','log'); view([125,7]); grid on title('ve 3D voi truc Z logarit'); xlabel('x'); zlabel('z log'); Chương trình thực hiẹn trong comman Window >> Bài 32 Đề bài 32 Vẽ hàm Z(x,y)=x.e-x2-y2 Chương trinh viết trong M-file như sau function Z=nice(X,Y); dx=0.5;dy=1.0; x=[-2:dx:2]; y=[-2:dy:2]; SV:Nguyễn Ngọc Tỉnh 30 Điện1A Trường ĐH Kinh tế kỹ thuật công nghiêp... Kinh tế kỹ thuật công nghiêp Đề bài 33: Dùng câu lệnh Interp1 để tính nội suy cho hàm sau: Y=sin(x), x thuộc khoảng [0:10] Hãy vẽ hàm đó bằng lệnh plot Chương trình trong M-file nh- sau x=0:10; y=sin(x); xi=0:0.25:10; yi=interp1(x,y,xi); plot(x,y,'o',xi,yi) Chương trình trong common window >> bai33 SV:Nguyễn Ngọc Tỉnh 32 Điện1A Trường ĐH Kinh tế kỹ thuật công nghiêp Đề bài 34: Cho mạch điên phi tuyến... mạch điện bằng Simulink Bài: 35: Minh họa chức năng của các khối Product (phép nhân), Gain (khuếch đại tín hiệu đầu vào), Math Function (trong khối này có một lượng khá lớn các hàm toán đã được chuẩn bị sẵn tại ô Function) Hãy mô phỏng phương trình sau: F(t)=80exp(-1/80t)sin(0.25+pi/3) Tham số của khối Sine wave Amplitude:1; Frequency (rad/sec):0.25, Phase:pi/3, Sample Time:0 Bài làm SV:Nguyễn Ngọc Tỉnh... trong comman Window >> a=linspace(0,10,100); >> b=linspace(0,6,100); >> [x,y]=meshgrid(a,b); >> z=sin(x).*cos(y); >> plot3(x,y,z) SV:Nguyễn Ngọc Tỉnh 28 Điện1A Trường ĐH Kinh tế kỹ thuật công nghiêp Đề bài 31: Vẽ 3D với trục Z tuyến tính và trục Z logarit, các hàm sau: X=cos(20.pi.t), y=sin(20.pi.t), z=t Với t=[0.01:0.05:0.09] Chương trình viết trong M-file như sau t=[0.01:0.05:0.99]; x=cos(20*pi*t); ...Trng H Kinh t k thut cụng nghiờp Chng I : Matlab c bn Bai 1: Lp chng trỡnh dang function gii phng trỡnh bc >> %vi du >> bt1 Nhap he so a=3... nh sau: SV:Nguyn Ngc Tnh in1A Trng H Kinh t k thut cụng nghiờp Kt qu chong trỡnh Commanwindows: Bài 5: Nhập ma trận A(m.n) tìm hình chữ nhật hình vuông ma trận chứa toàn số mà diện tích lớn Giải:... tu sai = dai=1 SV:Nguyn Ngc Tnh in1A Trng H Kinh t k thut cụng nghiờp rong= -1 sai=0 dai= -1 >> Bài 6: Cho mt qu cu nng cú lng l m=2kg c buc vo mt si dõy di l=1.5m mt ngi nm vo u dõy v quay trũn