matlab求解简谐振动问题x=dsolve('D2x=-k*x')结果为C15*exp((-k)^(1/2)*t) +