matlab教程(第3講-數(shù)組)

第3章Matlab數(shù)據(jù)及其運(yùn)算

Matlab數(shù)據(jù)類型非常豐富，除數(shù)值型、字符型等基本數(shù)據(jù)類型外，還有結(jié)構(gòu)體、單元等更為復(fù)雜的數(shù)據(jù)類型。各種數(shù)據(jù)類型都以矩陣形式存在，矩陣是Matlab最基本的數(shù)據(jù)對(duì)象，并且矩陣的運(yùn)算是定義在復(fù)數(shù)域上的。

2023/2/4ApplicationofMatlabLanguage2

2.1數(shù)值表示、變量及表達(dá)式數(shù)值的記述

Matlab的數(shù)只采用習(xí)慣的十進(jìn)制表示，可以帶小數(shù)點(diǎn)和負(fù)號(hào);其缺省的數(shù)據(jù)類型為雙精度浮點(diǎn)型（double）。例如：3-100.0011.3e101.256e-6變量命令規(guī)則變量名、函數(shù)名對(duì)字母的大小寫是敏感的。如myVar與myvar表示兩個(gè)不同的變量。變量名第一個(gè)字母必須是英文字母。變量名可以包含英文字母、下劃線和數(shù)字。變量名不能包含空格、標(biāo)點(diǎn)。變量名最多可包含63個(gè)字符（6.5及以后的版本）。

2023/2/4ApplicationofMatlabLanguage3

2.1數(shù)值表示、變量及表達(dá)式（續(xù)）Matlab預(yù)定義的變量變量名意義ans最近的計(jì)算結(jié)果的變量名epsMATLAB定義的正的極小值=2.2204e-16pi圓周率πinf∞值，無限大i或j虛數(shù)單元，sqrt(-1)NaN非數(shù)，0/0、∞/∞〖說明〗每當(dāng)MATLAB啟動(dòng)完成，這些變量就被產(chǎn)生。MATLAB中，被0除不會(huì)引起程序中斷，給出報(bào)警的同時(shí)用inf或NaN給出結(jié)果。用戶只能臨時(shí)覆蓋這些預(yù)定義變量的值，Clear或重啟MATLAB可恢復(fù)其值。

2023/2/4ApplicationofMatlabLanguage4

2.1數(shù)值表示、變量及表達(dá)式（續(xù)）運(yùn)算符和表達(dá)式運(yùn)算數(shù)學(xué)表達(dá)式MATLAB運(yùn)算符MATLAB表達(dá)式加a+b+a+b減a-b-a-b乘axb*a*b除a/b或a\b/或\a/b或a\b冪^a^b〖說明〗Matlab用“\”和”/”分別表示“左除”和“右除”。對(duì)標(biāo)量而言，兩者沒有區(qū)別。對(duì)矩陣產(chǎn)生不同影響。MATLAB表達(dá)式的書寫規(guī)則與“手寫方式”幾乎完全相同。表達(dá)式按與常規(guī)相同的優(yōu)先級(jí)自左至右執(zhí)行運(yùn)算。優(yōu)先級(jí)：指數(shù)運(yùn)算級(jí)別最高，乘除次之，加減最低。括號(hào)改變運(yùn)算的次序。

2023/2/4ApplicationofMatlabLanguage5

2.1數(shù)值表示、變量及表達(dá)式（續(xù)）復(fù)數(shù)及其運(yùn)算MATLAB中復(fù)數(shù)的表達(dá)：z=a+bi，其中a、b為實(shí)數(shù)。MATLAB把復(fù)數(shù)作為一個(gè)整體，象計(jì)算實(shí)數(shù)一樣計(jì)算復(fù)數(shù)?！纠?.3-1】復(fù)數(shù)z1=3+4i，z2=1+2i,z3=

計(jì)算z1=3+4*i,z2=1+2*i,z3=exp(i*pi/6),z=z1*z2/z3z_real=real(z),z_image=imag(z),z_angle=angle(z),z_length=abs(z),

2023/2/4ApplicationofMatlabLanguage6

2.2Matlab矩陣(數(shù)組)的表示數(shù)組的概念一維數(shù)組變量的創(chuàng)建二維數(shù)組變量的創(chuàng)建數(shù)組元素的標(biāo)識(shí)與尋訪數(shù)組運(yùn)算多維數(shù)組

