matlab欧拉方程求解微分方程并和ode45对比结果

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、参考论文

略

答1疑2咨3询4企q鹅q号：1762016542