function dy=yirong(t,y)
dy=zeros(5,1);
L=0.3e-3;
G1=0.42e-3;
G2=2.2e-3;
C1=7.8e-9;
a=-0.7;
b=0.5;
v1=(a+b*abs(y(4)))*y(3);
y(5)=v(1);
dy(1)=(y(2)-v1)/L;
dy(2)=(-y(1)-G2*y(2))/C1;
dy(3)=y(1)+G1*v1;
dy(4)=y(3);
end
function yirong_main
close all
clear all;
format long
tspan=0:5e-4:1;
Y0=[0.02 0 0 0 0];
na =1000;
[~,y]=ode45('yirong',tspan,Y0);
subplot(221)
plot(y(na:end,1),y(na:end,2),'r');
subplot(222)
plot(y(na:end,2),y(na:end,3),'r');
subplot(223)
plot(y(na:end,1),y(na:end,3),'r');
subplot(224)
plot(y(na:end,5),y(na:end,3),'r');
能不能帮我看看我仿真不出来是因为什么啊
不知道你这个问题是否已经解决, 如果还没有解决的话: