数字信号处理实验指导书-程序.docVIP

数字信号处理实验指导书 课程组 扬州大学信息工程学院 2013年10月 2 信号采样与重建.4 实验内容 实验内容一 MATLAB参考程序和仿真内容 %*******************************************************************% %f—余弦信号的频率 % M—基2 FFT幂次数 N=2^M为采样点数，这样取值是为了便于作基2的FFT分析 %2. 采样频率Fs %*******************************************************************% function samples(f,Fs,M) N=2^M; % fft点数=取样总点数 Ts=1/Fs; % 取样时间间隔 T=N*Ts; % 取样总时间=取样总点数*取样时间间隔 n=0:N-1; t=n*Ts; Xn=cos(2*f*pi*t); subplot(2,1,1); stem(t,Xn); axis([0 T 1.1*min(Xn) 1.1*max(Xn)]); xlabel(t --); ylabel(Xn); Xk=abs(fft(Xn,N)); subplot(2,1,2); stem(n,Xk); axis([0 N 1.1*min(Xk) 1.1*max(Xk)]); xlabel(frequency --); ylabel(!Xk!); %*******************************************************************% 实验内容二 MATLAB参考程序和仿真内容 % 信号取样与重构 % k = 1 临界取样 0k1 过取样 k1 欠取样 display(Please input the value of k); k = input(k = ); wm = 1; %信号带宽 Ts = k*pi/wm; %取样间隔 ws = 2*pi/Ts; %取样角频率 wc = ws/2; %为简便，滤波器截止频率wc=ws/2 n = -10:10; %定义取样点的数量 m = fix(length(n)/2); %单边取样周期数,fix函数的功能为取整 nTs = n*Ts; %计算每个取样点 dt = 0.05; t = -m*Ts:dt:m*Ts; %Sa(t)信号取样 f = sinc(nTs/pi); %计算Sa(t)在取样点处的函数值，得到取样信号 %Sa(t)信号重构 fa = f*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs*ones(1,length(t)))); %Sa(t)包络线 f2=sinc(t/pi); subplot(311); stem(nTs,f); %绘制取样信号 hold on plot(t,f2,r:); %绘制包络线 xlabel(t); title(对sa(t)信号进行取样); axis([-m*Ts m*Ts -0.5 1.2]); hold off subplot(312); h1=plot(t,fa); %绘制重构信号 hold on for i=-m*Ts:Ts:m*Ts ft = sinc(i/pi)*sinc(wc/pi*(t-i)); h2=plot(t,ft,m:); %绘制重构信号的各个分量 hold on yy = sinc(i/pi); yy1 = yy:yy; plot(i,yy1,o); end xlabel(t); title(重构信号); axis([-m*Ts m*Ts -0.5 1.2]); legend([h1,h2],重构信号,分信号); hold off subplot(313); %计算原信号和重构信号之间的误差 error=abs(fa-f2); plot(t,error); axis([-m*Ts m*Ts min(erro