计算机仿真第一章.docVIP

1.5 MATLAB的数值运算 1.5.1 矩阵运算 1. MATLAB环境下的矩阵 （1）矩阵赋值 其中一个常用的函数： [ n , m ] = size ( A ) 测取一个矩阵的维数。返回的两个参数n和m分别为A矩阵的行数、列数。 （2）矩阵的元素 MATLAB的矩阵元素可用任何表达式描述，也可以是复数。 B=[-1/3 1.3 sqrt(3) (1+2+3)*j ] （3）特殊矩阵 特殊矩阵可利用函数来建立。 单位矩阵函数 eye( ) 零矩阵函数 zeros( ) 1矩阵函数 ones( ) 随机元素矩阵函数 rand( ) 对角矩阵函数 diag( ) 伴随矩阵函数 compan( ) 上三角矩阵函数 triu( ) 下三角矩阵函数 tril( ) 2. 矩阵的基本运算 （1）矩阵的转置 B=A’ （2）加减运算 运算符“+” ，“—” （3）矩阵的乘法 两矩阵按乘法运算规则运算，运算符“*” 。 （4）矩阵的除法 矩阵的除法有两种运算符：左除“\” , 右除“/” 。 通常。 是的解，即A\B=inv(A)*B。同时B/A=B*inv(A) 例方程： ， X=A\B的解；解亦为 同样X=B/A是X*A=B的解，即B/A=B*inv(A) A B 执行X=B/A X= B*inv(A) 结果相同 X=[ 1 3 4 ]. （5）矩阵的乘方 运算符“^” B=A^3 （6）矩阵的翻转 fliplr(A) %左右翻转 flipup(A) %上下翻转 rot90(A) %旋转90度 （7）矩阵的超越函数 对矩阵的单个元素操作。 3. 矩阵的特殊运算 （1）矩阵行列式 格式： det(A) （2）矩阵求逆 inv(A) （3）矩阵的迹 定义： 指令格式： trace(A) （4）矩阵的秩 k=rank(A) （5）矩阵的三角分解 函数调用格式： [L , U]=lu(A) （6）矩阵的奇异值分解 函数调用格式： [U , S , V]=svd(A) 其中满足，A=USVT （7）矩阵的范数 N=norm(A , 选项) 例： 求1-范数 N=norm(A , 1) 求无穷范数 N=norm(A , inf ) （8）矩阵的特征值与特征向量 特征值与特征向量：对于一个矩阵A，如果存在一个非零向量x，且有一个标量满足 则称为矩阵A的一个特征值，而x为对应于特征值的特征向量。 矩阵的特征值与特征向量可由函数eig( )求出。 调用格式： [V , D]=eig(A) D为对角阵，对角线上的元素为矩阵A的特征值。 （9）矩阵的特征多项式、特征方程、特征根 矩阵A的特征多项式： 特征方程的根，为该矩阵的特征根。 求取特征多项式系数的函数为： p=poly(A) p为一个行向量，即。 求特征根的函数： V=roots(p) 1.5.2 向量运算 MATLAB向量运算符由矩阵运算符前面加一点“.”来表示。 点运算：对应元素的运算。 1. 向量的加减 2. 向量的乘法 3. 向量的除法 4. 向量的乘方 1.5.3 关系和逻辑运算 1. 关系运算 关系运算符见表1-8 表1-8关系运算符 关系操作符 意义 = = = = ~ = 小于 小于等于 大于 大于等于 等于 不等于 2. 逻辑运算 | ~ 与 或 非 1.5.4 多项式运算 1. 多项式的表示 数学中表示： MATLAB中表示： 注意：缺项补零。 多项式的系数向量 —— 多项式形式，函数为： f=poly2str( p , ’s’ ) p—系数向量，s—变量名。 2. 多项式的四则运算 （1）加减 阶次不同，首零填补。 （2）乘法 两个多项式p1，p2相乘的函数为： p=conv( p1 , p2 ) （3）除法 多项式除法的函数为： [p , r]=deconv( p1, p2 ) p—商多项式，r —余多项式。 3. 多项式的值及其导数 多项式的求值函数： f0=polyval( p , x0 ) x0—求值点 多项式的求导函数： dp=polyder( p ) 4. 多项式的求解 函数为： r=roots(p) p为多项式的系数向量，r为多项式的解。 5. 多项式的拟合 调用函数： [p , s]=polyfit(x , y , n) p_待拟合多项式； x , y_进行拟合的数据；