Octave常用語法

大體內容是轉的,但是繪圖部分做了補充和修改。

向量


  • 分號:分割行
  • 空格或逗號:分割列

創建與訪問

行向量

空格或逗號分割

>> v = [1 2 3] % 等效:v = [1, 2, 3]
v =
   1   2   3
>> v(2) % 只有一行,所以指定就是列
ans =  2

列向量

  • 分號分割
>> v = [1; 2; 3]
v =
   1
   2
   3
>> v(2) % 只有一列,所以指定就是行
ans =  2

矩陣


  • 與向量一樣,空格或逗號分割列,分號分割行
  • 冒號:代表所有,所有行或所有列

創建

>>   A = [1, 2; 3, 4] % 一行寫完
A =
   1   2
   3   4
>> A = [1, 2; % 分行寫
> 3, 4]
A =
   1   2
   3   4

創建技巧

  • a:c --- 從a到c
  • a:b:c --- 從a,間隔b,到c
  • 向量也可以用這樣的方法創建
>> A = [1:3; 4:6]
A =
   1   2   3
   4   5   6
>> A = [1:2:5; 2:2:6]
A =
   1   3   5
   2   4   6

特殊矩陣


單位矩陣

>> eye(3) % 單位矩陣
ans =
Diagonal Matrix
   1   0   0
   0   1   0
   0   0   1
>> flipud(eye(3))
ans =
Permutation Matrix
   0   0   1
   0   1   0
   1   0   0

轉置矩陣

A =
   1   2
   3   4
>> A'
ans =
   1   3
   2   4

逆矩陣

A =
   1   2
   3   4
>> pinv(A)
ans =
  -2.00000   1.00000
   1.50000  -0.50000
>> pinv(A) * A
ans =
   1.00000   0.00000
  -0.00000   1.00000

全1矩陣

>> ones(2, 3)
ans =
   1   1   1
   1   1   1

全0矩陣

>> 0*ones(2, 3)
ans =
   0   0   0
   0   0   0
>> zeros(2, 3)
ans =
   0   0   0
   0   0   0

隨機矩陣

>> rand(1, 3) % 0~1的隨機數,1行3列
ans =
   0.99291   0.65946   0.95102

高斯分布矩陣

>> randn(1, 3) % 高斯分布
ans =
   0.14646   2.02587   1.33266

訪問


  • 冒號:代表所有,所有行或所有列

訪問元素

A =
   1   2   3
   4   5   6
   7   8   9
>> A(2, 2) % 第二行,第二列,先定行再定列
ans =  5

訪問單行或單列

>> A(1, :) % 第一行,所有,元素
ans =
   1   2   3
>> A(:, 2) % 第二列,所有,元素
ans =
   2
   5
   8

訪問多行或多列

>> A(:, [1, 3]) % 第一列和第三列,所有,元素
ans =

   1   3
   4   6
   7   9

>> A([1, 3], :) % 第一行和第三行,所有,元素
ans =

   1   2   3
   7   8   9

連接

  • C = [A B],把B按列加到A上,生成C
  • C = [A; B],把B按行加到A上,生成C

添加元素

>> A = [1]
A =  1
>> A = [A, 2] % 把2,按列,添加到A上,再賦給A
A =
   1   2
>> A = [A; 3] % A有兩列,3只有一列,規模不匹配
error: vertical dimensions mismatch (1x2 vs 1x1)

添加行或列

>> A = [A; [3, 4]] % 將向量[3, 4],作為行,添加到A上,在賦給A
A =
   1   2
   3   4
>> A = [A, [5; 6]] % 將向量[5, 6],作為列,添加到A上,在賦給A
A =
   1   2   5
   3   4   6

矩陣連接

>> A = [1, 2; 3, 4]
A =
   1   2
   3   4
>> B = [5, 6; 7 ,8]
B =
   5   6
   7   8
>> [A; B] % 將B,作為行,添加到A上
ans =
   1   2
   3   4
   5   6
   7   8
所有的數據放入一個向量中
>> A
A =
   1   2
   3   4
>> A(:)
ans =
   1
   3
   2
   4
>> A(:)'
ans =
   1   3   2   4

賦值


  • 在訪問基礎上,給定相同規模的數據
A =
   1   2   3
   4   5   6
   7   8   9
>> A(3, 3) = 10 % 修改單個元素的值
A =
    1    2    3
    4    5    6
    7    8   10
