日志
- 分享 Collective Matlab codes of Duffing-Van der Pol Oscillator
- An example of a nonlinear dynamical system with numerous application in engineering is the single well-forced Duffing-Van der Pol oscillator: d^2y_1/dt^2-/Mu*(1-y_1^2)*d ...
- 分享 Cobweb of Logistic maps
- %Author: Thomas Lee function cobweb(x0,tol,a,b,n) %function cobweb(x0,tol,a,b,n) %x0=starting iteration value for x %tol=how small you want the difference to be between adjacent iteration values %a=lower x bound %b=upper x bound %n=number of x values you want to have between a and b % % EX: % %cob ...
- 分享 Basin of attraction using Newton's method
-
% Author: Thomas Lee % basin of attraction % using Newton's method to track initial values to % convergence of the roots for z^3 -1 % 3 blocks for each root, z1=1,z2,z3 as following x = -1:0.001:1; y = -1:0.001:1; = meshgrid(x,y); n = 300; z2 = -1./2. + i*sqrt(3.)/2.; z3 = -1./2. - i*sqr ...
- 分享 Bifurcation and MLE of Henon map
-
clear all d0=1e-12; A= ; for a=0:0.001:1.4 x1=0;y1=0; x2=0;y2=d0; lsum=0; for j=1:500 x3=a-x1*x1+0.3*y1; y1=x1; & ...
- 分享 Henon映射的分叉图的Matlab实现
-
% Author: Thomas Lee % E-mail: lixf1979@126.com % Corresponding: School of Mathematics, Physics and Software Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China %Plotting the bifurcation of the 2D Henon Map b=0.3; niter=2000; x=0; y=0; A= ; %for a=0 ...
- 分享 Lozi混沌映射系统的Matlab实现
- % Author: Thomas Lee % E-mail: lixf1979@126.com % Corresponding: School of Mathematics, Physics and Software Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China % if you want to get more information, please refer to one of the publishedarticles of the author: %李险峰 ...
- 分享 带有周期激励的非自治动力系统的Poincare截面程序
- %具体的讨论可以见我的个人空间！ %这是duffing方程的Poincare截面的程序，经过试验和调试是没有错误的！ %这是定义的duffing方程的式子！ %要调用poincare.m! function dy=duffing(t,y); global det global w global r dy= ; %Duffing 方程的Poincare截面程序！ %经过试验和调试是没有错误的！ %要调用duffing.m! func ...
- 分享 单摆方程的Poincare截面Matlab实现
- % Author: Thomas Lee % E-mail: lixf1979@126.com % Corresponding: School of Mathematics, Physics and Software Engineering, Lanzhou Jiaotong University, Lanzhou 730070, China function ydot=pend(t,y,n) F=1.15; a=0.5; ydot= ; clear; %set innitial conditions as y and dy/dt y0= ; ic=1; n=2/3; %inte ...
