数学实验 正文.pptVIP

第一篇 软件实验 实验1 矩阵及其基本运算 一、实验目的 熟悉Matlab软件环境 学会用Matlab软件进行矩阵的基本运算。 二、实验内容 （2）常用运算指令 1、常数矩阵输入 任何矩阵（向量），我们可以直接按行方式输入每个元素： 1)同一行中的元素用逗号（，）或者用空格符来分隔，且空格个数不 限； 2)不同的行用分号（；）分隔； 3)所有元素处于一方括号（[ ]）内； 4)当矩阵是多维（三维以上），且方括号内的元素是维数较低的矩阵 时，会有多重的方括号。 例如： Time = [11 12 1 2 3 4] ； X_Data = [2.32 3.43；4.37] ； Null_M = [ ]； 2、符号矩阵的生成 在MATLAB中输入符号向量或者矩阵的方法和输入数值类型的向量或者矩阵在形式上很相像，只不过要用到符号矩阵定义函数sym，或者是用到符号定义函数syms，先定义一些必要的符号变量，再像定义普通矩阵一样输入符号矩阵。 MATLAB还提供了一个将数值型转化成符号型的命令，即sym。 Digit_M = [1/3 sqrt（2）；exp（0.23） log（29）] Syms_M = sym（Digit_M） 注意：矩阵是用分数形式还是浮点形式表示的，将矩阵转化成符号矩阵后，都将以最接近原值的有理数形式表示或者是函数形式表示。 3、特殊矩阵的生成 1)全零阵函数 zeros 格式 B = zeros(n) %生成n×n全零阵 B = zeros(m,n) %生成m×n全零阵 2)单位阵函数 eye 格式 Y = eye(n) %生成n×n单位阵 Y = eye(m,n) %生成m×n单位阵 3)全1阵函数 ones 格式 Y = ones(n) %生成n×n全1阵 Y = ones(m,n) %生成m×n全1阵 4) (0,1)内均匀分布随机矩阵函数 rand 格式 Y = rand(n) %生成n×n随机矩阵 Y = rand(m,n) %生成m×n随机矩阵 5)对角矩阵函数 blkdiag 格式 out = blkdiag(a,b,c,d,…) %产生以a,b,c,d,…为对角线元素的矩阵 6)Magic(魔方)矩阵函数 magic 格式 M = magic(n) %产生n 阶魔方矩阵 4、矩阵运算 1)加、减、乘法运算符：＋ 、－ 、* 运算规则： 对应元素相加、减，乘即按线性代数中矩阵的 “十” 、“一” 、 “ * ”运算进行。 2)数乘 例：a=2*X 3)点乘函数 dot 格式 C = dot(A,B) %若A、B为向量，则返回向量A与B的 点积，A与B长度相同；若为矩阵，则A 与B有相同的维数。 C = dot(A,B,dim) %在dim维数中给出A与B的点积 4)向量叉乘函数 cross 格式 C = cross(A,B) %若A、B为向量，则返回A与B的叉乘。 5)除法运算 左除（\）和右除（/） 一般情况下，x=a\b是方程a*x =b的解，而x=b/a是方程x*a=b的解。 6)矩阵乘方运算符：^ 格式 A^P 注意：A必须是方阵。 7)矩阵转置运算符：′ 运算规则： 若A为复数矩阵，则A转置后的元素由A对应元素的共轭复数构成。 若仅希望转置，则用如下命令：A.′。 8)方阵的行列式函数 det 格式 d = det(X) %返回方阵X的多项式的值 9)逆函数 inv 格式 Y=inv(X) %求方阵X的逆矩阵。若X为奇异阵或近似奇异阵，将给出警告信息。 5、特征值与二次型 特征值与特征向量的求法 设A为n阶方阵，如果数K和n维列向量x使得关系式Ax=Kx成立，则称K为方阵A的特征值，非零向量x称为A对应于特征值K的特征向量。 函数 eig 格式 d = eig(A) （%求矩阵A的特征值d，以向量形式存放d。 ） d = eig(A,B) （ %A、B为方阵，求广义特征值d，以向量形式 存放d。 ） 6、秩与线性相关性 矩阵A的秩是矩阵A中最高阶非零子式的阶数；向 量组的秩通常由该向量组构成的矩阵来计算。 秩函数 rank 格式 k = rank(A