Gauss列主元消去法、QR(MATLAB)
更新时间:2023-03-20 23:53:01 阅读量: 实用文档 文档下载
例:用Gauss 列主元消去法、QR 方法求解如下方程组:
12342212141
312.4201123230x x x x ?????? ? ? ?- ? ? ?= ? ? ?-- ? ? ???????
1. 1)Gauss 列主元法源程序:
function x=Gauss(A,b)
[m,n]=size(A);
if m~=n
error('矩阵不是方阵')
return
end
B=[A,b];
n=length(A);
for j=1:n-1
q=[zeros(j-1,1);B(j:n,j)];
[c,r]=max(abs(q)); %c 为列主元,r 为所在行
if r~=j
temp=B(j,:); %交换两行
B(j,:)=B(r,:);
B(r,:)=temp;
end
for i=j+1:n
B(i,:)=B(i,:)-B(j,:)*(B(i,j)/c);
end
end
x(n)=B(n,n+1)/B(n,n);
for i=n-1:-1:1
for j=i:n-1
B(i,n+1)=B(i,n+1)-B(i,j+1)*x(j+1);
end
x(i)=B(i,n+1)/B(i,i);
end
2)在命令窗口输入A ,b,得到x 的近似解:
>> A=[2,2,1,2;4,1,3,-1;-4,-2,0,1;2,3,2,3];
>> b=[1;2;1;0];
>> x=Gauss(A,b)
x =
1.5417
-2.7500
0.0833
1.6667
2. 1)QR方法源程序:
function [Q,R,X]=qrfj(A,b)
[m,n]=size(A);
if m error('A不符合规则') return end R=A; Q=eye(n); for i=1:n-1 a=R(i:n,i); e=[1;zeros(n-i,1)]; w=a-norm(a)*e; Hw=eye(n-i+1)-(2/(w'*w))*(w*w'); H=blkdiag(eye(i-1),Hw); %块对角矩阵 Q=Q*H; R=H*R; end Y=Q'*b; X(n)=Y(n)/R(n,n); for i=n-1:-1:1 for j=i:n-1 Y(i)=Y(i)-R(i,j+1)*X(j+1); end X(i)=Y(i)/R(i,i); end 2)在命令窗口输入A,b,得到x的近似解: >> A=[2,2,1,2;4,1,3,-1;-4,-2,0,1;2,3,2,3]; >> b=[1;2;1;0]; >> [Q,R,X]=qrfj(A,b) Q = 0.3162 0.3705 -0.0284 -0.8729 0.6325 -0.4940 0.5966 -0.0000 -0.6325 0.0823 0.7386 -0.2182 0.3162 0.7822 0.3125 0.4364 R = 6.3246 3.4785 2.8460 0.3162 0.0000 2.4290 0.4529 3.6641 0.0000 0.0000 2.3864 1.0227 -0.0000 -0.0000 -0.0000 -0.6547 X = 1.5417 - 2.7500 0.0833 1.6667 (资料素材和资料部分来自网络,供参考。可复制、编制,期待你的好评与关注)
正在阅读:
二年级语文数学易错题03-23
会计基础试题(二)(1)01-04
2015国家公务员考试时事热点:“基层看作风”之二 - 网埭港村治脏记09-12
ArcSWAT使用手册整理06-23
2011年注册城市规划师考试原理真题03-13
天然药物活性筛选考卷01-27
北京语言大学澳大利亚留学预科03-09
PPT复习题02-27
2016高考生物新课标3卷解析(二)11-07
- 教学能力大赛决赛获奖-教学实施报告-(完整图文版)
- 互联网+数据中心行业分析报告
- 2017上海杨浦区高三一模数学试题及答案
- 招商部差旅接待管理制度(4-25)
- 学生游玩安全注意事项
- 学生信息管理系统(文档模板供参考)
- 叉车门架有限元分析及系统设计
- 2014帮助残疾人志愿者服务情况记录
- 叶绿体中色素的提取和分离实验
- 中国食物成分表2020年最新权威完整改进版
- 推动国土资源领域生态文明建设
- 给水管道冲洗和消毒记录
- 计算机软件专业自我评价
- 高中数学必修1-5知识点归纳
- 2018-2022年中国第五代移动通信技术(5G)产业深度分析及发展前景研究报告发展趋势(目录)
- 生产车间巡查制度
- 2018版中国光热发电行业深度研究报告目录
- (通用)2019年中考数学总复习 第一章 第四节 数的开方与二次根式课件
- 2017_2018学年高中语文第二单元第4课说数课件粤教版
- 上市新药Lumateperone(卢美哌隆)合成检索总结报告
- 消去法
- MATLAB
- Gauss
- QR
- 人教版小学数学三年级上册《8分数的初步认识:认识几分之一》赛课导学案_10
- 八年级物理上册全册全套试卷测试卷 (word版,含解析)
- 郭硕鸿《电动力学》课后答案
- 超声辅助提取杨梅渣膳食纤维工艺研究
- 四年级四则运算练习题带中括号
- 汽车内外塑料件项目可行性研究报告评审方案设计(2013年发改委标准案例范文)
- 对高中美术鉴赏的再认识与课堂气氛的调动
- 中班绘本说课稿范文
- 人教版数学七年级下册8.1《二元一次方程组(1)》 教案
- 小学新教师培训心得体会1000字.doc
- 公司日常备用金管理规定
- 数学知识点广东省广州市番禺区六校教育教学联合体_七年级数学10月统考试题(A卷)新人教版【含解析】
- 基于产业链整合的广告公司并购研究
- 最新公司共管账户暂行办法资料
- 江西省金溪县第一中学下册圆周运动达标检测卷(Word版 含解析)
- 新北京版二下数学:《探索规律》教案
- 水泥搅拌桩首件施工总结
- 心电监护仪的日常维护及保养事项
- 提请批准延长羁押期限意见书范文
- 信用社(银行)信用户信用村信用镇评定办法