第一章基礎(chǔ)準(zhǔn)備及入門

本章有兩個目的：一是講述MATLAB正常運行所必須具備的基礎(chǔ)條件；二是簡明系統(tǒng)

地介紹高度集成的Desktop操作桌面的功能和使用方法。

本章的前兩節(jié)分別講述：MATLAB的正確安裝方法和MATLAB環(huán)境的啟動。因為指

令窗是MATLAB最重要的操作界面，所以本章用第1.3、1.4兩節(jié)以最簡單通俗的敘述、算

例講述指令窗的基本操作方法和規(guī)則。這部分內(nèi)容幾乎對MATLAB各種版本都適用。

MATLAB6.X不同于其前版本的最突出之處是：向用戶提供前所未有的、成系列的交互式

工作界面。了解、熟悉和掌握這些交互界面的基本功能和操作方法，將使新老用戶能事半

功倍地利用MATLAB去完成各種學(xué)習(xí)和研究。為此，本章特設(shè)幾節(jié)用于專門介紹最常用的

交互界面：歷史指令窗、當(dāng)前目錄瀏覽器、工作空間瀏覽器、內(nèi)存數(shù)組編輯器、交互界面分

類目錄窗、M文件編輯/調(diào)試器、及幫助導(dǎo)航/瀏覽器。

本章是根據(jù)MATLAB6,5版編寫的，但大部分內(nèi)容也適用于其他6.x版。

1.1MATLAB的安裝和內(nèi)容選擇

ProductList

1.Selectdirectorywhereproductsvillbeinstalled.Spaceavailable：

|DMATLAB6P5

7523H

2.StlectinsitlltUonoption，

Sp?c?required

'?Installproductsanddocuntnlttion.forproduct(1),

「［Installproductsonly.HoreOptionsifany

Installdocumentationonly1115H

3.Selectlanguageofdocunentation.

required

'?Englishonly.

doCVDentation,if

「EnglishandJapanese,ifavailable.

any：

433H

4.Selectproductsand/ordocunent&tion.

GMATLAB6.5

GSinulink5.0

BAerospaceBlockset1.01Totdsptc*

GCDMAReferenceBlockset1.1

rConnunicationsBlockset

rConnunicationsToolbox2.11548M

圖1.1-1

1.2Desktop操作桌面的啟動

1.2.1MATLAB的啟動

1.2.2Desktop操作桌面簡介

一操作桌面的缺省外貌

圖1.2-1

二通用操作界面

1.3CommandWindow運行入門

1.3.1CommandWindow指令窗簡介

1.3.2最簡單的計算器使用法

【例1.321］求［12+2x(7—4)］+3?的算術(shù)運算結(jié)果。

(1)用鍵盤在MATLAB指令窗中輸入以下內(nèi)容

?(12+2*(7-4))/3A2

(2)在上述表達(dá)式輸入完成后，按【Enter】鍵，該就指令被執(zhí)行。

(3)在指令執(zhí)行后，MATLAB指令窗中將顯示以下結(jié)果。

ans

2

ri23/

【例1.322】簡單矩陣A=456*的輸入步驟。

0C

<789f

(1)在鍵盤上輸入下列內(nèi)容

A=［1,2,3;4,5,6;7,8,9］

(2)按［Enter］鍵，指令被執(zhí)行。

(3)在指令執(zhí)行后，MATLAB指令窗中將顯示以下結(jié)果:

A=

123

456

789

【例132-3】矩陣的分行輸入。

A=[l,2,3

4,5,6

7,8,9]

A=

123

456

789

【例1.3.2-4]指令的續(xù)行輸入

S=l-l/2+l/3-l/4+

1/5-1/6+1/7-1/8

S=

0.6345

1.3.3數(shù)值、變量和表達(dá)式

一數(shù)值的記述

二變量命名規(guī)則

三MATLAB默認(rèn)的預(yù)定義變量

四運算符和表達(dá)式

