实验一低通采样定理和内插与抽取实现.docVIP

实验一 低通采样定理和内插与抽取实现 一、实验目的 用Matlab 编程实现自然采样与平顶采样过程，根据实验结果给出二者的结论；掌握利用MATLAB实现连续信号采样、频谱分析和采样信号恢复的方法。 二、实验原理 1．抽样定理 若是带限信号，带宽为, 经采样后的频谱就是将的频谱 在频率轴上以采样频率为间隔进行周期延拓。因此，当时，不会发生频率混叠；而当 时将发生频率混叠。 2．信号重建 经采样后得到信号经理想低通则可得到重建信号，即：=* 其中：== 所以： =*=* = 上式表明，连续信号可以展开成抽样函数的无穷级数。 利用MATLAB中的来表示，有 ，所以可以得到在MATLAB中信号由重建的表达式如下： = 我们选取信号=作为被采样信号，当采样频率=2时，称为临界采样。我们取理想低通的截止频率=。下面程序实现对信号=的采样及由该采样信号恢复重建： ，试以以下采样频率对信号采样： (a) ; (b) ; (c)， 求x(t)信号原信号和采样信号频谱，及用采样信号重建原信号x’(t)时序图。假定重构滤波器为fs/2的理想低通滤波器，具有带通增益为Ts=1/fs. 四、 实验步骤 五、 实验结果分析结论 显然，以10000Hz采样，采样频谱发生了混叠。这是因为该采样频率不满足奈奎斯特采样定理：fs=2fh。 六、实验程序 %通信系统建模与仿真 实验一 带通采样定理和内插与抽取实现 %主要参数：ts采样间隔 fs采样频率 clear all clc fs=20000; %采样频率 ts=1/fs; %采样间隔 T=0.05; %截取时间 df=1/T; %频域间隔 t=0:ts:T-ts; %时间序列 x=0; for m=0:99 %原始信号 x=x+(m+1)*cos(2*pi*(100*m+50)*t); end figure(1) subplot(211) plot(t,x) %画原始信号波形图 title(原始信号波形);xlabel(时间t);ylabel(幅度); [X,f]=fftseq(x,ts); %傅里叶变换 subplot(212) plot(f,fftshift(abs(X)));title(原始频谱),xlabel(f),ylabel(X) fs1=10000; ts1=1/fs1; T1=0.05; df1=1/T1; t1=0:ts1:T1-ts1; x1=0; for m=0:99 x1=x1+(m+1)*cos(2*pi*(100*m+50)*t1); end figure(2) subplot(211) plot(t1,x1) title(10000Hz抽样波形);xlabel(时间t);ylabel(幅度); [X1,f1]=fftseq(x1,ts1); subplot(212) plot(f1,fftshift(abs(X1)));title(以10000Hz采样频谱),xlabel(f),ylabel(X) fs3=30000; %抽样频率 ts3=1/fs3; T3=0.05; df3=1/T3; t3=0:ts3:T3-ts3; x3=0; for m=0:99 x3=x3+(m+1)*cos(2*pi*(100*m+50)*t3); end figure(3) subplot(211) plot(t3,x3) title(30000Hz抽样波形);xlabel(时间t);ylabel(幅度); [X3,f3]=fftseq(x3,ts3); subplot(212) plot(f3,fftshift(abs(X3)));title(以30000Hz采样频谱),xlabel(f),ylabel(X) fp=fs/2-500;fs1=fs/2; [n,f0,m0,w]=firpmord([fp fs1],[1 0],[0.0001 0.001],fs); b=firpm(n,f0,m0,w); figure(4) freqz(b,1,1024,fs) y=conv(x,b); figure(5) subplot(211) plot(t,y(n:length(y)-1)); title(重建原始信号时域);xlabel(时间t);ylabel(幅度); [X,f]=fftseq(y(n:length(y)-1),ts); subplot(212) plot(f,fftshift(abs(X))) title(