第三章MATLAB在高等数学中的应用.pptVIP

第3章 MATLAB在高等数学中的应用 ;3.1 矩阵分析 ;格式：n＝norm(A，p) 功能：norm函数可计算几种不同类型的矩阵范数，根据p的不同可得到不同的范数。 ;⑶ p=‘fro’时，计算矩阵A的Frobenius范数，即sqrt(sum(diag(A’*A)))。;⑶ norm(A，inf) 返回值max(abs(A))。;格式：X=det(A) 功能：计算方阵A的行列式值。 ; 写成矩阵形式可表示为：AX＝B 或 XA＝B。其中系数矩阵A的阶数为m×n。在MATLAB中，引入矩阵除法求解。; %创建线性方程组的系数矩阵和向量;x = 1.0000 -1.0000 2.0000 ;4．矩阵的分解 ;格式二：[L,U,P]=lu(A) 功能：返回上三角矩阵U，真正下三角矩阵L，及一个置换矩阵P(用来表示排列规则的矩阵)，满足L*U=P*A；如果P为单位矩阵，满足A=L*U。; a=[1:3;4:6;4,2,6] a = 1 2 3 4 5 6 4 2 6;U = 4.0000 5.0000 6.0000 0 -3.0000 0 0 0 1.5000;U = 4.0000 5.0000 6.0000 0 -3.0000 0 0 0 1.5000 P = 0 1 0 0 0 1 1 0 0;格式二：[Q,R,E]=qr(A) 功能：产生一个置换矩阵E，一个上三角矩阵R(其对角线元素降序排列)和一个归一化矩阵Q，满足A*E=Q*R；;5．奇异值分解 ;格式二：S=svd(A) 功能：返回奇异值组成的向量。;格式二：[V,D]=eig(A) 功能：返回矩阵V和D。其中对角阵D的对角元素为A的特征值，V的列向量是相应的特征向量，使得A*V=V*D。;7．矩阵的幂次运算: A^p;(4) 矩阵的元素幂、按矩阵元素的幂 利用运算符“A.^p”实现矩阵的元素幂或按矩阵元素的幂运算。 ;8．矩阵结构形式的提取与变换 ; x=fliplr(a) x = 3 2 1 6 5 4;格式二 ：X=reshape(A,m,n,p,...) 或 X=reshape(A,[m,n,p,...]) 功能：从A中形成多维阵列(m×n×p×...)。; x=reshape(a,3,2) x = 1 5 4 3 2 6;(5) 对角矩阵和矩阵的对角化函数diag( );格式二：X=diag(A) 功能：当A为n元向量时，等同于k=0时的X=diag(A,k)，即产生A的元素处于主对角线的对角方阵。当A为矩阵时，X=diag(A)相当于k=0。; x=diag(a,1) x = 0 1 0 0 0 0 2 0 0 0 0 3 0 0 0 0; c1=diag(b) c1 = 8 5 2;(6) 取矩阵的左下三角部分函数tril( );格式一：X=triu(A,k) 功能：得到矩阵A的第k条对角线及其以上的元素；当k=0时表示主对角线，k0表示主对角线之上，k0表示主对角线以下。 格式二：X=triu(A) 功能：得到矩阵A的右上三角阵。;3.2 多项式运算;3．多项式相乘运算 格式：w=conv(u,v) 功能：返回u和v两向量的卷积，也就是u和v代表的两多项式的乘积。;1．多项式求导函数polyder 格式一：k=polyder(p) 功能：返回多项式p的一阶导数。 格式二：k=polyder(u,v) 功能：返回多项式u与v乘积的导数。 格式三：[q,d]=polyer(u,v) 功能：返回多项式商u/v的导数, 返回的格式为：q为分子, d为分母。; w=[1 6 20 5] w = 1 6 20 5 polyder(w) ans = 3 12 20; [q,d]=polyder(a,b) q = 2 12 6 d