五復(fù)數(shù)和復(fù)數(shù)矩陣

133-1］4=1+21/3=2e；z=33

【例復(fù)數(shù)=3+4Z,Z2表達(dá)'及計算

Z3

(1)

zl=3+4i

zl=

3.0000+4.0000i

(2)

z2=1+2*i

z3=2*exp(i*pi/6)

z=zl*z2/z3

z2=

1.0000+2.0000i

z3=

1.7321+l.OOOOi

z=

0.3349+5.580H

【例133-2】復(fù)數(shù)矩陣的生成及運算

A=[l,3;2z4]-[5z8;6z9]*i

B=[l+5iz2+6i;3+8*i,4+9*i]

C=A*B

A

1.0000-5.0000i3.0000-8.0000i

2.0000-6.0000i4.0000-9.0000i

B

1.0000+5.0000i2.0000+6.0000i

3.0000+8.0000i4.0000+9.0000i

C

1.0e+002火

0.99001.1600-0.0900i

1.1600+0.0900i1.3700

【例133-3】求上例復(fù)數(shù)矩陣C的實部、虛部、模和相角。

C_real=real(C)

C_imag=imag(C)

C_magnitude=abs(C)

C_phase=angle(C)*180/pi*以度為單位計算相角

Creal=

99116

116137

Cimag=

0-9

90

C_magnitude

99.0000116.3486

116.3486137.0000

C_phase=

0-4.4365

4.43650

【例1.3.3-4]用MATLAB計算V^8能得到-2嗎？

(1)

a=-8;

r=aA(1/3)

r=

1.0000+1.7321i

(2)

m=[0,lz2];

R=abs(a)A(1/3);

Theta=(angle(a)+2*pi*m)/3;

rrr=R*exp(i*Theta)

rrr=

1.0000+1.732H-2.0000+O.OOOOi1.0000-1.7321i

(3)

t=0:pi/20:2*pi;x=R*sin(t);y=R*cos(t);

