fdma頻分多址通信課程設(shè)計(jì)

1、<p><b>　　信息科學(xué)與工程學(xué)院</b></p><p><b>　　課程設(shè)計(jì)報(bào)告書</b></p><p>　　課程名稱： 數(shù)字信號(hào)處理課程設(shè)計(jì) </p><p>　　班 級(jí)： 電子信息工程2007級(jí) 6 班 </p><p>　　學(xué) 號(hào)：

2、 </p><p>　　姓 名： </p><p>　　指導(dǎo)教師： </p><p>　　二○一○年 01 月 </p><p><b>　　一、課程設(shè)計(jì)目的</b></p>

3、<p>　　綜合運(yùn)用本課程的理論知識(shí)進(jìn)行頻譜分析以及濾波器設(shè)計(jì)，通過理論推導(dǎo)得出相應(yīng)結(jié)論，并利用MATLAB作為編程工具進(jìn)行計(jì)算機(jī)實(shí)現(xiàn)，從而復(fù)習(xí)鞏固了課堂所學(xué)的理論知識(shí)，提高了對(duì)所學(xué)知識(shí)的綜合應(yīng)用能力，并從實(shí)踐上初步實(shí)現(xiàn)了對(duì)數(shù)字信號(hào)的處理。</p><p><b>　　課程設(shè)計(jì)內(nèi)容及要求</b></p><p><b>　　設(shè)計(jì)內(nèi)容：</b&

4、gt;</p><p>　　在Matlab 環(huán)境中，利用編程方法對(duì)FDMA通信模型進(jìn)行仿真研究。</p><p><b>　　設(shè)計(jì)要求：</b></p><p>　?。?）Matlab支持麥克風(fēng)，可直接進(jìn)行聲音的錄制，要求錄制3路不同人的語音信號(hào)，并對(duì)錄制的信號(hào)進(jìn)行采樣；畫出采樣后語音信號(hào)的時(shí)域波形和頻譜圖。</p><p&

5、gt;　?。?）將各路語音信號(hào)分別與各自的不同的高頻載波信號(hào)相乘，由于各高頻載波信號(hào)將各語音信號(hào)頻譜移到不同頻段，復(fù)用信號(hào)頻譜為各信號(hào)頻譜的疊加，因此，只需傳輸該復(fù)用信號(hào)便可在同一信道上實(shí)現(xiàn)各路語音信號(hào)的同時(shí)傳輸。畫出復(fù)用信號(hào)的頻譜圖。</p><p>　?。?）傳輸完成后，通過選擇合適的帶通濾波器，即可獲得各個(gè)已調(diào)信號(hào)；給定濾波器的性能指標(biāo)，設(shè)計(jì)合適數(shù)字濾波器，并畫出帶通濾波器的頻率響應(yīng)。</p>

6、<p>　?。?）再進(jìn)行解調(diào)，即將各個(gè)已調(diào)信號(hào)分別乘以各自的高頻載波信號(hào)，這樣，原始低頻信號(hào)被移到低頻段。畫出解調(diào)后3路信號(hào)各自的頻譜圖。</p><p>　　（5）最后通過選擇合適的低通濾波器恢復(fù)出各原始語音信號(hào)，從而實(shí)現(xiàn)FDMA通信傳輸。畫出低通濾波器的頻率響應(yīng)，恢復(fù)信號(hào)的時(shí)域波形和頻譜，并對(duì)濾波前后的信號(hào)進(jìn)行對(duì)比，分析信號(hào)的變化?；胤耪Z音信號(hào)。</p><p><b

7、>　　詳細(xì)設(shè)計(jì)過程</b></p><p>　　錄音模塊：首先，Matlab中支持麥克風(fēng)，因此可直接錄制三段語音信號(hào)，并保存和播放。由時(shí)域采樣定理知采樣頻率要大于信號(hào)高頻率的兩倍，而且調(diào)制時(shí)要將音頻調(diào)到不同的頻段，所以選取采樣頻率時(shí)，可選高一點(diǎn)，便于調(diào)制?？蛇xfs為44100Hz, 為CD音質(zhì)。</p><p>　　調(diào)制與解調(diào)模塊：首先設(shè)計(jì)載波頻率，各個(gè)載頻的間隔既要大于

