matlab中ode45的例子如何运行
来源:学生作业帮 编辑:神马作文网作业帮 分类:综合作业 时间:2024/10/05 14:22:43
matlab中ode45的例子如何运行
程序是:function dy = rigid(t,y)
dy = zeros(3,1); % a column vector
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,Y] = ode45(@rigid,[0 12],[0 1 1],options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')
程序是:function dy = rigid(t,y)
dy = zeros(3,1); % a column vector
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,Y] = ode45(@rigid,[0 12],[0 1 1],options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')
1、第一步将下面的文件命令保存为rigid.m的文件.
function dy = rigid(t,y)
dy = zeros(3,1); % a column vector
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);
2、在command窗口输入
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,Y] = ode45(@rigid,[0 12],[0 1 1],options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')
function dy = rigid(t,y)
dy = zeros(3,1); % a column vector
dy(1) = y(2) * y(3);
dy(2) = -y(1) * y(3);
dy(3) = -0.51 * y(1) * y(2);
2、在command窗口输入
options = odeset('RelTol',1e-4,'AbsTol',[1e-4 1e-4 1e-5]);
[T,Y] = ode45(@rigid,[0 12],[0 1 1],options);
plot(T,Y(:,1),'-',T,Y(:,2),'-.',T,Y(:,3),'.')
matlab ode45
ode45的问题,请求matlab大神帮助
matlab中如何用ode45 求出微分方程某个点的函数值?
我在编matlab程序的过程中,需要用ode45解一个一阶微分方程组,
请问如何实现matlab解一阶微分方程组(用其自身命令如ode45等怎么写能运行出结果)
matlab ode45 解方程
MATLAB 中的ODE45
大神求解动力学微分方程,用matlab中ode45编程实现!
matlab ode45 求方程(1+x^2)y''=2xy',y(0)=1,y'(0)=3 的数值解,我的运行不行,
用matlab的ode45求解如下微分方程的程序
关于Matlab的一个简单问题,ode45函数的
matlab ode45解线性微分方程时的参数