数字信号大作业(2)

title('y(n)的时域特性曲线')

subplot(2,1,2);stem(n2,hbn,'.') %绘制hb(n)的时域特性曲线 xlabel('n');ylabel('hb(n)') title('hb(n)的时域特性曲线')

figure(2) %用于分析y(n)的频域特性 w=linspace(-pi,pi,1000); %限制w的范围 Y=yn*exp(-1i*n1'*w); %对y(n)做DTFT变换 subplot(2,1,1);

plot(w/pi,abs(Y)); %绘制y(n)的幅度曲线 xlabel('\\omega/π');ylabel('幅度') title('DTFT[y(n)]的幅度');

Hb=hbn*exp(-1i*n2'*w); %对hb(n)做DTFT变换%用于分析hb(n)的频域特性 subplot(2,1,2);

plot(w/pi,abs(Hb)); %绘制hb(n)的幅度曲线 xlabel('\\omega/π');ylabel('幅度') title('DTFT[hb(n)]的幅度');

for k=1:2 N=10;

xcn=ones(1,N); han=ones(1,N); figure(2+k)

zn=conv(xcn,han); %观察系统ha(n)对信号xc(n)的响应特性 n3=0:length(zn)-1;

w=linspace(-pi,pi,1000); %限制w的范围 Z=zn*exp(-1i*n3'*w); %对z(n)做DTFT变换

plot(w/pi,abs(Z)); %绘制z(n)的幅度曲线 xlabel('\\omega/π');ylabel('幅度') title('DTFT[zn]的幅度'); End

xb(n)和系统hb(n)的时域和频域特性：

y(n)时域及频域特性：

N=10时。Y(n)的时域频域特性曲线

N=5时。Y(n)的时域频域特性曲线

③ 时域离散信号、 系统和系统响应分析。 实验代码：

clear all; A=1;a=0.4; m=2.0734; N=50;T=1; n1=0:N-1;

w=linspace(-pi*2,pi*2,1000); %设置w的范围

xn=A*exp(-a*n1*T).*sin(m*n1*T); hbn=[1,2.5,2.5,1];

yn=conv(xn,hbn); %计算信号xb(n)通过系统hb(n)的响应y(n) n2=0:length(hbn)-1; n3=0:length(yn)-1;

X=xn*exp(-1i*n1'*w); %对x(n)做DTFT变换 H=hbn*exp(-1i*n2'*w); %对hb(n)做DTFT变换 Y1=yn*exp(-1i*n3'*w); %对y(n)做DTFT变换 Y2=X.*H; figure(1);

subplot(2,1,1);stem(n1,xn,'.') %绘制y(n)的时域特性曲线 xlabel('n');ylabel('x(n)') title('x(n)的时域特性曲线') subplot(2,1,2);

plot(w/pi,abs(X)); %绘制y(n)的幅度曲线 xlabel('\\omega/π');ylabel('幅度') title('DTFT[x(n)]的幅度'); figure(2);

subplot(2,1,1);stem(n2,hbn,'.') %绘制hb(n)的时域特性曲线 xlabel('n');ylabel('hb(n)') title('hb(n)的时域特性曲线') subplot(2,1,2);

plot(w/pi,abs(H)); %绘制y(n)的幅度曲线 xlabel('\\omega/π');ylabel('幅度') title('DTFT[hb(n)]的幅度');

figure;

subplot(2,1,1); %绘制z(n)的幅度和相位曲线 plot(w/pi,abs(Y1));

xlabel('\\omega/π');ylabel('幅度') title('DTFT[yn]的幅度');

subplot(2,1,2); %绘制z(n)的幅度和相位曲线