信号处理实验三报告.pdf

实验三抽样

3.3.1抽样引起的混叠

【实验内容】

f=500HZ，fs分别取100HZ，1kHZ，10kHZ，绘出x[n]及其DTFT。

0

【程序】

%文件exp331.mxlabel(n);ylabel(x1[n]);

%问题3.3.1实验内容subplot(2,3,4);%绘制x1频域解图

closeall;clearall;plot(W1,X1);

f0=500;gridon;

B=0.25*pi;%初相title(fs=100HZ采样信号的DTFT);

n=0:20;xlabel(W1);ylabel(X1);

fs1=100;%采样频率100HZ%绘制x2时域采样结果

t1=n/fs1;subplot(2,3,2);stem(n,x2);

x1=sin(2*pi*f0*t1+B);title(fs=1kHZ时域采样结果);

[X1,W1]=dtft(x1,1000);%对采样信号做xlabel(n);ylabel(x2[n]);

DTFTsubplot(2,3,5);%绘制x2频域解图

fs2=1000;%采样频率1kHZplot(W2,X2);gridon;

t2=n/fs2;title(fs=1kHZ采样信号的DTFT);

x2=sin(2*pi*f0*t2+B);xlabel(W2);ylabel(X2);

[X2,W2]=dtft(x2,1000);%对采样信号做%绘制x3时域采样结果

DTFTsubplot(2,3,3);stem(n,x3);

fs3=10000;%采样频率10kHZtitle(fs=10kHZ时域采样结果);

t3=n/fs3;xlabel(n);ylabel(x3[n]);

x3=sin(2*pi*f0*t3+B);%绘制x3频域解图

[X3,W3]=dtft(x3,1000);%对采样信号做subplot(2,3,6);plot(W3,X3);gridon;

DTFTtitle(fs=10kHZ采样信号的DTFT);

hf=figure;%绘制x1时域采样结果xlabel(W3);ylabel(X3);

subplot(2,3,1);stem(n,x1);saveas(hf,exp331.fig);

title(fs=100HZ时域采样结果);

【结果及分析】

如图1.1所示，前三幅图分别为欠采样、临界采样和完全采样的结果，后三幅图为他们分别

对应的DTFT。

1/10

fs=100HZ时域采样结果fs=1kHZ时域采样结果fs=10kHZ时域采样结果

0.811

0.60.5