matlab解方程实验,MATLAB实验一解线性方程组的直接法

本次实验报告详细介绍了使用MATLAB进行线性方程组解法的实践,包括列主元高斯消去法和内部函数解法。实验比较了不同方法的解向量和条件数,分析了影响解的精度因素,并探讨了希尔伯特矩阵条件数与解分量有效位数的关系。通过实验,讨论了MATLAB内置函数与消去法的适用场景差异。
摘要由CSDN通过智能技术生成

41528d3028836879cd698677c3999917.gifMATLAB实验一解线性方程组的直接法

实 验 报 告 课程名称 数值分析 实验项目 解线性方程组的直接法 专业班级 姓 名 学 号 指导教师 成 绩 日 期 月 日 一. 实验目的 1、 掌握程序的录入和matlab的使用和操作; 2、 了解影响线性方程组解的精度的因素——方法与问题的性态。 3、 学会Matlab提供的“\”的求解线性方程组。 二. 实验要求 1、按照题目要求完成实验内容; 2、写出相应的Matlab 程序; 3、给出实验结果(可以用表格展示实验结果); 4、分析和讨论实验结果并提出可能的优化实验。 5、写出实验报告。 三. 实验步骤 1、用分解及列主元高斯消去法解线性方程组 a), 输出中系数分解的矩阵和,解向量和;用列主元法的行交换次序解向量和求;比较两种方法所得结果。 2、用列主高斯消元法解线性方程组。 (1)、 (2)、 分别输出,解向量,(1)中的条件数。分析比较(1)、(2)的计算结果 3、线性方程组的和分别为 , 则解. 用MATLAB内部函数求和的所有特征值和. 若令 , 求解,输出向量和,从理论结果和实际计算两方面分析线性方程组解的相对误差以及的相对误差的关系。 4、 希尔伯特矩阵,其中, (1)分别对计算,分析条件数作为的函数如何变化。(2)令,计算,然后用高斯消去法解线性方程组求出,计算剩余向量以及。分析当增加时解分量的有效位数如何随变化,它与条件数有何关系?当多大时连一位有效数字也没有了? 将每种情形的两个结果进行表格对比,如: n=6时: GAUSS列主消去法求得的 的有效数字 四、实验结果 五、讨论分析 (对上述算例的计算结果进行比较分析,主要说清matlab的算符与消去法的适用范围不同,自己补充) 六、改进实验建议 (自己补充) 1.列主元的高斯消去法 利用列主元的高斯消去法matlab程序源代码: 首先建立一个gauss.m的文件,用来实现列主元的消去方法。 function x=gauss(A,b) %高斯列主元消去法,要求系数矩阵非奇异的, % n = size(A,1); if abs(det(A))> a=[10 -7 0 1;-3 2.099999 6 2;5 -1 5 -1;2 1 0 2] a = 10.0000 -7.0000 0 1.0000 -3.0000 2.1000 6.0000 2.0000 5.0000 -1.0000 5.0000 -1.0000 2.0000 1.0000 0 2.0000 >> [l,u]=lu(a) l = 1.0000 0 0 0 -0.3000 -0.0000 1.0000 0 0.5000 1.0000 0 0 0.2000 0.9600 -0.8000 1.0000 u = 10.0000 -7.0000 0 1.0000 0 2.5000 5.0000 -1.5000 0 0 6.0000 2.3000 0 0 0 5.0800 >> b=[8 5.900001 5 1] b = 8.0000 5.9000 5.0000 1.0000 >> y=l\b y = 8.0000 1.0000 8.3000 5.0800 >> x1=U\x x1 = 0.0000

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值