>> A(1, :) = [0, 0, 0] % 修改一行的值
A =
    0    0    0
    4    5    6
    7    8   10
>> A(2:3, 2:3)
ans =
    5    6
    8   10
>> A(2:3, 2:3) = [0, 0; 0, 0] % 修改指定矩陣的值
A =
   0   0   0
   4   0   0
   7   0   0

運算


加,減

>> A = [1 1; 1 1]
A =
   1   1
   1   1
>> B = [2 2; 2 2]
B =
   2   2
   2   2
>> A + B
ans =
   3   3
   3   3
>> A - B
ans =
  -1  -1
  -1  -1
>> A - 1
ans =
  0  0
  0  0

相乘

>> A = [1 2; 3 4]
A =

   1   2
   3   4

>> B = [5 6; 7 8]
B =

   5   6
   7   8

>> A * B
ans =

   19   22
   43   50

>> -A % -1 * A
ans =

  -1  -2
  -3  -4

點運算


對應元素運算

  • 維度相同:對應元素相乘
  • 行維度相同:每行對應元素相乘
  • 列維度相同:每列對應元素相乘
  • A .* B = B .* A
  • 乘 *;除 /; 平方 ^;等

維度相同

A =
   1   1
   1   1
B =
   2   3
   2   3
>> A .* B
ans =
   2   3
   2   3

行維度相同

A =
   1   1
   1   1
B =
   5   6
>> A .* B
ans =
   5   6
   5   6

列維度相同

A =
   1   1
   1   1
B =
   5
   6
>> A .* B
ans =
   5   5
   6   6

點除

  • 矩陣乘以常數,A * 2,除可以是,A / 2
  • 反過來,2 * A沒問題,2 / A 就不行,要用2 ./ A
A =
   1   2
   3   4
>> 1 ./ A
ans =
   1.00000   0.50000
   0.33333   0.25000

邏輯


  • 每一個元素做比較,標記0或1
  • &gt, <, ==, !=(或~=), &&, ||
A =
   1   2
   3   4
>> A > 2
ans =
   0   0
   1   1

位運算


  • 或 | ,函數xor
  • 且 &
  • 非 ~
  • 異或 ^

控制語句


if

i =  1
>> if i == 1
>     disp(1)
>  elseif i == 2
>     disp(2)
>  else
>     disp(3)
>  end

for

>> for i = 1:3 % 從1到3
>     disp(i)
>  end
 1
 2
 3

while

>> while i <= 3
>     disp(i)
>     i = i + 1
>  end
 1
 2
 3

break,continue

  • 與C、C++、Java沒有區別

函數


  • size:獲取矩陣維度
  • length:獲取最大維度
  • who:變量列表
  • whos:變量詳情
  • clear 變量名:刪除指定變量
  • clear:刪除所有變量
  • find:返回符合條件元素的下標
  • log:log以e為底
  • exp:e的多少次方
  • abs:絕對值
  • floor:向下取整
  • ceil:向上取整
  • sum:求和
  • prop:求積

size

  • 獲取矩陣維度
A =

   1   2
   3   4
   5   6
>> aSize = size(A)
aSize =
   3   2
>> size(A, 1) % 3行
ans =  3
>> size(A, 2) % 2列
ans =  2

length

  • 獲取最大維度
A =
   1   2
   3   4
   5   6
>> length(A) % 輸出最大維度
ans =  3

who

變量列表

>> who % 現在都有哪里變量
Variables in the current scope:

A          aSize      ans        featuresX  w

whos

變量詳情

>> whos % 變量詳情
Variables in the current scope:

   Attr Name           Size                     Bytes  Class
   ==== ====           ====                     =====  =====
        A              3x2                         48  double
        aSize          1x2                         16  double
        ans            1x20                        20  char
        featuresX     27x2                        432  double
        w              1x10000                  80000  double

Total is 10082 elements using 80516 bytes

clear

刪除變量

>> clear featuresX
>> who
Variables in the current scope:

A      aSize  ans    w
>> clear % 清除全部變量
>> who % 一個變量都沒有了

find

返回符合條件元素的下標

A =

   5   6
   7   8

>> [r, c] = find(A > 6) % 符合要求的是:第二行第一個和第二行第二個
r =

   2
   2

c =

   1
   2

sum

求和

>> A = [1, 2; 3, 4]
A =

   1   2
   3   4

