应用高等数学第1章107.pptVIP

* 第一节 函数及其图形 10.7 数学实验 一、矩阵及行列式的运算 二、解方程组 三、线性规划求解 一、矩阵及行列式的计算 1。矩阵的生成 为了得到矩阵 可输入： A=[1 2 3；4 5 6；7 8 9] 在MATLAB中输入矩阵，以“[”开始，以“]”结束，元素间用“,”或空格分隔，换行用“;”或回车。 MATLAB还可以用一些看起来简单却十分有用的函数生成矩阵，称为工具函数，如表10-4 表10-4 将列矩阵M生成对角矩阵 diag(M，n) 生成m行n列随机矩阵 rand(m，n) 以M对角元素生成列矩阵 diag(M) 生成n阶单位矩阵 eye(n) 提取x的下三角矩阵 tril(x，n) 生成m行n列全1矩阵 ones(m，n) 提取x的上三角矩阵 triu(x，n) 生成m行n列零矩阵 zeros(m，n) 生成矩阵 函数 生成矩阵 函数 例如要生成8行5列的零矩阵，则可输入 zeros(8,5) 要生成5行5列的单位矩阵，则可输入 eye(5) 2。矩阵的运算 设k为任意实数，A，B为满足矩阵运算条件的矩阵，在MATLAB中，规定了矩阵的如下运算： A±B 矩阵的加减（A，B为同型矩阵） k*A 数乘矩阵 A*B 矩阵的乘法（A的列数等于B的行数） A.*B 两矩阵对应元素相乘（A，B为同型矩阵） A./B 两矩阵对应元素相除（A，B为同型矩阵） A.’ 矩阵的转置 A’ 矩阵的共轭转置 inv(A)或A^(-1) 矩阵A的逆 A^k 矩阵A的k次幂 A.^k A中每个元素k次幂 rank(A) 计算矩阵A的秩 A=[1 2 3;4 5 6;7 8 9]; A^2 ans = 30 36 42 66 81 96 102 126 150 A.^2 ans = 1 4 9 16 25 36 49 64 81 例如 3。行列式的计算 用det(A) 计算矩阵A对应行列式的值，例如 A=[5 2 3；6 5 4；7 8 9]； det(A) ans = 52 二、 解线性方程组 用X=A\B(左除)得到非齐次线性方程组AX=B（B≠0）的一个特解，用函数null(A，’r’)得到线性方程组AX=0的有理基础解系． 练习 解线性方程组 输入 format rat A=[2 -4 5 3 ；3 -6 4 2 ；4 -8 17 11 ]； B=[7 7 21]； AB=[A B]； RA=rank(A) RA = 2 RAB=rank(AB) RAB = 2 因为RA=RAB4，即系数矩阵的秩等于增广矩阵的秩小于未知数的个数，所以方程组有无穷多解，再输入 X0=A\B X0 = 0 -1/2 1 0 Y=null(A，r) Y = 2 2/7 1 0 0 -5/7 0 1 syms k1 k2 %申明k1，k2为符号量 X=X0+k1*Y(:，1)+k2*Y(:，2) X = [ 2*k1+2/7*k2] [ -1/2+k1] [ 1-5/7*k2] [ k2] 即为方程组的全部解． 三、线性规划求解 MATLAB提供了linprog()函数求线性规划问题中使目标函数最小的优化方案，其一般格式是： [X fval]=linprog(f，A，B，Aeq，Beq，LB，UB) 这里f是由目标函数系数构成的向量．A，B分别是约束条件中不等式组的系数矩阵和常数矩阵．Aeq，Beq分别是约束条件中方程组的系数矩阵和常数矩阵．LB、UB分别是决策变量的下界和上界．返回值X是目标函数取得最小值时决策变量的一组取值，fval是优化结束后得到的目标函数值． 注意，这个函数是计算满足目标函数取得最小值的一组变量的值，如果要求的是使目