目前市面上的Julia教程很多,但大部分教程都是从0开始,对我这种有其他语言编程基础的十分不友好,我真的不关注你的Unicode支持有多棒好吗,类似这种东西在教程里大讲特讲对我来说属实没必要。另一部分则是过于高深,各种高级特性满天飞,本想着先拿来写点东西,看到奇怪的语法当场被劝退。恰好最近在学数值分析,借这个机会刚好练习一下Julia的使用。我的使用仅关注功能的实现,在语法上仅限能跑通,对其他细节则不去关注。
本文的代码可以在这里看到
AppliedNumericalAnalysisgithub.com先介绍一下本文用到的Julia语法
- Julia的索引从1开始,与Python、C++不一样
- Array即我们需要的矩阵
x
3. size用来获取矩阵的大小
size
4. 矩阵的索引方法和matlab、python类似
x
5. 用.*可以实现对应元素相乘(和matlab一样)
6. Julia的map和python类似,第一个参数为表达式,第二个参数为对象
高斯消去法
高斯消去法是用来求解线性方程组的经典解法,其基本思路是对增广矩阵做行变换化为上三角矩阵后再回代求得方程组的解。以下均假设方程组有解。
首先定义函数的签名
function
其中,X是n阶方阵,b是n维列向量。
为了便于运算,我们将方阵和列向量合并为增广矩阵
A
接下来进行消去运算
for
以如下线性方程组为例
![2c54135965b019280ffddfd4f8773576.png](https://i-blog.csdnimg.cn/blog_migrate/79d3e7983a7fa17cd684cbc6d39664b2.png)
上述代码模拟的即是消去的过程
![ab3f2b16893756d97789f33733551042.png](https://i-blog.csdnimg.cn/blog_migrate/42d684c952fc6a9abb1ee8fa18653398.png)
在得到变换后的增广矩阵后即可进行回代
b
为了遵循课本的省内存的需求,将b作为结果返回。
高斯主元素法
高斯主元素法,又名列主元消去法,主要是为了解决在工程计算中由于计算机精度不足时除数过小所造成的精度丢失问题。其主要思想是,在每次消元时选取待选元素中最大的值作为当前目标。具体的思路可以在这里看到。
主元素法的基本思路和普通的消去法一致,仅多了选择最大列主元一项
target
参考文献
- 王开荣. 应用数值分析[M]. 高等教育出版社, 2010.