plot(x,yz'b:*),grid

holdon

11

plot(rrr(1)A.'z?MarkerSize',50,Color*z*r*)

1

plot(rrr([2,3])z'o,zMarkerSize',15,*Color*z*b*)

axis([-3,3,-3,3]),axissquare

holdoff

圖1.3-2

1.3.4計算結(jié)果的圖形表示

1_L

【例1.341】畫出衰減振蕩曲線y=/3sin3/及其它的包絡(luò)線方=/3?！旱娜≈捣秶?/p>

[0,44]。(圖1.3-3)

t=0:pi/50:4*pi;

y0=exp(-t/3);

y=exp(-t/3).*sin(3*t);

,,!ff

plot(t,y,-rztzyOz':bztz-yOz:b)

grid

si./+9)

【例1.3.4-21畫出z=—J,所表示的三維曲面(圖1.3-4)。的取值范圍是

[-8,8]?

clear;x=-8:0.5:8;

y=x';

X=ones(size(y))*x;

Y=y*ones(size(x));

R=sqrt(X.A2+Y.A2)+eps;%<5>

Z=sin(R)./R;%<6>

surf(X,Y,Z);%

colormap(cool)%

xlabel(*x*),ylabel('y'),zlabel('z*)

1.3-4

1.4CommandWindow操作要旨

1.4.1指令窗顯示方式的操作

一缺省顯示方式

二顯示方式的設(shè)置

1.4.2數(shù)值計算結(jié)果的顯示格式

1.4.3指令行中的標(biāo)點符號

1.4.4指令窗的常用控制指令

1.4.5指令窗中指令行的編輯

【例1.4.5-1］指令行操作過程示例。

1.5CommandHistory和實錄指令diary

1.5.1CommandHistory歷史指令窗簡介

JCommandHistoryn回區(qū)

FileEditViewWebWindowHelp

A(：)=-4：5喀用“全元素”法向AMMH

S0=2*(rtndQU,M)<0.5)-1

disp(ii*),disp(jj')

R(s)=lll

andfstate',23);

S0=2*(rand(Ns,M)<D.5)-1

->

囪1.5-1

1.5.2歷史指令行的再運行

【例152-1】再運行圖1.5-2所示歷史指令窗中的三行指令。

1.5.3指令窗實錄指令diary

1.6CurrentDirectory＞路徑設(shè)置器和文件管理

1.6.1CurrentDirectory當(dāng)前目錄瀏覽器簡介

圖1.6-1

一用戶目錄和當(dāng)前目錄設(shè)置

二借助當(dāng)前目錄瀏覽器獲取M和MAT文件信息

【例161-1］從圖1.2-1所示MATLAB缺省桌面開始，敘述引出圖1.6-1所示面貌的當(dāng)前目

錄瀏覽器的操作過程。

1.6.2MATLAB的搜索路徑

1.6.3MATLAB搜索路徑的擴展和修改

一何時需要修改搜索路徑

二利用設(shè)置路徑對話框修改搜索路徑

圖1.6-2

三利用指令path設(shè)置路徑

1.7WorkspaceBrowser和ArrayEditor

1.7.1WorkspaceBrowser工作空間瀏覽器簡介

1.7.2現(xiàn)場菜單用于內(nèi)存變量的查閱和刪除

一內(nèi)存變量查閱、刪除的指令操作法

【例1.7.2-I］在指令窗中運用who,whos查閱MATLAB內(nèi)存變量。

who

Yourvariablesare:

ABnumberDRXYZZ

BCDDXYx

whos

NameSizeBytesClass

A2x2230cellarray

Blxl264structarray

Bnumberlxl8doublearray

C2x2408symobject

D1x24chararray

DD2x28chararray

R33x338712doublearray

X33x338712doublearray

XYZ33x33x326136doublearray

Y33x338712doublearray

Z33x338712doublearray

X1x33264doublearray

y33x1264doublearray

Grandtotalis7722elementsusing62434bytes

【例1.722］在指令窗中運用clear指令可以刪除內(nèi)存中的變量。

clearBnumber

who

Yourvariablesare:

ABCDDDRXXYZYZxy

二內(nèi)存變量查閱和刪除的現(xiàn)場菜單操作法

1x18doublearray

旺/oublearray_____________

HROpen...k>ublearray

曲Thetaaplot

surf

2-DGraphics

3-DGraphicsmesh

Special2-DGraphics?surf

Special25-DGraphics?

Copy

Delete

ClearWorkspace

Rename

S1.7-2

【例1.723】通過“工作空間瀏覽器”的運作，采用圖形顯示內(nèi)存變量Z。

圖1.7-3

【例1.7.2-4］通過“工作空間瀏覽器”刪除內(nèi)存變量。

1.7.3ArrayEditor數(shù)組編輯器和大數(shù)組的輸入

1.7.4數(shù)據(jù)文件的存取

一存取數(shù)據(jù)文件的指令操作法

二通過內(nèi)存變量瀏覽器實現(xiàn)數(shù)據(jù)文件的存取

(1)產(chǎn)生保存全部內(nèi)存變量的數(shù)據(jù)文件的操作方法

閔1./一D

【例1.741］數(shù)據(jù)的存取。(假定內(nèi)存中已經(jīng)存在變量X,Y,Z)

(1)

mkdir('c:\\*my_dir');

cdc:\my_dir

savesafXYZ

dir

...saf.mat

(2)

clear

loadsafZ

who

Yourvariablesare:

Z

1.8LaunchPad交互界面分類目錄窗

^LaunchPad區(qū)|

圖1.8-1

1.9Editor/Debugger和腳本編寫初步

1.9.1Editor/DebuggerM文件編輯調(diào)試器簡介

一編輯調(diào)試器的開啟

小D:\MATLAB6p5\work\tryl.m匚］叵區(qū)］

FileEditViewTextDebugBreakpoints￥ebWindowHelp

口￥川電c?Mf.桓6Osty?

1

2制1.3.4-1trvl?

3%這是筒單的演示程序

4

5t=0：pi/50:4*pi,碇義目燹量取值喇B

6v0=.xp(-t/3);li+H與自變量*8?的力數(shù)犯

7V=?xp(-t/3).?sin(3*t);Mi十篁與自芟量相應(yīng)的做犯

plot(t,V,*-r*.t.vO.,K.t.-vO/:b')迎不同苗@、線型繪制曲線

scriptLn2COM4

囹1.y-1

二編輯器使用中的若干注意事項

1.9.2M腳本文件編寫初步

【例192-1］編寫解算例134-1題目的M腳本文件，并運行之。

操作步驟：

1.10幫助系統(tǒng)

l.io.l幫助方式概述

一“純文本”幫助

【例1.10.1-1]在指令窗中運行help的示例。

(1)

helphelp

HELPOn-linehelp,displaytextatcommandline.

HELP,byitself,listsallprimaryhelptopics.Eachprimarytopic

correspondstoadirectorynameontheMATLABPATH.

(2)

help

HELPtopics:

matlab\general-Generalpurposecommands.

matlab\ops-Operatorsandspecialcharacters.

matlab\lang-Programminglanguageconstructs.

matlab\elmat-Elementarymatricesandmatrixmanipulation.

matlab\elfun-Elementarymathfunctions.

Formorehelpondirectory/topic,type"helptopic".

(3)

helpelmat

Elementarymatricesandmatrixmanipulation.

Elementarymatrices.

zeros-Zerosarray.

ones-Onesarray.

eye-Identitymatrix.

(4)

helpeye

EYEIdentitymatrix.

EYE(N)istheN-by-Nidentitymatrix.

EYE(M,N)orEYE([M,N])isanM-by-Nmatrixwith1'son

thediagonalandzeroselsewhere.

EYE(SIZE(A))isthesamesizeasA.

SeealsoONES,ZEROS,RAND,RANDN.

【例1.10.1-2]在指令窗中，運用lookfor找Hl行(M函數(shù)文件的第一注釋行)

lookforfourier

FFTDiscreteFouriertransform.

FFT2Two-dimensionaldiscreteFourierTransform.

FFTNN-dimensionaldiscreteFourierTransform.

IFFTInversediscreteFouriertransform.

IFFT2Two-dimensionalinversediscreteFouriertransform.

IFFTNN-dimensionalinversediscreteFouriertransform.

XFOURIERGraphicsdemoofFourierseriesexpansion.

MOT563_FFTDiscreteFouriertransform.

MOT563_IFFTInversediscreteFouriertransform.

MOT566_FFTDiscreteFouriertransform.

MOT566_IFFTInversediscreteFouriertransform.

DFTMTXDiscreteFouriertransformmatrix.

INSTDFFTInversenon-standard1-DfastFouriertransform.

NSTDFFTNon-standard1-DfastFouriertransform.

FFTQuantizedFastFourierTransform.

FOURIERFourierintegraltransform.

IFOURIERInverseFourierintegraltransform.

二“導(dǎo)航/瀏覽器交互界面”幫助

三PDF幫助

四其他幫助

1.10.2HelpNavigator/Browser幫助導(dǎo)航/瀏覽器簡介

一Contents幫助文件目錄窗

【例1.10.2-1］通過鼠標(biāo)操作獲得如圖1.10-1所示的界面。

-Index幫助索引窗

【例1.10.2-21利用Idex搜索fourier這條術(shù)語。（注意把本例與例1.10.1-2,例1.10.2-3比

較。）

Productfilter:?AllOSelected

ContentsMexSearchDemosFavorites

Searchindexfor:

fourier

FProduct

A

fourierSymbolicMathToolbox

FourieranalysisMATLAB-Mathematics

basicfunctionWaveletToolbox

conceptsMATU\B-Mathematics

introductionWaveletToolbox

short-timeanalysis(S...WaveletToolboxV

<>

圖1.10-2

三Search搜索窗

【例1.10.2-3】利用"Search”窗搜索詞匯fourier。（注意把本例與例1.10.1-2、例1.1022

比較。）

Productfilter:⑥AllOSelected

Contents

Searchtype:

Searchfor:

TitleSection

Installer

MATLABDirectoiyStructureMATUkBInstallation

MATLAB

DataAnalysisandFourierTransformsFunctions-ByCate；

FourierAnalysisandtheFastFourierTra...DataAnalysisandE

fflMATU\BFunctionsv

<>

195pagescontaintheword:fourier

圖1.10-3

四Favorites書簽窗

第二章數(shù)值數(shù)組及其運算

數(shù)值數(shù)組(NumericArray)和數(shù)組運算(ArrayOperations)始終是MATLAB的核心內(nèi)

容。自MATLAB5.X版起，由于其“面向?qū)ο蟆钡奶卣?，這種數(shù)值數(shù)組(以下簡稱為數(shù)組)

成為了MATALB最重要的一種內(nèi)建數(shù)據(jù)類型(Built-inDataType),而數(shù)組運算就是定義

在這種數(shù)據(jù)結(jié)構(gòu)上的方法(Method)。

本章系統(tǒng)闡述：一、二維數(shù)值數(shù)組的創(chuàng)建、尋訪；數(shù)組運算和矩陣運算的區(qū)別；實現(xiàn)數(shù)

組運算的基本函數(shù)；多項式的表達(dá)、創(chuàng)建和操作；常用標(biāo)準(zhǔn)數(shù)組生成函數(shù)和數(shù)組構(gòu)作技法；

高維數(shù)組的創(chuàng)建、尋訪和操作；非數(shù)NaN、“空”數(shù)組概念和應(yīng)用；關(guān)系和邏輯操作。

順便指出：(1)本章所涉內(nèi)容和方法，不僅使用于數(shù)值數(shù)組，而且也將部分地延伸使

用于在其他數(shù)據(jù)結(jié)構(gòu)中。(2)MATLAB5.X和6.x版在本章內(nèi)容上的差異極微。(3)

MATLAB6.5版新增的兩種邏輯操作，在第2.13.2節(jié)給予介紹。

2.1引導(dǎo)

【例2.1-1】繪制函數(shù)y=xe-在0WxW1時的曲線。

x=0:0.1:1

y=x.*exp(-x)

1111

plot(xzy),xlabel('x),ylabel(y),title(y=x*exp(-x)')

x=

Columns1through7

00.10000.20000.30000.40000.50000.6000

Columns8through11

0.70000.80000.90001.0000

y=

Columns1through7

00.09050.16370.22220.26810.30330.3293

Columns8through11

0.34760.35950.36590.3679

圖2.1-1

2.2一維數(shù)組的創(chuàng)建和尋訪

2.2.1一維數(shù)組的創(chuàng)建

2.2.2一維數(shù)組的子數(shù)組尋訪和賦值

【例222-1】子數(shù)組的尋訪(Address)o

rand('state1,0)

x=rand(lz5)

X=

0.95010.23110.60680.48600.8913

x(3)

ans=

0.6068

x([l25])

ans=

0.95010.23110.8913

x(l:3)

ans=

0.95010.23110.6068

x(3:end)%

ans=

0.60680.48600.8913

x(3:-l:l)%

ans=

0.60680.23110.9501

x(find(x>0.5))

ans=

0.95010.60680.8913

x([12344321])

ans=

Columns1through7

0.95010.23110.60680.48600.48600.60680.2311

Column8

0.9501

【例222-2］子數(shù)組的賦值(Assign)o

x(3)=0

x=

0.95010.231100.48600.8913

x([l4])=[11]

x=

1.00000.231101.00000.8913

2.3二維數(shù)組的創(chuàng)建

2.3.1直接輸入法

【例2.3.1J】在MATLAB環(huán)境下，用下面三條指令創(chuàng)建二維數(shù)組C。

a=2.7358;b=33/79;

C=[l,2*a+i*b,b*sqrt(a);sin(pi/4)za+5*b,3.5+i]

C=

1.00005.4716+0.4177i0.6909

0.70714.82443.5000+l.OOOOi

2

【例2.3.1-2］復(fù)數(shù)數(shù)組的另一種輸入方式。

M_r=[l,2,3;4,5,6],M_i=[ll,12,13;14,15,16]

CN=M_r+i*M_i

M_r=

123

456

Mi=

111213

141516

CN=

1.0000+11.0000i2.0000+12.0000i3.0000+13.0000i

4.0000+14.0000i5.0000+15.0000i6.0000+16.0000i

2.3.2利用M文件創(chuàng)建和保存數(shù)組

【例2.3.2-1]創(chuàng)建和保存數(shù)組AM的MyMatrix.m文件。

(1)

%MyMatrix.mCreationandpreservationofmatrixAM

AM=[101,102,103,104,105,106,107,108,109;...

201,202,203,204,205,206,207,208,209;...

301,302,303,304,305,306,307,308,309];

(2)

(3)

2.4二維數(shù)組元素的標(biāo)識

2.4.1“全下標(biāo)”標(biāo)識

2.4.2“單下標(biāo)”標(biāo)識

2.4.3“邏輯1”標(biāo)識

丫一4-2024/

【例243-1】找出數(shù)組A8中所有絕對值大于3的元素。

<-3-1135f

A=zeros(2,5);

A(：)=-4:5

L=abs(A)>3

islogical(L)

X=A(L)

A=

-4-2024

-3-1135

10001

00001

ans=

1

X=

-4

4

5

【例243-2】演示邏輯數(shù)組與一般雙精度數(shù)值數(shù)組的關(guān)系和區(qū)別。（本例在例243-1基礎(chǔ)

3

上進(jìn)行)。

