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(:, ...

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