文章目录

简介
基本
符运算
- 算术运算符
- 逻辑运算符
- 赋值运算符
- 格式化输出
- 向量和矩阵
- 其他
移动数据
计算
绘图
基础语法

老师推荐人工智能最简单入手的是octave，用它写出来之后再转用python、c++会更方便。

简介

GNU Octave is a high-level language primarily intended for numerical computations. It is
typically used for such problems as solving linear and nonlinear equations, numerical linear
algebra, statistical analysis, and for performing other numerical experiments. It may also
be used as a batch-oriented language for automated data processing.

GNU Octave是一种高级语言，主要用于数值计算。它通常用于求解线性和非线性方程、数值线性代数、统计分析等问题，以及进行其他数值实验。它也可以用作自动化数据处理的面向批处理的语言。

这只是做个简单入门，更多的功能可以看一下官方文档：octave.pdf

基本

注释 %
- 类似于C++的//和python的#
ans 用于存储输出结果
如果你不喜欢输入行前缀，可以自定义修改PS1('你想改成的内容');
- 这个不是永久的，你电脑重启之后他又恢复默认了。
help 帮助
如果你忘了某一个语句的用法了，可以用这个查询。

比如输入help eye就会显示eye的用法。

5.多个语句同时执行用逗号分开即可。
c octave:226> a=4,b=5,c=6 a = 4 b = 5 c = 6
6. quit，exit可以关闭Octave（如果你不想点右上角的×的话）
7. addpath('路径')添加路径

符运算

算术运算符

加减乘除乘方

octave:1> 1+5
ans = 6
octave:2> 1*5
ans = 5
octave:3> 1/5
ans = 0.2000
octave:4> 1-5
ans = -4
octave:5> 5^2
ans = 25

逻辑运算符

等于 ==
不等于~=
或 ||
与 &&
异或 xor(a,b)
大于等于>=
大于 >

ans = 0
>>1==1
ans = 1
>>1 ~=2
ans = 1
>>1&&2
ans = 1
>>1&&0
ans = 0
>>1||0
ans = 1
>>xor(1,1)
ans = 0
>>xor(1,2)
ans = 0
>>xor(0,1)
ans = 1
>>1>=2
ans = 0
>>1<2
ans = 1

赋值运算符

=赋值运算符
- 语句后边加分号不直接打印赋值语句
- 语句后边不加分号直接打印赋值语句

>>a=1   %a=1后边没分号，所以输入回车之后下一行会直接打印a=1
a = 1
>>a     %再输入a回车看看是否赋值成功
a = 1
>>b=2;  %b=2后边加上分号了，所以下一行不会直接打印b的赋值
>>b     %输入b回车，查看给b的赋值
b = 2   >>c=pi
c = 3.1416

格式化输出

disp(x)：打印x的值
disp(sprintf('%0.2f',c))：格式化输出保留两位小数，这个和C语言一样
format long/short：对整数无影响，对小数输出格式有影响

>>disp(a)
1
>>disp(c)
3.1416
>>disp(sprintf('%0.2f',c))
3.14>>format long
>>a
a = 1
>>c
c = 3.141592653589793
>>format short
>>a
a = 1
>>c
c = 3.1416

向量和矩阵

声明矩阵官方原话写的是：

Vectors and matrices are the basic building blocks for numerical analysis. To create a new matrix and store it in a variable so that you can refer to it later, type the command

octave:1> A = [ 1, 1, 2; 3, 5, 8; 13, 21, 34 ]

Octave will respond by printing the matrix in neatly aligned columns.

>>martix = [1,2,3;4,5,6;7,8,9]
martix =1   2   34   5   67   8   9

当然还有很多其他写法，虽然不规范也能声明矩阵，这里就不列举了。

声明行向量，各个元素之间不用加符号，或者加逗号。

martix =1   22   33   4>>vector=[1 2 3]
vector =1   2   3>>vector = [4;5;6]
vector =456

其他特殊写法：

变量 = 起始:终止 ，生成一个向量，并且步长为1

变量 = 起始:步长:终止 ，就会生成一个向量

> v = 1:5
v =1   2   3   4   5> m = 1:0.5:1
m = 1
octave:3> m = 1:0.5:5
m =1.0000    1.5000    2.0000    2.5000    3.0000    3.5000    4.0000    4.5000    5.0000

