Matlab实验6.pptVIP

* 实验六：连续时间信号的数字处理 一、实验目的： 1、加深对信号时域采样和频域采样过程和效果的理解； 2、复习模拟低通滤波器的设计。 二、实验内容： 1、时域的采样过程和混叠效果； 2、频域中的采样效果； 3、模拟低通滤波器的设计； 4、语音信号的量化 1、 时域的采样过程和混叠效果 clf; T = 0.1;f = 13; n = (0:T:1); xs = cos(2*pi*f*n); t = linspace(-0.5,1.5,500); ya = sinc((1/T)*t(:,ones(size(n))) - (1/T)*n(:,ones(size(t))))*xs; plot(n,xs,o,t,ya);grid; xlabel(Time, msec);ylabel(Amplitude); title(Reconstructed continuous-time signal y_{a}(t)); axis([0 1 -1.2 1.2]) 三、实验例程 2、.频域中的混叠效果 clf; t = 0:0.005:10; xa = 2*t.*exp(-t); subplot(2,2,1) plot(t,xa);grid xlabel(Time, msec);ylabel(Amplitude); title(Continuous-time signal x_{a}(t)); wa = 0:10/511:10; ha = freqs(2,[1 2 1],wa); plot(wa/(2*pi),abs(ha));grid; xlabel(Frequency, kHz);ylabel(Amplitude); title(|X_{a}(j\Omega)|); axis([0 5/pi 0 2]); subplot(2,2,3) T = 1; n = 0:T:10; xs = 2*n.*exp(-n); k = 0:length(n)-1; stem(k,xs);grid; xlabel(Time index n);ylabel(Amplitude); title(Discrete-time signal x[n]); subplot(2,2,4) wd = 0:pi/255:pi; hd = freqz(xs,1,wd); plot(wd/(T*pi), T*abs(hd));grid; xlabel(Frequency, kHz);ylabel(Amplitude); title(|X(e^{j\omega})|); axis([0 1/T 0 2]) 3、 模拟低通滤波器的设计（巴特沃兹） clf; Fp = 3500;Fs = 4500; Wp = 2*pi*Fp; Ws = 2*pi*Fs; [N, Wn] = buttord(Wp, Ws, 0.5, 30,s); [b,a] = butter(N, Wn, s); wa = 0:(3*Ws)/511:3*Ws; h = freqs(b,a,wa); plot(wa/(2*pi), 20*log10(abs(h)));grid xlabel(Frequency, Hz);ylabel(Gain, dB); title(Gain response); axis([0 3*Fs -60 5]); *