五对角矩阵LU分解求解Ax=B =* function [x]=FivediagLU(A,a,b,c,C,B)%A,a,b,c,C依序为5从左到右5个对角线元素向量(其中不足n在对应元素补上0)，B为方程组右端向量 n=length(b);%初始化：L,l,u,r,R对应A,a,b,c,C；L=zeros(n,1);l=zeros(n,1);u=zeros(n,1);r=zeros(n,1);...
1)因位我们是针对五对角矩阵进行求解,由于这样的矩阵是稀疏的,观察题目1,,中给出的例子也可以看出其分解后的矩阵L和U都是稀疏矩阵,因此参考教材第 ,,21页的lin_band2.m程序,按如下方法对五对角矩阵A做LU分解。l121,,首先A的LU分解有如下形式:
第五章 特征值和特征向量、矩阵的对角化 扩展例题及求解
解析 解:矩阵的特征方程为 。 ………2分 解得。 ………4分 对由,解得基础解系。 ………6分 对由,解得基础解系………8分 由于正好正交,所以两两正交。 再将单位化,得。 故所求得正交矩阵且。………10分反馈 收藏
4 ?? 4??5function [x]=FivediagLU(A,a,b,c,C,B) %A,a,b,c,C依序为 5 从左到右 5 个对角线元素向量 (其中不足 n 在对应元素补上 0),B 为方程组右端向量n=length(b); %初始化: L,l,u,r,R对应 A,a,b,c,C;L=zeros(n,1);l=zeros(n,1);u=zeros(n,1);r=zeros(n,1);...
一些论著中也将这两种方法统称为多水平方法三、算法设计及分析1)因位我们是针对五对角矩阵进行求解,由于这样的矩阵是稀疏的,观察题目中给出的例子也可以看出其分解后的矩阵L和U都是稀疏矩阵,因此参考教材第21页的lin_band2.m程序,按如下方法对五对角矩阵A做LU分解。首先A的LU分解有如下形式: 11......
第五章 特征值和特征向量、矩阵的对角化 扩展例题及求解