function solve_odes
clear all;clc
y0=[ ];    % [%s] [%Mg]在t=0时刻的数值，即其初值，待填入
tspan=; %时间t的范围，待填入
[t,y]=ode45(@my_odes],[2 0])
figure(1),plot(t,y(:,1));     % [%s]对时间t作图
figure(2),plot(t,y(:,2);      % [%Mg]对时间t作图

function dydt=my_odes
beta= ;    %β的数值
Ip= ;    Ip的数值
......   %以下依次把已经参数的数值填入。

dydt(1)=        ;  % d[%s]/dt的表达式
dydt(2)=        ;  % d[%Mg]/dt的表达式

5楼: Originally posted by anxiaoxin2 at 2015-03-06 08:37:00
谢谢版主。。。这个不需要m文件之类的吗？

8楼: Originally posted by anxiaoxin2 at 2015-03-17 17:06:06
你好，按你提供的输入到matlab中了，显示
dydt(1)=(6*b*Ip*U1*ip*Pm*Km*(*exp(-Tmix/T)-CaOe))/(Ds*Wm*Ps)+(6*Q*irt*Pm*Kb*(*exp(-Tmix/T)-Mge))/Wm/Db+0.4167*k*Pm**exp(-Tmix/T)**exp(-Tmix/T)+A*Pz*Ks*Ls*(*ex ...

14楼: Originally posted by anxiaoxin2 at 2015-03-18 14:09:03
你好，我在以下表达式中,该如何表示？直接写y1,y2吗？dydt(1)=        ;  % d/dt的表达式
dydt(2)=        ;  % d/dt的表达式...

16楼: Originally posted by anxiaoxin2 at 2015-03-18 15:26:43
版主，这个还是运行不了啊。。。
这是代码，你帮我看看呗
function solve_odes
clear all;clc
y0=;    %  在t=0时刻的数值，即其初值，待填入
tspan=; %时间t的范围，待填入
=ode45(‘my_odes’, tspan, y0) ...

function solve_odes
clear all;clc
y0=[0.03,0.0001];    % [%S] [%Mg]在t=0时刻的数值，即其初值，待填入
tspan=[0:1:100]; %时间t的范围，待填入
[t,y]=ode15s(@my_odes, tspan, y0);
figure(1),plot(t,y(:,1));     % [%s]对时间t作图
figure(2),plot(t,y(:,2));      % [%Mg]对时间t作图

function dydt=my_odes(t,y)
b= 0.36392;   %β脱硫剂穿透比
Ip=0.233;  %脱硫剂喷粉速率
U1=0.4;  %CaO在脱硫剂中的质量分数
U3=0.6;  %Mg在脱硫剂中的质量分数
ip=25.85;  %粉剂在铁水中停留时间
irt=12.45;  %Mg气泡在铁水中的停留时间
Pm=7138;  %铁水密度
Ps=1380;  %脱硫剂密度        
Pz=3500;  %顶渣密度
Wm=210000;  %铁水重量
Wso=5000;  %初始渣量
Ds=0.001;  %脱硫剂径粒
Db=9.75E-05;  %Mg气泡直径
Km=1.42313E-05;  %铁水中硫的传质系数
Kb=0.001031861;  %硫在Mg气泡上的传质系数
Ks=4.19217E-09;  %硫在渣中的传质系数
k=5.20879E-05;  %表观反应常数
Tmix=240;  %混匀时间
Q=0.02777778;  %喷吹气量
A=8.24;  %渣金反应面积
Ls=169.6743082;  %硫的分配系数
B=0.9;  %Mg在铁液中的溶解比例
n=0.4;  %粉剂进入渣层的有效系数
SCaOe=5.71697E-05;  %CaO脱硫平衡时硫含量
SMge=1.42026E-05;  %Mg脱硫平衡时硫含量
S4=0.5;  %初始渣中硫含量
y1=y(1);
y2=y(2);

S3=0.1;

S2=(S4*Wso+(S3-y1*exp(-Tmix/t))*Wm)/(Wso+Ip*t*n);



dydt(1)=(6*b*Ip*U1*ip*Pm*Km*(y1*exp(-Tmix/t)-SCaOe))/(Ds*Wm*Ps)+(6*Q*irt*Pm*Kb*(y1*exp(-Tmix/t)-SMge))/Wm/Db+0.4167*k*Pm*y1*exp(-Tmix/t)*y2*exp(-Tmix/t)+A*Pz*Ks*Ls*(y1*exp(-Tmix/t)-S2/Ls)/Wm;
dydt(2)=100*B*U3*Ip/Wm-0.75*(0.4167*k*Pm*y1*exp(-Tmix/t)*y2*exp(-Tmix/t));

dydt=dydt';

17楼: Originally posted by 月只蓝 at 2015-03-18 16:29:34
S3没有定义，由于不知道S3的数值，以下代码中S3=0.3；程序已经可以跑通

function solve_odes
clear all;clc
y0=;    %  在t=0时刻的数值，即其初值，待填入
tspan=; %时间t的范围，待填入
=ode15s(@my_odes ...

19楼: Originally posted by anxiaoxin2 at 2015-03-19 08:09:30
版主，早上好啊。。。以上代码是怎么运行的？具体操作步骤能说下吗？我这还是不会运行。。。...