→ g1990219:[t,y]=(@SBR,....)應該不用加func吧? 03/18 14:05
→ pipi17g:3Q.. 03/19 00:27
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