2023/2/4ApplicationofMatlabLanguage7

2.2.1.數(shù)組(array)的概念數(shù)組定義：按行(row)和列(column)順序排列的實(shí)數(shù)或復(fù)數(shù)的有序集，被稱為數(shù)組。數(shù)組中的任何一個(gè)數(shù)都被稱為這個(gè)數(shù)組的元素，由其所在的行和列標(biāo)識(shí)，這個(gè)標(biāo)識(shí)也稱為數(shù)組元素的下標(biāo)或索引。Matlab將標(biāo)量視為1×1的數(shù)組。對(duì)m行、n列的2維數(shù)組a：計(jì)為m×n的數(shù)組a；*行標(biāo)識(shí)、列標(biāo)識(shí)均從1開始；行標(biāo)識(shí)從上到下遞增；列標(biāo)識(shí)從左到右遞增。12345222232425332333435442434445a(3,4)=34rowisfirsta=

2023/2/4ApplicationofMatlabLanguage8

數(shù)組的分類一維數(shù)組，也稱為向量(vector)。行向量(rowvector)、列向量(columnvector)。二維數(shù)組(矩陣matrix)。多維數(shù)組。有效矩陣：每行元素的個(gè)數(shù)必須相同，每列元素的個(gè)數(shù)也必須相同。2.2.1.數(shù)組(array)的概念

2023/2/4ApplicationofMatlabLanguage9

2.2.1、數(shù)組(array)的概念（續(xù)）行向量列向量a(2,1)=3a(1,2)=2b(3)=3c(2)=2

2023/2/4ApplicationofMatlabLanguage10

2.2.2、創(chuàng)建一維數(shù)組變量第一種方法：使用方括號(hào)“[]”操作符【例2-1】創(chuàng)建數(shù)組(行向量)a=[13pi3+5i]>>a=[13pi3+5*i]%ora=[1,3,pi,3+5*i]

a=1.00003.00003.14163.0000+5.0000i所有的向量元素必須在操作符“[]”之內(nèi)；向量元素間用空格或英文的逗點(diǎn)“,”分開。第二種方法：使用冒號(hào)“:”操作符【例2-2】創(chuàng)建以1~10順序排列整數(shù)為元素的行向量b。>>b=1:10b=12345678910

2023/2/4ApplicationofMatlabLanguage11

2.2.2、創(chuàng)建一維數(shù)組變量（續(xù)）【例2-3】鍵入并執(zhí)行c=1:2:10和d=1:2:9>>c=1:2:10c=13579>>d=1:2:9d=13579利用冒號(hào)“:”操作符創(chuàng)建行向量的基本語法格式：x=Start:Increment:EndStart表示新向量x的第一個(gè)元素；新向量x的最后一個(gè)元素不能大于End；Increment可正可負(fù)，若負(fù)，則必須Start>End；若正，則必須Start<End，否則創(chuàng)建的為空向量。若Increment=1,則可簡(jiǎn)寫為：x=Start:End。

2023/2/4ApplicationofMatlabLanguage12

2.2.2、創(chuàng)建一維數(shù)組變量（續(xù)）第三種方法：利用函數(shù)linspace函數(shù)linspace的基本語法x=linspace(x1,x2,n)該函數(shù)生成一個(gè)由n個(gè)元素組成的行向量；x1為其第一個(gè)元素；x2為其最后一個(gè)元素；x1、x2之間元素的間隔=(x2-x1)/(n-1)。如果忽略參數(shù)n，則系統(tǒng)默認(rèn)生成100個(gè)元素的行向量?！纠?-4】鍵入并執(zhí)行x=linspace(1,2,5)x=1.00001.25001.50001.75002.0000同學(xué)們可以在實(shí)驗(yàn)時(shí)察看x=linspace(1,2)執(zhí)行結(jié)果。

2023/2/4ApplicationofMatlabLanguage13

2.2.2、創(chuàng)建一維數(shù)組變量（續(xù)）第四種方法：利用函數(shù)logspace通過實(shí)驗(yàn)認(rèn)識(shí)該函數(shù)的功能。列向量的創(chuàng)建使用方括號(hào)“[]”操作符，使用分號(hào)“;”分割行?！纠?-5】鍵入并執(zhí)行x=[1;2;3]X=123使用冒號(hào)操作符【例2-6】鍵入并執(zhí)行x=(1:3)’%“’”表示矩陣的轉(zhuǎn)置