8、2倍的聲音的最高截止頻率,而且最高的載頻與采樣頻率也必須大于2倍的聲音的最高截止頻率。因?yàn)閒s為44100Hz，所以選擇載波頻率時(shí)，最好不超過fs/2，否則Ω+ω 容易與ω+fs-Ω相混淆（ω為載頻，Ω為聲音最高截止頻率，f為采樣頻率）。綜合考慮三個(gè)載頻分別取為4000Hz,11000Hz,18000Hz。然后將三路調(diào)制信號(hào)相加得到復(fù)用信號(hào)，并對(duì)其進(jìn)行頻譜分析。接著經(jīng)過帶通濾波出三路信號(hào)后，再解調(diào)并分析各路信號(hào)的頻譜。</p>

9、;<p>　　濾波器模塊：為達(dá)到較好效果，采用切比雪夫2型濾波器，并且取Rp=0.5;Rs=40; 三個(gè)帶通濾波器分別要濾出三路信號(hào)，其通頻帶要依據(jù)載頻而定，可以濾出上邊頻，也可以濾出下邊頻。綜合最后的試聽效果，分別同頻帶取4300~7400Hz，7600~11000Hz，14600~18000Hz。由于三個(gè)信號(hào)都是人的聲音信號(hào)，有著相同的頻率范圍，所以可以共用一個(gè)低通濾波器，且低通濾波器的通帶范圍為0~3400Hz。&l

10、t;/p><p><b>　　代碼及圖形</b></p><p>　　%(1)獲取錄音文件</p><p>　　fs=44100; %聲音的采樣頻率為44.1khz</p><p>　　duration=3; %錄音的時(shí)間<

11、/p><p>　　fprintf('按任意鍵開始錄音1:\n');</p><p><b>　　pause</b></p><p>　　fprintf('錄音中……\n');</p><p>　　sd1=wavrecord(duration*fs,fs); %duration*fs每次

12、獲得總的采樣數(shù)為132300,保存聲音</p><p>　　fprintf('放音中……\n'); %文件名為s1,以下類同.</p><p>　　wavplay(sd1,fs);</p><p>　　fprintf('錄音1播放完畢。\n');</p><p>　　wavwrite(sd1

13、,fs,'sound1.wav'); %將錄音文件保存為wav格式的聲音文件，下同</p><p>　　fprintf('按任意鍵開始錄音2:\n');</p><p><b>　　pause</b></p><p>　　fprintf('錄音中……\n');</p><

14、p>　　sd2=wavrecord(duration*fs,fs);</p><p>　　fprintf('放音中……\n');</p><p>　　wavplay(sd2,fs);</p><p>　　fprintf('錄音2播放完畢。\n');</p><p>　　wavwrite(sd2,fs,&#

15、39;sound2.wav');</p><p>　　fprintf('按任意鍵開始錄音3:\n');</p><p><b>　　pause</b></p><p>　　fprintf('錄音中……\n');</p><p>　　sd3=wavrecord(duration*fs

16、,fs);</p><p>　　fprintf('放音中……\n');</p><p>　　wavplay(sd3,fs);</p><p>　　fprintf('錄音3播放完畢。\n');</p><p>　　wavwrite(sd3,fs,'sound3.wav');</p>

17、<p>　　%(2)聲音樣本的時(shí)域和頻域分析</p><p>　　fs=44100; %聲音的采樣頻率為44.1khz</p><p>　　duration=3; </p><p>　　t=0:duration*fs-1; %總的采樣數(shù)&

18、lt;/p><p>　　[s1,fs]=wavread('sound1.wav'); %打開保存的錄音文件</p><p>　　[s2,fs]=wavread('sound2.wav');</p><p>　　[s3,fs]=wavread('sound3.wav');</p>

