洛伦茨方程的运动图像画法，可以控制播放速度

最新推荐文章于 2023-07-13 21:34:35 发布

chuanhoner

最新推荐文章于 2023-07-13 21:34:35 发布

阅读量815

点赞数

文章标签： matlab

本文链接：https://blog.csdn.net/chuanhoner/article/details/121661843

版权

本文介绍了如何在MATLAB中使用洛伦茨方程创建动态运动图像，并通过修改comet3函数内部代码实现播放速度的控制，以达到更好的视觉效果。

摘要由CSDN通过智能技术生成

首次建立m文件，带入洛伦茨方程的参数：

有两种建法：

1：

function dy=lorenz(t,y)
    dy=zeros(3,1); %建立一个三位列向量
    dy(1)=10*(-y(1)+y(2));
    dy(2)=28*y(1)-y(2)-y(1)*y(3);
    dy(3)=y(1)*y(2)-8*y(3)/3;

2：

function dy_dt=lorenz(t,y)
%y=[x,y,z]储存三个变量的矩阵
p=10;
r=28;
b=8/3;

dx_dt=p*(y(2)-y(1));
dy_dt=r*y(1)-y(2)-y(1)*y(3);
dz_dt=y(1)*y(2)-b*y(3);

dy_dt=[dx_dt;dy_dt;dz_dt];

我选用第二种，并配合rungekutta4阶精度来计算数据源：

function [t,y]=rungekutta4(f,tspan,y0,n)
%四阶Runge-Kutta算法
%f:dy/dt=f(t,y);求解区间tspan:[t0,tf];初值y0=y(t0)=[x0,y0,z0];n采样点数
h=(tspan(2)-tspan(1))/n;
t=tspan(1):h:tspan(2);
y=zeros(length(y0),length(t));
y(:,1)=y0;
for m=2:length(t)
    k1=f(t(m-1),y(:,m-1));
    k2=f(t(m-1)+h/2,y(:,m-1)+k1*h/2);
    k3=f(t(m-1)+h/2,y(:,m-1)+k2*h/2);
    k4=