MATLAB 基础培训教程

MATLAB基础教程

第一章基础准备及入门3

1.1MATLAB的安装和内容选择3

1.2Desktop操作桌面的启动3

1.3CommandWindow运行入4

1.4CommandWindow操作要旨8

1.5CommandHistory和实录指令diary8

1.6CurrentDirectory路径设置器和文件管理9

1.7WorkspaceBrowser和ArrayEditor10

1.8LaunchPad交互界面分类目录窗13

1.9Editor/Debugger和脚本编写初步14

1.10帮助系统14

第二章数值数组及其运算18

1.11引导18

1.12一维数组的创建和寻访19

1.13二维数组的创建20

1.14二维数组元素的标识20

1.15二维数组的子数组寻访和赋值21

1.16执行数组运算的常用函数22

1.17数组运算和矩阵运算22

1.18多项式的表达方式及其操作23

1.19标准数组生成函数和数组操作函数24

1.20数组构作技法综合26

1.21高维数组29

1.22“非数”和“空”数组33

1.23关系操作和逻辑操作36

第三章字符串、元胞和构架数组38

1.24字符串数组39

1.25元胞数组43

1.26构架数组46

1.27关于数据类型的归纳性说明53

第四章数值计算54

1.28引言54

1.29LU分解和恰定方程组的解55

1.30矩阵特征值和矩阵函数56

1.31奇异值分解59

1.32函数的数值导数和切平面60

1.33函数的零点62

1.34函数极值点65

1.35数值积分66

1.36随机数据的统计描述71

1.37多项式拟合和非线性最小二乘74

1.38插值和样条79

1.39样条函数及其应用80

1.40Fourier分析82

1.41常微分方程87

1.42稀疏矩阵94

第五章符号计算96

1.43符号对象和符号表达式96

1.44符号表达式和符号函数的操作99

1.45符号微积分103

1.46符号积分变换106

1.47符号代数方程的求解108

1.48符号微分方程的求解110

1.49利用MAPLE的深层符号计算资源111

1.50可视化数学分析界面113

第一章基础准备及入门

本章有两个目的：一是讲述MATLAB正常运行所必须具备的基础条件；二是简明系统

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

本章的前两节分别讲述：MATLAB的正确安装方法和MATLAB环境的启动。因为指

令窗是MATLAB最重要的操作界面，所以本章用第1.3、1.4两节以最简单通俗的叙述、算

例讲述指令窗的基本操作方法和规则。这部分内容几乎对MATLAB各种版本都适用。

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

式工作界面。了解、熟悉和掌握这些交互界面的基本功能和操作方法，将使新老用户能事

半功倍地利用MATLAB去完成各种学习和研究。为此，本章特设几节用于专门介绍最常用

的交互界面：历史指令窗、当前目录浏览器、工作空间浏览器、内存数组编辑器、交互界

面分类目录窗、M文件编辑/调试器、及帮助导航/浏览器。

本章是根据MATLAB6.5版编写的，但大部分内容也适用于其他6.x版。

1.1MATLAB的安装和内容选择

ProductList

1.Selectdirectorywhereproductsvillbeinstalled.Spaceavailable：

|D'MATLAB6P5

7523H

2.S«l«ctinstallationoption*.

Spasrt<nir«d

•Installproductsanddocun«ntttionforproductG),

「•InstallproductsonlyHoreOptions...ifany

「Installdocumentationonly.1115M

3.Selectlanguageofdocunentation.

required

GEnglishonly.

docnnentation,if

^EnglishandJapanese,ifavailable.

any-

433H

4.Selectproductsand/ordocunent&tion.

7MATLAB6.5

QSinulink5.0

BAerospaceBlockset1.01Toldsp«c«

GCDMAReferenceBLockset1.1r«qu>r«d

rIConnunxcationsBlockset

FConnunicationsToolbox211548M

图1.1-1

1.2Desktop操作桌面的启动

1.2.1MATLAB的启动

1.2.2Desktop操作桌面简介

一操作桌面的缺省外貌

图1.2-1

二通用操作界面

1.3CommandWindow运行入门

1.3.1CommandWindow指令窗简介

图1.3-1

1.3.2最简单的计算器使用法

【例1.321］求［12+2x(7—4)］+32的算术运算结果。

(1)用键盘在MATLAB指令窗中输入以下内容

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

