11楼: Originally posted by 月只蓝 at 2015-03-18 10:32:46
比如t=0时，初值为1 2；那么  y0=;...

4楼: Originally posted by 月只蓝 at 2015-03-05 21:00:42
function solve_odes
clear all;clc
y0=;    %  在t=0时刻的数值，即其初值，待填入
tspan=; %时间t的范围，待填入
=ode45(@my_odes],)
figure(1),plot(t,y(:,1));     % 对时间t作图
figure(2),plot(t,y(:, ...

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

15楼: Originally posted by 月只蓝 at 2015-03-18 14:35:28
,可直接写成y1 y2
代码按如下

function dydt=my_odes(t,y)
y1=y(1);
y2=y(2);
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 ...

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
版主，早上好啊。。。以上代码是怎么运行的？具体操作步骤能说下吗？我这还是不会运行。。。...