bài tập matlab (file đính kèm .m)

Nguyen Thi Loan D4H3 Bài 1.1 Hãy sử dụng lệnh Matlab để thực phép tính sau: a, 25,4+17*(34/4,2)-2,5 >> 25.4+17*(34/4.2)-2.5 KẾT QUẢ: ans = 160.5190 b, cos(5,3)+sin(3,7) >> cos(5.3)+sin(3.7) KẾT QUẢ : ans = 0.0245 c, sin2(12,4)-cos1/2(22,4) >> (sin(12.4)^2)-(cos(22.4))^1/2 KẾT QUẢ : ans = 0.4862 d, e2,5+tg(21,7) >> exp(2.5)+tan(21.7) KẾT QUẢ : ans = 11.8828 e, S = (2,7+j3,2)-(2+j1,5) >> S=(2.7+3.2i)-(2+1.5i) KẾT QUẢ : Nguyen Thi Loan D4H3 S= 0.7000 + 1.7000i f) Hãy tính giá trị biểu thức, biểu thị kết hình với n chữ số sau dấu phẩy đông/tĩnh(e/f ) với liệu lấy theo chữ họ, tên đệm tên người giải bảng 1.1 Họ tên: Nguyễn Thị Loan >> a=22;b=13;c=43;d=24; xi=[0.5 0.82 0.97 1.73 2.14 3.74]; yi=((cos(a^3+xi.^2)).^2-log(b)+d)./(sqrt(c^2+2*xi.^2)); fprintf('yi=%10.8e dv\n',yi) KẾT QUẢ: yi=4.99153581e-001 dv yi=4.99666272e-001 dv yi=5.03827107e-001 dv yi=5.04593819e-001 dv yi=5.13889510e-001 dv yi=5.10955546e-001 dv Bài 1.2 a, Thực thao tác cần thiết để hiển thị tên người giải lên hình: >> disp('Nguyen Thi Loan') KẾT QUẢ: Nguyen Thi Loan b, Sử dụng lệnh Matlab xác định giá trị dòng điện chạy mạch: >> U= 116.2+110.4i; Z=4.25+5.34i; disp('dong dien chay mach la') I=U/Z Nguyen Thi Loan D4H3 KẾT QUẢ: dong dien chay mach la I= 23.2596 - 3.2485i c, Biểu thị kết tính toán dòng điện dạng Imjα với Im module α góc pha: >> disp('module dong dien la:') Im=abs(I) disp('goc pha dong dien la:') theta=angle(I) disp('bieu thi lai gia tri dong dien:') I=Im*exp(j*theta) KẾT QUẢ: module dong dien la: Im = 23.4853 goc pha dong dien la: theta = -0.1388 bieu thi lai gia tri dong dien: I= 23.2596 - 3.2485i Bài 1.3 Điện áp định mức mạng điện U kV, điện thành phần trở tác dụng Ω, công suất truyền tải đường dây S=P+j*Q Hãy áp dụng lệnh Matlab để xác định tổn thất điện áp dụng lệnh fprintf để biểu thị kết hình với n chữ số sau dấu phẩy tĩnh/động(t/đ), biết thời gian tổn thất công suất cực đại năm τ h Các liệu tính toán lấy bảng 1.2 với họ tên người giải Nguyen Thi Loan D4H3 >> P=45.33; Q=36.2; U=10; R=5.43; t=4620; deltaA=((P^2+Q^2)/(U^2))*R*t; disp('ton that dien la:') fprintf('delta A=%10.2f\n',deltaA) KẾT QUẢ: ton that dien la: delta A= 844226.53 Bài 1.4 a, Chuyển Z1 dạng số mũ : >> Z1= 4+3i; Zm=abs(Z1); theta=angle(Z1); Z1=Zm*exp(j*theta);[Zm theta] KẾT QUẢ: ans = 5.0000 0.6435 b, Chuyển Z2 dạng đại số: >> Z2=2.71*exp(j*(pi/12)) KẾT QUẢ: Z2 = 2.6177 + 0.7014i Nguyen Thi Loan D4H3 c, Hãy tính giá trị biểu thức ghi lại kết dạng mũ, đối số (argument) biểu thị khoảng –pi +pi >> Z2=2.71*exp(j*(pi/12)); Z1=5*exp(j*0.6435); Z3=1.82*exp(j*(-1.2)); Z4=sqrt(3)-2i; C=Z1^2*Z2/Z3+Z4; Cm=abs(C); theta=angle(C); [Cm theta] KẾT QUẢ: ans = 34.8797 2.7828 Bài 2.1 Lập trình giải lại 1.3 >> function dA=tonthatdiennang (P,Q,U,R,t0) dA=((P^2+Q^2)/U^2)*R*t0; %Bieu thuc tinh ton that dien nang fprintf('%7.2f\n') lưu vào m.file Tại cửa sổ ta cần gõ sau: >> dA=tonthatdiennang(45.33,36.2,10,5.43,4620) KẾT QUẢ : dA = 8.4423e+005 Bài 2.2 Hãy xây dựng hàm xác định điện trở mạch gồm n nhánh song song áp dụng hàm vừa xây dựng để tính điện trở tương đương mạch gồm nhánh song song với điện trở lấy theo chữ đầu họ tên người giải bảng 2.1: Nguyen Thi Loan D4H3 R1=4.78; R2=5.35; R3=6,73; R4=4,67; R5=3,78 >> %R=[R1 R2 R3 ]; function Rss=machss2(R) g=1./R; Rss=1/(sum(g)) fprintf('%7.3f\n') lưu vào m.file Tại sổ ta cần gõ sau: >> Rss=machss2([4.78 5.35 6.73 4.67 3.78]) KẾT QUẢ: Rss = 0.9771 Bài 2.3 Hãy hiển thị hàm f1(x), cho bảng 2.2 với chữ đầu tên người giải, xác định giá trị hàm với x=0,3 x= 2,5 >> f1=inline('x^4+3*x^3-6*x^2+1.7*x^1/2','x') y1=f1(0.3) y2=f1(2.5) KẾT QUẢ: f1 = Inline function: f1(x) = x^4+3*x^3-6*x^2+1.7*x^1/2 y1 = -0.1959 y2 = 50.5625 Bài 2.4 Nguyen Thi Loan D4H3 Biểu thị hàm véc tơ F=[f2,f3], cho bảng 2.2 với chữ đầu họ tên người giải, dạng inline objects xác định giá trị hàm ứng với x1 x2 >> F3=inline('[x1^5+x2^3-3 4*x1^2*x2-x1*x2+log(x1)]','x1','x2') y3=F3(1.5,2.4) KẾT QUẢ: F3 = Inline function: F3(x1,x2) = [x1^5+x2^3-3 4*x1^2*x2-x1*x2+log(x1)] y3 = 18.4177 18.4055 Bài 2.5 Hãy sử dụng vòng lặp while cho toán xác định giá trị lớn n mà có tổng S=1^3+2^3+…+n^3 nhỏ 250 >> n=1; S=0; while S+(n+1)^3> A=rand(6,4) KẾT QUẢ : A= 0.3724 0.0527 0.4177 0.6981 0.1981 0.7379 0.9831 0.6665 0.4897 0.2691 0.3015 0.1781 0.3395 0.4228 0.7011 0.1280 0.9516 0.5479 0.6663 0.9991 0.9203 0.9427 0.5391 0.1711 b) Ma trận B có số cột gấp đôi số cột ma trận gốc >> B=rand(6,8) KẾT QUẢ : B= 0.6787 0.2769 0.4387 0.7094 0.9597 0.8909 0.8143 0.7577 0.0462 0.3816 0.7547 0.3404 0.9593 0.2435 0.8308 0.7431 0.0971 0.7655 0.2760 0.5853 0.3922 0.8235 0.7952 0.6797 0.3517 0.5472 0.9293 0.5853 0.2238 0.1386 0.3500 0.5497 0.6555 0.6948 0.1869 0.6551 0.7513 0.1493 0.1966 0.9172 0.1712 0.3171 0.4898 0.1626 0.2551 0.2575 0.2511 0.2858 Nguyen Thi Loan D4H3 c) Hãy thiết lập ma trận C gồm phần tử từ dòng thứ dến dòng thứ từ cột thứ dến cột thứ ma trận gốc >> C=A(1:2,1:3) KẾT QUẢ : C= 0.3724 0.0527 0.4177 0.1981 0.7379 0.9831 d)Hãy thiết lập ma trận D kích thước ma trận A >>D=rand(6,4) KẾT QUẢ : D= 0.0326 0.4607 0.1909 0.3846 0.5612 0.9816 0.4283 0.5830 0.8819 0.1564 0.4820 0.2518 0.6692 0.8555 0.1206 0.2904 0.1904 0.6448 0.5895 0.6171 0.3689 0.3763 0.2262 0.2653 e, Hãy thiết lập ma trận ngẫu nhiên E kích thước ma trận D >> E=rand(6,4) KẾT QUẢ: E= 0.8244 0.9063 0.9827 0.4253 0.5985 0.8797 0.3127 0.4709 Nguyen Thi Loan D4H3 0.7302 0.8178 0.1615 0.6959 0.3439 0.2607 0.1788 0.6999 0.5841 0.5944 0.4229 0.6385 0.1078 0.0225 0.0942 0.0336 f, Hãy áp dụng lệnh cat để thiết lập ma trận ba chiều F với ma trận liên kết A, D E >> cat(3,A,D,E) KẾT QUẢ: ans(:,:,1) = 0.3724 0.0527 0.4177 0.6981 0.1981 0.7379 0.9831 0.6665 0.4897 0.2691 0.3015 0.1781 0.3395 0.4228 0.7011 0.1280 0.9516 0.5479 0.6663 0.9991 0.9203 0.9427 0.5391 0.1711 0.0326 0.4607 0.1909 0.3846 ans(:,:,2) = 0.5612 0.9816 0.4283 0.5830 0.8819 0.1564 0.4820 0.2518 0.6692 0.8555 0.1206 0.2904 0.1904 0.6448 0.5895 0.6171 0.3689 0.3763 0.2262 0.2653 Nguyen Thi Loan D4H3 xichmax1 = 11.5005 xichmax2 = 14.4194 xichmay = 135.4012 1.4 Độ lệch trung bình bình phương đại lượng x1, x2 y: >> A=std(x1,1)% Do lech trung binh binh phuong cua dai luong x1 B=std(x2,1)% Do lech trung binh binh phuong cua dai luong x2 C=std(y,1)% Do lech trung binh binh phuong cua dai luong y KẾT QUẢ: A= 10.9103 B= 13.6794 C= 128.4528 1.5 Hệ số tương quan giá trị quan sát đại lượng Rx1_y Rx2_y: >> Rx1_y=corrcoef(x1,y)% He so tuong quan cua cac gia tri quan sat x1 va y Rx2_y=corrcoef(x2,y)% He so tuong quan cua cac gia tri quan sat x2 va y KẾT QUẢ: Rx1_y = 1.0000 0.9938 0.9938 1.0000 Rx2_y = Nguyen Thi Loan D4H3 1.0000 0.9895 0.9895 1.0000 Xây dựng hàm hồi quy thực nghiệm đánh giá sai số với hàm: 2.1 Bậc hai : y=f(x1) ; y=f(x2) >> x1=[3.13 7.84 13.06 15.78 18.16 23.34 25.92 27.63 33.53 40.21]; x2=[8.33 10.75 12.86 23.85 28.34 33.23 35.67 39.44 44.34 49.45]; y=[74.3 99.05 136.2 191.9 235.2 278.6 315.7 341.8 421.8 478.6]; a=polyfit(x1,y,2) b=polyfit(x2,y,2) ssb=(Y-y).^2; er=sqrt(sum(ssb)./size(x))*100/mean(Y); ert=er(1) KẾT QUẢ : a= 0.0358 10.1681 25.3768 b= 0.1010 3.6945 49.8446 ert = 96.9661 2.2 Tuyến tính : y=f(x1) ; y=f(x2) : >> c=polyfit(x1,y,1) d=polyfit(x2,y,1) er=sqrt(sum(ssb)./size(x))*100/mean(Y); ssb=(Y-y).^2; ert=er(1) KẾT QỦA : Nguyen Thi Loan D4H3 c= 11.7010 13.2323 d= 9.2918 -8.6735 ert = 96.9661 2.3 Hàm hồi quy bội : y=f(x1,x2) : >> x1=[3.13 7.84 13.06 15.78 18.16 23.34 25.92 27.63 33.53 40.21]'; x2=[8.33 10.75 12.86 23.85 28.34 33.23 35.67 39.44 44.34 49.45]'; y=[74.3 99.05 136.2 191.9 235.2 278.6 315.7 341.8 421.8 478.6]'; XX=[ones(size(x1)) x1 x2]; a=XX\y Y=XX*a; ssb=(Y-y).^2; er=sqrt(sum(ssb)./size(x))*100/mean(Y); ert=er(1) KẾT QUẢ : a= 2.2332 7.1792 3.6793 ert = 4.7227 Bài 4.3 Xác định giá trị nội suy ma trận y ứng với giá trị x >> x=[2:7]'; Nguyen Thi Loan D4H3 y=[x x.^1.5 x.^2.8 x.^4]; xi=[2.5 3.5 6.2]; yi=interp1(x,y,xi) KẾT QUẢ : yi = 1.0e+003 * 0.0025 0.0040 0.0143 0.0485 0.0035 0.0066 0.0351 0.1685 0.0062 0.0155 0.1672 1.5170 Bài 4.4 Xây dựng hàm y khoảng biến thiên x hàm nội suy theo liệu bảng 4.2 đánh giá sai số so với hàm gốc >> x=-2:2; y=2.5*sin(x); xx=-2:0.5:2; yy=2.5*sin(xx); yn=interp1(x,y,xx,'nearest'); y1=interp1(x,y,xx,'linear'); yc=interp1(x,y,xx,'cubic'); ys=interp1(x,y,xx,'spline'); plot(x,y,'o',xx,yy,xx,yn,xx,y1,xx,yc,xx,ys) grid; legend('data','function','bac thang','tuyen tinh hoa','lap phuong mo phong','lap phuong',2); sbt=max(abs(yn-yy)); Nguyen Thi Loan D4H3 stt=max(abs(y1-yy)); slpm=max(abs(yc-yy)); slp=max(abs(ys-yy)); disp('sbt stt slpm slp') fprintf('%g'), disp([sbt,stt,slpm,slp]) KẾT QUẢ: sbt stt slpm slp 2.5364 2.5364 2.5364 2.8143 Bài 5.2 Vẽ đồ thị hàm số y khoảng biến thiên x (xem liệu bảng 5.2 Dữ liệu lấy ứng với chữ đầu theo họ người giải): Họ Nguyễn: x=0:2pi; y=2.5*x^2+3.7*x+6; a) Với trục tọa độ chia theo tỷ lệ log số 10; >> x=0:2*pi; y=2.5*x.^2+3.7*x+6; loglog(x,y), grid xlabel('x,log10'); ylabel('y,log10') title('Do thi y=f(x) theo lenh loglog(x,y)') Nguyen Thi Loan D4H3 Do thi y=f(x) theo lenh loglog(x,y) y,log10 10 10 10 10 x,log10 b, Với trục y chia theo tỉ lệ log số 10: >> x=0:2*pi; y=2.5*x.^2+3.7*x+6; semilogy(x,y), grid xlabel('x'); ylabel('y,log10') title('Do thi y=f(x) theo lenh semilogy(x,y)') Nguyen Thi Loan D4H3 Do thi y=f(x) theo lenh semilogy(x,y) 10 y,log10 10 10 10 x c, Với trục x chia theo tỷ lệ log số 10: >> x=0:2*pi; y=2.5*x.^2+3.7*x+6; semilogx(x,y), grid xlabel('x,log10'); ylabel('y') title('Do thi y=f(x) theo lenh semilogx(x,y)') Nguyen Thi Loan D4H3 Do thi y=f(x) theo lenh semilogx(x,y) 120 100 y 80 60 40 20 0 10 x,log10 Bài 5.3 Xây dựng đồ thị hàm d khoảng biến thiên phi(xem liệu bảng 5.1) hệ tọa độ cực( liệu lấy ứng với chữ đầu theo tên đệm người giải) >> phi=0:0.01:5*pi; d=exp(cos(phi))-2*cos(4*phi); polar(phi,d) Do thi ham d=exp(cos(phi))-2cos(4phi) toa cuc 90 120 60 Bài 5.4 150 30 180 210 330 240 300 270 Vẽ đồ thị hàm Z khoảng biến thiên x y (xem liệu bảng 5.2) với bề mặt đặc; bề mặt với đường mức; bề mặt suốt bề mặt thẳng đứng (thác nước) Nguyen Thi Loan D4H3 phần trang giấy.( Dữ liệu lấy ứng với chữ đầu theo tên người giải) L: z=sin(x^2+y^2); x= 2:5; y=2:5 >> [x,y]=meshgrid(2:5); z=sin(x^2+y^2); subplot(2,2,1);surf(x,y,z);grid on; title('do thi ham so voi be mat dac'); subplot(2,2,2);meshc(x,y,z);grid on; title('do thi ham so voi cac duong muc'); subplot(2,2,3);meshz(x,y,z);grid on; title('do thi ham so voi be mat dung'); subplot(2,2,4);mesh(x,y,z);grid on;hidden off; title('do thi ham so voi be mat suot'); thi ham so voi be mat dac thi ham so voi cac duong muc 1 0 -1 -1 6 2 2 thi ham so voi be mat dung 0 2 Bài 6.1 4 thi ham so voi be mat suot -1 6 4 -1 6 2 Nguyen Thi Loan D4H3 Tìm nghiệm phương trình dạng đa thức y=f(x) với liệu cho bảng 6.1 Lấy theo chữ đầu tên người giải >> p=[1 0.7 -1.5 0.4 -5]; x=roots(p) x= -1.2652 + 0.9939i -1.2652 - 0.9939i 1.3688 0.2307 + 1.1653i 0.2307 - 1.1653i Bài 6.2 a, Phương pháp đồ thị: >> x=-5:0.01:5; y=2.1*sin(x); plot(x,y), grid 2.5 1.5 0.5 -0.5 -1 -1.5 -2 -2.5 -5 -4 -3 -2 -1 >> ginput Khi xuất dấu chữ thập, ta đưa dấu đến điểm cắt đồ thị trục hoành, ấn chuột trái ấn enter, kết quả: Nguyen Thi Loan D4H3 ans = -3.1452 -0.0073 -0.0115 -0.0073 3.1452 -0.0073 b, Dùng lệnh fzero: >> x0=fzero('2.1*sin(x)',0) x1=fzero('2.1*sin(x)',1.0) x2=fzero('2.1*sin(x)',2.0) x3=fzero('2.1*sin(x)',3.0) x0 = x1 = 1.5485e-024 x2 = 3.1416 x3 = 3.1416 Bài 6.3 >> A=[2 10 3;11 0.7 6.6 0.8;5 -3.2 4;0 4.5 -4.2]; b=[6;-5;6;5]; x=A\b KẾT QUẢ: x= -1.4457 -1.4838 1.8294 Nguyen Thi Loan D4H3 -0.1668 Bài 6.5 Sao cách ko được?? >> syms x1 x2 f1=x1.^4+x2.^3-7; f2=4*x1.^2*x2-x1*x2+log(x2); [x1,x2]=solve(f1,f2) Bài 9.1 Xác định giá trị dòng điện chạy mạch phương pháp dòng nhánh: >> R=[12.2,15.6,9.4,13.2,2.8]; E=[120,117,0,0,0]; Z=[1,1,1,1,1;R(1),-R(2),0,0,0;0,R(2),-R(3),0,0;0,0,R(3),-R(4),0;0,0,0,R(4),-R(5)] U=[0;E(1)-E(2);E(2);0;0]; I=Z\U; disp('Ket qua la:') disp('I1 I2 I3 I4 I5') fprintf('%g'), disp([I(1),I(2),I(3),I(4),I(5)]) KẾT QUẢ: Z= 1.0000 1.0000 1.0000 1.0000 1.0000 12.2000 -15.6000 0 0 15.6000 -9.4000 0 0 9.4000 -13.2000 0 Ket qua la: 13.2000 -2.8000 Nguyen Thi Loan D4H3 I1 I2 I3 I4 I5 7.7627 5.8785 -2.6910 -1.9163 -9.0340 Bài 9.2 Giải mạch phương pháp dòng điện nhánh với số phức: >> Z=[9+6.5i,5.2+7i,6.2+5.4i,5.4+7i,5.3+7i]; E=[145+88i,170+125i,0,150+110i,0]; Z=[1,1,1,1,1;Z(1),-Z(2),0,0,0;0,Z(2),-Z(3),0,0;0,0,Z(3),-Z(4),0;0,0,0,Z(4),-Z(5)] U=[0;E(1)-E(2);E(2);0;0]; I=Z\U; disp('Ket qua la:') disp('I1 I2 I3 I4 I5') fprintf('%g'), disp([I(1),I(2),I(3),I(4),I(5)]) KẾT QUẢ: Z= 1.0000 1.0000 1.0000 9.0000 + 6.5000i -5.2000 - 7.0000i 1.0000 5.2000 + 7.0000i -6.2000 - 5.4000i 0 0 1.0000 0 0 6.2000 + 5.4000i -5.4000 - 7.0000i 0 5.4000 + 7.0000i -5.3000 - 7.0000i Ket qua la: I1 I2 I3 I4 I5 8.7919 - 1.1577i 15.3431 - 4.5526i -8.6483 + 0.8756i -7.7275 + 2.3741i -7.7591 + 2.4606i Bài 9.3 Xác định đại lượng mạch điện hình sin: >> Um=100;tetau=20; Nguyen Thi Loan D4H3 Z=5.15;gam=63; tetai=tetau-gam; teta=(tetau-tetai)*pi/180; Im=Um/Z; wt=0:0.05:2*pi; u=Um*cos(wt); i=Im*cos(wt+tetai*pi/180); p=u.*i; U=Um/sqrt(2); I=Im/sqrt(2); P=U*I*cos(teta) Q=U*I*sin(teta) S=P+j*Q Kết quả: P= 440.7675 Q= 865.0549 S= 4.4077e+002 +8.6505e+002i Vẽ đồ thị đại lượng: Bài 9.7 >> disp('Cac Gia tri hieu dung cua dong dien chay mach dien'); U=120*exp(j*5); C=225*10^-6;R=3.25; L=120*10^-3;f=50; Nguyen Thi Loan D4H3 omega=2*pi*f; Z1=-j*1/(omega*C); Z2=R; Z3=j*omega*L; Z=Z1+Z2*Z3/(Z2+Z3); I=U/Z Ir=I*Z3/(Z1+Z3) Il=I-Ir disp('I Ir Il') Cac Gia tri hieu dung cua dong dien chay mach dien I= 8.4127 + 0.4975i Ir = 13.4660 + 0.7963i Il = -5.0533 - 0.2988i ... dv yi=5.10955546e-001 dv Bài 1.2 a, Thực thao tác cần thiết để hiển thị tên người giải lên hình: >> disp('Nguyen Thi Loan') KẾT QUẢ: Nguyen Thi Loan b, Sử dụng lệnh Matlab xác định giá trị dòng... dong dien: I= 23.2596 - 3.2485i Bài 1.3 Điện áp định mức mạng điện U kV, điện thành phần trở tác dụng Ω, công suất truyền tải đường dây S=P+j*Q Hãy áp dụng lệnh Matlab để xác định tổn thất điện... C1 = 0.9134 0.9706 0.9572 0.9595 I= Bài 4.2 Nguyen Thi Loan D4H3 Số liệu thống kê đại lượng x y cho bảng 4.1 lấy theo chữ đầu tên người giải: Sử dụng hàm matlab để xác định đại lượng: 1.1 Giá