24小时热门版块排行榜

返回列表

本帖产生 1 个博学EPI ，点击这里进行查看

yl19931018

银虫 (初入文坛)

应助: 0 (幼儿园)
金币: 218
帖子: 11
在线: 14.3小时
虫号: 2840694
注册: 2013-12-01
性别: MM
专业: 电力电子学

[求助] 分数阶统一混沌系统的离散算法

分数阶统一混沌系统的离散算法

https://www.codeforge.cn/article/36487

回复此楼

» 猜你喜欢

1楼 2014-08-17 08:34:35

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

回帖置顶 ( 共有1个 )

feixiaolin

荣誉版主 (文坛精英)

专家经验: +518
博学EPI: 39
应助: 942 (博后)
贵宾: 1.275
金币: 3430
散金: 58785
红花: 532
沙发: 11
帖子: 24215
在线: 2601.8小时
虫号: 2139575
注册: 2012-11-21
专业: 光学信息获取与处理
管辖: 数学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
yl19931018: 金币+50, 博学EPI+1, ★★★很有帮助, 谢谢啊！评晚了，刚接触小木虫。 2014-08-22 17:01:25

%This is the program for fractional order Unified system, a simple but fast
%method in time domain.
clc
clear
clf
a=1; %parameter a
q1=0.99;q2=0.99;q3=0.99; %the fractional order
x10=3;x20=1;x30=5;
h=0.01;N=1000;
x1(N+1)=[0];x2(N+1)=[0];x3(N+1)=[0];
xx1(N+1)=[0];xx2(N+1)=[0];xx3(N+1)=[0];
M1=0;M2=0;M3=0;N1=0;N2=0;N3=0;

xx1(1)=x10+h^q1*(25*a+10)*(x20-x10)/(q1*gamma(q1));
xx2(1)=x20+h^q2*((28-35*a-x30)*x10+(29*a-1)*x20)/(q2*gamma(q2));
xx3(1)=x30+h^q3*(x10*x20-x30*(8+a)/3)/(q3*gamma(q3));

x1(1)=x10+h^q1*(25*a+10)*(xx2(1)-xx1(1)+q1*(x20-x10))/gamma(q1+2);
x2(1)=x20+h^q2*((28-35*a-xx3(1))*xx1(1)+(29*a-1)*xx2(1)+q2*((28-35*a-x30)*x10+(29*a-1)*x20))/gamma(q2+2);
x3(1)=x30+h^q3*(xx1(1)*xx2(1)-xx3(1)*(8+a)/3+q3*(x10*x20-x30*(8+a)/3))/gamma(q3+2);

for n=1:N
M1=(n^(q1+1)-(n-q1)*(n+1)^q1)*(x20-x10);
M2=(n^(q2+1)-(n-q2)*(n+1)^q2)*((28-35*a-x30)*x10+(29*a-1)*x20);
M3=(n^(q3+1)-(n-q3)*(n+1)^q3)*(x10*x20-x30*(8+a)/3);

N1=((n+1)^q1-n^q1)*(x20-x10);
N2=((n+1)^q2-n^q2)*((28-35*a-x30)*x10+(29*a-1)*x20);
N3=((n+1)^q3-n^q3)*(x10*x20-x30*(8+a)/3);

for j=1:n
      M1=M1+((n-j+2)^(q1+1)+(n-j)^(q1+1)-2*(n-j+1)^(q1+1))*(x2(j)-x1(j));
      M2=M2+((n-j+2)^(q2+1)+(n-j)^(q2+1)-2*(n-j+1)^(q2+1))*((28-35*a-x3(j))*x1(j)+(29*a-1)*x2(j));
      M3=M3+((n-j+2)^(q3+1)+(n-j)^(q3+1)-2*(n-j+1)^(q3+1))*(x1(j)*x2(j)-x3(j)*(8+a)/3);

      N1=N1+((n-j+1)^q1-(n-j)^q1)*(x2(j)-x1(j));
      N2=N2+((n-j+1)^q2-(n-j)^q2)*((28-35*a-x3(j))*x1(j)+(29*a-1)*x2(j));
      N3=N3+((n-j+1)^q3-(n-j)^q3)*(x1(j)*x2(j)-x3(j)*(8+a)/3);
end

xx1(n+1)=x10+(25*a+10)*h^q1*N1/(q1*gamma(q1));
xx2(n+1)=x20+h^q2*N2/(q2*gamma(q2));
xx3(n+1)=x30+h^q3*N3/(q3*gamma(q3));

x1(n+1)=x10+(25*a+10)*h^q1*(xx2(n+1)-xx1(n+1)+M1)/gamma(q1+2);
x2(n+1)=x20+h^q2*((28-35*a-xx3(n+1))*xx1(n+1)+(29*a-1)*xx2(n+1)+M2)/gamma(q2+2);
x3(n+1)=x30+h^q3*(xx1(n+1)*xx2(n+1)-xx3(n+1)*(8+a)/3+M3)/gamma(q3+2);

end
x1=x1';x2=x2';x3=x3';

plot3(x3(100:end),x1(100:end),x2(100:end),'k');
xlabel('\itz','fontsize',18);
ylabel('\itx','fontsize',18);
zlabel('\ity','fontsize',18);
figure
plot(x1(100:end),x2(100:end),'k');
xlabel('\itx','fontsize',18);
ylabel('\ity','fontsize',18);
figure
plot(x1(100:end),x3(100:end),'k');
xlabel('\itx','fontsize',18);
ylabel('\itz','fontsize',18);
figure
plot(x2(100:end),x3(100:end),'k');
xlabel('\ity','fontsize',18);
ylabel('\itz','fontsize',18);

赞一下

回复此楼

2楼2014-08-17 14:04:48

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 yl19931018 的主题更新

返回列表