資源描述：
《通信系統(tǒng)建模與仿真實驗報告.doc》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在應(yīng)用文檔-天天文庫。
1、實驗報告課程名稱通信系統(tǒng)建模與仿真實驗項目名稱通信系統(tǒng)建模與仿真軟件實驗實驗類型實驗學(xué)時班級20110822學(xué)號2011082227姓名楊俊良指導(dǎo)教師何忠秋實驗室名稱實驗時間實驗成績預(yù)習(xí)部分實驗過程表現(xiàn)實驗報告部分總成績教師簽字日期哈爾濱工程大學(xué)教務(wù)處制實驗一：低通采樣定理和內(nèi)插與抽取實現(xiàn)一、實驗?zāi)康挠肕atlab編程實現(xiàn)自然采樣與平頂采樣過程，根據(jù)實驗結(jié)果給出二者的結(jié)論；掌握利用MATLAB實現(xiàn)連續(xù)信號采樣、頻譜分析和采樣信號恢復(fù)的方法。二、實驗原理1．抽樣定理若是帶限信號，帶寬為,經(jīng)采樣后的頻譜就是將的頻譜在頻率軸上以采樣頻率為間隔進行周期延拓。因此，當(dāng)時，
2、不會發(fā)生頻率混疊；而當(dāng)<時將發(fā)生頻率混疊。2．信號重建經(jīng)采樣后得到信號經(jīng)理想低通則可得到重建信號，即：=*其中：==，所以：=*=*=上式表明，連續(xù)信號可以展開成抽樣函數(shù)的無窮級數(shù)。利用MATLAB中的來表示，有，所以可以得到在MATLAB中信號由重建的表達式如下：=我們選取信號=作為被采樣信號，當(dāng)采樣頻率=2時，稱為臨界采樣。我們?nèi)±硐氲屯ǖ慕刂诡l率=。下面程序?qū)崿F(xiàn)對信號=的采樣及由該采樣信號恢復(fù)重建：三、實驗內(nèi)容已知信號，試以以下采樣頻率對信號采樣：(a);(b);(c)，求x(t)信號原信號和采樣信號頻譜，及用采樣信號重建原信號x’(t)時序圖。假定重構(gòu)濾
3、波器為fs/2的理想低通濾波器，具有帶通增益為Ts=1/fs.四、實驗步驟1.設(shè)置采樣時間間隔ts=1/fs，以及時間范圍。2.輸入已知信號。3.對原始信號進行傅里葉變換。4.畫出原信號的頻譜。5.利用firpmord，firpm以及freqz，conv函數(shù)求出采樣重建信號。fp=fs/2-500;fs1=fs/2;[n,f0,m0,w]=firpmord([fpfs1],[10],[0.00010.001],fs);b=firpm(n,f0,m0,w);figure(4)freqz(b,1,1024,fs)y=conv(x,b);6.對重建信號進行傅里葉變換。
4、7.畫出重建信號的波形及頻譜。流程圖五、實驗結(jié)果分析結(jié)論心得體會：本次實驗要用到采樣定理，信號內(nèi)插和重建的一些知識，學(xué)會了如dyadup命令等一些新的指令和思考方法，感覺到受益匪淺。附程序：clearallclcfs=20000;ts=1/fs;T=0.05;df=1/T;t=0:ts:T-ts;x=0;form=0:99x=x+(m+1)*cos(2*pi*(100*m+50)*t);endfigure(1)subplot(211)plot(t,x)title('原始信號波形');xlabel('時間t');ylabel('幅度');y=dyadup(x,0)
5、;y=dyadup(y,0);[X,f]=fftseq(y,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;form=0:99x1=x1+(m+1)*cos(2*pi*(100*m+50)*t1);endfigure(2)subplot(211)plot(t1,x1)title('10000Hz抽樣波形');xlabel('時間t');yla
6、bel('幅度');y1=dyadup(x1,0);y1=dyadup(y1,0);[X1,f1]=fftseq(y1,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;form=0:99x3=x3+(m+1)*cos(2*pi*(100*m+50)*t3);endfigure(3)subplot(211)plo
7、t(t3,x3)title('30000Hz抽樣波形');xlabel('時間t');ylabel('幅度');y3=dyadup(x3,0);y3=dyadup(y3,0);[X3,f3]=fftseq(y3,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([fpfs1],[10],[0.00010.001],fs);b=firpm(n,f0,m0,
8、w);figure(4)
