GIẢI ĐỀ, ỨNG DỤNG MÁY TÍNH

GIẢI ĐỀ Câu 1.1 1.1.a x la vecto hang y la vecto cot z la vecto hang 1.1.b >>x=[3 5] >>z=[1 3] >>Tong =x+z >>Tich =x.*z 1.1.c >>x=[3 5] >>y=[2;1;3] >>Tong=x+y % khong the tinh >>Tich=y*z Câu 1.2 1.2.a x=linspace(1,10,10) la mang 10 phan tu gom:1 10 1.2.b >>format bank >>x=linspace(1,10,10) >>y=x.^2./(1+x) 1.2.c >>x=linspace(1,10,10) >>z=y-x Câu 2.1.a >>format bank >>A=[5*10^-6 -0.0009 0.0527 (0.0668-0.6)] >>roots(A) 2.1.b function f= amdo(x) global RH f=1-exp(-0.0003053*(25+134.136)*x^1.2164)-RH; ………………………………………………………… global RH RH =0.6; fsolve('amdo',0.1) 2.2 x=linspace(0,100,100); RH1=5*10^-6.*x.^3-0.0009.*x.^2+0.0527.*x+0.0668; RH2=1-exp(-0.003053*(25+134.136).*x.^1.2164); plot(x,RH1,x,RH2) grid xlabel('x') ylabel('RH') title('am can bang dau xanh') legend('RH1','RH2') câu 3.1 function NDR=cdau(T) Cn=4.1762-9.0864*10^-5.*T+5.4731*10^-6.*T.^2; Cp=2.0082+1.2089*10^-3.*T-1.1329*10^-6.*T.^2; Cc=1.5488+1.9625*10^-3.*T-5.9399*10^-6.*T.^2; NDR=0.4*Cn+0.2*Cp+0.4*Cc; …………………………………………………… m=input('nhap khoi luong thuc pham:'); Q=m*quad('cdau',30,150,'reltol',1e-6) disp('kj') 3.2 function C=nhietdungrieng(T) global xn xp xc Cn=4.1762-9.0864*10^-5.*T+5.4731*10^-6.*T.^2; Cp=2.0082+1.2089*10^-3.*T-1.1329*10^-6.*T.^2; Cc=1.5488+1.9625*10^-3.*T-5.9399*10^-6.*T.^2; C=xn*Cn+xp*Cp+xc*Cc; ………………………………………… global xn xp xc xn=input('nhap gia tri xn:'); xp=input('nhap gia tri xp:'); xc=input('nhap gia tri xc:'); m=input('nhap gia tri m:'); Q=m*quad('nhietdungrieng',30,150,'reltol',1e-6) Câu : Ứng dụng matlab truyền nhiệt- truyền khối Đề : Một phẳng đặt gần vật đen phát xạ nhiệt với nhiệt độ T R =1000 0C Mặt lại có nhiệt độ 15 0C Dùng phương pháp lặp để tính T1 Với g T2 + g * TR4 =3.24*10^-10 ; T1 = + g * T13 Giải TR=1273; T2=288; T1=300; g=3.24*10^-10; f=abs((T2+g*TR^4)/(1+g*T1^3)-T1); while f>0.01 f=abs((T2+g*TR^4)/(1+g*T1^3)-T1); T1=abs(T2+g*TR^4)/(1+g*T1^3) End Câu Đề :Giải phương trình y’=r*y dùng phương pháp Runge-Kutta vẽ đồ thị khoảng ≤ t ≤ 0.5 r = -1 với điều kiện ban đầu y0=2 Giải function dtdy =f54(t,y) global r dtdy =r*y; global r r =-1; [t,y]=ode45 ('f54',[0 0.5],2); ycx =2*exp(-1*t); plot(t,y,'o',t,ycx) xlabel('t') ylabel('y') [t,y] =ode45 ('f54',0.1,2) (co the bo) Câu Đề: Dùng vòng lặp for vẽ vòng tròn bán kính r = 1,2,3,4 cm đồ thị Giải function[x,y]=tron(r) for r=1:1:4 alpha=linspace(0,2*pi,100); x=r*cos(alpha); y=r*sin(alpha); plot(x,y,'.') axis('equal') title('vongtron') hold on end Câu Đề:Giải hệ phương trình tìm ethanol sản xuất trình lên men Giải function E=ethanol(R,Y) A=[0,0,-1,-2;2,-1,-2,-1;0,-2,0,-6;R,0,-1,0]; (A ,b phu thuoc hpt da cho) b=[6-6*(144/(180*Y));3-6*(144/(180*Y));7-12*(144/(180*Y));0]; x=A\b; E=x(4); Câu Đề: Thiết lập function file để vẽ đường tròn Giải function[x,y]=vongtron(r) r=input('nhap ban kinh:') theta=linspace(0,2*pi,100); x=r*cos(theta); y=r*sin(theta); plot(x,y,'o','color','b') axis('equal'); xlabel('truc x') ylabel('truc y') title('vong tron ban kinh r') grid on % hinh nen dang luoi câu Đề: Giải phương trình y’=r*y dùng phương pháp Euler vẽ đồ thị khoảng ≤ t ≤ 0.5 r = -10 với điều kiện ban đầu y0=2; ∆ t =0.02; y(t) = 2*e-10t Giải r=-10; h= 0.02; y(1)= 2; k=0; for time =[h:h:0.5] k=k+1; y(k+1)=y(k) + r*y(k)*h; end t=[0:h:0.5]; ycx=2*exp(-10*t); plot(t,y,'o',t,ycx) xlabel('t') ylabel('y') grid on % hinh nen dang luoi Câu 10 Đề: chưng cất dạng mẽ tách benzen khỏi hổn hợp dung dịch benzen/toluene Giải function f = chungcat(x) global L f = 100*(x/0.6)^0.625*((1-x)/0.4)^(-1.625)-L; …………………………………………………… global L L =70; fzero('chungcat', 0.4)