(TIỂU LUẬN) bài tập matlab và simulink

BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KĨ THUẬT TP.HCM KHOA CƠ KHÍ ĐỘNG LỰC BÀI TẬP MATLAB VÀ SIMULINK Giảng viên hướng dẫn: Nguyễn Mạnh Cường Sinh viên thực hiện: Nguyễn Ngọc Hoài MSSV: 19145380 Hồ Chí Minh, ngày 30 tháng năm 2022 Mục lục Bài 1: Bai 2: Câu a Câu b: Bài 3: Câu a Câu b Bài 4: Bài 5: Bài 6: 10 Bai 7: 11 Bài 8: 12 Câu a 12 Câu b 14 Bài 1: Giải phương trình bậc hai sau : f (x 2 )=( a+b ) x +(d −E) x+ d =0 Biện luận trường hợp phương trinh trên, với a,b,d,E nhập từ bàn phím Bài làm : clc close all a=input('nhap a=') b=input('nhap b=') d=input('nhap d=') E=input('nhap E=') x=0 y=(a+b)*x^2+(d-E)*x+d*d y=0 if (a+b)==0 disp('phuong trinh co nghiem') x=-d*d/(d-E) else delta=(d-E)*(d-E)-4*(a+b)*d*d if delta==0 disp('phuong tring co nghiem kep') x=-(d-E)/(2*(a+b)) else if(delta>0) disp('phuong trinh co nghiem phan biet') x1=(-(d-E)+sqrt(delta))/(2*(a+b)) x2=(-(d-E)-sqrt(delta))/(2*(a+b)) else disp('phuong trinh vo nghiem') end end end Bài 2: Cho hàm số: f (x1 )=5 x 1−6 x 1+3 f (x2 )=7 x 2−3 x 2−2 √x1−7 √e x a.Dùng lệnh plot , subplot để vẽ hàm số biết ∈[20 :80 ] x ∈[ 1:100 ], x b.Nếu x=5 tính S=[ f ( x1)+ f ( x2 )]−f (x2 ) + Bài làm Câu a clc close all x1=1:1:100 x2=20:1:80 y1=5*x1.^3-6*x1.^2+3*sqrt(x1)-7 y2=7*x2.^4-3*x2.^3-2*sqrt(exp(x2)) subplot(211) plot(x1,y1,'-b') title('ve thi') xlabel('truc x') ylabel('truc y') grid on subplot(212) plot(x2,y2,'-r') title('ve thi') xlabel('truc x') √f ( x1 ) ylabel('truc y') grid on Đồ thị thu Câu b: clc close all x1=5 x2=5 y1=5*x1.^3-6*x1.^2+3*sqrt(x1)-7 y2=7*x2.^4-3*x2.^3-2*sqrt(exp(x2)) S=(y1+y2)-(y2).^2+sqrt(y1) Kết : x1 = 5 x2 = y1 = 474.7082 y2 = 3.9756e+03 S= -1.5801e+07 Bài 3: Cho hàm số cos( x )−sin ( x ) y1=x + x −8 x +3 sin( x ) cos( x ) e y 2=tan g( x)4 +e( tan g( x )−cos ( x )) a.Dùng lệnh plot vẽ hàm −|tan g (x )3−ecot y ,y ang( x )2 +sin( x )4 +cos ( x)3| b.Dùng lệnh sobplot để vẽ y ,y với x ∈[1:100 ] Nhận xét khác giưa lệnh plot subplot Bài làm : Câu a clc clear close all x1=1:1:100 y1=x1.^4+x1.^3-8*x1.^2+3*sin(x1).*cos(x1).*exp((cos(x1)-sin(x1))) y2=tan(x1.^4)+exp((tan(x1)-cos(x1)))-abs(tan(x1.^3)-exp(cot(x1.^2)) +sin(x1.^4)+cos(x1.^3)) plot(x1,y1,'-',x1,y2,' ') title('do thi ham so') xlabel('gia tri x') ylabel('gia tri y') grid on Đồ thị thu dược: Câu b clc clear close all x=1:1:100 y1=x.^4+x.^3-8*x.^2+3*sin(x).*cos(x).*exp((cos(x)-sin(x))) y2=tan(x.^4)+exp((tan(x)-cos(x)))-abs(tan(x.^3)-exp(cot(x.^2))+sin(x.^4)+cos(x.^3)) subplot(2,1,1) plot(x,y1,'+') title('do thi thu nhat') xlabel('gia tri x') ylabel('gia tri y') grid on subplot(2,1,2) plot(x,y2,'-') title('do thi thu hai') xlabel('gia tri x') ylabel('gia tri y') grid on Đồ thị thu được: Nhận xét: - Lệnh plot vẽ đồ thị hàm số hệ trục tọa độ - Còn lệnh subplot vẽ đồ thị hàm số hệ trục tọa độ riêng biệt Bài 4: Cho x nghiệm phương trình a x 2+ − = b.x c S=ex− x +1 + x+ + x+ + x+ n Tính n+ Với n,a,b,c nhập từ bàn phím Bài làm: clc clear close all a=input('a=') b=input('b=') c=input('c=') n=input('n=') x=sym('x') y=a*x*x+b*x-c y=0 if(a==0) if(b==0) disp('phuong trinh vo nghiem va khong co gia tri s') else x=((-b)/(a)) x=x s=exp(x)-((x+1)/(2))+((x+2)/(3))+((x+3)/(4))+((x+n)/(n+1)) end else delta=b*b-4*a*(-c) if(delta0) disp("Phuong trinh co nghiem phan biet") x1=log((b+sqrt(delta))/(2*a)) x2=log((-b-sqrt(delta))/(2*a)) S1=sqrt((x1)^3-sin(x1))+tan(x1)^3-3*cos(x1)+2cot(x1) S2=sqrt((x2)^3-sin(x2))+tan(x2)^33*cos(x2)+2-cot(x2) else disp("phuong trinh vo nghiem") end end end Bai 7: Nhập số thực x, số nguyên n tính tổng sau : S x3 10 13 Bài làm : clc clear close all x=input('type x = ') n=input('type n = ') if n>=1 t=0;i=1 for i=1:n t=t+(-1)^i*(x+i)/(2.*i) end disp('Sum of S = ') s=3-(x).^3+sqrt(10+t) else disp('Retype n') end Bài 8: Cho hàm số: f (x1 ) 5.x13 6.x12 x1 f (x2 ) 7.x24 3.x23 e x2 x x2 a.Dùng lệnh plot,subplot vẽ đồ thị hai hàm biết 80] b.Nếu x=5 tính x [1 : 100], x [20 : S [ f (x ) f (x )] f (x )2 f (x ) 2 Bài làm: Câu a clc clear close all 14 x1=1:0.5:100 x2=20:0.5:80 y1=5*x1.^3-6.*x1.^2+3.*sqrt(x1)-7 subplot(2,1,1) plot(x1,y1,'-r','Linewidth',2) title('First Graph of a function') xlabel('X axle') ylabel('Y axle') grid on y2= 7*x2.^4-3*x2.^3-2.*sqrt(exp(x2-sqrt(x2)))+x2 subplot(2,1,2) plot(x2,y2,'-g','Linewidth',2) title('Second Graph of a function') xlabel('X axle') ylabel('Y axle') grid on Đồ thị thu được: 15 Câu b clc clear close all x1=5 x2=5 y1=5*x1.^3-6.*x1.^2+3.*sqrt(x1)-7 y2= 7*x2.^4-3*x2.^3-2.*sqrt(exp(x2-sqrt(x2)))+x2 disp('Value of S') S=y1+y2-(y2).^2+sqrt(y1) Kết : 16 x1 = x2 = y1 = 474.7082 y2 = 3.9970e+03 Value of S 17 S= -1.5972e+07 18 ...Mục lục Bài 1: Bai 2: Câu a Câu b: Bài 3: Câu a Câu b Bài 4: Bài 5:... Bài 5: Bài 6: 10 Bai 7: 11 Bài 8: 12 Câu a 12 Câu b 14 Bài 1: Giải phương trình bậc hai... thị hàm số hệ trục tọa độ riêng biệt Bài 4: Cho x nghiệm phương trình a x 2+ − = b.x c S=ex− x +1 + x+ + x+ + x+ n Tính n+ Với n,a,b,c nhập từ bàn phím Bài làm: clc clear close all a=input(''a='')