试验五矩阵运算与方程组求解.DOC

项目五 矩阵运算与方程组求解 实验3 线性方程组 实验目的 熟悉求解线性方程组的常用命令,能利用Mathematica命令各类求线性方程 组的解. 理解计算机求解的实用意义. 基本命令 1.命令NullSpace,给出齐次方程组的解空间的一个基. 2.命令LinearSolve,给出非齐次线性方程组的一个特解. 3.解一般方程或方程组的命令Solve见Mathematica入门. 实验举例 求齐次线性方程组的解空间 设为矩阵,为维列向量,则齐次线性方程组必定有解. 若矩阵的 秩等于,则只有零解;若矩阵的秩小于,则有非零解,且所有解构成一向量空间. 命令 NullSpace给出齐次线性方程组的解空间的一个基. 例3.1 (教材 例3.1) 求解线性方程组 输入 Clear[A]; A={{1,1,?2,?1},{3,?2,?1,2},{0,5,7,3},{2,?3,?5,?1}}; NullSpace[A] 则输出 {{?2,1,?2,3}} 说明该齐次线性方程组的解空间是一维向量空间,且向量(?2,1,?2,3)是解空间的基. 注:如果输出为空集{ },则表明解空间的基是一个空集,该方程组只有零解. 例3.2 求解线性方程组 输入 Clear[A]; A={{1,1,2,-1},{3,-2,-3,2},{0,5,7,3},{2,-3,-5,-1}}; Nullspace[A] 输出为 { } 因此解空间的基是一个空集,说明该线性方程组只有零解. 例3.3 (教材 例3.2) 向量组是否线性相关? 根据定义, 如果向量组线性相关, 则齐次线性方程组 有非零解. 输入 Clear[A,B]; A={{1,1,2,3},{1,?1,1,1},{1,3,4,5},{3,1,5,7}}; B=Transpose[A]; NullSpace[B] 输出为 {{?2,?1,0,1}} 说明向量组线性相关,且 非齐次线性方程组的特解 例3.4 (教材 例3.3) 求线性方程组 的特解. 输入 Clear[A,b]; A={{1,1,?2,?1},{3,?2,?1,2},{0,5,7,3},{2,?3,?5,?1}}; b={4,2,?2,4} LinearSolve[A,b] 输出为 {1,1,?1,0} 注: 命令LinearSolve只给出线性方程组的一个特解. 例3.5 求线性方程组 的特解. 输入 Clear[A,b]; A={{1,1,2,-1},{3,-2,-1,2},{0,5,7,3},{2,-3,-5,-1}}; b={4,2,2,4} Linearsolve[A,b] 输出为 Linearsolve::nosol:Linear equation encountered which has no solution. 说明该方程组无解. 例3.6 向量是否可以由向量 ，， 线性表示? 根据定义, 如果向量可以由向量组线性相关, 则非齐次线性方程组 有解. 输入 Clear[A,B,b]; A={{1,2,-3,1},{5,-5,12,11},{0,5,7,3},{1,-3,6,3}}; B=Transpose[A]; b={2,-1,3,4}; Linearsolve[B,b] 输出为 {,,0} 说明可以由线性表示,且 例3.7 (教材 例3.4) 求出通过平面上三点(0,7),(1,6)和(2,9)的二次多项式并画出其图形. 根据题设条件有 输入 Clear[x]; A={{0,0,1},{1,1,1},{4,2,1}} y={7,6,9} p=LinearSolve[A,y] Clear[a,b,c,r,s,t];{a,b,c}.{r,s,t} f[x_]=p.{x^2,x,1}; Plot[f[x],{x,0,2},GridLines?Automatic,PlotRange?All]; 则输出的值为 {2,?3,7} 并画出二次多项式的图形(略). 非齐次线性方程组的通解 用命令Solve求非齐次线性方程组的通解. 例3.8 求出通过平面上三点(0,0),(1,1),(-1,3)以及满足的4次多项式 解 设则有 输入 Clear[a,b,c,d,e];q[x_]=a*x^4+b*x^3+c*x^2+d*x+e; eqs=[q[0]= =0,q[1]= =1,q[-1]= =3,q’[-1]= =20,q’[1]= =9]; {A,y}=LinearEquationsToMatrices[eqs,{a,b,c,d}]; p=LinearSolve[A,y];f[x_]=p.{x^4,x^3,x^2,x,1}; Plot[f[x],{x,-1,1},GridLines-Automa