MATLAB求解3对角系数矩阵方程,实验5.3 用追赶法求解三对角方程组
实验5.3 用追赶法求解三对角方程组
实验目的
会使用Matlab 语言编程使用追赶法求解线性方程组。
实验原理
设系数矩阵为三对角矩阵
112223311100000000
000000
00
n n n n
n b c a b c a b A a b c a b ---?? ? ? ?=
? ? ?
? ??
?
则方程组Ax=f 称为三对角方程组。
设矩阵A 非奇异,A 有Crout 分解A=LU ,其中L 为下三角矩阵,U 为单位上三角矩阵,记
1122233110
00010
000001000
0000100,0000000000
0001n n n
n b L U γαβγββγβ--???? ?
?
? ? ? ?
?==
?
? ? ? ? ?
? ?
? ???
?
?
?
可先依次求出L ,U 中的元素后,令Ux=y ,先求解下三角方程组Ly=f 得出y ,再求解
上三角方程组Ux=y 。
事实上,求解三对角方程组的2追赶法将矩阵三角分解的计算与求解两个三角方程组的计算放在一起,使算法吏为紧凑。其计算公式为:
1111,
1111
,111
,2,3,,,1,2,,1i
i i i i i i i i
i i i i i n n
i i i i c f b y i n c a b a f y y x y i n n x y x βγββαβγγβαβγ--+?
===??
=??
?==-=
???
-?=??
=??=--?=-?? 对对
实验程序
MATLAB求解3对角系数矩阵方程,实验5.3 用追赶法求解三对角方程组相关推荐
- C 语言求解N阶线性矩阵方程Ax=b,用C语言求解N阶线性矩阵方程Ax=b的简单解法
首先,为了能够求解N阶线性方程组(N由用户输入),所以需要定义一个大于N维的数组a[dim+1][dim+1](dim为设定的最大维数,防止计算量溢出),当用户输入的阶数N超过设定值时提示重启程序重新 ...
- c语言程序 用追赶法求解方程组,编写用追赶法解三对角线性方程组的程序,并解下列方程组(3页)-原创力文档...
计算方法与实习上机实验(二) 实验名称: 编写用追赶法解三对角线性方程组的程序,并解下列方程组: (1) (2)Ax=b,其中 A10×10=, b10×1= 程序代码: #include using ...
- c语言程序 用追赶法求解方程组,编写用追赶法解三对角线性方程组的程序,并解下列方程组...
计算方法与实习上机实验(二) 实验名称: 编写用追赶法解三对角线性方程组的程序,并解下列方程组: 2x1 x2 5, x 2x x 12, 123(1) x 2x x 11,34 2 x3 2x4 1 ...
- Matlab实现——严格对角占优三对角方程组求解(高斯赛尔德Gauss-Seidel迭代、超松弛)
欢迎前往个人博客 驽马点滴 和视频空间 哔哩哔哩-<挨踢日志> 严格对角占优三对角方程组求解 对中等规模的n阶的(n<100)线性方程组,直接法的准确性和可靠性,所以常采用直接法 对 ...
- matlab trisys,Matlab实现——严格对角占优三对角方程组求解(高斯赛尔德Gauss-Seidel迭代、超松弛) | 学步园...
严格对角占优三对角方程组求解 对中等规模的n阶的(n<100)线性方程组,直接法的准确性和可靠性,所以常采用直接法 对于较高阶的方程组,特别是地于某些偏微分方程离散化后得到的大型稀疏方程组(系统 ...
- matlab编程数字信号,MATLAB--数字信号实验.doc
MATLAB--数字信号实验 前言 MATLAB 是一套功能强大的工程计算及数据处理软件,广泛应用于工业,电子,医疗和建筑等众多领域.它是一种面向对象的,交互式程序设计语言,其结构完整又优良的可移植性 ...
- 电力系统matlab实验报告,基于matlab语言计算电力系统暂态稳定仿真程序实验报告.docx...
基于matlab语言计算电力系统暂态稳定仿真程序实验报告 BeijingJiaotongUniversity 电力系统分析 暂态稳定分析实验 学院:电气工程学院 班级:xxxxxxxx 学号:xxxx ...
- 遍历所有点的最短路径matlab_运筹学实验8 最短路的求解
实验8 最短路的求解 成 绩 实验类型:◆验证性实验 ◇综合性实验 ◇设计性实验 实验目的:学会使用Matlab求解最短路. 实验内容:1.Floyd算法:2.利用Matlab编程实现 ...
- matlab波形相似系数,如何用matlab计算波形相似系数
. 低通采样定理实验 1.1 实验目的 1.了解数字信号处理系统的一般构成: 2.掌握奈奎斯特抽样定理. 1.2 实验仪器 1.YBLD智能综合信号源测试仪 1台 2.双踪示波器 1台 3.MCOM- ...
- 基于matlab的传热学虚拟实验开发,基于MATLAB的传热学课程虚拟实验软件的开发
215教育现代化·2018 年 12 月第 49 期 教育信息技术 基于 MATLAB 的传热学课程虚拟实验软件的开发 周永利,李友荣,石万元,张力元,杨晨,卞煜,王国强,李俊,包键 ( 重庆大学 低 ...
最新文章
- 最新版本的Silverlight Tools For Silverlight 2 RTW正式版
- 协议森林13 9527 (DNS协议)
- cocos2d-x 2x 菜鸟编塔防 02 添加敌人
- Jeecg平台扩展性不好的地方收集启动。
- 调用者保存寄存器和非调用者保存寄存器
- php echo nbsp,关于include里面的函数echo的问题
- 揭露动画(Reveal Effect)实现时的注意事项(附上bug-logcat)...
- c语言ole错误,错误 '8002801c' 访问 OLE 注册表的错误
- [转]centos7.2 下 nginx 开机启动
- 局部变量AND全局变量
- 计蒜客 第一场 搜狗输入法
- Ci522/CI523替代RC522手机NFC开锁方案13.56Mhz读写器芯片
- github 更纱黑体_更纱黑体v0.12.6
- 关于lora和lorawan所涉及的名词解释
- 《Linux设备驱动开发详解 A》一一3.1 Linux内核的发展与演变
- FPGA学习笔记(1)简单的时序逻辑电路——流水灯
- 手游平台系统搭建sdk服务端接口文档
- win7python怎么设置环境变量_如何在win7下设置python的环境变量
- PCM音频 I2S音频
- java:计算π的值
热门文章
- [转载野猪的一篇文章] 设计模式--简单工厂模式在unity3d里面的使用
- 看完《一九四二》的感受--也写给正在找工作时迷茫的你
- 伟大的民族英雄赵充国
- 法国电子与计算机信息工程学校排名,法国工程学院的十大排名情况
- 【安全】如何关闭Windows系统的137, 139, 445端口?
- android 视频连续播放,VideoView实现视频无缝连续播放
- 千万不要急着 返城上班!
- pdf文件拆分为单个pdf_如何根据文件内容拆分重命名移动一批PDF文件
- 【ctfshow】- web189
- MyEclipse出现弹窗提示“当前页面的脚本发生错误”