(1)

Num=[1,0,0,0,1;0,0,0,0,1];

N_L=Num==L

c_N=class(Num)

c_L=class(L)

N_L=

11111

11111

C_N=

double

c_L=

double

(2)

islogical(Num)

Y=A(Num)

ans=

0

???Indexintomatrixisnegativeorzero.Seereleasenotesonchanges

to

logicalindices.

2.5二維數(shù)組的子數(shù)組尋訪和賦值

【例2.5-1］不同賦值方式示例。

A=zeros(2,4)

A=

0000

0000

A(:)=l:8

A=

1357

2468

s=[235];

A(s)

Sa=[102030]1

A(s)=Sa

ans=

235

Sa.=

10

20

30

A=

12030

10468

A(:,[23])=ones(2)

A=

1117

10118

2.6執(zhí)行數(shù)組運算的常用函數(shù)

4

2.6.1函數(shù)數(shù)組運算規(guī)則的定義:

2.6.2執(zhí)行數(shù)組運算的常用函數(shù)

【例2.621】演示pow2的數(shù)組運算性質(zhì)。

A=[l:4;5:8]

A=

1234

5678

pow2(A)

ans=

24816

3264128256

2.7數(shù)組運算和矩陣運算

2.7.1數(shù)組運算和矩陣運算指令對照匯總

