北邮 MATLAB编程基础

1、Electronic Engineer School of BUPT尹霄丽第2章 MATLAB编程基础EE of BUPTMATLAB在信号与系统课程中的应用主要内容l2.1 变量l2.2 数组，向量与矩阵的创建和访问l2.3 部分常用运算符l2.4 基本数学函数l2.5 基本绘图函数l2.6 字符与字符串l2.7 MATLAB 程序设计EE of BUPTMATLAB在信号与系统课程中的应用2.1 2.1 变量变量n标识符是标志变量名、常量名、函数名和文件名的字符串的总称。n变量和常量的标识符长度不超过31(6.5以后版本为63个)个字符。n标识符中的第一个字符必须是英文字母n标识符可以包含

2、下划线、数字，但不能为空格符、标点 。n函数和文件名通常不超过8（？）个字符；n大小写敏感；n变量无需定义即可使用。n只有一种数据格式：双精度（64位）二进制， 对应于10进制16位有效数字和308次幂。EE of BUPTMATLAB在信号与系统课程中的应用部分特殊变量和常数nans 最近生成的无名结果最近生成的无名结果neps 计算机的零阈值计算机的零阈值npi 3.14159265358979n i 虚数单位虚数单位n j 虚数单位虚数单位nInf 无穷大，例如无穷大，例如1/0 nNaN Not-a-Number，例如，例如0/0ninputname 输入变量名输入变量名n注：注：在在

3、MATLAB的的IEEE算法规则中，被算法规则中，被0除是允许的。这不除是允许的。这不会导致任何程序的中断，只是给出警告信息然后使用名称会导致任何程序的中断，只是给出警告信息然后使用名称Inf或或NaN来记述。这些名称可在后面步骤得到合理运用。来记述。这些名称可在后面步骤得到合理运用。EE of BUPTMATLAB在信号与系统课程中的应用2.2 数组，向量与矩阵的创建和访问矩阵的创建和赋值：矩阵的值放在方括号矩阵的值放在方括号中；中；矩阵的行间用分号矩阵的行间用分号”;”分隔；分隔；矩阵的列间用空格或逗号矩阵的列间用空格或逗号”,”分隔；分隔；直接输入法直接输入法：例如：创建一个：例如：创建

4、一个33的矩阵，输入的矩阵，输入 ： a=1 2 3;4 5 6;7 8 9 屏幕输出屏幕输出a = 1 2 3 4 5 6 7 8 9步长生长法步长生长法：初值初值:步长步长:终值终值，例如：，例如：t=0:1:10EE of BUPTMATLAB在信号与系统课程中的应用说明n在在MATLABMATLAB中，上面三个概念中，上面三个概念在创建和显示的时候没有在创建和显示的时候没有任何区别任何区别，可分为一维，二维和高维数组。，可分为一维，二维和高维数组。n矩阵的元素用圆括号矩阵的元素用圆括号“()”()”中的数字也称中的数字也称下标下标来注明；来注明；注意注意向量的第一个元素的向量的第一个元

5、素的下标是下标是1 1，而不是，而不是0 0。n如果赋值元素的下标超出了原来矩阵的大小，矩阵的行如果赋值元素的下标超出了原来矩阵的大小，矩阵的行列会自动扩展。列会自动扩展。 n语句的结尾可用回车符或逗号，可显示运算结果；如果语句的结尾可用回车符或逗号，可显示运算结果；如果不希望显示结果，则以分号结尾。不希望显示结果，则以分号结尾。EE of BUPTMATLAB在信号与系统课程中的应用复数的表示nz1=8+10i;%直接按直角坐标的方式输入直接按直角坐标的方式输入nz2=8+10*i;%运算符构成的直角坐标的方式输入运算符构成的直角坐标的方式输入nz3=25*exp(i*pi/3); %运算符

