实验3离散信号与系统的离散频域分析Word文件下载.docx

实验3离散信号与系统的离散频域分析Word文件下载.docx

《实验3离散信号与系统的离散频域分析Word文件下载.docx》由会员分享，可在线阅读，更多相关《实验3离散信号与系统的离散频域分析Word文件下载.docx（3页珍藏版）》请在冰点文库上搜索。

16点序列，作DFT，画出幅度谱和相位谱，并与原序列的DFT进行比较。

3、将信号以8为周期扩展，得到长为16的两个周期，作DFT，画出幅度谱和相位谱，并与原序列的DFT进行比较。

上机操作

1.

>

n=0:

7;

x=（0.9*exp（j*pi/3））.^n;

x1=[xzeros（1,8）];

x2=[x1zeros（1,16）];

Xk=dft（x,8）;

Xk1=dft（x1,16）;

Xk2=dft（x2,32）;

w=0:

pi/200:

2*pi;

X=dtft（x,n,w）;

magX=abs（X）;

angX=angle（X）;

magXk=abs（Xk）;

angXk=angle（Xk）;

magXk1=abs（Xk1）;

angXk1=angle（Xk1）;

magXk2=abs（Xk2）;

angXk2=angle（Xk2）;

subplot（4,2,1）;

plot（w/pi,magX）;

xlabel（'

w/pi'

）;

ylabel（'

幅度|X|'

gridon;

subplot（4,2,2）;

plot（w/pi,angX）;

相位（rad/π）'

subplot（4,2,3）;

stem（0:

7,magXk）;

k'

幅度|X（k）|'

axis（[0,8,0,6]）;

subplot（4,2,4）;

7,angXk）;

axis（[0,8,-2,2]）;

subplot（4,2,5）;

15,magXk1）;

幅度|X1（k）|'

axis（[0,16,0,6]）;

subplot（4,2,6）;

15,angXk1）;

axis（[0,16,-2,2]）;

subplot（4,2,7）;

31,magXk2）;

幅度|X2（k）|'

axis（[0,32,0,6]）;

subplot（4,2,8）;

31,angXk2）;

axis（[0,32,-2,2]）;

fori=1:

8 %将原序列每两个相邻样本之间插入一个零值，扩展为长为16的序列y1（2*i-1:

2*i）=[x（i）0];

end

Yk1=dft（y1,16）;

Warning:

FunctioncalldftinvokesinexactmatchD:

\MATLAB7\work\DFT.M.

magYk1=abs（Yk1）;

angYk1=angle（Yk1）;

y2=[xx];

%将信号以8为周期扩展，得到长为16的两个周期

Yk2=dft（y2,16）;

magYk2=abs（Yk2）;

angYk2=angle（Yk2）;

subplot（3,2,1）;

subplot（3,2,2）;

subplot（3,2,3）;

15,magYk1）;

幅度|Y1（k）|'

subplot（3,2,4）;

15,angYk1）;

subplot（3,2,5）;

15,magYk2）;

幅度|Y2（k）|'

subplot（3,2,6）;

15,angYk2）;