分数阶混沌程序MATLAB.docVIP

下载本文档

126
0
约1.79千字
约 3页
2017-12-20 发布于江西
举报
版权申诉

分数阶混沌程序MATLAB.doc

1、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

分数阶混沌程序MATLAB

q1=0.9;q2=0.9;q3=0.9; h=0.01;N=4000; x0=2;y0=1;z0=3; M1=0;M2=0;M3=0; x(N+1)=[0];y(N+1)=[0];z(N+1)=[0]; x1(N+1)=[0];y1(N+1)=[0];z1(N+1)=[0]; x1(1)=x0+h^q1*35*(y0-x0)/(gamma(q1)*q1); y1(1)=y0+h^q2*(-7*x0-x0*z0+28*y0)/(gamma(q2)*q2); z1(1)=z0+h^q3*(x0*y0-3*z0)/(gamma(q3)*q3); x(1)=x0+35*h^q1*(y1(1)-x1(1)+q1*(y0-x0))/gamma(q1+2); y(1)=y0+h^q2*(-7*x1(1)-x1(1)*z1(1)+28*y1(1)+q2*(-7*x0-x0*z0+28*y0))/gamma(q2+2); z(1)=z0+h^q3*(x1(1)*y1(1)-3*z1(1)+q3*(x0*y0-3*z0)/gamma(q3+2)); for n=1:N M1=(n^(q1+1)-(n-q1)*(n+1)^q1)*35*(y0-x0); M2=(n^(q2+1)-(n-q2)*(n+1)^q2)*(-7*x0-x0*z0+28*y0); M3=(n^(q3+1)-(n-q3)*(n+1)^q3)*(x0*y0-3*z0); N1=((n+1)^q1-n^q1)*35*(y0-x0); N2=((n+1)^q2-n^q2)*(-7*x0-x0*z0+28*y0); N3=((n+1)^q3-n^q3)*(x0*y0-3*z0); for j=1:n M1=M1+((n-j+2)^(q1+1)+(n-j)^(q1+1)-2*(n-j+1)^(q1+1))*35*(y(j)-x(j)); M2=M2+((n-j+2)^(q2+1)+(n-j)^(q2+1)-2*(n-j+1)^(q2+1))*(-7*x(j)-x(j)*z(j)+28*y(j)); M3=M3+((n-j+2)^(q3+1)+(n-j)^(q3+1)-2*(n-j+1)^(q3+1))*(x(j)*y(j)-3*z(j)); N1=N1+((n-j+1)^q1-(n-j)^q1)*35*(y(j)-x(j)); N2=N2+((n-j+1)^q2-(n-j)^q2)*(-7*x(j)-x(j)*z(j)+28*y(j)); N3=N3+((n-j+1)^q3-(n-j)^q3)*(x(j)*y(j)-3*z(j)); end x1(n+1)=x0+h^q1*N1/(gamma(q1)*q1); y1(n+1)=y0+h^q2*N2/(gamma(q2)*q2); z1(n+1)=z0+h^q3*N3/(gamma(q3)*q3); x(n+1)=x0+h^q1*(35*(y1(n+1)-x1(n+1))+M1)/gamma(q1+2); y(n+1)=y0+h^q2*(-7*x1(n+1)-x1(n+1)*z1(n+1)+28*y1(n+1)+M2)/gamma(q2+2); z(n+1)=z0+h^q3*(x1(n+1)*y1(n+1)-3*z1(n+1)+M3)/gamma(q3+2); end %subplot(1,1,1),plot3(x,y,z);xlabel(x),ylabel(y);zlabel(z); subplot(2,2,1),plot(x,y);xlabel(x),ylabel(y); subplot(2,2,2),plot(x,z);xlabel(x),ylabel(z); subplot(2,2,3),plot(y,z);ylabel(y),ylabel(z); subplot(2,2,4),plot3(x,y,z);xlabel(x),ylabel(y),zlabel(z); grid on A=[max(x) max(y) max(z)]; B=[min(x) min(y) min(z)]; 电路仿真