【例271-1】兩種不同轉(zhuǎn)置的比較

clear;A=zeros(2,3);

A(:)=1:6;

A=A*(1+i)

A_A=A.T

A_M=A'

A=

1.0000+l.OOOOi3.0000+3.0000i5.0000+5.0000i

2.0000+2.0000i4.0000+4.0000i6.0000+6.0000i

A_A=

1.0000+l.OOOOi2.0000+2.0000i

3.0000+3.0000i4.0000+4.0000i

5.0000+5.0000i6.0000+6.0000i

A_M=

1.00001.000012.00002.0000i

3.00003.0000i4.00004.0000i

5.00005.0000i6.00006.0000i

2.8多項式的表達(dá)方式及其操作

2.8.1多項式的表達(dá)和創(chuàng)建

-多項式表達(dá)方式的約定

二多項式行向量的創(chuàng)建方法

【例281.2-1】求3階方陣A的特征多項式。

A=[ll1213;141516;171819];

PA=poly(A)

PPA=poly2str(PA,*s1)

PA=

1.0000-45.0000-18.00000.0000

PPA=

sA3-45sA2-18s+1.8303e-014

【例281.2-2】由給定根向量求多項式系數(shù)向量。

R=[-0.5,-0.3+0.4*i,-0.3-0.4*i];

