數(shù)字信號(hào)處理實(shí)驗(yàn)六時(shí)域采樣與信號(hào)的重建

1、實(shí)驗(yàn)?zāi)康模?了解用matlab語(yǔ)言進(jìn)行時(shí)域抽樣與信號(hào)重建的方法2. 進(jìn)一步加深對(duì)時(shí)域信號(hào)抽樣與恢復(fù)的基本原理的理解3. 掌握采樣頻率的確定方法和內(nèi)插公式的編程方法。二.實(shí)驗(yàn)內(nèi)容1認(rèn)真閱讀并輸入實(shí)驗(yàn)原理與方法中介紹的例子，觀察輸出波形曲線，理解每一 條語(yǔ)句的含義。2 已知一個(gè)連續(xù)時(shí)間信號(hào)f(t)=sinc(t)。取最高有限帶寬頻率fm=1hz。(1) 分別顯示原連續(xù)時(shí)間信號(hào)波形和 fm=fm、fm=2fm、fm=3fm三種情況下抽 樣信號(hào)的波形。實(shí)驗(yàn)程序：dt=0.1; f0=1; to=1/fo;fm=fo; tm=1/fm;t=-2:dt:2;f=sin c(t);subplot(4,1,1

2、),plot(t,f,'k');axis(mi n(t) max(t) 1.1*mi n(f) 1.1*max(f); title('原連續(xù)信號(hào)和抽樣信號(hào)');for i=1:3;fs=i*fm; ts=1/fs;n=-2:ts:2;f=s inc(n);subplot(4,1,i+1),stem( n,f,'filled','k');axis(mi n(n) max( n) 1.1*mi n(f) 1.1*max (f); end實(shí)驗(yàn)截圖:(2) 求解原連續(xù)信號(hào)波形和抽樣信號(hào)所對(duì)應(yīng)的幅度譜。 實(shí)驗(yàn)程序：dt=0.1;t=-4:d

3、t:4;n=le ngth(t);f=si nc(t);tm=1;fm=1/tm; wm=2*pi*fm;k=1:n;w1=k*wm/n;f1=f*exp(-j*t'*w1)*dt; subplot(4,1,1),plot(w1/(2*pi),abs(f1);gridaxis(0 max(4*fm) 1.1*mi n(f1) 1.1*max(f1); for i=1:3;if i<= 2 c=0 ,else c=0.2,e ndfs=(4-i+c)*fm; ts=1/fs;n=-4:ts:4;f=s inc(n);n=le ngth( n); wm=2*pi*fs;k=1:n;w

4、=k*wm/n;f=f*exp(-j* n'*w)*ts; subplot(4,1,5-i),plot(w/(2*pi),abs(f),'k');grid axis(o max(4*fm) 1.1*min(f) 1.1*max(f); end實(shí)驗(yàn)截圖:3.52.543.53.543251 b 6 4 2 cs cs ci cs00.511523.525o2 1 fl 6 4 2(3) 用時(shí)域卷積方法(內(nèi)插公式)重建信號(hào)實(shí)驗(yàn)程序：f0=1; t0=1/f0; dt=0.01;fm=f0; tm=1/fm;t=-4:dt:4*f0;x=sin c(t);subplot(4,1,1),plot(t,x);axis(mi n(t),max(t),1.1*mi n(x),1.1*max(x);title('用時(shí)域卷積重建抽樣信號(hào)')for i=1:3;fs=i*fm; ts=1/fs;n=-4:(4*fo)/ts;t1=-4:ts:4*fo;x1=s inc(n /fs);t_n=o nes(le ngth( n),1)*t1- n'*ts*o nes(1,le ngth(t1) xa=x1*si nc(fs*pi*t_n); subplot(4,1,i+1),pl