matlab——6(1)数值计算.pptVIP

第6章 数值计算(1) 高等光学仿真 指导教师：张卫华 content 矩阵分析 多项式计算 数据统计处理 数据插值 离散傅立叶变换 线性方程组求解 非线性方程求解 常微分方程数值解法 函数极值 数值积分 数值微分 例： 先建立5×5矩阵A，然后将A的第一行元素乘以1，第二行乘以2，…，第五行乘以5。 程序： A=[17,0,1,0,15;23,5,7,14,16;4,0,13,0,22;10,12,19,21,3 ;11,18,25,2,19]; D=diag(1:5); U=D*A %用D左乘A，对A的每行乘以一个指定常数 矩阵的转置 (’) 矩阵的旋转 rot90(A,k)将矩阵A逆时针旋转90o的k倍。 矩阵的左右翻转 fliplr(A)。 矩阵的上下翻转 flipud(A)。 1．矩阵的逆 inv(A)。 例2 用求逆矩阵的方法解线性方程组。Ax=b其解为：x=A-1b 方阵的行列式 det(A)。 在MATLAB中，计算矩阵A的特征值和特征向量的函数是eig(A)，常用的调用格式有3种： E=eig(A)：求矩阵A的全部特征值，构成向量E。 [V,D]=eig(A)：求矩阵A的全部特征值，构成对角阵D，并求A的特征向量构成V的列向量。 [V,D]=eig(A,’nobalance’)：与第2种格式类似，但第2种格式中