状态方程:(i',v')=[-r/l-1/l[i[1/l1/c0]*v]+0]*v分别用了ODE45算子求解和matlab控制工具箱(exp3ss.m)