CH2 基本運算tobemodified.docx

Matlab_2 基本運算u 變量和數(shù)據(jù)類型u 整型和浮點型數(shù)據(jù)u 數(shù)組和矩陣運算u 字符串、單元體和結(jié)構(gòu)體u 關(guān)系和邏輯運算u 位運算、時間運算 變量命名規(guī)則1) 以字母開頭，可以包含63個字母、數(shù)字和下劃線（字母區(qū)分大小寫）。2) 不能使用Matlab關(guān)鍵字。3) 避免與已經(jīng)定義的變量或函數(shù)同名。相關(guān)的函數(shù)：namelengthmax % 變量的最大長度 isvarname % 檢查命名合法性iskeyword % 查看關(guān)鍵字（熟記?。〆xist % 檢測命名是否已存在which % 查找命名的位置Matlab定義的常用變量：ans, pi, Inf, NaN, i, j, 容易被用作變量的Matlab函數(shù)/變量： i, j, size, length, any, all, % 用戶可以重置Matlab已定義的變量或函數(shù)名pi = 2.718 % 雖然沒有必要，但是合法10*pi變量類型根據(jù)變量的生存時間和可訪問特性，Matlab定義了三種類型的變量：1) local 局部變量（默認）2) global 全局變量3) persistent 永久變量（僅用于函數(shù)中）Matlab搜索路徑搜索路徑指的是Matlab所有文件所在目錄的列表。當Matlab讀取用戶輸入的變量名（或M文件名），Matlab按以下順序查找：1） 是不是Matlab工作區(qū)的變量，如果不是，執(zhí)行下一步。2） 是不是Matlab內(nèi)置函數(shù)，3） 是不是當前目錄下的M文件，4） 按Matlab搜索路徑順序查找，5） 沒有找到，提示錯誤信息。更改搜索路徑的方法：1) 用戶可以通過path組函數(shù)調(diào)用，對搜索路徑內(nèi)容進行添加、刪除、排序等操作。2) 更簡便的方法是在Matlab桌面File菜單Set Path選項中進行操作。3) 也可以在Current Folder窗口，按在線提示直接設(shè)置。% exist查看命名是否已經(jīng)存在exist % var, func, mfile% 使用幫助文檔是掌握Matlab的必要環(huán)節(jié)！help exist % 查看返回值的意義doc exist % 查看幫助文檔% which查看變量或文件位置which sinwhich -all sinhelp which % 查看更多的用法doc which基本數(shù)據(jù)類型Matlab定義了15種基本數(shù)據(jù)類型1) 默認的數(shù)據(jù)類型為double（雙精度浮點型）；2) 除了function_handle（函數(shù)句柄），所有數(shù)據(jù)類型均支持數(shù)組運算；3) cell（單元體）和struct（結(jié)構(gòu)體）為承載其它數(shù)據(jù)的容器；4) 函數(shù)class和whos查看數(shù)據(jù)類型。整數(shù)型數(shù)值數(shù)據(jù)1) 整型數(shù)據(jù)必須用相應(yīng)的函數(shù)說明int8, int16, int32, int64uint8, uint16, uint32, uint642) intmax, intmin查看數(shù)值范圍3) 注意整型數(shù)的運算規(guī)則% 整型數(shù)數(shù)據(jù)范圍intmax(int8) intmin(int8)intmin(uint8) % 整型數(shù)的運算a = int8(pi)b = int16(3.5)class(a) % 查看變量類型class(b)whos % 查看當前工作區(qū)內(nèi)所有的變量c = a*b % 不同類型整數(shù)不能直接運算d = a*pi % 注意整數(shù)和浮點數(shù)的運算結(jié)果whos浮點型數(shù)值數(shù)據(jù)1) 單、雙精度浮點數(shù)分別用single、double函數(shù)說明2) realmax, realmin查看數(shù)值范圍3) format改變數(shù)據(jù)顯示格式4) eps顯示數(shù)據(jù)精度% 浮點數(shù)據(jù)范圍realmax singlerealmin singlerealmaxrealmin % 數(shù)據(jù)顯示格式piformat longpi help format% eps：一個浮點數(shù)能夠取得的最小增量。eps eps(single) eps(double)eps(1)eps(10)eps(1e10) % 浮點數(shù)的計算精度sin(pi)tan(pi/2)exp(pi*i)% 大數(shù)和小數(shù)的運算a = 1 + 1e-16b = 1a = b % 關(guān)系運算符=判別兩個操作數(shù)是否相等% 運算次序format longa = 0.1+0.4-0.3b = 0.1-0.3+0.4if(a=b) disp( okay. ) % disp顯示字符串內(nèi)容else disp( this is a bug! )enda b% 避免用運算后浮點數(shù)結(jié)果是否相等作為判別條件。% 一般浮點數(shù)都不能準確表達。% 整數(shù)可以用浮點格式準確表達。% 用format hex查看整數(shù)的存儲% 浮點數(shù)代替整數(shù)的運算% flintmax 浮點數(shù)能表達的連續(xù)整數(shù)的最大值253% bitmax 浮點數(shù)能表達的最大整數(shù)253-1（新版本中已刪除）% fix 取整數(shù)部分% round 按舍入原則取整% rem 計算余數(shù)%format longm = 253n = flintmaxm = na = 1.5, b = 1.4999fix(a)round(a)rem(10, 3)ream(1, 0.3)% 驗證Fermat數(shù)： F(n) = 2(2n) + 1, n = 0,1,2% F(5)不是質(zhì)數(shù)n = 5N = 2(2n) + 1for i = 2:sqrt(N)if( rem(N, i)=0 ) disp(A factor is found!) break;end;end復(fù)數(shù)運算復(fù)數(shù)可以直接用i, j生成，也可以用complex函數(shù)常用函數(shù)：real, imag, angle, abs, isreal幾乎所有的Matlab運算符和數(shù)學函數(shù)都支持復(fù)數(shù)運算對復(fù)數(shù)矩陣，操作符 表示共軛轉(zhuǎn)置運算%a = 1b = 2z = 1 + 2iz = a + b*1j % 寫成*1j避免與下標j混淆z = complex(a, b)real(z)imag(z)angle(z)abs(z)isreal(z)isreal(z*z)isreal(z*z.) % 注意操作符.和.的差別isreal( exp(pi*1j) ) % 注意浮點計算精度的影響z22/zz2sqrt(z)exp(z)數(shù)組和數(shù)組運算矩陣（matrix）是長度為m*n的二維數(shù)值數(shù)組。向量（vector）是長度為n的一維數(shù)組，可看做1*n或n*1的矩陣。標量（scalar）指的是一個數(shù)字，可看做1*1的矩陣。向量的創(chuàng)建1) 方括號內(nèi)直接輸入2) 冒號操作符3) linspace和logspace函數(shù)% 直接創(chuàng)建a = 1 2 3 b = 1*2 3-4c = 1+i 2-3i% 注意：% 方括號內(nèi)，表達式最好寫成緊湊格式（不留空格）。% 空格號與逗號的功能相同。%下列表達易產(chǎn)生錯誤d = 1 +2 3- 4 % 加減號可以是一元，也可以是二元操作符e = 1 *i 2/ i% 冒號操作符a = 1:6 % 缺省的步長為1b = 0:0.1:0.5c = -1:3:10% 函數(shù)創(chuàng)建d = linspace(1, 6, 6)e = logspace(1, 3, 3) % first:step:end 指定步長step（缺省為step = 1）% linspace(first, end, n) 指點節(jié)點數(shù)n% logspace(first, end, n) 10first 10end% Matlab數(shù)學函數(shù)支持向量輸入，計算結(jié)果也是以向量返回tan(a) log10(e)一維數(shù)組尋址% 數(shù)組下標用 ( ) 表達a = 1:10a(1), a(3), a(10), a(end)%下標可以是向量a(1:5)a(5:-1:1)a(7:end)a(1 3 5 7)% 下標只能是正整數(shù)a(3.14) % xa(0) % x% 下標不能超出end a(12) % x數(shù)組方向：行向量與列向量Matlab的冒號操作符和linspace創(chuàng)建的均為行向量創(chuàng)建列向量：1） 分號操作符2） 轉(zhuǎn)置操作符（撇號）%a = 1 2 3b = 1; 2; 3 % 分號的換行功能c = a. % 轉(zhuǎn)置操作% 復(fù)數(shù)問題a = 1:3c = a + a*1i c % 表示復(fù)數(shù)共軛轉(zhuǎn)置，對非復(fù)數(shù)與. 功能相同 c. % . 表示非復(fù)數(shù)共軛轉(zhuǎn)置 二維數(shù)組的創(chuàng)建1) 直接創(chuàng)建2) 由向量或數(shù)組合成3) 利用標準數(shù)組% 直接創(chuàng)建a = 1 2 3; 4 5 6 % 分號的功能b = 1 2 3; 4 5 6 7 % xa = 1:4b = 2:5c = a; b % 注意數(shù)組的方向d = a., b.e = c, 3:6 % 標準數(shù)組zeros % 全0數(shù)組ones % 全1數(shù)組rand % （01）隨機數(shù)數(shù)組randn % 正態(tài)分布（0,1）隨機數(shù)組eye % 單位矩陣diag % 對角矩陣magic % 魔方矩陣% 注：前4個函數(shù)可以用于多維數(shù)組% Matlab中函數(shù)的輸入、輸出是可變的eye(5)eye(3, 4)eye(3, 4, 2) % x%diag( 1 3 5 7 )diag( 1:2:7 )%magic(3)magic(5)數(shù)組運算標量-數(shù)組運算 加減乘除冪等運算是對數(shù)組中所有元素進行。Matlab通過標量擴展，將標量轉(zhuǎn)換為相應(yīng)大小的數(shù)組。a = magic(3)b = a-3c = (a-3)*2d = 2a數(shù)組-數(shù)組運算 大小相同的數(shù)組可以進行元素-元素的運算： 符號 .* ./ . . 數(shù)組-數(shù)組運算遵循矩陣運算規(guī)則 符號 * / 1./a % 標量擴展a-1 % 求逆陣1/a % x數(shù)組操作數(shù)組生成后，可以通過特定的下標，完成提取、賦值、重排等多種功能。% 賦值a = ones(3)a(2,2) = 9a(:, 1) = 0 0 0.a(:, 1) = 0 % 標量擴展% 賦值語句中，當下標超出范圍時，Matlab將自動擴展，未定義的位置填充0x(2, 3) = 10 % 數(shù)組x是未創(chuàng)建的數(shù)組x(5, 5) = 100 % 合法的% !不推薦使用此類賦值方法% 提取b = a(1:2, 2:3) % 提取數(shù)組元素時，下標超出范圍不合法c = a(10, 10) % x%多維數(shù)組的單下標尋址（線性索引）d = magic(3)d(1), d(6)d(:) % Matlab多維數(shù)組的存儲：列優(yōu)先原則d(3:7)d(10) % x索引號不能超出總元素數(shù)d(10) = 0 % x 數(shù)組形狀不能改變d(end, 1) % 注意end的取值d(end)% 了解函數(shù)sub2ind, ind2sub的用法% 查找：函數(shù)find的使用a = magic(4)ind = find(a = 5)r, c = find(a=5) % 賦值號左邊 的用法% 多個返回值b = magic(4)ind = find(b3)r, c = find(b3)doc find %了解find更多的用法% 排序：函數(shù)sort的使用s = randperm(6) % 生成1n隨機排列的數(shù)組sort(s) % 元素排序，缺省為ascendsort(s, ascend)sort(s,descend)% 矩陣或多維數(shù)組t = magic(4)sort(t) % 缺省為列優(yōu)先、升序sort(t,1)sort(t,2)sort(t, 2, descend)st, index = sort(t, 2, descend) % 熟悉函數(shù)sum, mean, max, min的用法% 了解函數(shù)reshape和repmat 的用法數(shù)組大小函數(shù)size, length, numel, ndims%a = 1:4; 5:8size(a) % 返回一個行向量 行數(shù), 列數(shù), 頁數(shù) length(a) % 返回數(shù)組的最大長度max(size(a) numel(a) % 返回數(shù)組的元素總數(shù) ndims(a) % 數(shù)組維數(shù)length(size(a)向量的乘積-內(nèi)積、外積、混合積Matlab定義了函數(shù)dot和cross分別計算向量的內(nèi)積和外積向量u, v, w的混合積的定義為u.(v*w)% 注：函數(shù)cross只能對長度為3的向量進行運算%u = 1 2 3v = 1 3 5w = 2 1 1% 直接按定義計算dot(u, cross(v, w)% 另一方法det ( u; v; w ) % 函數(shù)det計算方陣的行列式det ( u v w )線性方程組-矩陣的乘法和除法線性方程組a*x = b的解為x = ab 或 a-1*b%a = 1 2 3; 4 5 6; 7 8 9b = 1 2 3x = aba*x - bx = inv(a)*b % inv計算逆陣, inv(a)等同于a-1a*x - b% 問題的分析inv(a)*a % 線性方程組有唯一解的條件？rank(a), rank( a,b) % rank計算矩陣的秩%a(3,3) = 0 rank(a), rank(a,b) x = abx = inv(a)*binv(a)*a 多維數(shù)組格式：Array(row, column, page, )第三維稱為頁，更高的維數(shù)無統(tǒng)一的名稱創(chuàng)建方法、元素運算與二維類似乘除運算無定義%a(:, :,1) = ones(3);a(:,:,2) =eye(3);a % 分頁顯示% 函數(shù)zeros, ones, rand, randn可以創(chuàng)建多維數(shù)組randn(2, 3, 4) ones(2, 3, 3, 3) % 多維數(shù)組的運算次序：列-行-頁b(:,:,1) = diag(1 10 100);b(:,:,2) = magic(3);bsum(b)sum(b, 1) % 1/2/3分別表示對列/行/頁進行求和sum(b, 2)sum(b, 3)sum(sum(sum(b)sum( b(:) )稀疏矩陣當矩陣的非0元素小于總元素數(shù)的1%，稱為稀疏矩陣。存貯大量的0元素消耗內(nèi)存資源，對0元素的運算也浪費大量機器時間。%a = eye(10)s = sparse(a) % 函數(shù)sparse將普通矩陣轉(zhuǎn)換為稀疏矩陣whoss(1,5) = 100 % 稀疏矩陣的操作與普通矩陣相同size(s), length(s), numel(s) c = s(2,2)b = full(s) % 函數(shù)full將稀疏矩陣轉(zhuǎn)換為普通矩陣d = full(c) % 以上僅僅是說明sparse用法，創(chuàng)建稀疏矩陣最好用speye, spdiags等函數(shù)s = speye(1000);a = eye(1000);whos% 創(chuàng)建稀疏數(shù)組必須用Matlab指定的函數(shù)aa(1000, 1000) = 1; % 盡管非0元素極少，aa仍然是普通矩陣whosss = sparse(aa)whos字符串字符串是以字符顯示顯示的ASCII數(shù)值型數(shù)組。Matlab中，字符串是用單引號括起來的一串字符；每個字符稱為一個元素，通常用兩個字節(jié)存儲。%s = Hellow, Matlab.d = double(s)char(d)% 內(nèi)容提取s(2)s(1:6)s(9:14) % 字符串中還有單引號，必須寫成兩個單引號It doesnt work unless% 字符串合成a = Hellow, b = Matlab.a, b % 函數(shù)disp顯示字符串內(nèi)容，不顯示變量名disp(s)disp(Very good)% 漢字和英文字符的處理方法相同h = 漢字處理d = double(h)char(d) whos字符串數(shù)組多行字符構(gòu)成字符串數(shù)組，單每行的元素必須相同。字符數(shù)少的串后面填充空格。% 直接輸入s = one two three four size(s) % 用函數(shù)char可以避免輸入空格t = char(one, two, three, four) size(t)zh = char( 中國, 安徽省合肥市, 中國科學技術(shù)大學, 230026) size(zh)% 內(nèi)容提取t(2,1)t(1,4)t(:, 3)t(3, 1:3) t(12) % 線性索引% 賦值和替換t(4,1) = F t(1,:) = ONE % x，字符數(shù)目必須相等t(:, 3) = C % 標量擴展格式化輸入輸出和格式轉(zhuǎn)換% 函數(shù)disp回顯字符串或數(shù)組，無返回參數(shù)disp(Matlab is matrix language.);disp( magic(3) );% 函數(shù)sprintf將數(shù)據(jù)格式化寫入到字符串str = sprintf(Equation:n%d x %.5f = %.2e, 1000, pi, 1000*pi);disp(str);% sscanf從字符串中讀取格式化數(shù)據(jù)% % 函數(shù)num2str/str2num數(shù)字/字符串轉(zhuǎn)換a = 3.14s = num2str(a)str2num(s)a = s% 函數(shù)dec2bin/bin2dec無符號十進制數(shù)/二進制數(shù)顯示的轉(zhuǎn)換a = 14bits = dec2bin(a,8) % 整數(shù)顯示為二進制8位字符串bin2dec(bits)bin2dec( 000111 ) % 二進制字符串轉(zhuǎn)換為整數(shù)% 例：判別一個十進制整數(shù)中是否有某個數(shù)字a = 31*137s = num2str(a)aa = str2num(s.)find(aa = 3) %單元數(shù)組和結(jié)構(gòu)體單元數(shù)組（cell array）和結(jié)構(gòu)體（structure）將不同類型的數(shù)據(jù)集成到一個變量中，方便相關(guān)數(shù)據(jù)的處理。單元數(shù)組和結(jié)構(gòu)僅僅是承載其它數(shù)據(jù)類型的容器，具體的運算主要針對其包含的數(shù)據(jù)。單元數(shù)組和結(jié)構(gòu)體的差別在于標識的方式：cell array通過數(shù)字（下標）標識structure通過名稱（字段）標識單元數(shù)組的創(chuàng)建1） 直接用賦值語句創(chuàng)建2） 函數(shù)cell創(chuàng)建% 直接創(chuàng)建的兩種格式：% （1）按單元索引，待賦給單元的值必須用a(1,1) = 1:10 ;a(1,2) = 1+i ;a(2,1) = this is a test. ;a(2,2) = magic(3) % （2）按內(nèi)容索引，將寫在等式左邊b1,1 = 1:10;b1,2 = 1+i;b2,1 = this is a test.;b2,2 = magic(3) % 兩種方法完全等效% 用于訪問單元的值，( )用于標識單元% 單元體比較簡單時，可一步直接寫出c = 1:10, 1+i; this is a test., magic(3) % 比較a, b, c是否等同isequal(a, b, c) % 所有收入?yún)?shù)均相等返回1% 單元的內(nèi)容很長，可能無法回顯a(1,1)a1,1a(1,2)a1,2% 回顯單元體全部內(nèi)容用celldisp函數(shù)，或用冒號操作符celldisp(a)a:, :% 函數(shù)cellplot圖示單元體結(jié)構(gòu)cellplot(a)cell 函數(shù)a = cell(3,4) % 生成3x4的空單元組數(shù)a(1,1) = this doesnt work % xa(1,1) = this does work 單元數(shù)組的處理% % 注意與普通的數(shù)值型數(shù)組的差異% 和( )的差異a = 1:10, 1+i; this is a test., magic(3) size(a)numel(a) size(a(1,1)size(a1,1) class(a(1,1)class(a1,1) % （）不能用于處理單元的內(nèi)容% 內(nèi)容提取b = a(1,1) % b是1x1單元數(shù)組c = a1,1 % c是1x10普通數(shù)組% 賦值和替換a(1,1) = eye(4) % xa(1,1) = eye(4) a1,1 = eye(4)字符串單元數(shù)組c = one; two; three; four s = char(one, two, three, four)whos% 函數(shù)cellstr將字符串數(shù)組轉(zhuǎn)換為字符串單元數(shù)組，其功能與char相反cellstr(s)char(c)isequal( cellstr(s), c )isequal( char(c), s)c1,1 = a long words(1,1) = a long word %結(jié)構(gòu)體的創(chuàng)建結(jié)構(gòu)體用字段（fields）來進行索引，用點號訪問字段中的數(shù)據(jù)。創(chuàng)建方法：1） 直接賦值創(chuàng)建2） 函數(shù)struct創(chuàng)建% 創(chuàng)建結(jié)構(gòu)體circlecircle.radius = 2.3;circle.center = 1, 2;circle.linestyle = -% 結(jié)構(gòu)體數(shù)組circle(2).linestyle = -.; %circle(2).radius = 32;circle(2).center = 2, 3% 函數(shù)structa1 = 2.3, 32 ; % 特定的字段值必須寫成單元數(shù)組形式a2 = 1, 2, 2, 3 ;a3 = -, -. ;C = struct( radius, a1, center, a2, linestyle, a3 ) % 注意struct函數(shù)的輸入?yún)?shù)格式% 比較兩種創(chuàng)建方法結(jié)果是否等同isequal( circle, C) 結(jié)構(gòu)體處理% 添加新字段C(1).color = red C.color % 未定義的內(nèi)容填充空數(shù)組C(2).color = blueC.color% 內(nèi)容提取r1 = C(1).radiusr2 = C(2).radius% 同時提取不同元素同字段的值 r1, r2 = C.radius % 注意輸出格式c1, c2 = C.center 關(guān)系和邏輯運算關(guān)系和邏輯運算主要用于真/假判斷，控制程序執(zhí)行。運算結(jié)果為True/False （1/0），占一個字節(jié)。關(guān)系和邏輯運算中，所有非0值視為True。% 函數(shù)logical將數(shù)組轉(zhuǎn)換為邏輯數(shù)組a = diag( 1:5 )b = logical(a)%a = eye(3)b = logical(a)whos % 查看a和b的數(shù)據(jù)類型和存儲關(guān)系運算符=關(guān)系運算主要用來比較兩個尺度相等的數(shù)組對應(yīng)元素大小，返回與原數(shù)組尺度相同的邏輯數(shù)組。對于數(shù)組和標量比較，按標量擴展進行。a = rand(3)b = rand(3)a b a 0.5% 關(guān)系表達式可以和數(shù)學表達式進行混合運算；運算中，邏輯值轉(zhuǎn)換為浮點數(shù)。b - (a0.5)% 例：除數(shù)為0問題1/0 % 結(jié)果為無窮大，Inf0/0 % 結(jié)果不確定，非數(shù)值，NaN (Not a Number)%x = -2:0.5:2sin(x)./xx = x + (x=0)*epssin(x)./x邏輯數(shù)尋址% 按邏輯數(shù)組的索引，根據(jù)其內(nèi)容（1/0）完成相應(yīng)操作。a = magic(3)index = logical( 1 0 1 0 1 ) % 生成邏輯數(shù)組b = a(index) % 提取a(index) = 0 % 賦值a = magic(4)ii = logical( eye(3) )a(ii)a(ii) = -1%a = 1, pi, 0, -2.3index = find(a) logical(a)index = find( a=1 )a=1邏輯運算符& 數(shù)組間元素與操作| 數(shù)組間元素或操作 數(shù)據(jù)元素反操作& 標量關(guān)系表達式的短路式(short-circuiting)與操作| 標量關(guān)系表達式的短路式或操作%a = randperm(5)b = randperm(5)c = abd = c% short-circuiting是指當已知邏輯運算結(jié)果時，就不再執(zhí)行后面的表達式。a = 1, b = 2a = 1 | b = 1a = 2 & b =2關(guān)系和邏輯函數(shù)xor(x, y) 邏輯異或any(x) 判別數(shù)組中 是否 含非零元素all(x) 判別數(shù)組中 是否 全為非零元素% a = 3 1 0 0 b = 2 0 4 0 xor(a, b)any(a)all(a)% 注意運算次序：列優(yōu)先原則a = 1 0 0; 2 0 0; 3 4 0any(a) % 缺?。喊疵苛刑幚韆ny(a, 1) any(a, 2) % 按每行處理% 函數(shù)max/min的使用b = magic(3)max( b, 6 ) % 比較矩陣相應(yīng)元素的大小max( b) % 查找每列元素中的最大值max( b, , 2 ) % 查找每行元素中的最大值mb = max(max(b)i, j = find(b = mb) % b(i, j) % i,j可能不是標量index = find( b = mb ) %b(index)% is-簇函數(shù)：% 檢驗特定條件是否成立，返回邏輯值。isvarname, iskeyword, isequalischar, iscell, isstruct, islogicalisinf, isnan, isemptydoc is* % 查看is-簇函數(shù)特殊變量：Inf, NaN, Inf, -Inf 計算結(jié)果超出能表達的范圍，無窮大。NaN 計算結(jié)果不確定，非數(shù)值 是指一維或多維的長度為0的數(shù)組變量。%a = exp(1000)log(a)a*00/0a = 1 NaN Inf -Inf NaN a.2log(a)% 關(guān)系運算檢測計算結(jié)果是否含NaN或Inf，通常是不可取的a = NaNa = Inf% 正確的方法是采用Matlab函數(shù)isnan(a)isinf(a)% 替換NaN為0i = find( isnan(a) )a(i) = zeros( size(i) ) % a( find(isnan(a) ) = 0;% 替換Inf為一個數(shù)字i = find( isinf(a) )a(i) = 1e10*ones( size(i) )空數(shù)組 運算% 空數(shù)組的維數(shù)和大小a = ones(3, 0)size(a)length(a) % length和numel不計空數(shù)組numel(a)% 空數(shù)組直接進行關(guān)系運算，結(jié)果不確定a = a = a = 1a = 1% 正確的方法是用函數(shù)isempty判別isempty(a)% 空數(shù)組用途1：函數(shù)返回值a = 1:5b = find(a10)isempty(b)% 空數(shù)組用途2：數(shù)值操作a = 1:7a(3:6) = % 刪除指定元素b = magic(4)b(:, 3) = % 刪除整列b(1:3, 2) = % x矩陣形狀不能改變c = zeros(1, 3)n = 2 % repeat for 0 & -1d = ones(1, n), c % 數(shù)組合成運算符和特殊符號help + % 查看Matlab所有符號的功能help precedence % 查看運算符優(yōu)先權(quán)注：1) 同級運算按從左到右次序進行2) 圓括號（）可以改變優(yōu)先權(quán)3) 注意&, |, &, |的優(yōu)先權(quán)順