MatLab实现追赶法
最近在学数值分析,导师要求编写书上的程序,我的内心belike:这对我一个在IT行业干过的半个前端人士不是嘎嘎乱杀嘛~~!!!!
但是实际上写程序不算什么,重点是里面的数学逻辑,我在网上找了一些别人写的追赶法代码,看了半天,觉得有些复杂了,还是决定自己写吧。下面是我的代码,希望可以帮到有需要的孩子们哦!
简单说一下我的代码思路:我是用了三个向量来存储三对角矩阵的三个对角元素的,这样可以节省内存空间,减少不必要的计算,进而增加程序效率;最后将方程组的解存放在原来方程组的右边向量D里。(其实最后只要输出D就行了,但是我想看看我的主对角线上元素有没有问题,所以额外输出的一个向量B)
运行结果如下:(个人感觉我的代码还算简洁,没有参考别书本上的代码,纯属自己靠感觉来的!要是有好的代码可以一起学习一下呀~~)