ones(a,b) 生成a行b列只有1的矩阵
- 除了ones还有zeros，生成的元素都是0，但是没有twos、threes……
- ```
  >ones(2,3)ans =1   1   11   1   1
```
rand(n) 形成n阶方阵，元素都是0-1之间随机数

rand(a,b) 生成a行b列的0-1之间随机数的矩阵

  >rand(2,3)ans =0.4046   0.4508   0.80210.6986   0.8620   0.2631

max(rand(n),rand(n))生成量个方阵，取其中较大值形成一个新方阵

octave:124> c = max(rand(4),rand(4))
c =0.8010   0.9926   0.7241   0.90530.9007   0.3026   0.9856   0.57100.7679   0.2630   0.5200   0.96150.7333   0.5113   0.6957   0.5683

randn(a,b) 也是生成a行b列随机数矩阵，但是符合高斯分布（正态分布）。返回一个具有正态分布随机元素的矩阵，这些随机元素的平均值为零，方差为1。
- 正态分布：

randi()

randi(max,m,n)：生成1~max的m行n列矩阵
randi([min,max],m,n)：生成min~max的m行n列矩阵

  >randi(5,3,3)ans =1   5   55   2   15   1   1

  >randi([5,10],5,5)ans =6    8    6    8    510    9   10    8    75    7    7    9   1010    8   10    8    810    9    7    6    7

eye(n) 生成n维单位矩阵

>eye(5)
ans =Diagonal Matrix1   0   0   0   00   1   0   0   00   0   1   0   00   0   0   1   00   0   0   0   1

matrixB = matrixA(x:y,v:w):将matrixA的x~y行v~w列形成一个新矩阵
- :表示所有元素
- a:b表示a到b的元素
- [a b c]表示abc三行或列
- 以上三种表示方法在矩阵别的方法中也能使用，比如赋值、打印

矩阵拼接

在矩阵右侧添加一列

octave:30> A
A =4    2    7    51   10    7   105   10    9    9octave:31> A = [A,[2;4;5]]
A =4    2    7    5    21   10    7   10    45   10    9    9    5

matrixC=[matrixA matrixB] 横向拼接
写[matrixA, matrixB]也一样，加不加逗号都行

octave:33> A = randi(5,2,3)
A =4   1   32   3   5octave:34> B = randi(5,2,3)
B =5   4   51   1   3octave:35> c = [A B]
c =4   1   3   5   4   52   3   5   1   1   3

matrixC=[matrixA; matrixB]纵向拼接

octave:33> A = randi(5,2,3)
A =4   1   32   3   5octave:34> B = randi(5,2,3)
B =5   4   51   1   3octave:36> d=[A;B]
d =4   1   32   3   55   4   51   1   3

matrix(:) 将矩阵中的所有元素放入一个列向量

octave:30> A
A =4    2    7    51   10    7   105   10    9    9octave:31> A(:)
ans =415210107795109245

其他

hist(变量)

可以画出变量的图像

如果你觉得这样太少了。那可以使用hist(变量，条纹数)，下图就是画出50条纹后的图像。
检测大小
- size()
- length()
下图可以看出，length显示矩阵的长度。size显示矩阵有几行几列。size(矩阵名称，1/2)其中1代表行参数，2代表列参数。
```
v =9    9    5    8    85   10    7    6    86    7    7    7    67   10    8    7    97   10    7    5    6>length (v)
ans = 5
>size(v)
ans =5   5
>size(v,1)
ans = 5
```

索引：

matrix(i,j) 显示第i行第j列的元素

octave:26> A = randi(10,3,4) %生成一个1~10之间随机数的 3×4 的矩阵
A =4    2    7    51   10    7   105   10    9    9
octave:27> A(2,3)            %打印这个矩阵2行3列的元素
ans = 7

matrix(:,[a,b,c...])

:：代表全部
[]：方括号中有哪几个就选中哪几个

  octave:26> A = randi(10,3,4)A =4    2    7    51   10    7   105   10    9    9octave:27> A(:,[1,2,4])   %显示A所有行的1,2,4列ans =4    2    51   10   105   10    9

显示现有变量
- who：显示目前所有变量
- whos：显示目前所有变量和详细信息

删除某个变量 clear

