Mathematica 线性代数 数学软件与 及数学实验 .ppt

;1、Join[list1,list2] 把表list1，list2连成一个表;写出分段函数 的Mathematica自定义函数形式,并画出其在[3,3]上的图形 解: 因此Mathematica 命令为: f[x_]:=If[x1, x+Sin[x], x*Cos[x]] Plot[f[x],{x,-3, 3}] ;Do[expr,{n}] 重复计算expr n次;; 矩阵的裁剪与拼接 从一个矩阵中取出若干行（列）构成新矩阵称为裁 剪，将几个矩阵接在一起称为拼接。 例如： ; （2）B[[ 4 ]] 表示截取B的位于第4行的元素, 即{ 1，-2，0，2}。 （3）Transpose [B][[4]] 表示截取B的第4列的元素， 即{ 0，5，3，2 }; （4）B[[ {1,2},{3,4}]] 表示截取B的位于第1行和 第2行、第3列和第4列的元素， 即; ? 特殊函数生成的矩阵（特殊矩阵的构造） 1． 零矩阵和全同矩阵 输入语句： Table[0,{m},{n}}] MatrixForm[%] ;A=Table[Random[Integer,{0,9}],{I,5},{j,5}] MatrixForm[%];5. 三对角矩阵 ;B=Array[Prime,8] DiagonalMatrix[B]//MatrixForm;特殊矩阵例-2;特殊矩阵例-3;特殊矩阵例-4; （1）?矩阵的输入格式： A={{a,b,c},{d,e,f},{g,h,i}}, （2）???? 求A的转置：Transpose[A] （3）???? 求A加B：A+B （4）???? 求A减B：A - B （5）???? 求数k乘以A：k*A （6）???? 求A乘以B： A.B (注：区分A B) （7）???? 求A的行列式：Det[A] ;（8）???? 求A的逆：Inverse[A] （9）???? 求A的特征值组：Eigenvaluse[A] （10）???? 求A的数字特征值 Eigenvaluse[N[A]] （11）???? 求A的特征向量组: Eigenvectors[A];已知三阶方阵 A1 = {{1, 2, 3}, {2, 3, 1}, {3, 1, 2}}; Det[A1] Transpose[A1] E1 = Inverse[A1] 注：容易验证矩阵A1同它的逆矩阵E1之间有 A1.E1 = {{1，0，0}，{0，1，0}，{0，0，1}};Clear[A] A={{1,2,3},{2,2,1},{3,4,3}} MatrixForm[%] D3=IdentityMatrix[3] MatrixForm[%] AE=Transpose[Join[Transpose[A],D3]]; MatrixForm[%] EANi=RowReduce[AE] MatrixForm[%] EANi[[{1,2,3},{4,5,6}]] MatrixForm[%];是否线性无关？;如何求矩阵的极大线性无关向量组呢？ ;三、线性方程组的求解;语 句;矩阵形式： A X= 0;例2;例3;四、Mathematica 在线性规划中的应用;例 某工厂有甲、乙、丙、丁四个车间，生产A、B、C、D、E、F六种产品，根据车床性能和过去的生产情况，得知生产单位产品所需车间的工作小时数，每个车间每月工作小时的上限，以及产品的价格如下表所示： 问各种产品每月应生产多少？才能使这个工厂总值达到最大？;解： 问题的数学模型为： max z =0.4 x1 +0.28 x2 +0.32 x3 +0.72 x4 +0.64 x5 + 0.6x6 s.t. 0.02x1 + 0.05 x4 700 0.02x2 + 0.05 x5 100 0.03 x3