我只添加一个单独的答案来解释如何使用dsolve,然后如何策划解决方案。
你不使用一个循环来解决它。没有必要改变t在一个循环中,自解,如果能够找到一个dsolve,将t的函数。不管怎样,这些多个调用dsolve保存没有导致任何变量能够那么阴谋。
我可能会尝试:
信谊x (t)
Dx = diff (x, t);
索尔= dsolve (diff (x, t, 2) = = 0.002 * cos(型)sin (x) x (0) = = 0, Dx (0) = = 0)
警告:不能找到明确的解决方案。
在dsolve >(第190行)
索尔=
(空符号)
然而MATLAB放弃,无法找到一个解决方案。这可能意味着它只是需要一点帮助。例如,移位的余弦值转换成一对术语使用身份,因此;
因为(型)= cos (x) * cos (t) + sin (x) *罪(t)
然而这也失败了。
索尔= dsolve (diff (x, t, 2) = = 0.002 * (cos (x) * cos (t) + sin (x) * sin (t)) sin (x), x (0) = = 0, Dx (0) = = 0)
警告:不能找到明确的解决方案。
在dsolve >(第190行)
索尔=
(空符号)
不幸的是,它很容易编写一个微分方程,缺乏一个解决方案,至少dsolve可以处理。dsolve设法找到解决方案,例如在一个更简单的问题,然后从fplot情节可以获得。
索尔= dsolve (diff (x, t) = = sin (x) x (0) = = 1)
索尔=
2 *:exp (t +日志(tan (1/2))))
fplot(溶胶,[0,3 *π])
缺乏一个解决方案,你最好使用一个数值解算器。埃米尔显示怎么做了。
