




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、数值代数实验数值线性代数实验一一、实验名称:矩阵的 LU分解.二、实验目的:用不选主元的LU分解和列主元LU分解求解线性方程组 Ax=b,并比较这两种方法.三、实验内容与要求(1)用所熟悉的计算机语言将不选主元和列主元程序求解下面的84阶方程组r 6 18 (i 1:r28 61/3 , 8G1861f留8(j/ ,tS4/LU分解编成通用的子程序, 然后用编写的15 L5 *1515 14 7将计算结果与方程组的精确解进行比较,并就此谈谈你对Gauss消去法的看法.(2)写出追赶法求解三对角方程组的过程,并编写程序求该实验中的方程组Gauss消去法:用消去法解方程组的基本思想是用逐次消去未知
2、数的方法把原来方程组Ax=b化为与其等价的三角方程组,而求解三角方程组就容易了。换句话说,上 述过程就是用行的初等变换将原方程组系数矩阵化为简单形式,从而将求解原方程组的问题转化为求解简单方程组的问题。利用Gauss消去法对线性方程组 Ax=b进行求解。用MATLASt立m文件DelGauss.m,程序如下:function x=DelGauss(a,b)n,m=size(a);nb=length(b);det=1;x=zeros(n,1);for k=1:n-1for i=k+1:nif a(k,k)=0returnendm=a(i,k)/a(k,k);for j=k+1:na(i,j)=a
3、(i,j)-m*a(k,j);endb(i)=b(i)-m*b(k);enddet=det*a(k,k);enddet=det*a(n,n);for k=n:-1:1for j=k+1:nb(k)=b(k)-a(k,j)*x(j);endx(k)=b(k)/a(k,k);End在matlab中输入如下:结果如下:方程组的精确解为X1=X2=- =X84=1.0000,与Gauss消去法求得的解差距很大,所 得结果不够准确,计算简单但其消元过程有时不能进行到底而使求解出现解失真 的情况。数值线性代数实验二一、实验名称:实对称正定矩阵的A的Cholesky分解.二、实验目的:用平方根法和改进的平方
4、根方法求解线性方程组Ax=b三、实验内容与要求 用所熟悉的计算机语言将Cholesky分解和改进的Cholesky分解编成通用的子程序,然后用编写的程序求解对称正定方程组Ax=b,其中(1) b随机的选取,系数矩阵为100阶矩阵/ 1 I 1() 1 110 1+ - , + 110 1I L0 I 11口 /(2)系数矩阵为40阶Hilbert矩阵,即系数矩阵A的第i行第j列元素为 斯J = 壮向量b的第i个分量为八=出口(3)用实验一的程序求解这两个方程组,并比较所有的计算结果,然后评价各个方法的 优劣。平方根法:平方根法就是利用对称正定矩阵的三角分解而得到的求解对称正定方程组的一种有效方
5、法。平方根法递推公式可以证明对于对称正定矩阵A,可以唯一地分解成A=LLT,其中L是非奇异下三角形矩阵。模型二:利用平方根法对线性方程组Ax=b进行求解。用MATLABt立m文件pingfg.m,程序如下:function x=pingfg(A,b)%Cholesky 分解n,n=size(A); L=zeros(n,n);%实际上不用为 L申请空间,使用 A即可L(1,1)=sqrt(A(1,1);for k=2:nL(k,1)=A(k,1)/L(1,1);endfor k=2:n-1L(k,k)=sqrt(A(k,k)-sum(L(k,1:k-1).A2); for i=k+1:nL(i,
6、k)=(A(i,k)-sum(L(i,1:k-1).*L(k,1:k-1)/L(k,k); end endL(n,n)=sqrt(A(n,n)-sum(L(n,1:n-1).A2);%犁下三角方程组 Ly=by=zeros(n,1); for k=1:nj=1:k-1;y(k)=(b(k)-L(k,j)*y(j)/L(k,k);end%犁上三角方程组L'x=yx=zeros(n,1); U=L' for k=n:-1:1j=k+1:n;x(k)=(y(k)-U(k,j)*x(j)/U(k,k); End模型三:利用改进的平方根法对线性方程组 Ax=b进行求解%!改进平方法求解
7、Ax=b%1为n*n矩阵%M n*n的主对角矩阵%L的主对角元素均为 1用MATLABt立m文件ave.m,程序如下:function x=ave(A,b,n)L=zeros(n,n);D=diag(n,0);S=L*D;for i=1:nL(i,i)=1;end%佥证A是否为对称正定矩阵%人的特征值小于 0或A非对称时,输出 wrongfor i=1:nfor j=1:nif (eig(A)<=0)|(A(i,j)=A(j,i)disp( 'wrong' ); break ; end end endHilbert 矩阵用MATLABt立m文件Hil.m,程序如下: fu
8、nction b=Hil() for k=1:40for m=1:40s=0;t=s+1/(k+m-1);s=t; end b(k,1)=s;end在matlab中输入如下:在输入:输出为:问题3:Gauss消去法所得的结果与平方根法和改进的平方根法求得的结果差距很 大,而且Gauss消去法所得的结果大部分为零,显然平方根法和改进的平方根法 求得的结果与方程的精确解比Gauss消去法的更接近,更准确。但不管是哪一类 算法都只能在预定的计算步骤内或给定的精度内得到近似解,有一定的误差。数值线性代数实验三一、实验名称:矩阵 A的QR分解二、实验目的:应用改进的 Gram Schmidt方法和 Ho
9、useholder变换的方法计算矩阵 A 的 QR 分解.其中 A (aj) Rmn(m n), rank A = n三、实验内容与要求输入:A的各列a1,a2,冏防,aj" 1, ,n)输出:Q的各列元素(存放在A的相应位置上)以及R的元素rji 1, , n, j i, , n)数值线性代数实验四一、实验名称:用迭代法求解方程组及超松弛迭代和最佳松弛因子的确定二、实验目的: 应用Jacobi迭代法、Gauss Seidel迭代法和 超松弛迭代方法求解 线性方程组 Ax=b,并选择不同的松弛因子,观察松弛因子对松弛迭代法计算效果的影响.三、实验内容与要求dy dy a将常微分方程d
10、x2 dx (0 a 1 )离散化得到差分方程Ax=b,取y(0) 0,y(1) 11a 0.5,n 100,h ,应用Jacobi迭代法、Gauss Seidel迭代法和超松弛迭代 n方法求解线性方程组,分别取 1,0.1,0.01,0.0001 ,用SO越代法计算对应的数值解,并与精确解进行比较.写出这三种迭代法求解线性方程组的步骤,并对计算结果进行分析.四、实验原理将0,1区间n等分,方程离散化得差分方程2(h)yi 1 (2h)yiyi 1 ah差分方程对应的系数矩阵和右端项分别为(2 h) h(2 h) hA(2 h) OO O其中 y0y(0), yn y(i).SOR1代法的迭代矩阵和常数项分别为B (D L) 1 U (1)D, f (DSO砒代法的迭彳t公式:yk 1 B yk f五、实验过程ah2yoah2bMhah22, 一(2 h) (n 1) (n 1)ah ( h)yn (n 1) 1L) 1b利用Matlab语言实现SOR1代公式的计算, 分别对于 1,0.1,0.01,0.0001 ,松弛因子在(1, 2)范围内,先按步长 0.1搜索
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 个人贷款协议的书写要点
- 2025年《健康管理师》考试练习题及答案
- 线上教育培训协议模板
- 社交媒体影响者营销策略-洞察与解读
- 2025年福建省漳州市圆山劳务派遣服务有限公司招聘若干人考前自测高频考点模拟试题完整答案详解
- 2025年事业单位招聘考试综合类面试真题模拟试卷:情景模拟
- 跨境贸易合作协议执行保证函6篇
- 菏泽地理会考试题及答案
- 多版本回归测试技术-洞察与解读
- 2025国考无锡市监狱管理岗位行测必刷题及答案
- 2024-2025学年第一学期九年级期中考试英语试卷
- 人工智能 教育行业营销策略方案
- 医疗咨询行业营销策略方案
- 人教版八年级数学上册重难考点专题04因式分解(知识串讲+13大考点)特训(原卷版+解析)
- 左心发育不良综合征超声
- 招标文件范本三篇
- 幼儿园小班语言课件:《小狗找朋友》
- 南方主要虫害图谱-稻飞虱课件讲解
- 临床技术操作规范重症医学分册资料
- (正式版)JBT 106-2024 阀门的标志和涂装
- 心脏搭桥手术病历
评论
0/150
提交评论