|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
本帖最后由 牛小贱 于 2015-1-25 16:48 编辑
- function xdot=fjyqxyctbdj(t,x,options,b)
- global b;
- a=30;
- %b=25;
- %state functionxdot
- xdot=[-0.67*x(1)+x(3)*x(2)-0.27;-x(2)-x(3)*x(1)+a*x(3)-0.42;b*(x(2)-x(3))+0.98*x(1)*x(2)-1.20];
复制代码
运行程序:
- clear
- options=odeset;
- t0=8000;tf=8500;
- x0=[0.05;0.02;0.05];
- [t,x]=ode45('fjyqxyctbdj',[t0,tf],x0,[]);
- for b=10:5:40.0;
- [P,F]=periodogram(x(:,1),[],[],10)
- figure
- plot(F,P);
- plot3(b,F,P);
- %hold on;
- end
复制代码 就是以b变量,画三维图,为什么不行?
|
|