《MATLAB實踐教程》課件第2章

2.1矩陣運算的操作步驟

2.2創(chuàng)建矩陣

2.3矩陣變換

2.4矩陣算術運算

2.5陣列算術運算

2.6陣列關系運算

2.7陣列邏輯運算

2.8測定矩陣大小

2.9編程中的初始化命令

第2單元MATLAB矩陣運算

矩陣運算的操作步驟如下：

(1)點擊目錄下拉列表框或直接在目錄文本框內鍵入文件夾全路徑，選定工作文件夾，如圖2-1所示。不選則缺省文件夾是work。2.1矩陣運算的操作步驟圖2-1選定工作文件夾

(2)點擊工具條上的New按鈕或File_New_M-file菜單命令，啟動(打開)Editor窗口，如圖2-2所示。圖2-2點擊File菜單上的New_M-File命令(3)在Editor窗口輸入用戶的矩陣運算程序，如圖2-3所示。圖2-3在Editor窗口為矩陣運算編程

(4)點擊Editor窗口的Save按鈕或執(zhí)行File_Save菜單命令，給程序命名并存盤(如myjuzhen01.m)，如圖2-4所示。圖2-4給所編程序命名并存盤

(5)激活CommandWindow窗口，在命令編輯區(qū)鍵入所存程序的文件名(不含擴展名)，回車執(zhí)行并觀察程序運行結果(如鍵入myjuzhen01)，如圖2-5所示。

(6)在CommandWindow中直接輸入程序的各個語句，回車執(zhí)行并觀察程序運行結果。將此結果與在Editor窗口編寫的程序的運行結果進行比較。圖2-5在CommandWindow鍵入程序名并回車2.2.1冒號法創(chuàng)建矩陣

冒號法又稱做from:step:to方式，“from”、“step”、“to”分別表示開始值、步長和結束值。操作步驟如下：

(1)在Editor窗口輸入下述命令并自行命名存盤：

low=1;step=0.5;up=5;

x=low:step:up

y=1:0.5:5

z=1:step/5:5

u=low*pi/2:pi/10:up*pi/52.2創(chuàng)建矩陣

(2)激活CommandWindow窗口，在命令編輯區(qū)鍵入所存程序的文件名，運行結果如下：

x=

1.00001.50002.00002.50003.00003.50004.00004.50005.0000

y=

1.00001.50002.00002.50003.00003.50004.00004.50005.0000

z=

Columns1through9

1.00001.10001.20001.30001.40001.50001.60001.70001.8000Columns10through18

1.90002.00002.10002.20002.30002.40002.50002.60002.7000

Columns19through27

2.80002.90003.00003.10003.20003.30003.40003.50003.6000

Columns28through36

3.70003.80003.90004.00004.10004.20004.30004.40004.5000

Columns37through41

4.60004.70004.80004.90005.0000

u=

1.57081.88502.19912.51332.82743.14162.2.2方括弧法創(chuàng)建矩陣

矩陣元素直接排列在方括弧內，每行的元素使用空格或者逗號隔開，行與行之間使用分號隔開。大的矩陣可以分行輸入，也可用回車鍵代替分號。矩陣元素可以是純數字，也可以是變量。操作步驟如下：

(1)在Editor窗口輸入下述命令并自行命名存盤：

x=[11.522.533.544.55]

y=[12345;678910;1112131415]

z=[1234;5678;9101112]

u=[yz;x]

(2)激活CommandWindow窗口，在命令編輯區(qū)鍵入所存程序的文件名，部分運行結果如下：

x=

1.00001.50002.00002.50003.00003.50004.00004.50005.0000

y=

12345

678910

1112131415

z=

1234

5678

91011122.2.3下標法創(chuàng)建矩陣

下標法創(chuàng)建矩陣的操作步驟如下：

(1)在Editor窗口輸入下述命令并自行命名存盤：

x=[12345;678910;1112131415]

z1=x([13],[135])

z2=x(:,2:4)

z3=x([23],:)

z4=x(1:2,3:5)

z5=x(11)

z6=x(3,5)

x(2:3,:)=[]

(2)激活CommandWindow窗口，在命令編輯區(qū)鍵入所存程序的文件名，運行結果如下：

x=

12345

678910

1112131415

z1=

135

111315z2=

234

789

121314

z3=

678910

1112131415

z4=

345

8910

z5=9

z6=15

x=123452.2.4函數法創(chuàng)建矩陣

函數法創(chuàng)建矩陣的步驟如下：

(1)在Editor窗口輸入下述命令并自行命名存盤：

x1=zeros(3,5),x2=zeros(5)

y1=ones(3,5),y2=ones(5)

z1=eye(3,5),z2=eye(5)

