MATLAB微分方程组求解全部代码

发布网友 发布时间：2022-03-25 19:58

我来回答

共2个回答

热心网友 时间：2022-03-25 21:28

w=100*pi;R=1;L=0.01;t=0:0.000001:2;Us=220.*sqrt(2).*sin(100.*pi.*t).*((1/6)<100.*t-fix(100.*t)<1);i=dsolve('Di=(Us-R*i)/L','i(0)=0','t')i = (Us - Us*exp(-(R*t)/L))/Ri=(220.*sqrt(2).*sin(100.*pi.*t).*((1/6)<100.*t-fix(100.*t)<1)-220.*sqrt(2).*sin(100.*pi.*t).*((1/6)<100.*t-fix(100.*t)<1).*exp(-(1.*t)/0.01))/1;plot(t,i)

热心网友 时间：2022-03-25 22:46

任意一个微分方程都可以吗？如：
syms y(t)
mu=20;
V = odeToVectorField(diff(y, 2)-mu*(1 - y^2)*diff(y)+y,y(1)==2,y(2)==0);
M = matlabFunction(V,'vars', {'t','Y'});
sol = ode45(M,[0 200],[2 0]);
x = linspace(0,100,1000);
y = deval(sol,x,1);
y1=deval(sol,x,2);
subplot(2,1,1)
plot(x,y);
grid on
subplot(2,1,2)
plot(x,y1);
grid on追问这个表达是什么，或者求的是什么？

追答这就是vanderpol方程的解程序