6、构成的极坐标的方式输运算符构成的极坐标的方式输入入nA=z1,z2,z3;%矩阵形式矩阵形式EE of BUPTMATLAB在信号与系统课程中的应用n访问矩阵的第访问矩阵的第3行第行第2列的元素：列的元素： a(3,2) ans= 8n访问矩阵某行访问矩阵某行(列）的所有元素：列）的所有元素：a(3,:) a(:,2)这两个符号分别是这两个符号分别是1*3行向量和行向量和3*1列向量列向量n确定矩阵的维度：确定矩阵的维度：nLENGTH Length of vector.nSIZE Size of array. (2) 矩阵下标及矩阵元素的访问b=1 2 3; 4 5 6length(b)le

7、ngth(b)n,m=size(b)EE of BUPTMATLAB在信号与系统课程中的应用由小矩阵生成大矩阵nf1=1 2 3;nf2=f1 f1 f1;nf3=f1;f1;f1;nf4=f3 f3;f3 fliplr(f3);EE of BUPTMATLAB在信号与系统课程中的应用(3) 创建矩阵的常用函数nzeros: 全全0矩阵矩阵nones: 全全1矩阵矩阵nmagic 魔方矩阵魔方矩阵nrand: 0 1之间的之间的均匀分布均匀分布的伪随机数的伪随机数nrandn: 均值为均值为0，方差为，方差为1的的正态分布正态分布n% hist distributionneye: 单位矩阵单位

8、矩阵nlinspace 线性分布线性分布nlogspace 对数分布对数分布 （以（以10为底）为底）n:EE of BUPTMATLAB在信号与系统课程中的应用nZEROS Zeros array.n ZEROS(N) is an N-by-N matrix of zeros.n ZEROS(M,N) or ZEROS(M,N) is an M-by-N matrix of zeros.nLINSPACE Linearly spaced vector.n LINSPACE(X1, X2) generates a row vector of 100 linearly equally space

9、d points between X1 and X2.n LINSPACE(X1, X2, N) generates N points between X1 and X2. For N fliplr(a)ans = 3 2 1diag(1 2 3)ans = 1 0 0 0 2 0 0 0 3a=1 2 3;4 5 6;7 8 9a = 1 2 3 4 5 6 7 8 9 diag(a)ans = 1 5 9EE of BUPTMATLAB在信号与系统课程中的应用举例a=1 2 3;4 5 6;7 8 9a = 1 2 3 4 5 6 7 8 9tril(a)ans = 1 0 0 4 5 0

10、 7 8 9triu(a)ans = 1 2 3 0 5 6 0 0 9a=1 2 3 4 5 6 7 8 9；reshape(a,3,3)ans = 1 4 7 2 5 8 3 6 9EE of BUPTMATLAB在信号与系统课程中的应用2.3 2.3 部分常用运算符部分常用运算符(1) 算术运算符算术运算符n 加加n 减减n * 乘乘(包括标量乘，矩阵乘，标量与矩阵乘包括标量乘，矩阵乘，标量与矩阵乘)n / 除除n标量除标量除n矩阵除标量矩阵除标量n右除：右除：X=B*inv(D)=B/D % inv() 为矩阵求逆运算为矩阵求逆运算n 左除：左除：X=inv(D)*B=DB n 矩阵求

11、幂（矩阵必须为方阵）矩阵求幂（矩阵必须为方阵） EE of BUPTMATLAB在信号与系统课程中的应用元素群运算元素群运算(array operation)(array operation)n.* 数组相乘数组相乘n ./ 数组相除数组相除n. Left array dividen . 数组求幂数组求幂注意：注意：.* ./ 表示两个表示两个同维同维数组中的对应元素做乘，除数组中的对应元素做乘，除 ； . 表示对数组的每个元素求幂表示对数组的每个元素求幂EE of BUPTMATLAB在信号与系统课程中的应用举例a=1 2 3;4 5 6;7 8 9b=ones(3,3)c=a*bd=a.*

12、bc = 6 6 6 15 15 15 24 24 24d = 1 2 3 4 5 6 7 8 9a.*a和a.2相同EE of BUPTMATLAB在信号与系统课程中的应用(2) 关系运算符n 小于小于n 大于大于n = 大于等于大于等于n = 等于等于n = 不等于不等于运算法则：若关系式成立，结果为运算法则：若关系式成立，结果为1； 若关系式不成立，结果为若关系式不成立，结果为0。EE of BUPTMATLAB在信号与系统课程中的应用举例a=1 2 3;4 5 6;7 8 9b=a5b = 0 0 0 0 0 1 1 1 1EE of BUPTMATLAB在信号与系统课程中的应用(3)

13、 逻辑运算符n& 与与 n | 或或n 非非运算法则：运算法则： 若逻辑真，结果为若逻辑真，结果为1； 若逻辑假，结果为若逻辑假，结果为0。a=1 2 3;4 5 6;7 8 9b=a5c=a8d=b&ce=b|cf=bEE of BUPTMATLAB在信号与系统课程中的应用(4) 其他常见符号n = 变量赋值变量赋值n % 注释符注释符n 共轭转置符共轭转置符n : 冒号运算符冒号运算符n n:s:m 产生产生nm,步长为步长为s的序列，的序列，s可以为可以为 正或负或者小数，默认值为正或负或者小数，默认值为1。EE of BUPTMATLAB在信号与系统课程中的应用2.4

14、基本数学函数n三角函数：三角函数：nsin cos tan sinh cosh sech (双曲正割双曲正割)nasin acos atant=0:0.1:2*pi;x1=sin(t);plot(t,x1)y=asin(x1)figure;plot(t,y)EE of BUPTMATLAB在信号与系统课程中的应用指数函数nEXP: Exponential. 以以e为底的指数为底的指数nPOWER: . Array power. C = POWER(A,B) is called for the syntax A . B when A or B is an object. nSQRT：Square

15、 root. nPOW2: Base 2 power and scale floating point number. nNEXTPOW2: Next higher power of 2.a=1 2 3exp(a)power(a,2)power(a,a)pow2(a)nextpow2(15 )EE of BUPTMATLAB在信号与系统课程中的应用对数函数nlog :自然对数自然对数 %不要当最以不要当最以10为底的对数为底的对数nlog10: 以以10为底的对数为底的对数 nlog2:以以2为底的对数为底的对数log(exp(1)log10(1e6) %1e6 10的6次方EE of BUP

16、TMATLAB在信号与系统课程中的应用复数运算nabs:取模取模 nangle:辐角辐角 nreal:取实部取实部 nimag:取虚部取虚部nconj:取共轭取共轭 nisreal：True for real array. nunwrap: unwrap phase angle. ncplxpair: sort numbers into complex conjugate pairs.EE of BUPTMATLAB在信号与系统课程中的应用举例x=1+i 1+2i 1-i 1-2icplxpair(x)ans = Columns 1 through 3 1.0000 - 2.0000i 1.0

17、000 + 2.0000i 1.0000 - 1.0000i Column 4 1.0000 + 1.0000iEE of BUPTMATLAB在信号与系统课程中的应用举例nt=0:0.1:7nx=2*atan(tan(t)nplot(t,x)nhold onnphase=unwrap(x)nplot(t,phase,r:)01234567-5051015t originalunwrapedEE of BUPTMATLAB在信号与系统课程中的应用取整函数nround: Round towards nearest integer. 四舍五入四舍五入 nfix: Round towards zer

18、o. 保留整数部分保留整数部分 nfloor: Round towards minus infinity. nceil: Round towards plus infinity. nsign: Signum function. nrem: Remainder（余数）（余数） after division. nmod: Modulus （模数（模数) after division.EE of BUPTMATLAB在信号与系统课程中的应用2.5 基本绘图函数nMATLAB提供了强大的图形绘制功能。在大多数情提供了强大的图形绘制功能。在大多数情况下，用户只需要指定绘图的方式，提供绘图数据，况下，用户

19、只需要指定绘图的方式，提供绘图数据，利用利用MATLAB提供的丰富的二维，三维图形函数，提供的丰富的二维，三维图形函数，就可以绘制出所需的图形。就可以绘制出所需的图形。n(1) 绘制二维连续函数绘制二维连续函数 MATLAB中最常用的绘图函数是中最常用的绘图函数是plot，plot的命令格的命令格式有以下几种：式有以下几种：EE of BUPTMATLAB在信号与系统课程中的应用n(1) plot(y) 当当y为一向量时，以为一向量时，以y的的序号序号作为作为X轴，按向量轴，按向量y的值绘制曲线。的值绘制曲线。n(2) plot(x,y) x,y均为向量时，以均为向量时，以x向量作为向量作为X

20、轴，向量轴，向量y作为作为Y轴绘制曲线。轴绘制曲线。n注意：注意：x和和y中元素的个数必须相同中元素的个数必须相同！EE of BUPTMATLAB在信号与系统课程中的应用help plotPLOT Linear plot. PLOT(X,Y) plots vector Y versus vector X. If X or Y is a matrix, then the vector is plotted versus the rows or columns of the matrix, whichever line up. If X is a scalar and Y is a vector

21、, length(Y) disconnected points are plotted.PLOT(Y) plots the columns of Y versus their index. If Y is complex, PLOT(Y) is equivalent to PLOT(real(Y),imag(Y). In all other uses of PLOT, the imaginary part is ignored.EE of BUPTMATLAB在信号与系统课程中的应用Various line types, plot symbols and colors may be obtai

22、ned with PLOT(X,Y,S) where S is a character string made from one element from any or all the following 3 columns:EE of BUPTMATLAB在信号与系统课程中的应用设置 b blue . point - solid g green o circle : dotted r red x x-mark -. dashdot c cyan + plus - dashed m magenta * star (none) no line y yellow s square k black

23、d diamond v triangle (down) triangle (up) triangle (right) p pentagram 五边形五边形 h hexagram 六边形六边形默认的颜色是blue，多条曲线按照顺序循环使用。EE of BUPTMATLAB在信号与系统课程中的应用PLOT(X1,Y1,S1,X2,Y2,S2,X3,Y3,S3,.) combines the plots defined by the (X,Y,S) triples, where the Xs and Ys are vectors or matrices and the Ss are strings.

24、PLOT returns a column vector of handles to line series objects, one handle per plotted line.See also plottools, semilogx, semilogy, loglog, plotyy, plot3, grid, title, xlabel, ylabel, axis, axes, hold, legend, subplot, scatter.EE of BUPTMATLAB在信号与系统课程中的应用例子n=1:20;x1=0:8;x=x1*ones(1,20)+rand(9,20);h=

25、plot(n,x)legend(1,2,3,4,5,6,7,8,9);set(h,linewidth,2);set(h(1),linewidth,4);EE of BUPTMATLAB在信号与系统课程中的应用图形024681012141618200123456789 123456789EE of BUPTMATLAB在信号与系统课程中的应用(2)绘制二维离散序列n在在MATLAB用用stem命令实现离散序列的绘制。命令实现离散序列的绘制。 stem命令的格式有以下几种命令的格式有以下几种:(1)stem(y)以以x=1,2,3（缺省）（缺省）为各点数据的为各点数据的x坐标，以坐标，以y向量的各

26、个对应元素为向量的各个对应元素为y坐标，在（坐标，在（X,Y）坐标）坐标面画一个空心小圆圈。面画一个空心小圆圈。(2)stem(x,y,filled) 以以x向量的各个元素为向量的各个元素为x坐标，以坐标，以y向量的各个对向量的各个对应元素为应元素为y坐标，在（坐标，在（x,y）坐标面画一个实心）坐标面画一个实心小圆圈，并连接一条线段到小圆圈，并连接一条线段到X轴。轴。EE of BUPTMATLAB在信号与系统课程中的应用 举例% exm1_2.m, used to show the usage of stemn=0:20;x1=sin(pi/4*n);x2=exp(-0.5*n);h1=s

27、tem(n,x1);set(h1,linewidth,2);hold onh2=stem(n,x2,filled,r);set(h2,linewidth,2);legend(sin,exp);hold offEE of BUPTMATLAB在信号与系统课程中的应用图形02468101214161820-1-0.8-0.6-0.4-0.200.20.40.60.81 sinexpEE of BUPTMATLAB在信号与系统课程中的应用其他二维图形命令nstairs 阶梯图阶梯图nbar 条形图条形图nerrorbar 误差条形图误差条形图nhist 直方图直方图nfill(t,y,颜色标注符颜色

28、标注符) 在曲线和坐标轴之间的封闭区在曲线和坐标轴之间的封闭区填以指定的颜色填以指定的颜色npolar 极坐标绘图极坐标绘图nloglog 双对数双对数X-Y坐标绘图坐标绘图nsemilogx 半对数半对数X坐标绘图坐标绘图nsemilogy 半对数半对数Y坐标绘图坐标绘图nplotyy 用左、右两种用左、右两种y坐标绘图坐标绘图EE of BUPTMATLAB在信号与系统课程中的应用stairs: 小台阶，相当于零阶保持小台阶，相当于零阶保持x=0:0.25:10;stairs(x,sin(x);012345678910-1-0.8-0.6-0.4-0.200.20.40.60.81x=1

29、1.1 1.3 1.4 1.6y=x+1stairs(x,y)11.11.21.31.41.51.61.71.822.12.22.32.42.52.62.72.8EE of BUPTMATLAB在信号与系统课程中的应用hist: 可以用来表示分布可以用来表示分布请调出一个班级的数据，找出最高分，最低分，计算请调出一个班级的数据，找出最高分，最低分，计算平均分和画出分布图。平均分和画出分布图。作业：找一个班级的数据，内容同上作业：找一个班级的数据，内容同上1220283644526068768405101520253035hist(score)EE of BUPTMATLAB在信号与系统课程中的

30、应用semilogxt=0:0.2:1000;x=power(1.1,t);subplot(2,1,1);plot(x);subplot(2,1,2);semilogy(x);010002000300040005000600000.511.522.5x 104101000200030004000500060001001050SEMILOGY(.) is the same as PLOT(.), except a logarithmic (base 10) scale is used for the Y-axis.EE of BUPTMATLAB在信号与系统课程中的应用plotyyx1=0:0.

31、1:2*pi;x2=x1;y2=x2.2;plotyy(x1,y1,x2,y2)01234567-1-0.500.5101234567010203040EE of BUPTMATLAB在信号与系统课程中的应用二维图形函数库中其他的主要函数n坐标控制坐标控制naxis 控制坐标轴比例和外观控制坐标轴比例和外观nhold 保持当前图形保持当前图形nsubplot 在平铺位置建立图形轴系在平铺位置建立图形轴系n图形注释图形注释ntitle 标出图形名称标出图形名称nxlabel X轴标注轴标注nylabel Y轴标注轴标注ntext 在图上标文字在图上标文字ngtext 用鼠标定位文字用鼠标定位文字

32、nlegend 标注图例标注图例ngrid 图上加坐标网格图上加坐标网格EE of BUPTMATLAB在信号与系统课程中的应用2.6 字符与字符串字符与字符串n数个字符可以构成字符串数个字符可以构成字符串(Strings)，一个字符串是被视，一个字符串是被视为一个为一个行向量行向量(Row vector)进行存储，而这个字符串中进行存储，而这个字符串中的每一个字符（含空格符），是以其的每一个字符（含空格符），是以其ASCII(The American Standard Code for Information Interchange )码的形式存放于此向量中的每一个元素码的形式存放于此向量中

33、的每一个元素(Elment)，只是，只是其外显形式仍然是可读的字符。其外显形式仍然是可读的字符。n例如例如string=Beijing University of Posts and Telecomunnicationsn多个字符串变量并排，可得到一个新的字符串变量多个字符串变量并排，可得到一个新的字符串变量nstr1=you ;nstr2=and menstr3=str1 str2EE of BUPTMATLAB在信号与系统课程中的应用单引号和类型检测单引号和类型检测nMATLAB用用单引号单引号（Single Quote）来界定字）来界定字符串变量，而符串变量，而C语言则是用语言则是用“双

34、引号双引号”（Double Quote）来界定字符串变量，这是两者的不同之处。）来界定字符串变量，这是两者的不同之处。n要输入含有单引号的字符串，可重复单引号的使用，要输入含有单引号的字符串，可重复单引号的使用，例如例如nstring=Ive got a booknlength(string)nclass(string)nischarnisa(a,char)EE of BUPTMATLAB在信号与系统课程中的应用变量转换变量转换n若要查看字符串变量的存储内容（即若要查看字符串变量的存储内容（即ASCII），可），可使用使用double命令命令ndouble(北邮北邮)nstring_ascii

35、=double(string)n若要将字符串变量储存的若要将字符串变量储存的ASCII码转换成字符串形码转换成字符串形式，可用式，可用char命令命令nchar(string_ascii);nchar(97)nMATLAB中每个字符占用两个字节，汉字也占用两中每个字符占用两个字节，汉字也占用两个字节（中文的个字节（中文的ASCII 码大于码大于128）EE of BUPTMATLAB在信号与系统课程中的应用字符串命令的执行字符串命令的执行n执行某一特定字符串，可以使用执行某一特定字符串，可以使用eval命令，其效果就如命令，其效果就如同直接在同直接在MATLAB命令窗口内输入此一特定字符串命令

36、窗口内输入此一特定字符串nstr=x=1 2 3; y=x.2neval(str)neval命令特别适用于在命令特别适用于在for-loop内内自动产生变量自动产生变量，例如，例如nfor ii=3:6n str=x int2str(ii) =rand( int2str(ii) )n eval(str);nendnfeval()，它可以用字符串代表函数名称，并进行对此函，它可以用字符串代表函数名称，并进行对此函数的调用。数的调用。n feval(sin,1:4)EE of BUPTMATLAB在信号与系统课程中的应用使用一个变量来存储多个字符串使用一个变量来存储多个字符串n两种方法：两种方法：

37、n使用二维字符数组（使用二维字符数组（Two dimensional character arrays），但必须先确认每个字符串（即每一行）的），但必须先确认每个字符串（即每一行）的长度相同长度相同，否则就必须在短字符串结尾补上空格符。，否则就必须在短字符串结尾补上空格符。ndepartments=ee ;cs ;econn使用使用char命令命令ndepartments=char(ee,cs,econ)n在从二维字符数组抽取出字符串时，切记要使用在从二维字符数组抽取出字符串时，切记要使用deblank命令来删除尾部的空格符。命令来删除尾部的空格符。ndeblank(departments(1

38、,:)ncellstr转成异质数组转成异质数组EE of BUPTMATLAB在信号与系统课程中的应用字符串的比较、寻找、代换、分解与结合字符串的比较、寻找、代换、分解与结合nstrcmpnCompare strings.nstrcmp(char(97),a)nstrncmpnCompare first n characters of strings.nstrrepnSTRREP Replace string with another.n S = STRREP(S1,S2,S3) replaces all occurrences of the string S2 in string S1 wi

39、th the string S3. The new string is returned.nstrrep(my mother is a kind lady,my,your)nstrtok：Find token in string.nstrcat：Concatenate strings.EE of BUPTMATLAB在信号与系统课程中的应用字符串、数值与数组的转换字符串、数值与数组的转换nint2str()nnum2str()ndec2hex(): 10进制到进制到16进制进制nhex2num()nhex2dec()nbin2dec()ndec2bin()：Convert decimal in

40、teger to a binary string.nbase2dec()：BASE2DEC(S,B) converts the string number S of base B into its decimal (base 10) equivalent. nmat2str()：将矩阵转换为字符串：将矩阵转换为字符串EE of BUPTMATLAB在信号与系统课程中的应用举例举例nnum_1=input(Please input the first number: );nnum_2=input(Please input the second number: );nall=num_1+num_2

41、;nstr1=num2str(num_1);nstr2=num2str(num_2);nstr_all=num2str(all);nstring=str1 + str2 = str_all;ndisp(string);EE of BUPTMATLAB在信号与系统课程中的应用2.7 MATLAB 2.7 MATLAB 程序设计程序设计n MATLAB的工作方式有二种的工作方式有二种n交互式的交互式的指令操作指令操作方式。即用户在命令窗口方式。即用户在命令窗口中输入命令并按下回车键后，系统执行该指中输入命令并按下回车键后，系统执行该指令并立即给出运算结果。令并立即给出运算结果。nM文件文件的编程方

42、式。的编程方式。M文件是由文件是由 MATLAB语句构成的文件语句构成的文件,且文件名必须以且文件名必须以.m为扩展为扩展名，如名，如example.m。用户可以用任何文件编。用户可以用任何文件编辑器来对辑器来对M文件进行编辑。文件进行编辑。 EE of BUPTMATLAB在信号与系统课程中的应用2.7.1 程序流程控制程序流程控制n顺序结构：顺序结构：MATLAB 从上到下依次执行各语句，从上到下依次执行各语句，该结构最简单。该结构最简单。n循环结构循环结构n(1)for-end 循环循环n(2) while-end 循环循环n分支结构分支结构nif 语句语句n多条件语句多条件语句EE o

43、f BUPTMATLAB在信号与系统课程中的应用(1)for-end 循环用于循环次数事先确定的，格式为用于循环次数事先确定的，格式为for ii=n:s:m 语句体语句体end s 为步长，可以为正数，负数或小数。为步长，可以为正数，负数或小数。clcfor ii=1:5 kk=2*ii; %disp(kk); disp(strcat(kk=,int2str(kk); pause;endEE of BUPTMATLAB在信号与系统课程中的应用举例n编程生成一个编程生成一个6阶矩阵，使其主对角线上元素皆为阶矩阵，使其主对角线上元素皆为1，与，与主对角线相邻元素皆为主对角线相邻元素皆为2，其余元

44、素皆为，其余元素皆为0。000000000000000022222222220000111111a a=eye(6);for ii=1:5 a(ii,ii+1)=2; a(ii+1,ii)=2;enddiag(1 1 1 1 1 1)+diag(2 2 2 2 2,1)+diag(2 2 2 2 2,-1)EE of BUPTMATLAB在信号与系统课程中的应用用于循环次数不能事先确定的，用于循环次数不能事先确定的，格式为格式为 while 表达式表达式 语句体语句体 end 只有表达式为真，就执行语句体，只有表达式为真，就执行语句体，表达式为假，终止该循环。表达式为假，终止该循环。kk=-3;while kk0 y=1; disp(if x0, then y=1);else y=0; disp(else, y=0);end分支结构EE of BUPTMATLAB在信号与系统课程中的应用多条件语句格式为格式为if 表达式表达式1 语句体语句体1elseif 表达式表达式2 语句体语句体2.el