5

P=poly(R)

PR=real(P)

PPR=poly2str(PRZ*x*)

P=

1.00001.10000.55000.1250

PR=

1.00001.10000.55000.1250

PPR=

xA3+1.1xA2+0.55x+0.125

2.8.2多項式運算函數(shù)

(s~+2)(s+4)(s4-1)

【例2.8.2-1】求的“商”及“余”多項式。

5+5+1

pl=conv([lz0,2]zconv([lz4]z[1,1]));

p2=[l011];

[q,r]=deconv(pl,p2);

cq=，商多項式為1;余多項式為,;

1

disp([cqzpoly2str(q,*s)])zdisp([cr,poly2str(r,'s')])

商多項式為s+5

余多項式為5sA2+4s+3

【例2.822】兩種多項式求值指令的差別。

S=pascal(4)

P=poly(S);

T

PP=poly2str(Pz's)

PA=polyval(P,S)

PM=polyvalm(P,S)

s=

1111

1234

13610

141020

pp=

A

s人4-29s人3+72s2-29s+1

PA=

1.0e+004*

0.00160.00160.00160.0016

0.00160.0015-0.0140-0.0563

0.0016-0.0140-0.2549-1.2089

0.0016-0.0563-1.2089-4.3779

PM=

1.0e-010★

0.00160.00330.00900.0205

0.00450.01010.02860.0697

0.00950.02100.06530.1596

0.01630.03870.12260.3019

【例2.8.2?3】部分分式展開。

a=[lz3z4z2,7z2];

b=[3,2,5,4,6];

[r,szk]=residue(b,a)

r=

1.1274+1.1513i

1.1274-1.1513i

-0.0232-0.0722i

-0.0232+0.0722i

0.7916

s=

6

-1.7680+1.2673i

-1.7680-1.2673i

0.4176+1.1130i

0.4176-1.1130i

-0.2991

k=

[]

2.9標(biāo)準(zhǔn)數(shù)組生成函數(shù)和數(shù)組操作函數(shù)

2.9.1標(biāo)準(zhǔn)數(shù)組生成函數(shù)

【例291-1]標(biāo)準(zhǔn)數(shù)組產(chǎn)生的演示。

ones(1,2)

ans=

11

ones(2)

ans=

11

11

randn(*state,z0)

randn(2,3)

ans=

-0.43260.1253-1.1465

-1.66560.28771.1909

D=eye(3)

D=

100

010

001

diag(D)

ans=

1

1

1

diag(diag(D))

ans=

100

010

001

repmat(D,1,3)

ans=

Columns1through8

10010010

01001001

00100100

Column9

0

0

1

7

2.9.2數(shù)組操作函數(shù)

【例2.9.2-1]diag與reshape的使用演示。

a=-4:4

A=reshape(a,3,3)

a=

Columns1through8

-4-3-2-10123

Column9

4

A=

-4-12

-303

-214

al=diag(A,1)

al=

-1

3

Al=diag(al,-1)

Al=

000

-100

030

【例292-2】數(shù)組轉(zhuǎn)置、對稱交換和旋轉(zhuǎn)操作后果的對照比較。

A

A=

-4-12