u1=linspace(1,10,5),u2=linspace(1,10)

v1=logspace(1,10,5),v2=logspace(1,10)

w1=rand(3,5),w2=rand(5)

s1=randn(3,5),s2=randn(5)

(2)激活CommandWindow窗口，在命令編輯區(qū)鍵入所存程序的文件名，部分運行結果如下：

x1=

00000

00000

00000

x2=

00000

00000

00000

00000

00000y1=

11111

11111

11111

y2=

11111

11111

11111

11111

11111

z1=

10000

01000

00100z2=

10000

01000

00100

00010

00001

u1=

1.00003.25005.50007.750010.0000

v1=

1.0e+010*

0.00000.00000.00000.00561.0000w1=

0.69790.85370.89980.81800.2897

0.37840.59360.82160.66020.3412

0.86000.49660.64490.34200.5341

w2=

0.72710.70270.79480.97970.1365

0.30930.54660.95680.27140.0118

0.83850.44490.52260.25230.8939

0.56810.69460.88010.87570.1991

0.37040.62130.17300.73730.2987

s1=

-0.6436-0.0195-0.31790.42820.5779

0.3803-0.04821.09500.89560.0403

-1.00910.0000-1.87400.73100.6771s2=

0.5689-0.23400.62320.23790.3899

-0.25560.11840.7990-1.00780.0880

-0.37750.31480.9409-0.7420-0.6355

-0.29591.4435-0.99211.0823-0.5596

-1.4751-0.35100.2120-0.13150.4437

Linspace(p1,p2,p3)函數和logspace(p1,p2,p3)函數的三個輸入參數分別表示開始值、結束值和生成的數據個數。省略參數p3，則Linspace缺省產生100個數據，logspace缺省產生50個數據。

MATLAB提供了許多生成矩陣的函數，表2-1列舉出其中的一部分。表2-1創(chuàng)建矩陣的MATLAB函數2.2.5混合法創(chuàng)建矩陣

混合使用冒號法、方括弧法、下標法、函數法及算式，可創(chuàng)建更為復雜的矩陣。操作方法如下：

(1)在Editor窗口輸入下述命令并自行命名存盤：

a=1;b=2;c=7;

x=[a:b/2:c;zeros(3),ones(3,4);eye(4),rand(4,3)]

x(2,:)=[]

(2)激活CommandWindow窗口，在命令編輯區(qū)鍵入所存程序的文件名，運行結果如下：x=

1.00002.00003.00004.00005.00006.00007.0000

0001.00001.00001.00001.0000

0001.00001.00001.00001.0000

0001.00001.00001.00001.0000

1.00000000.17300.87570.8939

01.0000000.97970.73730.1991

001.000000.27140.13650.2987

0001.00000.25230.01180.6614

x=

1.00002.00003.00004.00005.00006.00007.0000

0001.00001.00001.00001.0000

0001.00001.00001.00001.0000

1.00000000.17300.87570.8939

01.0000000.97970.73730.1991

001.000000.27140.13650.2987

0001.00000.25230.01180.66142.3.1函數變換

矩陣的正弦、余弦、對數等運算，定義在矩陣的單個元素上，即對矩陣的每個元素分別進行正弦、余弦、對數等運算，此類面向矩陣的函數運算稱做函數變換。

MATLAB基本數學函數見表2-2。2.3矩陣變換表2-2MATLAB基本數學函數2.3.2幾何變換

操作方法如下：

(1)在Editor窗口輸入下述命令并自行命名存盤：

a=[12345;678910;1112131415]

x1=flipud(a)

x2=fliplr(a)

x3=rot90(a)

x4=sort(a)

x5=reshape(a,5,3)

(2)激活CommandWindow窗口，在命令編輯區(qū)鍵入所存程序的文件名，運行結果如下：a=

12345

678910

1112131415

x1=

1112131415

678910

12345

x2=

54321

109876

1514131211x3=

51015

4914

3813

2712

1611

x4=

12345

678910

1112131415

x5=

1129

6314

1185

21310

7415矩陣算術運算符有+(加)、-(減)、*(乘)、\(左除)、/(右除)、^(冪乘)、’(矩陣轉置)等。矩陣算術運算遵循線性代數規(guī)則，兩矩陣左除行數相同，右除列數相同。操作步驟如下：

2.4矩陣算術運算

(1)在Editor窗口輸入下述命令并自行命名存盤：

A=[12310;54611;97815]

B=[3250;4307;5912]

C=[123;546;978]

D=[24;17;35]

X=[A+B,A-B]

Y=C*B

Z=[C\D,B/A]

U=C^3

V=A’+B’

(2)激活CommandWindow窗口，在命令編輯區(qū)鍵入所存程序的文件名，運行結果如下：

