Question

如下MATLAB程序有什么问题吗？

二阶电路如图所示，已知电压源us=10v,R1=4 ,R2=0.2 ,C1时开关闭合，试绘出电容电压uc的波形。电路的状态方程表示如下(见式（3）): (3)Matlab程序：M文件为：d.mFunction yd=g(t,y) %函数定义U=10;R1=4;R2=0.2;C=1;L=1; %赋值yd=[-(1/r1*c)]*y(1)-(1/c)*y(2)+(1/(r1*c))*u (1/l)*y(1)-(R2/L)*Y(2)]; %描述方程 调用函数并作图（见图5）：tspan=[0,40]; %定义求数值解得时间区间y0=[0;0]; %变量的初始值[t,YY]=ode45('d',tspan,y0); %解方程plot(t,YY(:,1)); %作图xlabel('t'),ylabel('uc'); %描述x轴，y轴





Answer 1

题主提供的MATLAB程序有问题，请按下列代码对照一下。应该这样来改

tspan=[0,40];%定义求数值解得时间区间

y0=[0;0];%变量的初始值

[t,YY]=ode45(@g,tspan,y0);%解方程

plot(t,YY(:,1)); %作图

xlabel('t'),ylabel('u_c'); %描述x轴，y轴

function yd=g(t,y) 

U=10;R1=4;R2=0.2;C=1;L=1; %赋值

yd=[-(1/R1*C)*y(1)-(1/C)*y(2)+(1/(R1*C))*U;(1/L)*y(1)-(R2/L)*y(2)]; %函数定义

end

运行结果