19、<p>　　figure(1) %圖一為三個(gè)聲音樣本的時(shí)域波形</p><p>　　subplot(311)</p><p>　　plot(t,s1);xlabel('單位：s');ylabel('幅度');</p><p>　　title('三個(gè)聲音樣本的時(shí)

20、域波形');</p><p>　　subplot(312)</p><p>　　plot(t,s2);xlabel('單位：s');ylabel('幅度');</p><p>　　subplot(313)</p><p>　　plot(t,s3);xlabel('單位：s');ylabe

21、l('幅度');</p><p>　　figure(2) %圖二為三個(gè)聲音樣本的頻譜分析</p><p>　　subplot(311)</p><p>　　stem(t,abs(fft(s1)),'.');xlabel('單位:Hz');ylabel('幅度

22、'); %fft對(duì)聲音信號(hào)進(jìn)行快速傅里葉變換</p><p>　　title('三個(gè)聲音樣本的頻譜分析');</p><p>　　subplot(312)</p><p>　　stem(t,abs(fft(s2)),'.');xlabel('單位:Hz');ylabel('幅度');&l

23、t;/p><p>　　subplot(313)</p><p>　　stem(t,abs(fft(s3)),'.');xlabel('單位:Hz');ylabel('幅度');</p><p>　　%(3)調(diào)制,將三個(gè)聲音信號(hào)用高頻載波</p><p>　　x1=4*s1'.*cos(2*p

24、i*4000*t/fs);</p><p>　　x2=4*s2'.*cos(2*pi*11000*t/fs);</p><p>　　x3=4*s3'.*cos(2*pi*18000*t/fs);</p><p>　　s=x1+x2+x3;%復(fù)用信號(hào)頻譜為各信號(hào)頻譜的疊加</p><p>　　figure(3)

25、%圖三為復(fù)用信號(hào)的頻譜分析</p><p>　　stem(t,abs(fft(s)),'.');xlabel('單位:Hz');ylabel('幅度');</p><p>　　title('復(fù)用信號(hào)的頻譜分析');</p><p>　　%(4)帶通濾波器設(shè)計(jì)</p>&

26、lt;p>　　Rp=0.5;Rs=40; %用切比雪夫2型設(shè)計(jì)帶通濾波器1； </p><p>　　Wp1=[4000 8000]/22050; %fs/2=22050 </p><p>　　Ws1=[3800 8500]/22050;</p&

27、gt;<p>　　[n1,Wn1]=cheb2ord(Wp1,Ws1,Rp,Rs);</p><p>　　[b1,a1]=cheby2(n1,Rs,Wn1);</p><p>　　[h1,w1]=freqz(b1,a1);</p><p>　　mag1=abs(h1);</p><p>　　db1=20*log10((mag1+

28、eps)/max(mag1));</p><p>　　Wp2=[9000 13000]/22050; %用切比雪夫2型設(shè)計(jì)帶通濾波器2； </p><p>　　Ws2=[8000 14000]/22050;</p><p>　　[n2,Wn2]=cheb2ord(Wp2,Ws2,Rp,Rs);</p><p

29、>　　[b2,a2]=cheby2(n2,Rs,Wn2);</p><p>　　[h2,w2]=freqz(b2,a2);</p><p>　　mag2=abs(h2);</p><p>　　db2=20*log10((mag2+eps)/max(mag2));</p><p>　　Wp3=[14500 18500]/22050;

30、 %用切比雪夫2型設(shè)計(jì)帶通濾波器3； </p><p>　　Ws3=[14000 19000]/22050;</p><p>　　[n3,Wn3]=cheb2ord(Wp3,Ws3,Rp,Rs);</p><p>　　[b3,a3]=cheby2(n3,Rs,Wn3);</p>

31、;<p>　　[h3,w3]=freqz(b3,a3);</p><p>　　mag3=abs(h3);</p><p>　　db3=20*log10((mag3+eps)/max(mag3));</p><p>　　figure(4);</p><p>　　subplot(3,1,1);</p><p>

32、;　　plot(w1/pi,db1);axis([0 1 -50 20]);xlabel('w/pi');ylabel('20lg|H(ejw)|');</p><p>　　title('用切比雪夫2型設(shè)計(jì)三個(gè)帶通濾波器');</p><p>　　subplot(3,1,2);</p><p>　　plot(w2/pi

33、,db2);axis([0 1 -50 20]);xlabel('w/pi');ylabel('20lg|H(ejw)|');</p><p>　　subplot(3,1,3);</p><p>　　plot(w3/pi,db3);axis([0 1 -50 20]);xlabel('w/pi');ylabel('20lg|H(ejw

34、)|');</p><p>　　y1=filter(b1,a1,s); %濾出三路未解調(diào)信號(hào)</p><p>　　y2=filter(b1,a1,s);</p><p>　　y3=filter(b1,a1,s);</p><p><b>　　%(5)解調(diào)</b><

35、;/p><p><b>　　fs=44100</b></p><p>　　y01=y1.*cos(2*pi*4000*t/fs); %各個(gè)已調(diào)信號(hào)分別乘以各自</p><p>　　y02=y2.*cos(2*pi*11000*t/fs);%的高頻載波信號(hào)</p><p>　　y03=y3.*cos(2*

36、pi*18000*t/fs);</p><p>　　figure(5) %圖五為解調(diào)后3路信號(hào)各自的頻譜圖</p><p>　　subplot(311)</p><p>　　stem(t,abs(fft(y01)),'.');xlabel('單位:Hz');ylabel('幅度');</p&

37、gt;<p>　　title('解調(diào)后3路信號(hào)各自的頻譜圖');</p><p>　　subplot(312)</p><p>　　stem(t,abs(fft(y02)),'.');xlabel('單位:Hz');ylabel('幅度');</p><p>　　subplot(313)&

38、lt;/p><p>　　stem(t,abs(fft(y03)),'.');xlabel('單位:Hz');ylabel('幅度');</p><p><b>　　%(6)低通濾波</b></p><p>　　Rp=0.5; %低通濾波器參數(shù)選擇</p><p&

39、gt;<b>　　Rs=40;</b></p><p>　　Wp1=3400/(22050);</p><p>　　Ws1=4000/(22050);</p><p>　　[n1,Wn1]=cheb2ord(Wp1,Ws1,Rp,Rs); %采用切比雪夫2型(cheby2)帶通濾波器</p><p>　　[b1,a

40、1]=cheby2(n1,Rs,Wn1);</p><p>　　[h1,w1]=freqz(b1,a1);</p><p>　　mag1=abs(h1);</p><p>　　db1=20*log10((mag1+eps)/max(mag1));</p><p>　　figure(6) %圖六為低通濾波器的頻率響應(yīng)<

41、/p><p>　　plot(w1/pi,db1);axis([0 1 -50 20]);xlabel('w/pi');ylabel('20lg|H(ejw)|');</p><p>　　title('低通濾波器的頻率響應(yīng)');</p><p>　　%（7）回復(fù)信號(hào)的時(shí)域波形和頻譜分析</p><p>

42、;　　yy1=filter(b1,a1,y01);</p><p>　　yy2=filter(b1,a1,y02);</p><p>　　yy3=filter(b1,a1,y03);</p><p>　　figure(7) %圖七為恢復(fù)信號(hào)的時(shí)域波形</p><p>　　subplot(311)</p>&

43、lt;p>　　plot(t,yy1);xlabel('單位：s');ylabel('幅度');</p><p>　　subplot(312)</p><p>　　plot(t,yy2);xlabel('單位：s');ylabel('幅度');</p><p>　　subplot(313)</

44、p><p>　　plot(t,yy3);xlabel('單位：s');ylabel('幅度'); </p><p>　　title('恢復(fù)信號(hào)的時(shí)域波形');</p><p>　　figure(8) %圖八為恢復(fù)信號(hào)的頻譜分析</p><p>　　subplot(311)&l

45、t;/p><p>　　stem(t,abs(fft(yy1)));xlabel('單位:Hz');ylabel('幅度');</p><p>　　subplot(312)</p><p>　　stem(t,abs(fft(yy2)));xlabel('單位:Hz');ylabel('幅度');</p&

46、gt;<p>　　subplot(313)</p><p>　　stem(t,abs(fft(yy3)));xlabel('單位:Hz');ylabel('幅度');</p><p>　　title('恢復(fù)信號(hào)的頻譜分析');</p><p>　　wavplay(yy1,fs); %恢復(fù)

47、聲音信號(hào)的再現(xiàn)</p><p>　　wavplay(yy2,fs);</p><p>　　wavplay(yy3,fs);</p><p><b>　　結(jié)果分析與體會(huì)</b></p><p><b>　　結(jié)果分析：</b></p><p>　　經(jīng)過不斷的修改調(diào)試，在MATLA

48、B上仿真頻分多址通信技術(shù)取得了較好的效果。錄音的聲音在經(jīng)過調(diào)試和解調(diào)后的信號(hào)與原來相比較為接近。我覺得仿真的成功關(guān)鍵在于載波頻率的選擇以及帶通和低通濾波器的參數(shù)設(shè)置。（詳細(xì)設(shè)計(jì)中已有介紹，不在贅述）</p><p>　　另外在低通濾波階段，得到的恢復(fù)信號(hào)與原始信號(hào)基本一致，但在t=0附近有所失真，這是由于頻譜混疊所致，各信號(hào)頻譜混疊部分均為高頻部分，恢復(fù)信號(hào)在附近的波峰變化最快。即為頻率最高的區(qū)域，引起高頻部分失

49、真，這是因?yàn)殇浺羝陂g引入頻率高于語音信號(hào)的噪聲，所以如果在完全無噪音的環(huán)境中進(jìn)行錄音，可得無失真的恢復(fù)信號(hào)。最后仿真結(jié)果分析表明，信號(hào)在頻分復(fù)用時(shí)還存在著頻間干擾的問題。對(duì)此，采用了適當(dāng)加大采樣頻率的方法，在較大程度打夯使該問題得以解決。至于完全消除頻譜間的干擾，還有待進(jìn)一步研究的研究與完善。</p><p><b>　　體會(huì)： </b></p><p>　　這次