A=

12310

54611

97815

B=

3250

4307

5912

C=

123

546

978D=

24

17

35

X=

44810-2?0-210

97618?1?1?64

1416917?4-2?713

Y=

2635820

61763140

951115365Z=

-10.3333-1.77223.3376-1.2516

21.8667?1.2469-3.4938?2.1852

?-3.1333?0.3921-4.7497?3.4061

U=

544473612

120210391341

182315712024

V=

4914

4716

869

101817陣列算術運算符有+(加)、-(減)、.*(陣列乘)、.\(陣列左除)、./(陣列右除)、.^(陣列冪乘)、.’(陣列轉置)等。陣列算術運算遵循元素對元素的規(guī)則，參與運算的兩個陣列或維數相同或其中一個是標量。操作步驟如下：2.5陣列算術運算

(1)在Editor窗口輸入下述命令并自行命名存盤：

A=[12310;54611;97815]

B=[3250;4307;5912]

X=[A+B,A-B]

Y=[A.*B,A.^2]

Z=[A./B,A.\B]

U=[A’zeros(4,3);zeros(4,3)B’]

(2)激活CommandWindow窗口，在命令編輯區(qū)鍵入所存程序的文件名，運行結果如下：A=

12310

54611

97815

B=

3250

4307

5912

X=

44810-20-210

97618?11?64

14169174-2713Y=

34150149100

2012077251636121

4563830814964225

Warning:Dividebyzero.

Z=

0.33331.00000.6000Inf3.00001.00001.66670

1.25001.3333Inf1.57140.80000.750000.6364

1.80000.77788.00007.50000.55561.28570.12500.1333U=

159000

247000

368000

101115000

000345

000239

000501

000072陣列關系運算符有==(等于)、<(小于)、<=(小于等于)、>(大于)、>=(大于等于)、～=(不等于)等。陣列關系運算遵循元素對元素的規(guī)則，參與運算的兩個陣列(包括算式產生的陣列)或維數相同或其中一個是標量。操作步驟如下：2.6陣列關系運算

(1)在Editor窗口輸入下述命令并自行命名存盤：

A=[12310;54611;97815]

B=[3250;4307;5912]

X1=A>B,

X2=A>=B

X3=A~=B

X4=A>5

X5=(A.*B)>(A./B)

X6=(A.*B)~=(A./B)

X7=A.*sin(B)<B.*cos(A)

(2)激活CommandWindow窗口，在命令編輯區(qū)鍵入所存程序的文件名，部分運行結果如下：A=

12310

54611

97815

B=

3250

4307

5912

X1=

0001

1111

1011X2=

0101

1111

1011

X3=

1011

1111

1111

Warning:Dividebyzero.

X5=

1110

1101

1101Warning:Dividebyzero.

X6=

1111

1111

1101

X7=

1000

1000

1100

(3)在Editor窗口輸入下述命令并自行命名存盤：

x=-3:0.1:3;y=2*x.^2-5*x+15;

maxy=y(1);miny=y(1);n=length(y);

maxp=[x(1)y(1)];minp=[x(1)y(1)];

fori=1:n

ify(i)>maxy

maxp=[x(i)y(i)];maxy=y(i);

elseify(i)<=miny

minp=[x(i)y(i)];miny=y(i);

end

end

maxp=maxp,minp=minp

plot(x,y,maxp(1),maxp(2),'*',minp(1),minp(2),'p')

(4)激活CommandWindow窗口，在命令編輯區(qū)鍵入所存程序的文件名，運行結果如下：

maxp=

-348

minp=

1.300011.8800

程序運行結果如圖2-6所示。圖2-6二次曲線及其搜索范圍內的最大點和最小點邏輯運算符有&(與)、|?(或)、～(非)、xor(異或)等。陣列邏輯運算遵循元素對元素的規(guī)則，參與運算的兩個陣列(包括算式產生的陣列)或維數相同或其中一個是標量。操作步驟如下：2.7陣列邏輯運算

(1)在Editor窗口輸入下述命令并自行命名存盤：

X=[12310;54611;97815]

Y=[3250;4307;5912]

Z1=[X&Y,X&5]

Z2=[X|YX|5]

Z3=[X>2|Y<5X>2&Y<5]

Z4=[~X~Y~(X&Y)]

Z5=[xor(X,Y)xor(X,5)xor(X>2,Y<5)]

(2)激活CommandWindow窗口，在命令編輯區(qū)鍵入所存程序的文件名，運行結果如下：X=

12310

54611

97815

Y=

3250

4307

5912

Z1=

11101111

11011111

11111111Z2=

11111111

11111111

11111111

Z3=

11110001

11