| 24小時(shí)熱門(mén)版塊排行榜 |
| 2 | 1/1 | 返回列表 |
| 查看: 387 | 回復(fù): 1 | ||||
| 本帖產(chǎn)生 1 個(gè) 博學(xué)EPI ,點(diǎn)擊這里進(jìn)行查看 | ||||
yl19931018銀蟲(chóng) (初入文壇)
|
[求助]
分?jǐn)?shù)階統(tǒng)一混沌系統(tǒng)的離散算法
|
|||
|
分?jǐn)?shù)階統(tǒng)一混沌系統(tǒng)的離散算法 https://www.codeforge.cn/article/36487 |
榮譽(yù)版主 (文壇精英)
![]() |
專家經(jīng)驗(yàn): +518 |
|
%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 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 08工科 320總分 求調(diào)劑 +6 | 梨花珞晚風(fēng) 2026-03-17 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 二本跨考鄭大材料306英一數(shù)二 +3 | z1z2z3879 2026-03-17 | 3/150 |
|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +12 | 墨墨漠 2026-03-18 | 13/650 |
|
|
[考研]
|
.6lL 2026-03-18 | 8/400 |
|
|
[考研] 296求調(diào)劑 +6 | www_q 2026-03-18 | 10/500 |
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +5 | tcxiaoxx 2026-03-20 | 5/250 |
|
|
[考研] 265求調(diào)劑 +12 | 梁梁校校 2026-03-19 | 13/650 |
|
|
[考研] 319求調(diào)劑 +3 | 小力氣珂珂 2026-03-20 | 3/150 |
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考研] 0856調(diào)劑,是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
[考研] 085410人工智能專碩317求調(diào)劑(0854都可以) +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
[考研] 081700化工學(xué)碩調(diào)劑 +3 | 【1】 2026-03-16 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +10 | 妮妮ninicgb 2026-03-15 | 14/700 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 268求調(diào)劑 +6 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 6/300 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|