>> sum(A) % 等效:sum(A, 1)
ans =

   4   6
   
>> sum(A, 2)
ans =

   3
   7

prod
求積

>> A
A =

   1   2
   3   4

>> prod(A) % 等效:prod(A, 1)
ans =

   3   8

>> prod(A, 2)
ans =

    2
   12

max

  • 向量
>> a = [1 5 2 3.3]
a =
   1.0000   5.0000   2.0000   3.3000
>> max(a)
ans =  5
>> [val, ind] = max(a)
val =  5
ind =  2
  • 矩陣比較
    兩個矩陣比較每個元素,保留大的
>> A = rand(3)
A =
   0.2620788   0.6346345   0.4659161
   0.0880455   0.1258945   0.0079559
   0.0296765   0.7917592   0.4321800
>> B = rand(3)
B =
   0.039237   0.672424   0.214649
   0.491320   0.362929   0.197626
   0.821090   0.675265   0.698960
>> max(A, B)
ans =
   0.26208   0.67242   0.46592
   0.49132   0.36293   0.19763
   0.82109   0.79176   0.69896
  • 行與列的最大值
    參數二:與之比較的矩陣參數三:按行還是按列
>> A = [1 2; 3 4]
A =
   1   2
   3   4
>> max(A) % 每列最大值
ans =
   3   4
>> max(A, [], 1) % 每列最大值
ans =
   3   4
>> max(A, [], 2) % 每行最大值
ans =
   2
   4

randperm

生成亂序序列

>> A = [2, 3, 4 ,5 ,6]
A =

   2   3   4   5   6
>> rand_indices = randperm(length(A))
rand_indices =

   1   2   4   5   3
>> A(:, rand_indices(1:3))
ans =

   5   6   2

自定義函數


創建步驟


  • 創建存放函數的文件(文件名要與函數名相同,.m結尾)
square.m
  • 創建函數
% 一個返回值
function 返回值 = 函數名(參數列表)
    函數體
end
% 多個返回值
function [返回值1, 返回值2] = 函數名(參數列表)
    函數體
end
function y = square(x)
  y = x^2;
end
function [y1, y2] = squareAndCube(x)
  y1 = x^2;
  y2 = x^3;
end

調用


進入到存放函數文件的目錄
添加存放函數文件的目錄,到,搜索路徑

>> cd Downloads/ % 進入存放函數文件的目錄
>> square(2)
ans =  4

>> addpath('~/Downloads') % 添加到搜索路徑
>> cd ..
>> square(2)
ans =  4

數據

路徑

與Linux命令一樣

>> pwd % 當前路徑
ans = /Users/XXX

>> cd ~/Developer

>> pwd
ans = /Users/XXX/Developer

>> ls % 列出目錄文件和文件夾

>> cd .. % 回到上級目錄

載入與保存

>> load featuresX.dat % 載入數據,變量名為featuresX
>> data = load('ex1data1.txt'); % 變量名為data

>> save featuresX.mat featuresX; % 將featuresX中的數據保存到featuresX文件中
>> save hello.txt v -ascii % 將變量v中得數據已ascii編碼方式保存到hello.txt中

其他


命令

  • 簡化命令行:PS1('>> ');
  • 關閉圖表:close或close all
  • 命令行清屏:clc

注釋

%后為注釋

打印


自動打印

a = 3 % 會打印數據
a = 3; % 不會打印數據,分號可以防止輸出

格式化

disp(pi) % 輸出:3.1416
disp(sprintf('pi is %.2f', pi)) % c語言風格

默認格式

format long % 修改默認打印格式

幫助


help size % 查看size函數的幫助文檔
help help % 查看help函數的幫助文檔

繪圖


plot

屬性

  • LineWidth:線寬
  • MarkerFaceColor:標記顏色
  • MarkerSize:標記大小
t = [0:0.1:0.98];
y1 = sin(2*pi*4*t);
plot(t,y1);
y2 = cos(2*pi*4t);

hold on ;  %  使用 hold on 可以將接下來繪制的函數放在一張圖里

plot(t, y2, 'r'); % 第三個參數,表示顏色

xlabel('time');
ylabel('value');  %  標記出x和y軸

legend('sin','cos') %  在圖像右上角,用來區分多個函數的圖例,依次表示兩個曲線所表示的內容

tiitle('my plot') %  在圖像的正上方給圖像加上標題

print -dpng 'myPlot.png'  % 輸出圖像,保存在'pwd'目錄下

