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';

2楼: Originally posted by CelestialCYJ at 2015-03-05 16:11:26
你不给数值怎么算

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的表达式