PI-LNK409 源代码

function phi = dAngle(a,b) d = a-b; phi = mod(d,2*pi); if phi > pi phi = 2*pi-phi; end

[t,y]=nlbound('c7fun4','c7fun3',[0,pi/2],[-1,1],1e-8); plot(t,y); set(gca,'xlim',[0,pi/2]); y0=tan(t-pi/4); norm(y(:,1)-y0) …

缺省线型所绘的图形 set(0,'DefaultAxesLineStyleOrder',{'-o',':s','--+'}) set(0,'DefaultAxesColorOrder',[0.4,0.4,0.4]) x = 0:pi/10:2…

[t,y]=nlbound(@c3fun4,@c3fun3,[0,pi/2],[-1,1],1e-8); plot(t,y); set(gca,'xlim',[0,pi/2]); y0=tan(t-pi/4); norm(y(:,1)-y0) n…

Poincare截面 　　在相空间中适当（要有利于观察系统的运动特征和变化，如截面不能与轨线相切，更不能包含轨线）选取一截面，在此截面上某一对共轭变量如x1和x.1取固定值，称此截面为Poincare截面，相空间的连续轨迹与Poincare截面的交点成为截…

% Poincare_section[绘制庞加莱截面图] betaa=0.25; F=1.093; v=2/3; Poin=inline(['[x(2);',... '-2*betaa*x(2)-x(2).^2.*sin(x(…

% Poincare_section[绘制庞加莱截面图] [t,x]=ode45(@Poincare,[0,2800],[0,1.5,0]); x(:,2)=mod(x(:,2),2*pi)-pi; phi0=pi*2/3; % 选择phi=2…

%function ff=zonghemoxing1(t) clear all; clc; close all; format long; a=pi/2; h=0.01; t=1; w=pi/3; %T=2*pi/w; bt=rem(t,h…

%function ff=zonghemoxing1(t) clear all; clc; close all; format long; a=pi/4; h=1.000000000000000e-004; t=1; w=2*pi*1000/6…

clear all; clc; close all; format long; a=pi/2; h=0.001; t=5; w=2*pi*100; bt=rem(t,h); bbt=t-bt; n1=bbt/h; k=0; fo…