我的代码是
dsolve(Eq_sub_diff_all_value(1,1),VC1(0)==V0_C1)
dsolve(Eq_sub_diff_all_value(1,2),VC2(0)==V0_C2)
Vsolve=dsolve(Eq_sub_diff_all_value,VC1(0)==V0_C1,VC2(0)==V0_C2);
其中Eq_sub_diff_all_value(1,1)为V_Dsc2 - V_D1 - VC1(t) - C_C1r_c1diff(VC1(t), t) == 0
Eq_sub_diff_all_value(1,2)为V_D1 - V_Dsc1 + Vdc + VC2(t) - r_s3*(I_IR1 - C_C2diff(VC2(t), t)) + C_C2r_c2diff(VC2(t), t) == 0
前两行是对两个一阶微分方程单独求解
解分别为V_Dsc2 - V_D1 + exp(-t/(C_C1r_c1))(V_D1 - V_Dsc2)
V_Dsc1 - V_D1 - Vdc + I_IR1r_s3 + exp(-t/(C_C2*(r_c2 + r_s3)))(V_D1 - V_Dsc1 + Vdc - I_IR1r_s3)
都是正确的
最后一行是合起来求解
解分别为-exp(-t/(C_C1r_c1))(V_Dsc2 - V_D1 + exp(t/(C_C1r_c1))(V_D1 - V_Dsc2))
-exp(-t/(C_C2r_c2 + C_C2r_s3))(V_Dsc1 - V_D1 - Vdc + I_IR1r_s3 + exp(t/(C_C2r_c2 + C_C2r_s3))(V_D1 - V_Dsc1 + Vdc - I_IR1r_s3))
比正确解多乘了一个-exp，这是为什么？