矩阵分解与线性方程组求解
更新时间:2024-04-03 17:27:01 阅读量: 综合文库 文档下载
一、 用列主元素高斯削去法求解下述线性方程组:
?x1?13x2?2x3?34x4?13?2x?6x?7x?10x??22?1234 ??10x?x?5x?9x?141234????3x1?5x2?15x4??36程序:
function x=gaussa(a)
m=size(a); n=m(1); x=zeros(n,1); for k=1:n-1
[c,i]=max(abs(a(k:n,k))); q=i+k-1; if q~=k
d=a(q,:);a(q,:)=a(k,:);a(k,:)=d end
for i=k+1:n
a(i,:)=a(i,:)-a(k,:)*a(i,k)/a(k,k) end end
for j=n:-1:1
x(j)=(a(j,n+1)-a(j,j+1:n)*x(j+1:n))/a(j,j) end
执行过程:
>> a=[1 13 -2 -34 13;2 6 -7 -10 -22;-10 -1 5 9 14; -3 -5 0 15 -36] a =
-10 -1 5 9 14 2 6 -7 -10 -22 1 13 -2 -34 13 -3 -5 0 15 -36 >> gaussa(a) a =
-10.0000 -1.0000 5.0000 9.0000 14.0000 0 5.8000 -6.0000 -8.2000 -19.2000 1.0000 13.0000 -2.0000 -34.0000 13.0000 -3.0000 -5.0000 0 15.0000 -36.0000 a =
-10.0000 -1.0000 5.0000 9.0000 14.0000 0 5.8000 -6.0000 -8.2000 -19.2000 0 12.9000 -1.5000 -33.1000 14.4000 -3.0000 -5.0000 0 15.0000 -36.0000 a =
-10.0000 -1.0000 5.0000 9.0000 14.0000 0 5.8000 -6.0000 -8.2000 -19.2000 0 12.9000 -1.5000 -33.1000 14.4000 0 -4.7000 -1.5000 12.3000 -40.2000
a =
-10.0000 -1.0000 5.0000 9.0000 14.0000 0 12.9000 -1.5000 -33.1000 14.4000 0 5.8000 -6.0000 -8.2000 -19.2000 0 -4.7000 -1.5000 12.3000 -40.2000 a =
-10.0000 -1.0000 5.0000 9.0000 14.0000 0 12.9000 -1.5000 -33.1000 14.4000 0 0.0000 -5.3256 6.6822 -25.6744 0 a =
-10.0000 0 0 0 a =
-10.0000 0 0 0 x =
0 0 0 10.7786 x =
0 0 18.3452 10.7786 x =
0 30.9062 18.3452 10.7786 x =
14.3827 30.9062 18.3452 10.7786 ans =
14.3827 30.9062 18.3452
-4.7000 -1.0000 12.9000 0.0000 0 -1.0000 12.9000 0.0000 -0.0000 -1.5000 12.3000 -40.2000 5.0000 9.0000 14.0000 -1.5000 -33.1000 14.4000 -5.3256 6.6822 -25.6744 -2.0465 0.2403 -34.9535 5.0000 9.0000 14.0000 -1.5000 -33.1000 14.4000 -5.3256 6.6822 -25.6744 0 -2.3275 -25.0873 10.7786
二、用矩阵A的杜丽特尔三角分解A=LU,求解方程组:
?157010??x1??????618159??x2??010287??x?=
3?????50635??x????4??8????6??4? ???2???程序:
function x=lua(a,b)
m=size(a);n=m(1);x=zeros(n,1);y=zeros(n,1); for k=1:n
a(k,k:n) = a(k,k:n)-a(k,1:k-1)*a(1:k-1,k:n);
a(k+1:n,k)=(a(k+1:n,k)-a(k+1:n,1:k-1)*a(1:k-1,k))/a(k,k); end
U = triu(a,0)
L =eye(n)+ tril(a,-1) for i=1:n
y(i) = b(i)-L(i,1:i-1)*y(1:i-1); end
for i=n:-1:1
x(i)=(y(i)-U(i,i+1:n)*x(i+1:n))/U(i,i); end 执行过程:
>> a=[15 7 0 10;6 18 15 9;0 10 28 7;5 0 6 35] a =
15 7 0 10 6 18 15 9 0 10 28 7 5 0 6 35 >> b=[8;6;4;2] b = 8 6 4 2 >> lua(a,b) U =
15.0000 7.0000 0 10.0000 0 15.2000 15.0000 5.0000 0 0 18.1316 3.7105 0 0 0 30.7351 L =
1.0000 0 0 0 0.4000 1.0000 0 0
0 0.6579 1.0000 0 0.3333 -0.1535 0.4579 1.0000 ans =
0.5264 0.0718 0.1272 -0.0399
三、用乔列斯基分解计算下述线性方程组:
00??x1??5??4?10??????0??x2??8???14?10?0?14?10? ?x?=?16? ???3???0?14?1??x4??24??0?????000?14????x5??36?程序:
function x=choleskey(a,b)
m=size(a);n=m(1);x=zeros(n,1);y=zeros(n,1); G=zeros(m); for j=1:n
G(j,j)=sqrt(a(j,j)-G(j,1:j-1)*G(j,1:j-1)')
G(j+1:n,j)=(a(j+1:n,j)-G(j+1:n,1:j-1)*G(j,1:j-1)')/G(j,j) end
for i=1:n
y(i)=(b(i)-G(i,1:i-1)*y(1:i-1))/G(i,i); end
for i=n:-1:1
x(i)=(y(i)-G(i+1:n,i)'*x(i+1:n))/G(i,i); end 运行过程:
>> a=[4 -1 0 0 0;-1 4 -1 0 0;0 -1 4 -1 0;0 0 -1 4 -1;0 0 0 -1 4] a =
4 -1 0 0 0 -1 4 -1 0 0 0 -1 4 -1 0 0 0 -1 4 -1 0 0 0 -1 4 >> b=[5 ;8;16;24;36] b = 5 8 16 24 36
>> choleskey(a,b) G =
2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 G =
2.0000 0 0 0 0 -0.5000 0 0 0 G =
2.0000 -0.5000 0 0 0 G =
2.0000 -0.5000 0 0 0 G =
2.0000 -0.5000 0 0 0 G =
2.0000 -0.5000 0 0 0 G =
2.0000 -0.5000 0 0 0 G =
0 0 0 0 0 1.9365 0 0 0 0 1.9365 -0.5164 0 0 0 1.9365 -0.5164 0 0 0 1.9365 -0.5164 0 0 0 1.9365 -0.5164 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1.9322 0 0 0 0 1.9322 -0.5175 0 0 0 1.9322 -0.5175 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1.9319 2.0000 0 0 0 0 -0.5000 1.9365 0 0 0 0 -0.5164 1.9322 0 0 0 0 -0.5175 1.9319 0 0 0 0 -0.5176 0 G =
2.0000 0 0 0 0 -0.5000 1.9365 0 0 0 0 -0.5164 1.9322 0 0 0 0 G =
2.0000 -0.5000 0 0 0 ans =
2.3910 4.5641 7.8654 10.8974 11.7244
0 0 0 1.9365 -0.5164 0 0 -0.5175 0 0 0 1.9322 -0.5175 0 1.9319 0 -0.5176 1.9319 0 0 0 0 0 0 1.9319 0 -0.5176 1.9319
2.0000 0 0 0 0 -0.5000 1.9365 0 0 0 0 -0.5164 1.9322 0 0 0 0 -0.5175 1.9319 0 0 0 0 -0.5176 0 G =
2.0000 0 0 0 0 -0.5000 1.9365 0 0 0 0 -0.5164 1.9322 0 0 0 0 G =
2.0000 -0.5000 0 0 0 ans =
2.3910 4.5641 7.8654 10.8974 11.7244
0 0 0 1.9365 -0.5164 0 0 -0.5175 0 0 0 1.9322 -0.5175 0 1.9319 0 -0.5176 1.9319 0 0 0 0 0 0 1.9319 0 -0.5176 1.9319
正在阅读:
矩阵分解与线性方程组求解04-03
可追溯的产品内部批号编制规则和实施细则06-29
专题17 论述类文本阅读(教学案)-2017年高考语文二轮复习精品资料(解析版)Word版含解析高品质版11-22
试论公安指挥中心职能03-15
市人民法院党组书记院长述职述廉报告09-07
区住房和城乡建设局2021年工作总结及下一年工作计划08-02
大学生增强国家安全忧患意识的重要性论文05-22
上海市普陀区2013学年度高三年级第一次质量调研语文试卷06-06
自动化数电试卷A及答案03-13
- 多层物业服务方案
- (审判实务)习惯法与少数民族地区民间纠纷解决问题(孙 潋)
- 人教版新课标六年级下册语文全册教案
- 词语打卡
- photoshop实习报告
- 钢结构设计原理综合测试2
- 2014年期末练习题
- 高中数学中的逆向思维解题方法探讨
- 名师原创 全国通用2014-2015学年高二寒假作业 政治(一)Word版
- 北航《建筑结构检测鉴定与加固》在线作业三
- XX县卫生监督所工程建设项目可行性研究报告
- 小学四年级观察作文经典评语
- 浅谈110KV变电站电气一次设计-程泉焱(1)
- 安全员考试题库
- 国家电网公司变电运维管理规定(试行)
- 义务教育课程标准稿征求意见提纲
- 教学秘书面试技巧
- 钢结构工程施工组织设计
- 水利工程概论论文
- 09届九年级数学第四次模拟试卷
- 方程组
- 求解
- 矩阵
- 线性
- 分解
- 辽东山区软枣猕猴桃商品化栽培
- 政审思想品德鉴定表
- 省、市、县官方网站新闻投稿电子邮箱
- 水温自动控制系统实验报告课案
- 最实用的医学影像学题库
- 题库:终身学习与职业发展
- 项目管理目标责任书(最终版)
- 南昌县关于土地上房屋征收及补偿办法
- 诗歌朗诵比赛主持串词
- 管线保护专项施工方案
- 人文-2014医师定期考核题库40
- 生物制剂临床使用管理办法
- (精美排版)北京市文化创意产业提升规划(2014-2020年)(全文
- 工艺管道单线图英文翻译
- 小学生信息素养培养探微
- 中海成本控制的四大重点
- 2015-2016学年度人教版五年级数学下册教案(全册)
- 村卫生室各项规章制度册
- 安徽理工大学2014-2015学年度共青团工作先进集体和先进个人名单
- 合肥市司法局办公室 第64期 2008年12月9日