untitled2.m

来自「研究了利用混沌相变进行弱信号检测的理论及仿真试验( 对基于)振子初值敏感性检测弱」· M 代码 · 共 27 行

clc;close all;
tic;
a=0:0.02:0.3;
ts=linspace(0,12000,200);
p0=0.05*0.95.^[0:28];
for k=1:29;eval(['P',num2str(k),'=[];']);end;
for k=1:length(a);
    [t,p]=ode23('haha',ts,p0,[],a(k));
    for h=1:29;
        E=p(:,h);
        eval(['P',num2str(h),'=[P',num2str(h),',E];']);
    end
end
figure;mesh(a,ts,P2);% here input Pi
xlabel('\ita','fontsize',14);
ylabel('\itt','fontsize',14);
figure;mesh(a,ts,P3);% here input Pi
xlabel('\ita','fontsize',14);
ylabel('\itt','fontsize',14);
figure;mesh(a,ts,P4);% here input Pi
xlabel('\ita','fontsize',14);
ylabel('\itt','fontsize',14);
figure;mesh(a,ts,P5);% here input Pi
xlabel('\ita','fontsize',14);
ylabel('\itt','fontsize',14);
t=toc