Matlab dasturini kimyoviy injineringda qo'llash

Matlab dasturini kimyoviy injineringda qо'llash Reja Matlab dasturini kimyoviy injineringda qо'llash usuli. Matlab dasturi yordamida masalalar echish. Ishdan maqsad: aloqa operatsiyalari, mantiqiy operatsiyalar va shartli operatorlar bilan tanishish, ulardan tarmoqlangan hisoblashlarda foydalanish malakasiga ega bulish. MATLABda ikki turdagi fayllar muhim ahamiyatga ega - kengaytirilgan mat . va m. Birinchisi binar (qo'sh) fayllar bo'lib, ularda o'zgaruvchilarning qiymatlari saqlanadi, ikkinchisi fayl o'zida tekst fayllarini saqlab, ularda funksiya sistemasi va buyruqlarni aniqlaydigan tashqi programmalar mavjud. Funksiya va buyruqlarning katta qismi aynan shunga tegishli. Ko'p darchali to'g'irlovchi -redaktor bo'sh oynada m-faylni redaktorlash uchun buyruqlar qatoridan Edit buyrug'ini yoki buyruqlar menyusidan File New M-fiLe tanlanadi. MATLAB DASTURINI KIMYOVIY INJINERINGDA QО'LLASH Davriy distilyatsiya kolonkasida suv-etanol qorishmasini distillyatsiyalash 1-Mashq. 9X-Y plot Past bosimda Lto V=0.9; xd=0.996; PhiL=1; xf = 0.63; xw =0.008; X-Y plot Yuqori bosimda Lto V=0.9; xd=0.99; PhiL =0.5; xf =0.79; xw =0.63 2-Mashq. Oddiy majmuini Reaktori mustaqil o'zgaruvchilar bir funktsiyasi sifatida yozilgan tenglamalar[A-1, 124-bet]. function Reactor=Reac(t,x) global Temperature Qgen %Units, kg, mol,m^3 y s %Lambda==M; Mu==s Differensial tenglamalar uchun reaktorning (i:1) masalasini yordamida bajaramiz: Reactor=Reac(w,x) Natijalarni olish uchun Subplot operatori yordamida aniqlaymiz. Subplot(n, m, x) Bu erda n-sonlar qatorining qiymati, m-sonlarning ustun qiymati, x-manzili.[A-1 143-bet] 3-Mashq. M file: ReactorSO2: [a,b]=ode15s('Reac',[0 20],[0,750,202650]); plot(a,b(:,1)) xlabel('W (kg)') ylabel('X') Figure plot(a,b(:,2)) xlabel('W') ylabel('T (K)') Figure plot(a,b(:,3)) xlabel('W') ylabel('P (Pa)') M file: Reac: function Reactor=Reac(w,x) X=x(1); T=x(2); Presion=x(3); k=9.8692e-3*exp(-176008(1.8*(T-273.15)+491.67) -110.1*log((1.8*(T-273.15)+491.67))+912.8); Kp=0.0031415*exp(42311(1.987*(1.8*(T-273.15)+491.67)) -11.24); Ta=700; Pto=202650; ySO2o=0.10; yO2o=0.11; yN2o=0.79; PhiO2=yO2oySO2o; PSO2o=Pto*ySO2o; Fto=0.02153; Fao=Fto*ySO2o; epsilon=-0.05; G=0.433; ph=0.45; rhoo=0.866; Dp=0.00457; visc=3.72e-5; Ac=0.0016; rhob=542; To=770; sumCp=300.85-0.0402*T+0.0018*T^2-9.071e-8*T^3; dCp=-21.535+0.0789*T-7.112*10^(-5)*T^2+2.447e-8*T^3; deltaHr=-98480-21.535*(T-298)+0.0395*(T^2-298^2)-2.371*10^ (-5)*(T^3-298^3)+6.117*10^(-9)*(T^4-298^4); Ucorreg=9.3; if X, =, ...