50、課程設(shè)計(jì)中，我不僅復(fù)習(xí)鞏固了課堂所學(xué)的理論知識(shí)，提高了對(duì)所學(xué)知識(shí)的綜合應(yīng)用能力，并從實(shí)踐上初步實(shí)現(xiàn)了對(duì)數(shù)字信號(hào)的處理。而且通過這次較綜合的課程設(shè)計(jì)，對(duì)以前在課本中問題沒有弄懂的問題，我都有了更深入的理解。比如時(shí)域采樣定理、濾波器參數(shù)設(shè)計(jì)等。</p><p>　　總的來說，通過數(shù)字信號(hào)處理的這次課程設(shè)計(jì)，讓我對(duì)《信號(hào)與系統(tǒng)》這門課程有了深刻的認(rèn)識(shí)，收獲很大。上學(xué)期我們專注于信號(hào)的學(xué)習(xí)，涉及到的復(fù)雜的數(shù)學(xué)計(jì)算分析一度

51、讓我失去學(xué)習(xí)《信號(hào)與系統(tǒng)》的信心和興趣。下學(xué)期的系統(tǒng)方面或者說數(shù)字信號(hào)處理，以及通信相關(guān)知識(shí)的學(xué)習(xí)，我重新認(rèn)識(shí)了信號(hào)處理這門課的地位，更重要的是對(duì)于將來學(xué)習(xí)通信技術(shù)我也增強(qiáng)了很多信心！</p><p><b>　　六、附錄或參考資料</b></p><p>　　頻分多址接入模型設(shè)計(jì)及MATLAB仿真計(jì)算</p><p>　　作者：陳慧慧,鄭賓