作者pipi17g (ray78326)
看板NTHU_ENGI12
标题[认真] 化工程序控制第一次MATLAB作业请参考
时间Fri Mar 18 09:52:27 2011
Problem1
M-file
----------------------------------------------------------
clear
clc
x=[-5:0.01:5];
y1=x.^2;
y2=x-1;
plot(x, y1, 'b--', x, y2, 'm:')
title('x^2 and x-1')
xlabel('x')
ylabel('f(x)')
legend('x^2','x-1')
----------------------------------------------------------
Problem2
微分方程式主程式funcSBR写法
-----------------------------------------------------------
function dy=SBR(t,y)
dy= zeros(3,1); % a column vector
V0 = 5; %单位L
v0 =0.05 %单位L/sec
V= V0 + v0*t;
Cao =0.05; %单位mol/L
Cbo =0.025; %单位mol/L
V0 = 5; %单位L
v0 =0.05 %单位L/sec
K =2.2 %单位L/sec*mol
% Assign value into dy(n) for function output
dy(1) =-K*y(1)*y(2)-(v0/V)*y(1); %y(1)= Ca
dy(2)= -K*y(1)*y(2)+ (v0/V)*( Cbo- y(2)); %y(2)= Cb
dy(3)= K*y(1)*y(2)-(v0/V)* y(3); %y(3)= Cc
end
-------------------------------------------------------------
M-file(解ODE45写法)
clear
clc
Cao =0.05; %单位mol/L
Cbo =0.025; %单位mol/L
[t,y] = ode45(@funcSBR,[0 500],[Cao Cbo 0]);
plot(t,y(:,1),'r.',t,y(:,2),'k-.',t,y(:,3),'b:');
xlabel('时间');
ylabel('浓度');
legend('Ca','Cb','Cc')
----------------------------------------------------------------
以上是个人拙见......
请工院化工的学弟妹....给予指教...感谢
--
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 114.37.161.117
1F：→ g1990219:[t,y]=(@SBR,....)应该不用加func吧? 03/18 14:05
2F：→ pipi17g:3Q.. 03/19 00:27