ode，从符号到数字代码

发布的罗兰美国舒尔，2010年4月8日

19次(过去30天 0喜欢| 2的评论

在最近的一篇文章中求解常微分方程，读者Jason想知道是否有一种方法来消除使用的手工代数步骤符号数学工具箱．目前还没有一个固定的方法来做它，所以结果是一个MATLAB程序，但你可以在MuPAD，符号数学工具箱的符号引擎。现在，我想向您展示如何在MuPAD找到封闭形式的解决方案的情况下生成ODE解决方案。

内容

定义符号变量

让我们首先定义我们需要描述简谐振子方程的变量，使用与上面链接的相同的符号。

信谊xt米k

解决系统

使用函数dsolve，我们就能找到解决办法。

索尔= dsolve (' m * D2x + k * x (t)”，“x(0) = 0”，“Dx(0) = 1”，“t”）;

重新格式化适合张贴的解决方案显示。

csol = char (sol);disp (csol (1:45) disp (csol(46:结束)

(m * exp ((t * (- k * m) ^ (1/2)) / m)) / (2 * (- k * m) ^ (1/2)) - m / (2 * exp ((t * (- k * m) ^ (1/2)) / m) * (- k * m) ^ (1/2))

从MuPAD转换解决方案

使用matlabFunction，我接下来生成一个函数句柄的MATLAB版本的解决方案。

f = matlabFunction (sol);

再次，重新格式化表达式，以便在post中显示得更漂亮。

cf = func2str (f);disp (cf(一)disp ([' 'cf (66)):)

@ (k、m t) m。* exp (t。*√(m - k . *)) / m)。* 1.0. /√(m - k . *) * (1.0. / 2.0) - m。* exp (- (t。* sqrt (m - k . *)) / m)。* 1.0. /√(m - k . *) * (1.0. / 2.0)