(2)在上述表达式输入完成后，按【Enter】键，该就指令被执行。

(3)在指令执行后，MATLAB指令窗中将显示以下结果。

ans=

2

123

【例1.3.2-2】简单矩阵/1=456的输入步骤。

789

(1)在键盘上输入下列内容

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

(2)按［Enter］键，指令被执行。

(3)在指令执行后，MATLAB指令窗中将显示以下结果:

A=

123

456

789

【例132-3】矩阵的分行输入。

A=[l,2,3

4,5,6

7,8,9]

A=

123

456

789

【例1.324】指令的续行输入

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

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

S=

0.6345

1.3.3数值、变量和表达式

一数值的记述

二变量命名规则

三MATLAB默认的预定义变量

四运算符和表达式

五复数和复数矩阵

【例1.33-1］复数Zi=3+4，,3=1+万,5=2eN表达,及计算z=4

4

(1)

zl=3+4i

zl=

3.0000+4.OOOOi

(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

【例1.3.3-2]复数矩阵的生成及运算

A=[lf3;2,4]-[5z8;6z9]*i

B=[l+5i,2+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

1.0000+5.0000i2.0000+6.0000i

3.0000+8.0000i4.0000+9.0000i

.0e+002★

0.99001.1600-0.0900i

1.1600+0.0900i1.3700

【例133-3】求上例复数矩阵C的实部、虚部、模和相角。

C_real=real(C)

C_imag=imag(C)

C_magnitude=abs(C)

C_phase=angle(C)*180/pi%以度为单位计算相角

C_real=

99116

116137

C_imag=

0-9

90

C_magnitude

99.0000116.3486

116.3486137.0000

C_phase=

0-4.4365

4.43650

【例1.3.3-4】用MATLAB计算超二W能得到-2吗？

(1)

a=-8;

r=aA(1/3)

r=

1.0000+1.73211

2

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.732H

(3)

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

1

plot(xzyz'b:),grid

holdon

1111

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

1111

plot(rrr([2,3])z'o',MarkerSize,,15,Color,1b)

axis([-3Z3,-3,3])zaxissquare

holdoff

图L3-2

1.3.4计算结果的图形表示

【例1.341】画出衰减振荡曲线y=J§sin3f及其它的包络线汽f的取值范围是

[0,4力。(图1.3-3)

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

y0=exp(-t/3);

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

plot(t,yJ:b11,-yOJ:b')

grid

图1.3-3

【例1.3.4-2】画出」=也史，2，所表示的三维曲面(图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('*'),ylabel(*y'),zlabel('z')

1.4CommandWindow操作要旨

1.4.1指令窗显示方式的操作

一缺省显示方式

二显示方式的设置

1.4.2数值计算结果的显示格式

1.4.3指令行中的标点符号

1.4.4指令窗的常用控制指令

1.4.5指令窗中指令行的编辑

【例145-1】指令行操作过程示例。

1.5CommandHistory和实录指令diary

1.5.1CommandHistory历史指令窗简介

图1.5-1

1.5.2历史指令行的再运行

【例1.521】再运行图1.5-2所示历史指令窗中的三行指令。

图1.5-2

1.5.3指令窗实录指令diary

1.6CurrentDirectory＞路径设置器和文件管理

1.6.1CurrentDirectory当前目录浏览器简介

图1.6-1

一用户目录和当前目录设置

二借助当前目录浏览器获取M和MAT文件信息

【例1.6.1-1］从图1.2-1所示MATLAB缺省桌面开始，叙述引出图1.6-1所示面貌的当前目

录浏览器的操作过程。

1.6.2MATLAB的搜索路径

1.6.3MATLAB搜索路径的扩展和修改

一何时需要修改搜索路径

二利用设置路径对话框修改搜索路径

图1.6-2

三利用指令path设置路径

1.7WorkspaceBrowser和ArrayEditor

1.7.1WorkspaceBrowser工作空间浏览器简介

图1.7-1

1.7.2现场菜单用于内存变量的查阅和删除

一内存变量查阅、删除的指令操作法

【例1.721］在指令窗中运用who,whos查阅MATLAB内存变量。

who

Yourvariablesare:

ABnumberDRXYZZy

BCDDXYX

whos

NameSizeBytesClass

A2x2230cellarray

Blxl264structarray

Bnumberlxl8doublearray

C2x2408symobject

D1x24chararray

DD2x28chararray

R33x338712doublearray

X33x338712doublearray

XYZ33x33x326136doublearray

Y33x338712doublearray

Z33x338712doublearray

X1x33264doublearray

y33x1264doublearray

Grandtotalis7722elementsusing62434bytes

【例172-2］在指令窗中运用clear指令可以删除内存中的变量。

clearBnumber

who

Yourvariablesare:

ABCDDDRXXYZYZxy

二内存变量查阅和删除的现场菜单操作法

I

SAlixi8doublearray

■

■

■：：E_______100x10080000doublearray

HROpen...oublearray

I4tiThetaaplot

surf

SelectAll

2-DGraphics

ImportData...3-DGraphicsmesh

SaveSelectionAs...Special2-DGraphics►surf

SaveWorkspaceAs...Special2.5-DGraphics►

Copy

Delete

ClearWorkspace

Rename

图1.7-2

【例1.723】通过“工作空间浏览器”的运作，采用图形显示内存变量Z。

图1.7-3

【例1.7.2-4】通过“工作空间浏览器”删除内存变量。

1.7.3ArrayEditor数组编辑器和大数组的输入

图1.7-4

1.7.4数据文件的存取

一存取数据文件的指令操作法

二通过内存变量浏览器实现数据文件的存取

(1)产生保存全部内存变量的数据文件的操作方法

图1.7-5

【例174-1］数据的存取。(假定内存中已经存在变量X,Y,Z)

(1)

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

cdc:\my_dir

savesafXYZ

dir

・・,saf.mat

(2)

clear

loadsafZ

who

Yourvariablesare:

Z

1.8LaunchPad交互界面分类目录窗

图1.8-1

1.9Editor/Debugger和脚本编写初步

1.9.1Editor/DebuggerM文件编辑调试器简介

一编辑调试器的开启

图1.9-1

二编辑器使用中的若干注意事项

1.9.2M脚本文件编写初步

【例1.9.2-1］编写解算例134-1题目的M脚本文件，并运行之。

操作步骤：

1.10帮助系统

1.10.1帮助方式概述

一“纯文本”帮助

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

(1)

helphelp

HELPOn-linehelp,displaytextatcommandline.

HELP,byitself,listsallprimaryhelptopics.Eachprimarytopic

correspondstoadirectorynameontheMATLABPATH.

(2)

help

HELPtopics:

matlab\general-Generalpurposecommands.

matlab\ops-Operatorsandspecialcharacters.

matlabMang-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.

1

EYE(MzN)orEYE([MzN])isanM-by-Nmatrixwith1son

thediagonalandzeroselsewhere.

EYE(SIZE(A))isthesamesizeasA.

SeealsoONES,ZEROS,RAND,RANDN.

【例1.10.1-2]在指令窗中，运用lookfor找Hl行(M函数文件的第一注释行)

lookforfourier

FFTDiscreteFouriertransform.

FFT2Two-dimensionaldiscreteFourierTransform.

FFTNN-dimensionaldiscreteFourierTransform.

IFFTInversediscreteFouriertransform.

IFFT2Two-dimensionalinversediscreteFouriertransform.

工FFTNN-dimensionalinversediscreteFouriertransform.

XFOURIERGraphicsdemoofFourierseriesexpansion.

MOT563_FFTDiscreteFouriertransform,

MOT563_IFFTInversediscreteFouriertransform.

MOT566_FFTDiscreteFouriertransform,

MOT566IFFTInversediscreteFouriertransform.

DFTMTXDiscreteFouriertransformmatrix.

INSTDFFTInversenon-standard1-DfastFouriertransform.

NSTDFFTNon-standard1-DfastFouriertransform.

FFTQuantizedFastFourierTransform.

FOURIERFourierintegraltransform.

IFOURIERInverseFourierintegraltransform.

二“导航/浏览器交互界面”帮助

三PDF帮助

四其他帮助

1.10.2HelpNavigator/Browser帮助导航/浏览器简介

图1.10-1

一Contents帮助文件目录窗

【例1.10.2-1】通过鼠标操作获得如图1.10-1所示的界面。

二Index帮助索引窗

【例1.10.2-2】利用Idex搜索fourier这条术语。（注意把本例与例例1.10.2-3比

较。）

Productfilter:⑥AllOSelected

ContentsIndexSearchDemosFavorites

Searchindexfor:

fourier

FProduct

人

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.10.2-2

比较。）

Productfilter:(g)AllOSelected

TitleSection

Installer

MATLABDirectoryStructureMATLABInstallation

MATLAB

DataAnalysisandFourierTransformsFunctions-ByCate；

FourierAnalysisandtheFastFourierTra...DataAnalysisandS

fftMATLABFunctionsv

195pagescontaintheword:fourier

图1.10-3

四Favorites书签窗

第二章数值数组及其运算

数值数组（NumericArray）和数组运算（ArrayOperations）始终是MATLAB的核心内

容。自MATLAB5.X版起，由于其“面向对象”的特征，这种数值数组（以下简称为数组）

成为了MATALB最重要的一种内建数据类型（Built-inDataType）,而数组运算就是定义

在这种数据结构上的方法（Method）。

本章系统阐述：、二维数值数组的创建、寻访；数组运算和矩阵运算的区别；实现

数组运算的基本函数；多项式的表达、创建和操作；常用标准数组生成函数和数组构作技

法；高维数组的创建、寻访和操作：非数NaN、“空”数组概念和应用；关系和逻辑操作。

顺便指出：（1）本章所涉内容和方法，不仅使用于数值数组，而且也将部分地延伸使

用于在其他数据结构中。（2）MATLAB5.X和6.x版在本章内容上的差异极微。（3）

MATLAB6.5版新增的两种逻辑操作，在第2.13.2节给予介绍。

1.11引导

【例2.1-1］绘制函数y=xe'x在0WxW1时的曲线。

x=0:0.1:1

y=x.*exp(-x)

11111

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

Columns1through7

00.09050.16370.22220.26810.30330.3293

Columns8through11

0.34760.35950.36590.3679

1.12一维数组的创建和寻访

1.12.1一维数组的创建

1.12.2一维数组的子数组寻访和赋值

【例2.2.2-1]子数组的寻访(Address)。

rand('state',o)

x=rand(1,5)

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］子数组的赋值(Assign)o

x(3)=0

x=

0.95010.231100.48600.8913

x([l4])=[11]

x=

1.00000.231101.00000.8913

1.13二维数组的创建

1.13.1直接输入法

【例2.3.1-1]在MATLAB环境下，用下面三条指令创建二维数组C。

a=2.7358;b=33/79;

C=[lz2*a+i*bzb*sqrt(a);sin(pi/4),a+5*bz3.5+i]

C=

1.00005.4716+0.4177i0.6909

0.70714.82443.5000+l.OOOOi

【例231-2]复数数组的另一种输入方式。

M_r=[1,2,3;4,5,6]12z13;14,15,16]

CN=M_r+i*M_i

M_r=

123

456

111213

141516

CN=

1.0000+11.OOOOi2.0000+12.0000i3.0000+13.0000i

4.0000+14.0000i5.0000+15.OOOOi6.0000+16.OOOOi

1.13.2利用M文件创建和保存数组

【例2.3.2-1]创建和保存数组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)

