|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
为什么没有分叉啊 下面是我的程序 是否什么地方有错 期待解答
function dy=dddd(t,y)
global ff
dy=[y(2);-0.1*y(2)-48.704*y(1)-25.35e4*y(1)^3+ff*cos(2*t)];
clc
clear all
global f ww ff
ww=2;
f=0.1:0.1:10;
www=ww;
for i=1:length(f)
ff=f(i);
period=1/www;
tspan=[0:period/100:100*period];
[t,y]=ode45('dddd',tspan,[1e-8 0]);
plot(f(i),y(1000:100:end,2),'k.','MarkerSize',2);
hold on;
end
|
|