MATLAB作业2参考 答案(2018).docVIP

MATLAB作业二参考答案 1、试求下面线性代数方程的解析解与数值解，并检验解的正确性。 【求解】求出A, [A;B] 两个矩阵的秩，可见二者相同，所以方程不是矛盾方程，应该有无穷多解。 A=[2,-9,3,-2,-1; 10,-1,10,5,0; 8,-2,-4,-6,3; -5,-6,-6,-8,-4]; B=[-1,-4,0; -3,-8,-4; 0,3,3; 9,-5,3]; [rank(A), rank([A B])] ans = 4 4 用下面的语句可以求出方程的解析解，并可以验证该解没有误差。 x0=null(sym(A)); x_analytical=sym(A)\B; syms a; x=a*[x0 x0 x0]+x_analytical x = [ a+967/1535, a-943/1535, a-159/1535] [ -1535/1524*a, -1535/1524*a, -1535/1524*a] [ -3659/1524*a-1807/1535,-3659/1524*a-257/1535,-3659/1524*a-141/1535] [ 1321/508*a+759/1535, 1321/508*a-56/1535, 1321/508*a-628/1535] [ -170/127*a-694/307