% 绘制x(n)的幅度谱
k=-250:250;
W=pi/125*k;
X=x*(exp(-j*pi/125)).^(n'*k);
magX=abs(X);
subplot(1,2,2),plot(W/pi,magX)
title('理想采样序列的幅度谱')
当把其中的T该为响应的值就可以得到其他的几种频率的图形。
作业二
程序和图形如下：
n=0:5;
x=[1,1,1,1,0,0]
h=[1,1,1,1,0,0]
a=0:10;
y=conv(x,h);
stem(a,y);
title('序列y(n)的值')
理论长度是n1+n2-1应该是4+4-1=7
图形中非0值也是7个，所以和理论值是一样的。
作业三
方法一：n=0:50;
A=444.128;
a=50*sqrt(2.0)*pi;
T=1/1000;
w0=50*sqrt(2.0)*pi;
x=A*exp(-a*n*T).*sin(w0*n*T)
h=[1,2.5,2.5,1];
y=conv(x,h);
a=0:53;
subplot(1,2,1),stem(a,y)
title('序列 y(n)')
k=-250:250;
W=pi/125*k;
Y=y*(exp(-j*pi/125)).^(a'*k);
magy=abs(Y);
subplot(1,2,2),plot(W/pi,magy)
title('理想采样序列的幅度谱')
方法二：
> n=0:50;
A=444.128;
a=50*sqrt(2.0)*pi;
T=1/1000; % T分别取为1/1000、1/300和1/200
w0=50*sqrt(2.0)*pi;
x=A*exp(-a*n*T).*sin(w0*n*T);
k=-250:250;
W=pi/125*k;
X=x*(exp(-j*pi/125)).^(n'*k);
h=[1,2.5,2.5,1,];
m=0:3;
H=h*(exp(-j*pi/125)).^(m'*k);
y=H.*X;
Y=abs(y);
plot(W/pi,Y)
图形如下：
俩种方法的对比：由图形可以看出来，俩种方法的结果是相同的，我们可以看出俩种虽然一个是对时域进行分析，而另一种则是从频域开始分析，分析的方法不同，但是结果却是相同的，所以时域和频域是紧密相连的，我们在今后的学习里绝对不可以把俩者分开。