第一部分 MAAB基础.docVIP

第一部分 MATLAB基础 一.运算符 矩阵运算：+ - * / \ 一般运算：.+ ._ .* ./ 二.矩阵 [1 2 3 4 5 ]; [1 2 3 4 5;1 2 3 4 5; 1 2 3 4 5]; 第一个数:间隔:最后一个数; 1:1:5-[1 2 3 4 5] Zeros(n,n) Ones() Eyes() 三.画图 连续：Plot(x,y); 离散：Stem(x,y); stairs Xlabel Ylabel Title() Axis() Figure() Subplot 四.函数 Sin,cos,tan… Exp,real,imag,abs,length 连续信号：Sign(符号信号),sinc(抽样信号),rectpuls(矩形脉冲信号),tripuls(三角波脉冲信号),square(周期性矩形脉冲),sawtooth(周期性三角波信号) 五.离散时间信号 1．正弦序列 离散正弦序列的MATLAB表示与连续信号类似，只不过是用stem函数而不是用plot函数来画出序列的波形。下面就是正弦序列的MATLAB源程序。程序运行结果如图1.19所示。 %正弦序列实现程序 k=0:39; fk=sin(pi/6*k); stem(k,fk) 图1.19，可用MATLAB中的数组幂运算（即点幂运算）c*来实现。下面为用MATLAB编写绘制离散时间实指数序列波形的函数。 function dszsu(c,a,k1,k2) %c：指数序列的幅度 %a：指数序列的底数 %k1：绘制序列的起始序号 %k2：绘制序列的终止序号 k=k1:k2; x=c*(a.^k); stem(k,x,filled) hold on plot([k1,k2],[0,0]) hold off 利用上述函数，实现实指数波形MATLAB程序如下（其中值分别为）。 %离散时间实指数序列实现程序 subplot 221; dszsu(1,5/4,0,20); xlabel(k); title(f1[k]); subplot 222 dszsu(1,3/4,0,20); xlabel(k); title(f2[k]); subplot 223; dszsu(1,-5/4,0,20); xlabel(k); title(f3[k]); subplot 224; dszsu(1,-3/4,0,20); xlabel(k); title(f4[k]); 程序运行结果如图1.20所示。如图可知，对于离散时间实指数序列，当的绝对值大于1时，序列为随时间发散的序列，当的绝对值小于1时，序列为随时间收敛的序列。同时可见，当的值小于零时，其波形在增长或衰减的同时，还交替地改变序列值的符号。 图1.20function[]=dxzsu(n1,n2,w) %n1:绘制波形的虚指数序列的起始时间序号 %n2：绘制波形的虚指数序列的终止时间序号 %w：虚指数序列的角频率 k=n1:n2; f=exp(i*w*k); Xr=real(f) Xi=imag(f) Xa=abs(f) Xn=angle(f) subplot(2,2,1), stem(k,Xr,filled),title(实部); subplot(2,2,3), stem(k,Xi,filled),title(虚部); subplot(2,2,2), stem(k,Xa,filled),title(模); subplot(2,2,4), stem(k,Xn,filled),title(相角); 利用上述函数，实现虚指数波形MATLAB程序如下（其中虚指数分别为） %离散时间虚指数实现程序 figure(1); dxzsu(0,20,pi/4); figure(2); dxzsu(0,20,2); 程序运行结果如图1.21（a）、（b）所示。由图可见，只有当虚指数序列的角频率满足为有理数时，信号的实部和虚部和相角都为周期序列，否则为非周期序列。 （a）波形 (b) 波形 图1.21 虚指数序列波形 对于复指数序列，其一般形式为 可以通过调用下面绘制复指数序列时域波形的MATLAB函数。 function dfzsu(n1,n2,r,w) %n1:绘制波形的虚指数序列的起始时间序号 %n2：绘制波形的虚指数序列的终止时间序号 %w：虚指数序列的角频率 %r: 指数序列的底数 k=n1:n2; f=(r*exp(i*w)).^k; Xr=real(f); Xi=imag(f); Xa=abs(f); Xn=angle(f); subplot(2,2,1),