1.14二维数组元素的标识

1.14.1“全下标”标识

1.14.2“单下标”标识

1.14.3“逻辑1”标识

-20

【例243-1】找出数组1=中所有绝对值大于3的元素。

-11

A=zeros(2,5);

A(:)=-4:5

L=abs(A)>3

islogical(L)

X=A(L)

A=

-4-2024

-3-1135

L=

10001

00001

ans=

1

X=

-4

4

5

【例2.4.3-2J演示逻辑数组与一般双精度数值数组的关系和区别。(本例在例2.4.3-1基础

上进行)。

(1)

Num=[lz0z0,0,l;0z0z0,0zl];

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.

1.15二维数组的子数组寻访和赋值

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

A=zeros(2,4)

A

0000

0000

A(:)=l:8

A=

1357

2468

s=[235];

A(s)

Sa=[102030]

A(s)=Sa

ans=

235

Sa=

10

20

30

A=

120307

10468

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

A=

1117

108

1.16执行数组运算的常用函数

1.16.1函数数组运算规则的定义:

1.16.2执行数组运算的常用函数

【例2.621】演示pow2的数组运算性质。

A=[l:4;5:8]

A=

1234

5678

pow2(A)

ans=

24816

3264128256

1.17数组运算和矩阵运算

1.17.1数组运算和矩阵运算指令对照汇总

