1、内容简介
matlab欧拉方程求解微分方程并和ode45对比结果
2、内容说明
略
3、仿真分析
clc
close all
clear
%% ode45方法
y0 = [8.5;2;1];%初始值
[T,Y] = ode45(@CW4_201718,[0 10],y0);
figure
plot(T,Y)
xlabel 时间/t
ylabel xxx
title ode45
%% 欧拉方法
Ni = 1000; % 总步数
h = 0.001; % 步长
y1(1:Ni,1:3) = 0;
y1(1,1) = 8.5;y1(1,2) = 2;y1(1,3) = 1; %初始值
for i = 2:Ni
y1(i,1:3) = y1(i-1,1:3)+h*func(y1(i-1,1:3))';
end
t1 = 0:h:(Ni-1)*h;
figure
plot(t1,y1)
xlabel 时间/t
ylabel xxx
title 欧拉方法
4、参考论文
略
