全文预览已结束
下载本文档
版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
数值分析实验报告 一、实验目的与要求1掌握高斯消去法的基本思路和迭代步骤;2培养编程与上机调试能力。二、实验内容1编写用高斯消元法解线性方程组的MATLAB程序,并求解下面的线性方程组,然后用逆矩阵解方程组的方法验证.(1) (2)2编写用列主元高斯消元法解线性方程组的MATLAB程序,并求解下面的线性方程组,然后用逆矩阵解方程组的方法验证.(1) (2)三MATLAB计算源程序1. 用高斯消元法解线性方程组的MATLAB程序 输入的量:系数矩阵和常系数向量; 输出的量:系数矩阵和增广矩阵的秩RA,RB, 方程组中未知量的个数n和有关方程组解及其解的信息.function RA,RB,n,X=gaus(A,b)B=A b; n=length(b); RA=rank(A); RB=rank(B);zhica=RB-RA;if zhica0,disp(请注意:因为RA=RB,所以此方程组无解.)returnendif RA=RB if RA=ndisp(请注意:因为RA=RB=n,所以此方程组有唯一解.) X=zeros(n,1); C=zeros(1,n+1); for p= 1:n-1for k=p+1:n m= B(k,p)/ B(p,p); B(k,p:n+1)= B(k,p:n+1)-m* B(p,p:n+1);endend b=B(1:n,n+1);A=B(1:n,1:n); X(n)=b(n)/A(n,n); for q=n-1:-1:1 X(q)=(b(q)-sum(A(q,q+1:n)*X(q+1:n)/A(q,q); endelse disp(请注意:因为RA=RB0,disp(请注意:因为RA=RB,所以此方程组无解.)returnendif RA=RB if RA=ndisp(请注意:因为RA=RB=n,所以此方程组有唯一解.) X=zeros(n,1); C=zeros(1,n+1); for p= 1:n-1Y,j=max(abs(B(p:n,p); C=B(p,:);B(p,:)= B(j+p-1,:); B(j+p-1,:)=C;for k=p+1:n m= B(k,p)/ B(p,p); B(k,p:n+1)= B(k,p:n+1)-m* B(p,p:n+1);endend b=B(1:n,n+1);A=B(1:n,1:n); X(n)=b(n)/A(n,n); for q=n-1:-1:1 X(q)=(b(q)-sum(A(q,q+1:n)*X(q+1:n)/A(q,q); endelse disp(请注意:因为RA=RBn,所以此方程组有无穷多解.)endend3 实验过程: 1(1)编写高斯消元法的MATLAB文件如下: clear; A=0.101 2.304 3.555;-1.347 3.712 4.623;-2.835 1.072 5.643; b=1.183;2.137;3.035; RA,RB,n,X =gaus (A,b) 运行结果为: 请注意:因为RA=RB=n,所以此方程组有唯一解. RA = 3 RB = 3 n = 3 X = -0.3982 0.0138 0.3351 (2)编写高斯消元法MATLAB文件如下: clear; A=5 2 1;2 8 -3;1 -3 -6; b=8;21;1; RA,RB,n,X =gaus (A,b) 运行结果为: 请注意:因为RA=RB=n,所以此方程组有唯一解. RA = 3 RB = 3 n = 3 X = 1 2 -1 在MATLAB中利用逆矩阵法检验结果: (1) 在command windows中直接运行命令: A=0.101 2.304 3.555;-1.347 3.712 4.623;-2.835 1.072 5.643; b=1.183;2.137;3.035;X=Ab 运行结果为: X = -0.3982 0.0138 0.3351 (2) 在command windows中直接运行命令: A=5 2 1;2 8 -3;1 -3 -6; b=8;21;1;X=Ab 运行结果为: X = 1 2 -1两小题所得结果相同,检验通过 2(1)编写列组高斯消元法MATLAB文件如下: clear; A=0.101 2.304 3.555;-1.347 3.712 4.623;-2.835 1.072 5.643; b=1.183;2.137;3.035; RA,RB,n,X =liezhu(A,b) 运行结果: 请注意:因为RA=RB=n,所以此方程组有唯一解. RA = 3 RB = 3 n = 3 X = -0.3982 0.0138 0.3351 (2)编写列组高斯消元法的MATLAB文件如下: clear; A=5 2 1;2 8 -3;1 -3 -6; b=8;21;1; RA,RB,n,X =liezhu(A,b) 运行结果为: 请注意:因为RA=RB=n,所以此方程组有唯一解. RA = 3 RB = 3 n = 3 X = 1 2 -1与题 1 中逆矩阵计算所得结果相同,检验通过 四.实验体会: 通过实验我掌握了消元法解方程的一些基本算法以及用matlab实现矩阵的几种基本计算。对MATLAB软件
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2024年四川住院医师-四川住院医师妇产科笔试参考题库含答案
- 2024年广东高州市鉴龙水务建设投资集团有限公司招聘笔试参考题库含答案解析
- 2024年3月中国农业银行厦门市分行春季招考聘用60人笔试参考题库附带答案详解
- 2024年安徽肥西县乡村振兴建设发展有限责任公司招聘笔试参考题库含答案解析
- 2024年03月福建省宁德市畲族歌舞艺术传承中心(宁德市畲族歌舞团)2024年公开招考2名紧缺急需及高层次人才笔试参考题库附带答案详解
- 2024年03月湖南省中方县县直企事业单位2024年第一轮引进高层次急需紧缺人才笔试参考题库附带答案详解
- 2024年卫生资格(中初级)-精神病学主治医师笔试参考题库含答案
- 2024年03月浙江省台州市人力资源和社会保障局下属事业单位2024年招考编制外劳动合同工笔试参考题库附带答案详解
- 2024年03月浙江丽水莲都区融媒体中心招考聘用编外用工5人笔试参考题库附带答案详解
- 2024年03月江苏常州国家高新区新北区卫生健康局下属事业单位招考聘用18人笔试参考题库附带答案详解
- 不锈钢楼梯扶手安装合同
- 医疗技术临床应用论证制度
- 2024年国家电网内蒙古东部电力有限公司招聘笔试参考题库含答案解析
- 智慧口岸规划建设方案
- 体育产业合作框架协议
- 爱国卫生健康教育知识
- 山东省济南市历城二中2024届高考物理倒计时模拟卷含解析
- JB T 7392-2006数字压力表行业标准
- 2024临床免疫学定性检验程序性能验证指南
- 幼儿园游戏活动材料投放与指导课件
- 合肥家装行业分析
评论
0/150
提交评论