close % 關閉圖像

figure(1); plot(t,y1); % 新開第一個窗口繪制y1
figure(2); plot(t,y1); % 新開第二個窗口繪制y2

% pos為向量中為1的行序號組成的向量,neg為向量中為0的行序號組成的向量 
plot(X(pos,1),X(pos,2),'k+','LineWidth',2,'MarkerSize',7); 
plot(X(neg,1),X(neg,2),'ko','markerFaceColor','y','MarkerSize',7);
Paste_Image.png
 Format arguments:
    linestyle
          '-'  Use solid lines (default).
          '--' Use dashed lines.`
          ':'  Use dotted lines.
          '-.' Use dash-dotted lines.
    markerstyle
          '+'  crosshair
          'o'  circle
          '*'  star
          '.'  point
          'x'  cross
          's'  square
          'd'  diamond
          '^'  upward-facing triangle
          'v'  downward-facing triangle
          '>'  right-facing triangle
          '<'  left-facing triangle
          'p'  pentagram
          'h'  hexagram
    color
          'k'  blacK
          'r'  Red
          'g'  Green
          'b'  Blue
          'm'  Magenta
          'c'  Cyan
          'w'  White 

常用的修改屬性有: 
 "linestyle", 
 "linewidth",
 "color",
 "marker",
 "markersize",
 "markeredgecolor",
 "markerfacecolor".

直方圖


>> w = -6 + sqrt(10)*randn(1, 10000);
>> hist(w)
>> hist(w, 30) % 30組
最后編輯于
?著作權歸作者所有,轉載或內容合作請聯系作者
平臺聲明:文章內容(如有圖片或視頻亦包括在內)由作者上傳并發布,文章內容僅代表作者本人觀點,簡書系信息發布平臺,僅提供信息存儲服務。
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市,隨后出現的幾起案子,更是在濱河造成了極大的恐慌,老刑警劉巖,帶你破解...
    沈念sama閱讀 229,619評論 6 539
  • 序言:濱河連續發生了三起死亡事件,死亡現場離奇詭異,居然都是意外死亡,警方通過查閱死者的電腦和手機,發現死者居然都...
    沈念sama閱讀 99,155評論 3 425
  • 文/潘曉璐 我一進店門,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人,你說我怎么就攤上這事。” “怎么了?”我有些...
    開封第一講書人閱讀 177,635評論 0 382
  • 文/不壞的土叔 我叫張陵,是天一觀的道長。 經常有香客問我,道長,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 63,539評論 1 316
  • 正文 為了忘掉前任,我火速辦了婚禮,結果婚禮上,老公的妹妹穿的比我還像新娘。我一直安慰自己,他們只是感情好,可當我...
    茶點故事閱讀 72,255評論 6 410
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著,像睡著了一般。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發上,一...
    開封第一講書人閱讀 55,646評論 1 326
  • 那天,我揣著相機與錄音,去河邊找鬼。 笑死,一個胖子當著我的面吹牛,可吹牛的內容都是我干的。 我是一名探鬼主播,決...
    沈念sama閱讀 43,655評論 3 444
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了?” 一聲冷哼從身側響起,我...
    開封第一講書人閱讀 42,838評論 0 289
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后,有當地人在樹林里發現了一具尸體,經...
    沈念sama閱讀 49,399評論 1 335
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內容為張勛視角 年9月15日...
    茶點故事閱讀 41,146評論 3 356
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發現自己被綠了。 大學時的朋友給我發了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 43,338評論 1 372
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖,靈堂內的尸體忽然破棺而出,到底是詐尸還是另有隱情,我是刑警寧澤,帶...
    沈念sama閱讀 38,893評論 5 363
  • 正文 年R本政府宣布,位于F島的核電站,受9級特大地震影響,放射性物質發生泄漏。R本人自食惡果不足惜,卻給世界環境...
    茶點故事閱讀 44,565評論 3 348
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧,春花似錦、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 34,983評論 0 28
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至,卻和暖如春,著一層夾襖步出監牢的瞬間,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 36,257評論 1 292
  • 我被黑心中介騙來泰國打工, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人。 一個月前我還...
    沈念sama閱讀 52,059評論 3 397
  • 正文 我出身青樓,卻偏偏與公主長得像,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子,可洞房花燭夜當晚...
    茶點故事閱讀 48,296評論 2 376

推薦閱讀更多精彩內容