2023/2/4ApplicationofMatlabLanguage14

2.2.2、創(chuàng)建一維數(shù)組變量

（續(xù)）創(chuàng)建數(shù)組變量的一般方法創(chuàng)建變量的賦值語句的一般格式var=expressionvar為變量名expression為MATLAB合法表達(dá)式可以是單獨(dú)的常數(shù)值或數(shù)值數(shù)組；也可以由常數(shù)值、其他變量（部分或全部）、數(shù)值數(shù)組和運(yùn)算符（+、-等）構(gòu)成。【例2-7】鍵入并執(zhí)行a=[01+6];b=[a67];c=[6a7];d=[6a7a];

2023/2/4ApplicationofMatlabLanguage15

2.2.2、創(chuàng)建一維數(shù)組變量

（續(xù)）一旦被創(chuàng)建，變量就被存儲(chǔ)在工作空間，可以通過“Workspace”窗口或在“CommandWindow”執(zhí)行“whos”命令察看。操作一：使用“Workspace”窗口察看變量操作二：使用“whos”命令察看變量

2023/2/4ApplicationofMatlabLanguage16

2.2.3、創(chuàng)建二維數(shù)組變量第一種方法：使用方括號(hào)“[]”操作符使用規(guī)則數(shù)組元素必須在“[]”內(nèi)鍵入；行與行之間須用分號(hào)“;”間隔，也可以在分行處用回車鍵間隔；行內(nèi)元素用空格或逗號(hào)“,”間隔?！纠?-1】鍵入并執(zhí)行a2=[123;456;789]a2=123456789【例3-2】鍵入并執(zhí)行a2=[1:3;4:6;7:9]%結(jié)果同上

2023/2/4ApplicationofMatlabLanguage17

2.2.3、創(chuàng)建二維數(shù)組變量（續(xù)）【例3-3】由向量構(gòu)成二維數(shù)組。>>a=[123];b=[234];>>c=[a;b];>>c1=[ab];第二種方法：函數(shù)方法函數(shù)ones(生成全1矩陣)、zeros(生成全0矩陣)、reshape“helpelmat”獲得基本的矩陣生成和操作函數(shù)列表【例3-4】創(chuàng)建全1的3x3數(shù)組。>>ones(3)【例3-5】創(chuàng)建全1的3x4數(shù)組。>>ones(3,4)

2023/2/4ApplicationofMatlabLanguage18

2.2.3、創(chuàng)建二維數(shù)組變量（續(xù)）【例3-5】reshape的使用演示>>a=-4:4a=-4-3-2-101234>>b=reshape(a,3,3)b=-4-12-303-214數(shù)組元素的排列順序，從上到下按列排列，先排第一列，然后第二列，…要求數(shù)組的元素總數(shù)不變。

2023/2/4ApplicationofMatlabLanguage19

4、數(shù)組元素的標(biāo)識(shí)與尋訪數(shù)組元素的標(biāo)識(shí)“全下標(biāo)（index）”標(biāo)識(shí)經(jīng)典數(shù)學(xué)教科書采用“全下標(biāo)”標(biāo)識(shí)法：每一維對(duì)應(yīng)一個(gè)下標(biāo)。如對(duì)于二維數(shù)組，用“行下標(biāo)和列下標(biāo)”標(biāo)識(shí)數(shù)組的元素，a(2,3)就表示二維數(shù)組a的“第2行第3列”的元素。對(duì)于一維數(shù)組，用一個(gè)下標(biāo)即可，b(2)表示一維數(shù)組b的第2個(gè)元素，無論b是行向量還是列向量?！皢蜗聵?biāo)”（linearindex）標(biāo)識(shí)所謂“單下標(biāo)”標(biāo)識(shí)就是用一個(gè)下標(biāo)來表明元素在數(shù)組的位置。對(duì)于二維數(shù)組，“單下標(biāo)”編號(hào)：設(shè)想把二維數(shù)組的所有列，按先后順序首尾相接排成“一維長列”，然后自上往下對(duì)元素位置執(zhí)行編號(hào)。兩種“下標(biāo)”標(biāo)識(shí)的變換：sub2ind、ind2sub

