求解下列线性规刘问题(要求分别用Matlab和Lingo编程) ,其中的矩阵A =放在Matlab数据文件data. mat中。
其中x0是给定的x(t)的初始值,xp0是任意给定的x(1)的初始值,fixed_:x0和fixed_xp0是与xp0同维数的列向量,其分量为1表示需要保留的初值,为0表示需要求解的初始值。若fixed_x0和fixed_xp0等于空矩阵[],表示允许所有的初值分量可以发生变化。分别用显式和隐式解法求下列微分方程的数值解
比较四种方法的计算精度,并体会显式格式与隐式格式的区别.
用 Lingo软件求解:
式中:c=[6,8,4,2]' ,Q是三对角线矩阵,主对角线上元索全为-1,两条次对角线上元素全为2。