摘要: 1、 直接写方程 diff(y, 2) == k*(1 - y^2)*diff(y) - y,借助 odeToVectorField 和matlabFunction 和ode45函数求解。优点是,构造的方程非常直接; syms y(t); [V] = odeToVectorField(diff(y, 阅读全文
posted @ 2022-05-12 11:44 bH1pJ 阅读(257) 评论(0) 推荐(0) 编辑
摘要: 范德波振子的李雅普诺夫指数 keng code %% 计算范德波振子的李雅普诺夫指数 clear all;close all;clc; Z=zeros(1,100); %保存结果; d0=1e-6; %蝴蝶煽动翅膀; ks = linspace(0,5,100); transient = 50; f 阅读全文
posted @ 2022-05-12 10:59 bH1pJ 阅读(82) 评论(0) 推荐(0) 编辑
摘要: how to calcuate the Lyapunov exponents https://www.youtube.com/watch?v=-xSNqJQRoo4&ab_channel=nptelhrd code如下 %% 计算范德波振子的李雅普诺夫指数 clear all;close all;c 阅读全文
posted @ 2022-05-12 10:59 bH1pJ 阅读(63) 评论(0) 推荐(0) 编辑