2023/2/4ApplicationofMatlabLanguage20

4、數(shù)組元素的標(biāo)識(shí)與尋訪（續(xù)）【例4-1】單下標(biāo)的使用>>a=zeros(2,5);>>a(:)=-4:5a=-4-2024-3-1135注意數(shù)組的排列順序。

2023/2/4ApplicationofMatlabLanguage21

4、數(shù)組元素的標(biāo)識(shí)與尋訪（續(xù)）元素與子數(shù)組的尋訪與賦值【例4-3】一維數(shù)組元素與子數(shù)組的尋訪與賦值>>a=linspace(1,10,5)a=1.00003.25005.50007.750010.0000>>a(3)%尋訪a的第3個(gè)元素ans=5.5000>>a([125])%尋訪a的第1、2、5個(gè)元素組成的子數(shù)組ans=1.00003.250010.0000

2023/2/4ApplicationofMatlabLanguage22

4、數(shù)組元素的標(biāo)識(shí)與尋訪（續(xù)）>>a(1:3)%尋訪前3個(gè)元素組成的子數(shù)組ans=1.00003.25005.5000>>a(3:-1:1)%由前3個(gè)元素倒序構(gòu)成的子數(shù)組ans=5.50003.25001.0000>>a(3:end)ans=5.50007.750010.0000>>a(3:end-1)ans=5.50007.7500%第3個(gè)及其后所有元素構(gòu)成的子數(shù)組函數(shù)end作為參數(shù)使用，返回最后一個(gè)元素的下標(biāo)

2023/2/4ApplicationofMatlabLanguage23

>>a([12355321])ans=1.00003.25005.500010.000010.00005.50003.25001.0000數(shù)組元素可以被任意重復(fù)訪問，構(gòu)成長度大于原數(shù)組的新數(shù)組。>>a(6)???Indexexceedsmatrixdimensions.下標(biāo)值超出了數(shù)組的維數(shù)，導(dǎo)致錯(cuò)誤>>a(2.1)???Subscriptindicesmusteitherberealpositiveintegersorlogicals.下標(biāo)值只能取正整數(shù)或邏輯值4、數(shù)組元素的標(biāo)識(shí)與尋訪（續(xù)）

2023/2/4ApplicationofMatlabLanguage24

4、數(shù)組元素的標(biāo)識(shí)與尋訪（續(xù)）>>a(3)=0%修改數(shù)組a的第3元素值為0a=1.00003.250007.750010.0000>>a([25])=[11]a=1.00001.000007.75001.0000可以修改指定數(shù)組元素的值一次可以修改多個(gè)數(shù)組元素的值要修改的數(shù)組元素的個(gè)數(shù)應(yīng)與送入數(shù)組的元素個(gè)數(shù)相同

2023/2/4ApplicationofMatlabLanguage25

4、數(shù)組元素的標(biāo)識(shí)與尋訪（續(xù)）【例4-3】二維數(shù)組元素與子數(shù)組的尋訪與賦值>>a_2=zeros(2,4)%創(chuàng)建2x4的全0數(shù)組a_2=00000000>>a_2(:)=1:8a_2=13572468>>a_2([258])%單下標(biāo)方式尋訪多個(gè)元素ans=258注意元素的排列順序

2023/2/4ApplicationofMatlabLanguage26

4、數(shù)組元素的標(biāo)識(shí)與尋訪（續(xù)）>>a_2([258])=[102030]a_2gt;>a_2(:,[23])=ones(2)%雙下標(biāo)方式尋訪并修改a_2=1117101130二維數(shù)組可以“單下標(biāo)”方式或“全下標(biāo)”方式訪問、賦值；“單下標(biāo)”方式賦值時(shí)，等號(hào)兩邊涉及的元素個(gè)數(shù)必須相等；“全下標(biāo)”方式賦值時(shí)，等號(hào)右邊數(shù)組的大小必須等于原數(shù)組中涉及元素構(gòu)成的子數(shù)組的大小。

2023/2/4ApplicationofMatlabLanguage27

