初始值不同,结果略有差异:function hhfor A=3:-0.2:2.2vin=-(A-2)/10-(1-A/6)*(A-2)/10^3;ain=A/10^2-A*(A-2)/(2*10^4);[T,Y]=ode45(@collapse,[10,0.01],[vin,ain]);hold onplot(T,Y);endfunction dy=collapse(t,y)dy=zeros(2,1);dy(1)=(y(2)*(t-y(1))-2/t)*(t-y(1))/((t-y(1))^2-1);dy(2)=(y(2)-2/t*(t-y(1)))*(t-y(1))/((t-y(1))^2-1);...
|