第五讲线性代数中的数值计算问题.pptVIP

第1页，共32页，星期日，2025年，2月5日【引例】求下列三阶线性代数方程组的近似解MATLAB程序为：A=[2-54;15-2;-124]b=[5;6;5]x=A\bx=

2.7674

1.1860

1.3488第2页，共32页，星期日，2025年，2月5日在MATLAB命令窗口，先输入下列命令构造系数矩阵A和右端向量b：A=[2-54;15-2;-124]A=2-5415-2-124b=[5;6;5]b=565然后只需输入命令x=A\b即可求得解x：x=A\bx=2.76741.18601.3488第3页，共32页，星期日，2025年，2月5日一、特殊矩阵的实现第4页，共32页，星期日，2025年，2月5日1.零矩阵:所有元素值为零的矩阵称为零矩阵。零矩阵可以用zeros函数实现。zeros是MATLAB内部函数，使用格式如下：zeros(m)：产生m阶零矩阵；zeros(m,n)：产生m*n阶零矩阵，当m=n时同上；zeros(size(A))：产生与矩阵A同样大小的零矩阵。一、特殊矩阵的实现常见的特殊矩阵有零矩阵、幺矩阵、单位矩阵、三角形矩阵等，这类特殊矩阵在线性代数中具有通用性；还有一类特殊矩阵在专门学科中有用，如有名的希尔伯特(Hilbert)矩阵、范德蒙(Vandermonde)矩阵等。第5页，共32页，星期日，2025年，2月5日2.幺矩阵:所有元素值为1的矩阵称为幺矩阵。幺矩阵可以用ones函数实现。它的调用格式与zeros函数一样。【例1】试用ones分别建立3*2阶幺矩阵、和与前例矩阵A同样大小的幺矩阵。用ones(3,2)建立一个3*2阶幺阵：ones(3,2)%一个3*2阶幺阵ans=111111一、特殊矩阵的实现第6页，共32页，星期日，2025年，2月5日3.单位矩阵:主对角线的元素值为1、其余元素值为0的矩阵称为单位矩阵。它可以用MATLAB内部函数eye建立，使用格式与zeros相同。4.数量矩阵:主对角线的元素值为一常数d、其余元素值为0的矩阵称为数量矩阵。显然，当d=1时，即为单位矩阵，故数量矩阵可以用eye(m)*d或eye(m,n)*d建立。一、特殊矩阵的实现第7页，共32页，星期日，2025年，2月5日一、特殊矩阵的实现5.对角阵:对角线的元素值为常数、其余元素值为0的矩阵称为对角阵。我们可以通过MATLAB内部函数diag，利用一个向量构成对角阵；或从矩阵中提取某对角线构成一个向量。使用格式为diag(V)和diag(V,k)两种。设V为具有m个元素的向量，diag(V)将产生一个m阶对角阵，其主对角线的元素值即为向量的元素值；diag(V,k)将产生一个n(n=m+|k|，k为一整数)阶对角阵，其第k条对角线的元素值即为向量的元素值。注意：当k＞0，则该对角线位于主对角线的上方第k条；当k＜0，该对角线位于主对角线的下方第|k|条；当k=0，则等同于diag(V)。用diag建立的对角阵必是方阵。第8页，共32页，星期日，2025年，2月5日一、特殊矩阵的实现【例2】已知向量v，试建立以向量v作为主对角线的对角阵A；建立分别以向量v作为主对角线两侧的对角线的对角阵B和C。MATLAB程序如下：v=[1;2;3];%建立一个已知的向量vA=diag(v)A=100020003B=diag(v,1)B=0100002000030000C=diag(v,-1)C=0000100002000030第9页，共32页，星期日，2025年，2月5日6.从矩阵中提取某对角线我们也可以用diag从矩阵中提取某对角线构成一个向量。设A为m*n阶矩阵，diag(A)将从矩阵A中提取其主对角线产生一个具有min(m,n)个元素的向量。diag(A,k)的功能是：当k