4、數(shù)組元素的標(biāo)識(shí)與尋訪（續(xù)）>>a_2(:,end)ans=730>>a_2(:,end-1)ans=11>>a_2(:,end:-1:3)ans=71301>>a_2(end,:)ans=101130>>a_2(end,[2:4])ans=1130>>a_2([46])=6:7a_2=1117106730>>a_2(end,[2:end-1])

Whatistheresult?

2023/2/4ApplicationofMatlabLanguage28

4、數(shù)組元素的標(biāo)識(shí)與尋訪（續(xù)）【例4-4】

size、length函數(shù)>>a=ones(4,6)*6>>m=size(a)>>len=length(a)>>b=1:5;>>length(b)>>c=b’>>length(c)

size函數(shù)返回變量的大小，即變量數(shù)組的行列數(shù)

length函數(shù)返回變量數(shù)組的最大維數(shù)

2023/2/4ApplicationofMatlabLanguage29

4、數(shù)組元素的標(biāo)識(shí)與尋訪（續(xù)）雙下標(biāo)到單下標(biāo)的轉(zhuǎn)換【例4-5】

sub2ind函數(shù)-雙下標(biāo)轉(zhuǎn)換為單下標(biāo)>>A=[172418;222714;461320];>>A(:,:,2)=A-10>>A(2,1,2)>>sub2ind(size(A),2,1,2)>>A(14)

2023/2/4ApplicationofMatlabLanguage30

4、數(shù)組元素的標(biāo)識(shí)與尋訪（續(xù)）單下標(biāo)到雙下標(biāo)的轉(zhuǎn)換【例4-6】

ind2sub函數(shù)-雙下標(biāo)轉(zhuǎn)換為單下標(biāo)>>b=zeros(3);>>b(:)=1:9>>IND=[3456]>>[I,J]=ind2sub(size(b),IND)

2023/2/4ApplicationofMatlabLanguage31

6、數(shù)組的算術(shù)運(yùn)算MATLAB數(shù)組支持線性代數(shù)中所有的矩陣運(yùn)算。建立特有的數(shù)組運(yùn)算符，如：“.*”、“./”等。

MATLAB數(shù)組運(yùn)算符列表運(yùn)算 運(yùn)算符 含義說明

加 +相應(yīng)元素相加減 - 相應(yīng)元素相減乘 * 矩陣乘法點(diǎn)乘 .* 相應(yīng)元素相乘冪 ^ 矩陣冪運(yùn)算點(diǎn)冪 .^ 相應(yīng)元素進(jìn)行冪運(yùn)算左除或右除 \或/ 矩陣左除或右除左點(diǎn)除或右點(diǎn)除.\或./ A的元素被B的對(duì)應(yīng)元素除

2023/2/4ApplicationofMatlabLanguage32

6、算術(shù)運(yùn)算（續(xù)）

【例5-1】數(shù)組加減法>>a=zeros(2,3);>>a(:)=1:6;>>b=a+2.5b=3.50005.50007.50004.50006.50008.5000>>c=b-ac=2.50002.50002.50002.50002.50002.5000

2023/2/4ApplicationofMatlabLanguage33

6、算術(shù)運(yùn)算（續(xù)）【例5-2】畫出y=1/(x+1)的函數(shù)曲線，x∈[0,100]。x=0:100;y=1./(x+1);plot(x,y);legend(‘y=1/(x+1)’);

2023/2/4ApplicationofMatlabLanguage34

6、算術(shù)運(yùn)算（續(xù)）【例5-2】生成一個(gè)信號(hào)：x=sin(2*pi*t)+sin(4*pi*t)t=[0:199]./100;%采樣時(shí)間點(diǎn)%生成信號(hào)x=sin(2*pi*t)+sin(4*pi*t);plot(t,x);legend(‘x=sin(2*pi*t)+sin(4*pi*t)’);

2023/2/4ApplicationofMatlabLanguage35

6、算術(shù)運(yùn)算（續(xù)）

2023/2/4ApplicationofMatlabLanguage36

6、算術(shù)運(yùn)算（續(xù)）【例5-2】點(diǎn)冪“.^”舉例>>a=1:6a=123456>>b=reshape(a,2,3)b=135246>>a=a.^2a=149162536>>b=b.^2b=192541636

2023/2/4ApplicationofMatlabLanguage37

7、關(guān)系運(yùn)算