【例271-1】两种不同转置的比较

clear;A=zeros(2,3);

A(:)=1:6;

A=A*(1+i)

A_A=A.1

A_M=A'

A=

1.0000+1.OOOOi3.0000+3.OOOOi5.0000+5.OOOOi

2.0000+2.OOOOi4.0000+4.OOOOi6.0000+6.OOOOi

A_A=

1.0000+1.OOOOi2.0000+2.OOOOi

3.0000+3.OOOOi4.0000+4.OOOOi

5.0000+5.0000i6.0000+6.0000i

A_M=

1.0000-l.OOOOi2.0000-2.0000i

3.0000-3.0000i4.0000-4.OOOOi

5.0000-5.OOOOi6.0000-6.OOOOi

1.18多项式的表达方式及其操作

1.18.1多项式的表达和创建

一多项式表达方式的约定

二多项式行向量的创建方法

【例2.8.1.2-1]求3阶方阵A的特征多项式。

A=[ll1213;141516;171819];

PA=poly(A)

PPA=poly2str(PA,*s*)

PA=

1.0000-45.0000-18.00000.0000

PPA=

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

【例2.8.1.2・2】山给定根向量求多项式系数向量。

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

P=poly(R)

PR=real(P)

PPR=poly2str(PR,*x1)

P=

