zxx4477 发表于 2014-3-11 21:28

用符号推到的长公式,怎么赋予具体数值,用于计算?

本帖最后由 zxx4477 于 2014-3-11 21:36 编辑

先写一个M文件,推一个拉格朗日方程,

有5个广义坐标q1-q5,同时定义其一阶导数为dq1-dq5, 还有二阶导数ddq1-ddq5。


然后得出系统动能T=1/2*'*M*;


T是一个很长的式子 后面还有 U=f(q........dq.....)


这个M文件得到一系列用 syms 的 q1~q5 dq1~dq5 ddq1~ddq5表示的 变量T、U、Q、等


接下来要在另外一个M文件中进行计算,用for i=0:h:T, 在每个i时刻有具体数值,


也就是要将上个M文件中的符号变量用到此M文件中做变量。




我的问题是:怎么把之前用符号变量写成的公式,转化成变量。求指导。
<div class="blockcode"><blockquote>syms q1 q2 q3 q4 q5
syms m1 m2
q=;
M=;
T=1/2*q.'*M*q;
%得到T =(m1*q1^2)/2 + (m1*q2^2)/2 + (m1*q3^2)/2 + q4*(q1/2 + q2/2 + q3/2 + (m2*q4)/2) + ...+(m2*q5)/2)



Matlab 学习交流
页: [1]
查看完整版本: 用符号推到的长公式,怎么赋予具体数值,用于计算?