分数阶相关程序Word文档下载推荐.docx
- 文档编号:22618227
- 上传时间:2023-02-04
- 格式:DOCX
- 页数:14
- 大小:718.23KB
分数阶相关程序Word文档下载推荐.docx
《分数阶相关程序Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《分数阶相关程序Word文档下载推荐.docx(14页珍藏版)》请在冰豆网上搜索。
y0=1;
z0=1.9;
u=2;
a=1;
step=0.002;
Alpha=0.7;
n=30000;
clear;
data45;
var1=(step^(Alpha))/(gamma(Alpha+2));
var2=(step^(Alpha))/(gamma(Alpha+1));
x_sub=[];
y_sub=[];
z_sub=[];
Gam=[];
Nu=[];
Gam=Alpha;
Nu=1;
x_sub=-u*x0+y0*(z0+a);
y_sub=-u*y0+x0*(z0-a);
z_sub=z0-x0*y0;
x_p=x0+var2*sum(Nu.*x_sub);
y_p=y0+var2*sum(Nu.*y_sub);
z_p=z0+var2*sum(Nu.*z_sub);
x
(1)=x0+var1*(-u*x_p+y_p*(z_p+a))+var1*sum(Gam.*x_sub);
y
(1)=y0+var1*(-u*y_p+x_p*(z_p-a))+var1*sum(Gam.*y_sub);
z
(1)=z0+var1*(z_p-x_p*y_p)+var1*sum(Gam.*z_sub);
fori=1:
n
x_sub=[x_sub;
-u*x(i)+y(i)*(z(i)+a)];
y_sub=[y_sub;
-u*y(i)+x(i)*(z(i)-a)];
z_sub=[z_sub;
z(i)-x(i)*y(i)];
Nu=[(i+1)^Alpha-i^Alpha;
Nu];
Gam
(1)=(i+1)^(Alpha+1)+(i-1)^(Alpha+1)-2*(i^(Alpha+1));
Gam=[i^(Alpha+1)-(i-Alpha)*((i+1)^Alpha);
Gam];
x_p=x0+var2*sum(Nu.*x_sub);
y_p=y0+var2*sum(Nu.*y_sub);
z_p=z0+var2*sum(Nu.*z_sub);
x(i+1)=x0+var1*(-u*x_p+y_p*(z_p+a))+var1*sum(Gam.*x_sub);
y(i+1)=y0+var1*(-u*y_p+x_p*(z_p-a))+var1*sum(Gam.*y_sub);
z(i+1)=z0+var1*(z_p-x_p*y_p)+var1*sum(Gam.*z_sub);
end
plot(x,z);
Xlabel('
x'
);
Ylabel('
z'
title('
q1=q2=q3=0.7'
)
图3-3
data455.m
x0=5;
y0=10;
z0=12;
step=0.01;
q1=0.9;
q2=1;
q3=1;
Fun455.m
functiont=fun455(x,y,z)
data455;
t(:
1)=-u*x+y*(z+a);
2)=-u*y+x*(z-a);
3)=z-x*y;
clc;
data455;
Gamm1=[];
Gamm2=[];
Nu2=[];
Gamm3=[];
Nu3=[];
Gamm1=q1;
Nu1=1;
Gamm2=q2;
Nu2=1;
Gamm3=q3;
Nu3=1;
t=fun455(x0,y0,z0);
x_sub=t(:
1);
y_sub=t(:
2);
z_sub=t(:
3);
x_p=x0+(step^(q1))/(gamma(q1+1))*sum(Nu1.*x_sub);
y_p=y0+(step^(q2))/(gamma(q2+1))*sum(Nu2.*y_sub);
z_p=z0+(step^(q3))/(gamma(q3+1))*sum(Nu3.*z_sub);
t=fun455(x_p,y_p,z_p);
x
(1)=x0+(step^(q1))/(gamma(q1+2))*(t(:
1)+sum(Gamm1.*x_sub));
y
(1)=y0+(step^(q2))/(gamma(q2+2))*(t(:
2)+sum(Gamm2.*y_sub));
z
(1)=z0+(step^(q3))/(gamma(q3+2))*(t(:
3)+sum(Gamm3.*z_sub));
t=fun455(x(i),y(i),z(i));
1)];
2)];
3)];
Nu1=[(i+1)^q1-i^q1;
Nu1];
Nu2=[(i+1)^q2-i^q2;
Nu2];
Nu3=[(i+1)^q3-i^q3;
Nu3];
Gamm1
(1)=(i+1)^(q1+1)+(i-1)^(q1+1)-2*(i^(q1+1));
Gamm1=[i^(q1+1)-(i-q1)*(i+1)^q1;
Gamm1];
Gamm2
(1)=(i+1)^(q2+1)+(i-1)^(q2+1)-2*(i^(q2+1));
Gamm2=[i^(q2+1)-(i-q2)*(i+1)^q2;
Gamm2];
Gamm3
(1)=(i+1)^(q3+1)+(i-1)^(q3+1)-2*(i^(q3+1));
Gamm3=[i^(q3+1)-(i-q3)*(i+1)^q3;
Gamm3];
x_p=x0+(step^(q1))/(gamma(q1+1))*sum(Nu1.*x_sub);
y_p=y0+(step^(q2))/(gamma(q2+1))*sum(Nu2.*y_sub);
z_p=z0+(step^(q3))/(gamma(q3+1))*sum(Nu3.*z_sub);
t=fun455(x_p,y_p,z_p);
x(i+1)=x0+(step^(q1))/(gamma(q1+2))*(t(:
y(i+1)=y0+(step^(q2))/(gamma(q2+2))*(t(:
z(i+1)=z0+(step^(q3))/(gamma(q3+2))*(t(:
q1=0.9,q2=q3=1'
n=10000
图3-4
Data455.m
n=4000;
q1=1;
q2=0.73;
q1=1,q2=0.73,q3=1'
N=4000
N=50000
图3-5
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 分数 相关 程序