find matlab 重复_微分方程数值解+ode45+MATLAB代码

最新推荐文章于 2023-08-26 22:37:24 发布

weixin_39945792

最新推荐文章于 2023-08-26 22:37:24 发布

阅读量121

点赞数

文章标签： find matlab 重复

本文链接：https://blog.csdn.net/weixin_39945792/article/details/111663103

版权

MATLAB代码

clear;

%% 赋值及调用
tspan = [-50,50];
x0 = [0;0;0]; %初值
a = 0.2;
b = 0.2;
c = 2.5;
%options = odeset('RelTol',1e-10,'AbsTol',[1e-10 1e-10 1e-10]);
[t,x] = ode45(@myODE45,tspan,x0,[],a,b,c); %[]表示options的内容

%% 选择重复部分，因为开始时的初值不一定是解，先找到一个解作为初值
%tol = 2; %选择强度
%len = size(x,1);
%panDuan = zeros(len,1);
%for i = 1:len
%    if norm(x(i,:)-x(end,:))<tol
%        panDuan(i) = 1;
%    end
%end
%[index] = find(panDuan == 1);
%plot3(x(index(1):end,1),x(index(1):end,2),x(index(1):end,3));

%% 将得到的点作为初值在循环
x0 = x(end,:);
[t,x] = ode45(@myODE45,tspan,x0,[],a,b,c);
plot3(x(:,1),x(:,2),x(:,3));

%% 直接绘图
%plot3(x(:,1),x(:,2),x(:,3));

%% 微分方程组
function [dx] = myODE45(t,x,a,b,c)
dx = zeros(3,1);
dx(1) = -(x(2)+x(3));
dx(2) = x(1)+a*x(2);
dx(3) = b+x(3)*(x(1)-c);
end