clear 变量 删除三个变量
`clear 删除所有变量

octave:12> a = 15   % 给a赋值15
a = 15
octave:13> a        % 打印a
a = 15              % 显示a的值
octave:14> clear a
octave:15> a        % 再次打印a
error: 'a' undefined near line 1, column 1    %报错

移动数据

用过其他命令行的应该能适应，octave也同样支持某些命令。

pwd：列出当前目录
cd xxx文件夹：进入xxx文件夹
ls：列出当前目录信息
加载文件：
- load 文件名
- load('文件名')
save 文件名变量; ：将变量保存到文件中，比如save data.txt num，把num变量保存到data.txt中（如果操作目录中没有这个文件会创建一个新文件，如果目录中有这个文件会重写）

下边给出两个文件的数据，featureX和priceY，自己搞两个文件存起来。
featureX:

priceY:

首先进入对应的操作目录。

用who查看当前有的变量。

打印依稀下X6price确实已经成功读入文件

用size查看一下，导入的数据是存在矩阵中的。

把矩阵X6feature的1~5行，第2列赋值给矩阵x，并将矩阵x保存到data.txt中

计算

除了最开始提到的算术运算符，还有以下运算。

注意一个特殊的.*，矩阵同位相乘，即矩阵相同位置上的元素相乘。

octave:37> A
A =4   1   32   3   5octave:38> B
B =5   4   51   1   3octave:39> A * B    %AB无法相乘
error: operator *: nonconformant arguments (op1 is 2x3, op2 is 2x3)
octave:40> A .* B   %AB可以点乘
ans =20    4   152    3   15

.^给矩阵每个元素做乘方

octave:43> A
A =4   1   32   3   5octave:44> A .^2
ans =16    1    94    9   25octave:45>

矩阵/数即可让矩阵的每个元素都除以该数字，但是如果你想让一个数除以矩阵的每个元素，结果仍旧是矩阵，那就需要./

octave:38> B
B =5   4   51   1   3
octave:50> 2 ./B
ans =0.4000   0.5000   0.40002.0000   2.0000   0.6667

如果想让矩阵每个元素都加或减去一个数，直接用+-即可。

octave:71> C
C =2  -1  -4   1-5  -2  -1  -11  -1  -5  -4octave:72> C+1
ans =3   0  -3   2-4  -1   0   02   0  -4  -3octave:73> C-2
ans =0  -3  -6  -1-7  -4  -3  -3-1  -3  -7  -6octave:74> 1+C
ans =3   0  -3   2-4  -1   0   02   0  -4  -3

log(n)这里的log底数是e，即 log ⁡ e n \log_en logen或 ln ⁡ n \ln n lnn

octave:57> log(2)
ans = 0.6931
octave:58> log(e)
ans = 1
octave:59> A
A =4   1   32   3   5octave:60> log(A)
ans =1.3863        0   1.09860.6931   1.0986   1.6094

exp(n)相当于执行 e n e^n en

octave:61> exp(1)
ans = 2.7183
octave:62> exp(2)
ans = 7.3891
octave:63> A
A =4   1   32   3   5octave:64> exp(A)
ans =54.5982     2.7183    20.08557.3891    20.0855   148.4132

abs()求绝对值

octave:66> C = randi([-5,2],3,4)    %随机生成一个3×4的矩阵，取值在-5\~2之间
C =2  -1  -4   1-5  -2  -1  -11  -1  -5  -4octave:67> abs(-1)
ans = 1
octave:68> abs(C)
ans =2   1   4   15   2   1   11   1   5   4

matrix'转置，矩阵名称后加英文单引号

octave:75> C
C =2  -1  -4   1-5  -2  -1  -11  -1  -5  -4octave:76> C'
ans =2  -5   1-1  -2  -1-4  -1  -51  -1  -4

flipud使矩阵垂直翻转。就是左边变右边，右边变左边。

result1 =17    0    0    0    00    5    0    0    00    0   13    0    00    0    0   21    00    0    0    0    9octave:140> flipud(result1)
ans =0    0    0    0    90    0    0   21    00    0   13    0    00    5    0    0    017    0    0    0    0

矩阵<数值大于啊等于啊都同理。会返回一个同纬度的矩阵，每个元素都是对应位置元素进行比较返回的布尔值

a =1   5   3   4octave:85> a<2
ans =1  0  0  0

magic(n)生成n阶幻方，所谓幻方就是行、列、对角线加起来都是相同的值

octave:95> magic(4)
ans =16    2    3   135   11   10    89    7    6   124   14   15    1octave:96> magic(3)
ans =8   1   63   5   74   9   2

find(公式)

找到向量中符合的数据并返回其索引

a =1   5   3   4octave:88> find(a<3)
ans = 1
octave:89> find(a>2)
ans =2   3   4

找到矩阵中符合的数据

C =9    7    3    5
4    2    5   10
4    1    9   10octave:114> [r,c]=find(C==1)
r = 3
c = 2
octave:115> [r,c]=find(C>7)
r =1323c =1344

sum()求和prod()求积

c =1344octave:117> sum(c)
ans = 12
octave:118> prod(c)
ans = 48

floor()向下取整ceil()向上取整

a =0.50001.50002.00002.0000octave:120> floor(a)
ans =0122octave:121> ceil(a)
ans =1222

max(矩阵)取最大值
如果是矩阵，那就显示每一列的最大值；如果是向量，那就只显示一个最大值

v=max(matrix)：将矩阵的最大值赋给v，不是硬性要求，v可以换别的名称
[v,i]=max(matrix)：将矩阵最大值及其位置赋值给v和i
max(max(matrix))：取得整个矩阵最大值
max(matrix(:))：取得整个矩阵最大值，这个是将矩阵转化成响亮之后再求最大值

c =2   1   4   15   2   1   11   1   5   4octave:78> val = max(c)
val =5   2   5   4octave:79> [val,ind]=max(c)
val =5   2   5   4ind =2   2   3   3

max(矩阵,[],1/2) 取得矩阵每一行或者列的最大值形成一个向量

1 取每一列
2 取每一行

c =8   9   7   99   3   9   57   2   5   97   5   6   5octave:129> max(c,[],1)
ans =9   9   9   9octave:130> max(c,[],2)
ans =9997

小练习：
生成一个幻方，然后测试其每行，每列，对角线之和都是同一个数

octave:131>  m = magic(5)     %生成一个5×5的幻方
m =17   24    1    8   1523    5    7   14   164    6   13   20   2210   12   19   21    311   18   25    2    9octave:132> sum(m,1)          %每列相加都是65
ans =65   65   65   65   65octave:133> sum(m,2)          %每行相加都是65
ans =6565656565octave:135> e=eye(5)          %生成一个单位矩阵
e =Diagonal Matrix1   0   0   0   00   1   0   0   00   0   1   0   00   0   0   1   00   0   0   0   1octave:136> result1 = e .* m   %单位矩阵同位相乘原矩阵，新矩阵只剩主对角线元素
result1 =17    0    0    0    00    5    0    0    00    0   13    0    00    0    0   21    00    0    0    0    9
octave:139> sum(result1(:))   %主对角线相加为65
ans = 65octave:141> result2 = m .* flipud(e)
result2 =0    0    0    0   150    0    0   14    00    0   13    0    00   12    0    0    011    0    0    0    0octave:142> sum(result2(:))   %副加为65
ans = 65

pinv()求逆矩阵

m =8   1   63   5   74   9   2octave:145> pinv(m)
ans =0.147222  -0.144444   0.063889-0.061111   0.022222   0.105556-0.019444   0.188889  -0.102778

绘图

plot(A,B)A是横坐标向量，B纵坐标

octave:173> t=[-1:0.01:1];
octave:174> p1 = sin(2*pi*t);
octave:175> plot(t,p1)

octave:176> p2 = cos(2*pi*t);
octave:177> plot(t,p2)

画完一张图以后想要画另一张就会重新绘制，用hold on可以保留，在原图基础上继续画，并且会自动给你切换颜色。

octave:177> plot(t,p2)
octave:178> hold on
octave:179> plot(t,p1)

给横轴纵轴添加标签。

octave:183> xlabel('time')
octave:184> ylabel('speed')

添加图例

legend('cos','sin','y')

添加标题

title('my demo')

保存文件

print -dpng '文件名'

保存之后你的当前目录下边就会多出一个图片。如过你不知道当前目录是啥，输入pwd查看当前目录的路径。

如果你想同时生成多个图片，那在前边加上figure(n);则会按顺序生成图片。

octave:194> figure(1);plot(t,p1)
octave:195> figure(2);plot(t,p2)

合并显示图像subplot(n,m,index)，将图像分成n行m列，取index位置绘图。

octave:197> subplot(2,3,4)
octave:198> plot(t,p1)
octave:199> subplot(2,3,2)
octave:200> plot(t,p2)

改变纵横轴取值范围axis([横轴起始终止纵轴起始终止])

octave:203> subplot(1,2,1)
octave:204> plot(t,p2)
octave:205> subplot(1,2,2)
octave:206> plot(t,p2)
octave:207> axis([0 1 -2 2])

下图两个是同一副图，右边是改变坐标轴范围之后的图像。

clf清空整幅图，close不用点右上角×关闭图像

矩阵可视化
imagesc()

a =8   1   63   5   74   9   2octave:217> imagesc (a)

说个题外话，作为一个曾经的生物学生，我好像突然get到了我们以前的热图是怎么出来的。

octave:218> load 6feature.txt   % 还记得这组数据吗，就是我上边给的那两组数据的第一组。
octave:219> imagesc (X6feature)

画出来以后图像长这样，恩我真的get到热图怎么画的了。

现在让我们生成一个热图（假的，不是热图，看起来像而已）

octave:220> r = rand(5,50);
octave:221> imagesc (r)

colorbar添加比色的图例。

octave:223> colorbar

colormap gray使其变为灰度图

 colormap gray

基础语法

循环
```
for var = expression
body
endfor
```
```
while (condition)
body
endwhile
```
- 其实上边的endwhile，endfor还有下边的endif可以直接写为end
- break、continue也可以用

判断

if (condition) then-body
endif

if (condition) then-body
elseif (condition) elseif-body
else else-body
endif

举个栗子

v =8    4    8    2    57    4    6    1    310    6    6    4   106    6    8    8   109    1   10    5    7> for i = 1:5
>     for j = 1:5
>         if v(i,j)<5
>             continue
>         else
>             disp(v(i,j))
>         endif
>     endfor
> endfor
8
8
5
7
6
10
6
6
10
6
6
8
8
10
9
10
5
7

函数

普无返回值的函数
```
function name body
endfunctio
```

带返回值的函数

function [ret-list] = name (arg-list) body
endfunction

在octave中写函数：

>function say_hi(name)
> str = ['hello' name];
> disp(str)
> endfunction
>say_hi ('Sian')
helloSian

让octave使用外部的函数。
- 在当前目录或者path中放你的函数文件
- 文件名就是函数名
- 后缀.m
举个栗子
现在写一个函数show_matrix：

存放在path中，并且命名为show_matrix.m

回到octave中声明一个矩阵，看看是否能用该函数成功打印
再整一个带返回值的

Octave简明教程相关推荐

CGIC简明教程(转摘)
CGIC简明教程本系列的目的是演示如何使用C语言的CGI库"CGIC"完成Web开发的各种要求. ********************************* 基础 ...
kangle web server源代码安装简明教程
kangle web server源代码安装简明教程 - kangle使用交流 - kangle软件是一款高性能web服务器,反向代理服务器,提供虚拟主机管理系统及代理服务器,web服务器架设 - ...
CentOs6.5中安装和配置vsftp简明教程
这篇文章主要介绍了CentOs6.5中安装和配置vsftp简明教程,需要的朋友可以参考下一.vsftp安装篇复制代码代码如下: # 安装vsftpd yum -y install vsftpd # ...
sqlalchemy mysql_SQLAlchemy简明教程
原文可见:SQLAlchemy简明教程 - Jiajun的编程随想 SQLAlchemy是Python中常用的一个ORM,SQLAlchemy分成三部分: ORM,就是我们用类来表示数据库schema ...
python tcp server_python scoket 编程 | tcp server client - 简明教程
TCP 和 UDP 的区别这两个协议都是传输层的协议,解决的问题,都是端口与端口的通信问题. TCP 每次建立通信,都需要三次握手,确定双方状态完毕,在发送数据.如果发送的数据出现了异常,TCP 也 ...
简明python教程pdf-python简明教程中文pdf
python简明教程中文pdf电子书是作者通过自己对计算机语言多年来的研究和学习得出的学习经验,对于python新手来说非常有用,值得大家下载学习. python简明教程中文pdf文章目录 1.介绍 ...
简明python教程在线-Python简明教程
Python简明教程在线阅读地址: https://bop.molun.net/ DocStrings 该文档字符串所约定的是一串多行字符串,其中第一行以某一大写字母开始,以句号结束.第二行为空行,后 ...
linux 防火墙iptables简明教程
前几天微魔部落再次遭受到个别别有用心的攻击者的攻击,顺便给自己充个电,复习了一下linux下常见的防火墙iptables的一些内容,但是无奈网上的很多教程都较为繁琐,本着简明化学习的目的,微魔为大家剔 ...
Python re(正则表达式）简明教程
Python re(正则表达式)简明教程(Python3) 举个栗子常用函数举例按字母顺序查询模块

Octave简明教程