




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
计算方法实验报告1【课题名称】用列主元高斯消去法和列主元三角分解法解线性方程【目的和意义】高斯消去法是一个古老的求解线性方程组的方法,但由它改进得到的选主元的高斯消去法则是目前计算机上常用的解低阶稠密矩阵方程组的有效方法。用高斯消去法解线性方程组的基本思想时用矩阵行的初等变换将系数矩阵A约化为具有简单形式的矩阵(上三角矩阵、单位矩阵等),而三角形方程组则可以直接回带求解用高斯消去法解线性方程组(其中A∈Rn×n)的计算量为:乘除法运算步骤为,加减运算步骤为。相比之下,传统的克莱姆法则则较为繁琐,如求解20阶线性方程组,克莱姆法则大约要次乘法,而用高斯消去法只需要3060次乘除法。在高斯消去法运算的过程中,如果出现abs(A(i,i))等于零或过小的情况,则会导致矩阵元素数量级严重增长和舍入误差的扩散,使得最后的计算结果不可靠,所以目前计算机上常用的解低阶稠密矩阵方程的快速有效的方法时列主元高斯消去法,从而使计算结果更加精确。2、列主元三角分解法高斯消去法的消去过程,实质上是将A分解为两个三角矩阵的乘积A=LU,并求解Ly=b的过程。回带过程就是求解上三角方程组Ux=y。所以在实际的运算中,矩阵L和U可以直接计算出,而不需要任何中间步骤,从而在计算过程中将高斯消去法的步骤进行了进一步的简略,大大提高了运算速度,这就是三角分解法采用选主元的方式与列主元高斯消去法一样,也是为了避免除数过小,从而保证了计算的精确度【计算公式】列主元高斯消去法设有线性方程组Ax=b,其中设A为非奇异矩阵。方程组的增广矩阵为第1步(k=1):首先在A的第一列中选取绝对值最大的元素,作为第一步的主元素:然后交换(A,b)的第1行与第l行元素,再进行消元计算。设列主元素消去法已经完成第1步到第k-1步的按列选主元,交换两行,消元计算得到与原方程组等价的方程组A(k)x=b(k)第k步计算如下:对于k=1,2,…,n-1(1)按列选主元:即确定t使(2)如果t≠k,则交换[A,b]第t行与第k行元素。(3)消元计算消元乘数mik满足:(4)回代求解列主元三角分解法对方程组的增广矩阵经过k-1步分解后,可变成如下形式:clearx;disp('AX=b的解x是')x=b;【调用函数解题】【列主元三角分解法Matlab主程序】clear;augm=input('Pleaseinputtheaugumentalmatric:');[n,t]=size(augm);fori=1:n%d控制列三角主元的层数b=augm(i,:);p=i;%4到10行是选列主元并交换forw=i+1:nifb(1,1)<=augm(w,i)b=augm(w,:);p=w;endendc=augm(i,:);augm(i,:)=b;augm(p,:)=c;ifp~=i%只有发生了换行才将这种效果输出augm,endforj=i:t%首先变换与augm(i,i)同一行的元素,其列指标从i到ts=0;fork=1:i-1s=s+augm(i,k)*augm(k,j);endaugm(i,j)=(augm(i,j)-s);endforI=i+1:n%再变换与augm(i,i)同一列的元素,其行指标I从i+1到n,列数为is=0;fork=1:i-1%下三角部分列数与第i层一致s=s+augm(I,k)*augm(k,i);endaugm(I,i)=(augm(I,i)-s)/augm(i,i);endaugm,endx(n)=augm(n,t)/augm(n,n);%回代fori=n-1:-1:1s=0;forj=n:-1:i+1s=s+x(j)*augm(i,j);endx(i)=(augm(i,t)-s)/augm(i,i),end【调用函数解题】输入[Ab]可以清晰地看到augm(A)的变化过程:结果输出:【列主元三角分解法程序流程图】【编程疑难】这是第一次用matlab编程,对matlab的语句还不是非常熟悉,因此在编程过程中,出现了许多错误提示。并且此次编程的两种方法对矩阵的运算也比较复杂。问题主要集中在循环控制中,循环次数多了一次或者缺少了
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 印花铰链行业深度研究分析报告(2024-2030版)
- 中国精准营销行业发展监测及投资战略规划研究报告
- 银制艺术品项目投资可行性研究分析报告(2024-2030版)
- 2025年中国酒饮冲调行业市场发展前景及发展趋势与投资战略研究报告
- 分公司企业管理办法
- 区长epc管理办法
- 劳务公司安全制度
- 融创考勤管理制度
- 会议室对外出租管理制度
- 奥体中心考勤管理制度
- 上海华东师大三附中2025年物理高一下期末调研模拟试题含解析
- 2025至2030科技金融行业市场发展分析及前景趋势与行业市场深度研究及发展前景投资可行性分析报告
- 养老机构管理课件
- 加密数字资产管理办法
- 仓库账务知识培训课件
- 哈三中2024-2025学年度高一下学期期末考试(选考)物理试题含答案
- 2026版高三一轮总复习(数学) 第六章 第5课时 数列的综合应用 课件
- XXX医院保安部合作单位安全培训计划
- 设计投标响应方案(3篇)
- 工程设计文档编制规范
- 精密数控机床制造项目运营管理方案
评论
0/150
提交评论