代码:
disp("二阶电路零输入")
syms t;
syms Uc;
U0=input('请输入0+时刻电压:');
I0=input('请输入0+时刻电流:');
R=input('请输入等效串联电阻:');
L=input('请输入电感:');
C=input('请输入电容:');
Uc=dsolve('L*C*D2Uc+R*C*DUc+Uc=0','Uc(0)=U0','DUc(0)=-I0/C','t');
ezplot(Uc)
想要由输入量控制最终Uc的表达式,但现在求得的解是一个纯代数表达式
Uc =
(exp(-(t*((-C*(- C*R^2 + 4*L))^(1/2) + C*R))/(2*C*L))*(U0*(-C*(- C*R^2 + 4*L))^(1/2) + 2*I0*L - C*R*U0))/(2*(-C*(4*L - C*R^2))^(1/2)) + (exp((t*((-C*(- C*R^2 + 4*L))^(1/2) - C*R))/(2*C*L))*(U0*(-C*(- C*R^2 + 4*L))^(1/2) - 2*I0*L + C*R*U0))/(2*(-C*(4*L - C*R^2))^(1/2))
导致后续画图也画不出来
有大佬指点指点吗?
disp("二阶电路零输入")
syms t;
syms Uc;
U0=input('请输入0+时刻电压:');
I0=input('请输入0+时刻电流:');
R=input('请输入等效串联电阻:');
L=input('请输入电感:');
C=input('请输入电容:');
Uc=dsolve('L*C*D2Uc+R*C*DUc+Uc=0','Uc(0)=U0','DUc(0)=-I0/C','t');
ezplot(Uc)
想要由输入量控制最终Uc的表达式,但现在求得的解是一个纯代数表达式
Uc =
(exp(-(t*((-C*(- C*R^2 + 4*L))^(1/2) + C*R))/(2*C*L))*(U0*(-C*(- C*R^2 + 4*L))^(1/2) + 2*I0*L - C*R*U0))/(2*(-C*(4*L - C*R^2))^(1/2)) + (exp((t*((-C*(- C*R^2 + 4*L))^(1/2) - C*R))/(2*C*L))*(U0*(-C*(- C*R^2 + 4*L))^(1/2) - 2*I0*L + C*R*U0))/(2*(-C*(4*L - C*R^2))^(1/2))
导致后续画图也画不出来
有大佬指点指点吗?