Octave快速入门（转载）
原文链接https://blog.csdn.net/iszhenyu/article/details/78712228
以下是我的补充：
PS1('>>') 将原来的octave:1>提示输入改成>>提示输入。
octave:75> PS1('>>')
>>randn(m,n) 函数生成一个m行n列的矩阵，矩阵的所有元素服从高斯分布，均值为0，标准差或者方差为1。
hist(变量）为对一个随机变量画出直方图。
octave:5> w=-6+sqrt(11)*(randn(1,10000));
octave:6> hist(w)矩阵/向量
exp(向量) 求以e为底数，以向量中的元素为指数的幂运算
octave:36> exp([1 2 3])
ans =
2.7183 7.3891 20.0855其他一些运算：
% a = [1 15 2 0.5]，求最大值
>> val = max(a)
val = 15
% 求最大值，并返回最大值的索引
>> [val, idx] = max(a)
val = 15
idx = 2
% 矩阵对应元素的逻辑运算
>> a <= 1
ans =
1 0 0 1
% 返回a中小于3的元素的索引
>> find(a < 3)
ans =
1 3 4
% 计算之和
>> sum(a)
ans = 18.500
% 计算乘积
>> prod(a)
ans = 15
% 向下取整（向下取最接近的整数，即直接舍掉小数部分）
>> floor(a)
ans =
1 15 2 0
% 向上取整（向上取最接近的整数，即小数部分变为加1）
>> ceil(a)
ans =
1 15 2 1
% 生成一个随机矩阵，矩阵元素的值位于0-1之间
>> rand(3)
ans =
0.458095 0.323431 0.648822
0.481643 0.789336 0.559604
0.078219 0.710996 0.797278
% 返回一个矩阵，这个矩阵的每个元素取两个随机矩阵对应元素中的最大值的结果
octave:56> max(rand(3),rand(3))
ans =
0.916676 0.875362 0.881009
0.409057 0.489740 0.069200
0.494241 0.633359 0.961897
octave:57> A
A =
1 2
3 4
5 6
% 得到每一列的最大值作为一个行向量，默认。
octave:58> max(A,[],1) % 同 max(A)
ans =
5 6
% 得到每一行的最大值作为一个列向量
octave:59> max(A,[],2)
ans =
2
4
6
% 求整个矩阵最大值的两种方法：
% 1.先求每一列的最大值作为一个行向量，再求最大值
octave:61> max(max(A))
ans = 6
% 2.先将A的所有元素变为一个列向量，再求最大值
octave:63> max(A(:))
ans = 6
% 矩阵按行上下对换
>> flipud(eye(4))
ans =
Permutation Matrix
0 0 0 1
0 0 1 0
0 1 0 0
1 0 0 0magic(n）生成一个n维的幻方（每行、列和对角线上的数字和都相等的矩阵）
octave:65> M=magic(8)
M =
64 2 3 61 60 6 7 57
9 55 54 12 13 51 50 16
17 47 46 20 21 43 42 24
40 26 27 37 36 30 31 33
32 34 35 29 28 38 39 25
41 23 22 44 45 19 18 48
49 15 14 52 53 11 10 56
8 58 59 5 4 62 63 1
octave:66> sum(M)
ans =
260 260 260 260 260 260 260 260
octave:67> sum(M,2)
ans =
260
260
260
260
260
260
260
260
% 以下是得到对角线和
octave:70> M .* eye(8)
ans =
64 0 0 0 0 0 0 0
0 55 0 0 0 0 0 0
0 0 46 0 0 0 0 0
0 0 0 37 0 0 0 0
0 0 0 0 28 0 0 0
0 0 0 0 0 19 0 0
0 0 0 0 0 0 10 0
0 0 0 0 0 0 0 1
octave:71> sum(M .* eye(8))
ans =
64 55 46 37 28 19 10 1
octave:72> sum(sum(M .* eye(8)))
ans = 260
% 以下是求另一条对角线的和
octave:73> flipud(eye(8)) % 使矩阵垂直翻转
ans =
Permutation Matrix
0 0 0 0 0 0 0 1
0 0 0 0 0 0 1 0
0 0 0 0 0 1 0 0
0 0 0 0 1 0 0 0
0 0 0 1 0 0 0 0
0 0 1 0 0 0 0 0
0 1 0 0 0 0 0 0
1 0 0 0 0 0 0 0
octave:74> sum(sum(M .* flipud(eye(8))))
ans = 260>>t=[0:0.01:0.98];
>>y1 = sin(2 * pi * 4 * t);
>>y2 = cos(2 * pi * 4 * t);
>>plot(t,y1);画出的sin函数：
>>subplot(1,2,1); %将图像分为一个1*2的格子，现在要使用第一个格子画图
>>plot(t,y1); % 在当前格子里画那个正弦函数 
>>subplot(1,2,2);
>>plot(t,y2);
>>axis([0.5 1 -1 1]); % 将横坐标范围改为0.5~1，纵坐标范围改为-1~1。注意因为此时在第二个格子，因此只会改第二个
clf; 清除当前绘图的整个图像
一种数据可视化方向
现在先建一个矩阵
>>A=magic(5);
>>A
A =
17 24 1 8 15
23 5 7 14 16
4 6 13 20 22
10 12 19 21 3
11 18 25 2 9
>>imagesc(A) %将矩阵A中的元素数值按大小转化为不同颜色，并在坐标轴对应位置处以这种颜色染色
>>imagesc(A),colorbar; % 同时运行这两个函数，colorbar函数显示一个颜色条
将颜色变成灰度的就可以。这样就数值越大的颜色越深，数值越小的颜色越浅。
函数
我们还是先看例子，然后再说明具体的写法：
>> function y = squareNum(x)
y = x^2;
end;
>> squareNum(3)
ans = 9
1
2
3
4
5在Octave中，定义一个函数需要使用function 关键字，然后紧跟在 function 后面的是函数的声明，包括返回值，函数名称和参数，之后换行来实现具体的函数功能。
Octave的函数不需要显示的返回语句，Octave会将函数第一行声明的返回值返回给调用方，因此，我们在函数体中只需将最终的计算结果赋给定义的返回值，比如上面例子中的y。
也可以将函数部分写在文件里，文件名命名为“函数名.m”，要注意的就是此.m文件必须在你当前Octave路径下。即以下部分为此文件内容：
function y = squareNum(x)
y = x^2;然后在Octave中就可以
还有一点需要说明的是，在Octave中，函数可以返回多个值：
>> function [y1, y2] = calVal(x)
y1 = x^2;
y2 = x^3;
end;
>> [a, b] = calVal(3)
a = 9
b = 27