【1】ode求解微分方程只能求解显式方程,因此需要把高阶微分方程化解为一阶方程组。
是几阶微分方程就通过代换法写出几个式子,然后化成等效的微分方程,用matlab编辑函数之后,再用ode求解。
这个式子,先写出[y1=y,y2=y'],然后得到等效“微分”方程[y1'=y2,y2'=(A/B)ty1],注意要写成显式微分方程的形式,最后的y1',y2'的“=”右边不能再有微分形式了,都要用y1,y2替代掉。
【2】sort函数的使用
sort(G,1)表示把数组每列升序排列,结果输出给a,b中没列的元素对应a中该位置的元素在G中对应列的相对位置,如a的第一列“2,3,4”,分别是G中第一列的“2,1,3”行,因此b的第一列即为“2,1,3”.sort(G,2)表明把G数组的每一行按照升序排列,其余解释与sort(G,1)相同。