你好同学,这个用符号工具箱直接给定初值计算即可:
syms I1(t) I2(t) I10 I20 eq = diff([I1;I2], t, 1)==[-24,40;-16,-40]*[I1;I2] + [6;4]; conds = [I1(0);I2(0)]==[I10;I20] [I1,I2] = dsolve(eq, conds)