Matlab提供了6種關(guān)系運(yùn)算符：<、>、<=、>=、==、~=（不等于）關(guān)系運(yùn)算符的運(yùn)算法則：1、當(dāng)兩個(gè)標(biāo)量進(jìn)行比較時(shí)，直接比較兩數(shù)大小。若關(guān)系成立，結(jié)果為1，否則為0。2、當(dāng)兩個(gè)維數(shù)相等的矩陣進(jìn)行比較時(shí)，其相應(yīng)位置的元素按標(biāo)量關(guān)系進(jìn)行比較，并給出結(jié)果，形成一個(gè)維數(shù)與原來相同的0、1矩陣。3、當(dāng)一個(gè)標(biāo)量與一個(gè)矩陣比較時(shí)，該標(biāo)量與矩陣的各元素進(jìn)行比較，結(jié)果形成一個(gè)與矩陣維數(shù)相等的0、1矩陣。

2023/2/4ApplicationofMatlabLanguage38

7、關(guān)系運(yùn)算（續(xù)）【例】建立5階方陣A，判斷其元素能否被3整除。A=[24,35,13,22,63;23,39,47,80,80;...90,41,80,29,10;45,57,85,62,21;37,19,31,88,76]A=2435132263233947808090418029105785622119318876P=rem(A,3)==0%被3除，求余P=1000101000100001100100000

2023/2/4ApplicationofMatlabLanguage39

8、邏輯運(yùn)算

Matlab提供了3種邏輯運(yùn)算符：&（與）、|（或）、~（非）邏輯運(yùn)算符的運(yùn)算法則：1、在邏輯運(yùn)算中，確認(rèn)非零元素為真（1），零元素為假（0）。

2、當(dāng)兩個(gè)維數(shù)相等的矩陣進(jìn)行比較時(shí)，其相應(yīng)位置的元素按標(biāo)量關(guān)系進(jìn)行比較，并給出結(jié)果，形成一個(gè)維數(shù)與原來相同的0、1矩陣；3、當(dāng)一個(gè)標(biāo)量與一個(gè)矩陣比較時(shí)，該標(biāo)量與矩陣的各元素進(jìn)行比較，結(jié)果形成一個(gè)與矩陣維數(shù)相等的0、1矩陣；4、算術(shù)運(yùn)算優(yōu)先級(jí)最高，邏輯運(yùn)算優(yōu)先級(jí)最低。

2023/2/4ApplicationofMatlabLanguage40

【例】在[0，3π]區(qū)間，求的值。要求消去負(fù)半波，即(π,2π)區(qū)間內(nèi)的函數(shù)值置零。x=0:pi/100:3*pi;y=sin(x);y1=(y>=0).*y;%消去負(fù)半波

2023/2/4ApplicationofMatlabLanguage41

【例】建立矩陣A，找出在[10，20]區(qū)間的元素的位置。

A=[4,15,-45,10,6;56,0,17,-45,0];find(A>=10&A<=20)％找到非零元素的位置A=415-4510656017-450ans=367

2023/2/4ApplicationofMatlabLanguage42

7.數(shù)據(jù)分析與統(tǒng)計(jì)最大值和最小值

MATLAB提供的求數(shù)據(jù)序列的最大值和最小值的函數(shù)分別為max和min，兩個(gè)函數(shù)的調(diào)用格式和操作過程類似。

1、求向量的最大值和最小值求一個(gè)向量X的最大值的函數(shù)有兩種調(diào)用格式，分別是：

(1)y=max(X)：返回向量X的最大值存入y，如果X中包含復(fù)數(shù)元素，則按模取最大值；

(2)[y,I]=max(X)：返回向量X的最大值存入y，最大值的序號(hào)存入I，如果X中包含復(fù)數(shù)元素，則按模取最大值。求向量X的最小值的函數(shù)是min(X)，用法和max(X)完全相同。

2023/2/4ApplicationofMatlabLanguage43

【例7-1】求向量的最大值>>x=[-43,72,9,16,23,47];>>y=max(x)%求向量x中的最大值y=72>>[y,l]=max(x)%求向量x中的最大值及其該元素的位置y=72l=2

2023/2/4ApplicationofMatlabLanguage44

9.數(shù)據(jù)分析與統(tǒng)計(jì)2.求矩陣的最大值和最小值求矩陣A的最大值的函數(shù)有3