1.00001.10000.55000.1250

PR=

1.00001.10000.55000.1250

PPR=

xA3+1.1xA2+0.55x+0.125

1.18.2多项式运算函数

【例2.821】求"2+2}5+4)($+1)的“商”及“余”多项式。

S+S+1

pl=conv([1,0,2],conv([1,4]z[1,1]));

p2=[l011];

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

商多项式为1;余多项式为1;

11

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

商多项式为s+5

余多项式为5sA2+4s+3

【例2.8.2・2】两种多项式求值指令的差别。

S=pascal(4)

P=poly(S);

PP=poly2str(P,1s*)

PA=polyval(P,S)

PM=polyvalm(P,S)

S=

1111

1234

13610

141020

pp

AA

s429s人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

1.0e-010*

0.00160.00330.00900.0205

0.00450.01010.02860.0697

0.00950.02100.06530.1596

0.01630.03870.12260.3019

【例2.823】部分分式展开。

a=[l,3,4,2,7,2];

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

[r,s,k]=residue(b,a)

r=

1.1274+1.1513i

1.1274-1.1513i

-0.0232-0.0722i

-0.0232+0.0722i

0.7916

s=

-1.7680+1.2673i

-1.7680-1.2673i

0.4176+1.1130i

0.4176-1.1130i

-0.2991

k=

[]

1.19标准数组生成函数和数组操作函数

1.19.1标准数组生成函数

【例291-1】标准数组产生的演示。

ones(1,2)

ans=

11

ones(2)

ans=

11

1

randn(*stateI0)

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,1z3)

ans=

Columns1through8

10010010

01001001

00100100

Column9

0

0

1

1.19.2数组操作函数

【例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】数组转置、对称交换和旋转操作后果的对照比较。

A

A=

-4-12

-303

-214

A.'

ans=

-4-3-2

-101

234

flipud(A)

ans=

-214

-303

-4-12

fliplr(A)

ans=

2-1-4

30-3

41-2

rot90(A)

ans=

234

-101

-4-3—2

【例2.9.2・3】演示Kronecker乘法不具备“可交换规律”。

B=eye(2)

C=reshape(1:4,2,2)

B=

10

01

C=

13

24

kron(B,C)

ans=00

1300

2413

0024

00

kron(C,B)

ans=30

1003

0140

2004

02

1.20数组构作技法综合

【例2.10-11数组的扩展。

(1)数组的赋值扩展法

A=reshape(1:9,3,3)

A=

147

258

369

A(5Z5)=111

A=147

25800

36900

00000

00000

0111

A(:,6)=222

A

(2)多次寻访扩展法

AA=A(:Z[1:6,1:6])

AA=

1470022214700222

2580022225800222

3690022236900222

0000022200000222

00001112220000111222

(3)合成扩展法

B=ones(2,6)

B=

111111

111111

r=A;B

ABr-

-147oo222

258oo222

369oo222

00