序列的产生和运算.docVIP

序列的产生和运算 序列的产生 正弦序列离散正弦序列的MATLAB表示与连续信号类似，只不过是用stem函数而不是用plot函数来画出序列的波形。下面是正弦序列的MATLAB运行结果如图1.1所示。 % Generation of the sinusoidal sequence =0:39; x=sin(pi/6*n); stem(n,x);stem(n,x)stem(x,filled); 图1.1 正弦序列波形“stem”用于绘制针状图将需要绘制的数据存放在一个数组中stem(x); %绘制数组x中的数据，从1开始 stem(n,x); %绘制数组x中的数据，与n对应 stem(n,x,filled); %填充 实指数序列 指数序列的一般形式为，可用MATLAB中的数组幂运算（即点幂运算）来实现。下面是正弦序列的运行结果如图1.所示% Generation of the exponential sequence n=0:20; x=2*(1.2.^n); stem(n,x);图1.2 指数序列的波形 更通用地，可编写绘制离散时间实指数序列波形的函数% Generation of the exponential sequence function dszsu(k,a,n1,n2) %k：指数序列的幅度 %a：指数序列的底数 %1：绘制序列的起始序号 %2：绘制序列的终止序号 =n1:n2; x=k*(a.^n); stem(n,x,filled) hold on plot([n1,n2],[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.所示。如图可知，对于离散时间实指数序列，当的绝对值大于1时，序列为随时间发散的序列，当的绝对值小于1时，序列为随时间收敛的序列。同时可见，当的值小于零时，其波形在增长或衰减的同时，还交替地改变序列值的符号。图1.序列指数序列 离散时间复指数序列的一般形式为function[]=dComplexExp(a,b,n1,n2) %k：复指数序列的比例因子 %a,b：复指数(a+bi) %n1：绘制波形的复指数序列的起始时间序号 %n2：绘制波形的复指数序列的终止时间序号 n=n1:n2; c=a+b*i; f=k*exp(c*n); Xre=real(f) Xim=imag(f) Xmag=abs(f) Xphs=angle(f) subplot(2,2,1), stem(n,Xre,filled),title(实部); subplot(2,2,3), stem(n,Xim,filled),title(虚部); subplot(2,2,2), stem(n,Xmag,filled),title(模); subplot(2,2,4), stem(n,Xphs,filled),title(相角); ） %离散时间虚指数实现程序 figure(1); dComplexExp(1,0,pi/4,0,20); figure(2); dComplexExp(1,0,2,0,20); figure(3); dComplexExp(1,0.1,pi/4,0,20); figure(4); dComplexExp(1,-0.1,2,0,20); 程序运行结果如图1.4所示。由图可见，只有为有理数时，信号的实部和虚部和相角都为周期序列，否则为非周期序列。（a）波形(b) 波形 （a）波形(b) 波形 图1.4 复指数序列波形 单位冲激序列 可以通过借助MATLAB中的零矩阵函数zeros表示。全零矩阵zeros(1,N)产生一个由N个零组成的列向量，对于有限区间的可以通过以下MATLAB程序表示% 单位抽样序列实现程序 k=-30:30; delta=[zeros(1,30),1,zeros(1,30)]; stem(k,delta) 程序运行结果如图1.所示。单位阶跃序列 可以通过借助MATLAB中的单位矩阵函数ones表示。单位矩阵ones(1,N)产生一个由N个1组成的列向量，对于有限区间的可以通过以下MATLAB程序表示 % 单位阶跃