MATLAB解一阶微分方程

最新推荐文章于 2024-09-02 09:59:06 发布

da_kao_la

最新推荐文章于 2024-09-02 09:59:06 发布

阅读量2.1w

点赞数 6

分类专栏： MATLAB 文章标签： MATLAB 一阶微分方程龙格-库塔法 ode45

本文链接：https://blog.csdn.net/da_kao_la/article/details/80258895

版权

MATLAB里ode45命令用4阶龙格-库塔解一阶微分方程（常系数/变系数）都能解。

首先把微分方程化为dx/dt = f(x)的标准形式写成函数放在一个文件里，再写主程序用ode45调用这个文件。举例如下：

函数文件：方程的定义（一定要用函数名命名）

% di/dt = L-1u - L-1(WL+R)i
function iprime = hw11_2(t,i)
L = [.0072, 0, 0, .1, .0054, 0;
    0, .007, 0, 0, 0, .0026;
    0, 0, .001, 0, 0, 0;
    1.5*.1, 0, 0, 2.5, .125, 0;
    1.5*.0054, 0, 0, .125, .0068, 0;
    0, 1.5*.0026, 0, 0, 0, .0016];
W = zeros(6,6);
W(1,2) = -100*pi;
W(2,1) = 100*pi;
R = diag([.002,.002,.002,.4,.015,.015]);
u = zeros(6,1);
u(4) = 400;
iprime = L\u - L\(W*L+R)*i;

主程序文件：调用ode45

% solve function hw11_2
theta0 = 0;
h = 1e-3;                       % step
i0 = zeros(1,6);
i0(4) = 400/.4;
ts = 0:h:1;                     % time span
[t,i] = ode45('hw11_2',ts,i0);  %