matlab ode45求解常微分方程模板
Posted studyer_domi
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了matlab ode45求解常微分方程模板相关的知识,希望对你有一定的参考价值。
1、内容简介
matlab ode45求解常微分方程模板
307-可以交流、咨询、答疑
2、内容说明
matlab ode45求解常微分方程模板
3、仿真分析
主函数:
clc
close all
clear
tspan = [0 100];
x0 = [0 -1 0 -1]';
[t,x]=ode45('fun',tspan,x0);
figure
plot(t,x(:,1))
xlabel 时间/s
ylabel u1
figure
plot(t,x(:,2))
xlabel 时间/s
ylabel du1/dt
figure
plot(t,x(:,3))
xlabel 时间/s
ylabel u2
figure
plot(t,x(:,4))
xlabel 时间/s
ylabel du2/dt
子函数
function y = fun(t,x)
%---------------
y=zeros(4,1);
y(1)=x(2);
y(2)=-x(2)+x(4)+100*sin(10*t)+150*(x(3)-x(1));y(2)=y(2)/2.25;
y(3)=x(4);
y(4)=x(2)-x(4)+150*(x(1)-x(3));y(4)=y(4)/2;
4、参考论文
略
以上是关于matlab ode45求解常微分方程模板的主要内容,如果未能解决你的问题,请参考以下文章