第二讲 矩阵的数值运算.pptVIP

第二讲 MATLAB的数值计算 —— matlab 具有出色的数值计算能力，占据世界上数值计算软件的主导地位 数值运算的功能 创建矩阵 矩阵运算 多项式运算 线性方程组 数值统计 线性插值 函数优化 微分方程的数值解 一、命令行的基本操作 1. 创建矩阵的方法 其中符号的作用: 逗号和分号的作用 逗号和分号可作为指令间的分隔符，matlab允许多条语句在同一行出现。 分号如果出现在指令后，屏幕上将不显示结果。 用matlab函数创建矩阵 空阵 [ ] — matlab允许输入空阵，当一项操作无结果时，返回空阵。 rand(1,N)—— 产生[0 1]上均匀分布的随机序列，长度为N randn(1,N)—— 产生均值为0，方差为 1的高斯随机序列，即白噪声序列，长度N eye —— 单位矩阵 zeros ——全部元素都为0的矩阵 ones ——全部元素都为1的矩阵 2. 矩阵的修改 ? 直接修改 可用?键找到所要修改的矩阵，用?键移动到要修改的矩阵元素上即可修改。 ? 指令修改 可以用A(?,?)= ? 来修改。 例： a=[1 2 0;3 0 5;7 8 9] a =1 2 0 3 0 5 7 8 9 a(3,3)=0 a =1 2 0 3 0 5 7 8 0 二、矩阵运算 1. 矩阵加、减（＋,－）运算 2. 矩阵乘（?）运算 3. 在matlab中有两种矩阵除运算，矩阵左除和矩阵右除，\和/ 。如果A矩阵是非奇异方阵，则A\B和B/A运算可以实现。A\B等效于A的逆左乘B矩阵，也就是inv(A)*B，而B/A等效于A矩阵的逆右乘B矩阵，也就是B*inv(A)。 4. 矩阵乘方—— a^n，a^p，p^a a ^ p —— a 自乘p次幂 例： a=[1,2,3;4,5,6;7,8,9];a^2 ans =30 36 42 66 81 96 102 126 150 例： a=[1,2,3;4,5,6;7,8,9];a^0.5 ans = 0.4498 + 0.7623i 0.5526 + 0.2068i 0.6555 -0.3487i 1.0185 + 0.0842i 1.2515 + 0.0228i 1.4844 - 0.0385i 1.5873 - 0.5940i 1.9503 - 0.1611i 2.3134 + 0.2717i 5. 矩阵的其它运算 inv(A) —— 矩阵求逆 det(A) —— 矩阵A所对应的行列式的值 eig (A)—— 矩阵的特征值及特征向量 diag (A) —— 提取A对角元素或生成对角矩阵 sqrtm(A) ——矩阵A的平方根 logm(A) ——矩阵A的自然对数 expm（A） —— A的矩阵指数函数 trace(A) —— 求矩阵的迹 rank（A） ——求矩阵的秩 6. 矩阵的一些特殊操作 矩阵的变维 例： a=[1 2 3 4;5 6 7 8;9 10 11 12] a = 1 2 3 4 5 6 7 8 9 10 11 12 b=reshape(a,4,3) % 把指定的矩阵改变形状，但是元素个数不变 b = 1 6 11 5 10 4 9 3 8 2 7 12 关系运算符的运算法则为： 当两个比较量是标量时，直接比较两数的大小若关系成立，关系表达式结果为1，否则为0。 当参与比较的量是两个维数相同的矩阵时，比较是对两矩阵相同位置的元素按标量关系运算规则逐个进行，并给出元素比较结果。最终的关系运算的结果是一个维数与原矩阵相同的矩阵，它的元素由0或1组成。 9. 数组运算 数组运算指元素对元素的算术运算，与通常 意义上的由符号表示的线性代数矩阵运算不同。 数组加减(.+,.-) a=[1 2 3;4 5 6;7