目前市面上的Julia教程很多,但大部分教程都是从0开始,对我这种有其他语言编程基础的十分不友好,我真的不关注你的Unicode支持有多棒好吗,类似这种东西在教程里大讲特讲对我来说属实没必要。另一部分则是过于高深,各种高级特性满天飞,本想着先拿来写点东西,看到奇怪的语法当场被劝退。恰好最近在学数值分析,借这个机会刚好练习一下Julia的使用。我的使用仅关注功能的实现,在语法上仅限能跑通,对其他细节则不去关注。
本文的代码可以在这里看到
AppliedNumericalAnalysisgithub.com先介绍一下本文用到的Julia语法
- Julia的索引从1开始,与Python、C++不一样
- Array即我们需要的矩阵
x = [1. 2. 3.; 2 7 5; 1 4 9] # 3 * 3的矩阵
b = [1; 6; -3.] # 长度为3的列向量
3. size用来获取矩阵的大小
size(x) # 返回二维矩阵的行数和列数
size(x, 1) # 返回矩阵的行数
size(x, 2) # 返回矩阵的列数
4. 矩阵的索引方法和matlab、python类似
x = [1 2 3; 4 5 6; 7 8 9]
x[1, 1] # 1
x[1, :] # 取第一行的值,返回一个列向量
x[end, end] # 9
x[:, 2] # 取第二列的值,返回一个列向量