版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請盡快對帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

本帖產(chǎn)生 1 個(gè) 博學(xué)EPI ，點(diǎn)擊這里進(jìn)行查看

yl19931018

銀蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 218
帖子: 11
在線: 14.3小時(shí)
蟲號(hào): 2840694
注冊: 2013-12-01
性別: MM
專業(yè): 電力電子學(xué)

[求助] 分?jǐn)?shù)階統(tǒng)一混沌系統(tǒng)的離散算法

分?jǐn)?shù)階統(tǒng)一混沌系統(tǒng)的離散算法

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

回復(fù)此樓

» 猜你喜歡

346求調(diào)劑[0856] 已經(jīng)有7人回復(fù)
一志愿山大07化學(xué) 332分四六級已過本科山東雙非求調(diào)劑！已經(jīng)有3人回復(fù)
310求調(diào)劑已經(jīng)有3人回復(fù)
機(jī)械專碩299求調(diào)劑至材料已經(jīng)有4人回復(fù)
070300化學(xué)319求調(diào)劑已經(jīng)有7人回復(fù)
08工科 320總分求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有12人回復(fù)
307求調(diào)劑已經(jīng)有3人回復(fù)
265求調(diào)劑已經(jīng)有3人回復(fù)
301求調(diào)劑已經(jīng)有10人回復(fù)

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

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

feixiaolin

榮譽(yù)版主 (文壇精英)

專家經(jīng)驗(yàn): +518
博學(xué)EPI: 39
應(yīng)助: 942 (博后)
貴賓: 1.275
金幣: 3430
散金: 58785
紅花: 532
沙發(fā): 11
帖子: 24215
在線: 2601.8小時(shí)
蟲號(hào): 2139575
注冊: 2012-11-21
專業(yè): 光學(xué)信息獲取與處理
管轄: 數(shù)學(xué)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
yl19931018: 金幣+50, 博學(xué)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);

贊一下

回復(fù)此樓

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