【正文】 計(jì)： \n39。)。 snr=10*log10((s*s39。)。 figure(3) stem(t,abs(fft(s)),39。)。xlabel(39。單位： Hz39。)。 pause figure(2) %圖二為三個(gè)聲音樣本的頻譜分析 subplot(311) stem(t,abs(fft(sd1)),39。單位： s39。xlabel(39。)。39。)。錄音 3 播放完畢。 fprintf(39。)。錄音 2 播放完畢。 fprintf(39。)。錄音 1 播放完畢。 %duration*fs 每次獲得總的采樣數(shù)為 132300，保存聲音文件名為 sd1 fprintf(39。)。在制作畢業(yè)設(shè)計(jì)過(guò)程中我曾經(jīng)向老師們和同學(xué)們請(qǐng)教過(guò)不少的問(wèn)題，老師們的熱情解答和同學(xué)們的熱心幫助才使我的畢業(yè)設(shè)計(jì)能較為順利的完成。同理也可以設(shè)計(jì)出另外兩路信號(hào)的恢復(fù)。獲得的解調(diào)后的 3 路信號(hào)的頻譜如下圖所示。若 x 是復(fù)數(shù)的的情況下，則加入的噪聲為復(fù)噪聲。運(yùn)用 MATLAB 中 plot(x,y)繪圖命令 ，以 y為 縱 坐標(biāo)，x為 橫 坐標(biāo)， 來(lái)繪制圖形。因?yàn)樾盘?hào)噪聲是疊加在信號(hào)上的，而且噪聲是永遠(yuǎn)存在與通信系統(tǒng)中的沒(méi)有信號(hào)傳輸時(shí)也有噪聲，所以稱噪聲為加性噪聲 ]1[ 。 調(diào)用形式為： R s )R p ,W s ,p,c h e b2 o r d ( WW n ][ n , ? （ 7） 通過(guò)上式求得濾波器的 n 和 Wn 后，在通過(guò) MATLAB 軟件中的 cheby2 函數(shù)的調(diào)用得到濾波器 傳遞函數(shù)的 分母系數(shù) a 和 濾波器傳遞函數(shù)的分子 系數(shù) b。當(dāng) Wp 和Ws為二元矢量時(shí)，為帶通或帶阻濾波器。采樣好的信號(hào)用 MATLAB 中的 wavwrite來(lái)保存信號(hào)。 由頻譜分析理論可知：要保證抽樣的頻譜與信號(hào)的頻譜相一致則要使抽樣的頻率不發(fā)生重疊。在通過(guò)不同的帶通濾波器將信號(hào)提取出來(lái)，在通過(guò)解調(diào)和低通濾波器把原始信號(hào)恢復(fù)出來(lái)。 在 FDMA 通信系統(tǒng)中，把傳輸頻帶劃分成互不重疊且較窄的子頻帶若干個(gè)。 頻分復(fù)用是將信號(hào)資源劃分為多個(gè)子頻帶，每個(gè)子頻帶占用不同的頻率。 MATLAB技術(shù)可以對(duì)非線性動(dòng)態(tài)系統(tǒng)建模和仿真為科研工程設(shè)計(jì)以及數(shù)值計(jì)算提供了解決方案 [3]。 本人簽名： 日期： 巢湖學(xué)院本科畢業(yè) 論文 (設(shè)計(jì) )使用授權(quán)說(shuō)明 本人完全了解巢湖學(xué)院有 關(guān)收集、保留和使用畢業(yè)論文 (設(shè)計(jì) )的規(guī)定，即：本科生在校期間進(jìn)行畢業(yè)論文 (設(shè)計(jì) )工作的知識(shí)產(chǎn)權(quán)單位屬巢湖學(xué)院。本人完全意識(shí)到本聲明的法律結(jié)果由本人承擔(dān)。缺點(diǎn)就是濾波器特性不夠理想和信道間有噪聲干擾 [2]。復(fù)用也就是在一條信道上同時(shí)傳輸幾路信號(hào)以提高資源 [2]。雖然混頻本身不是線性設(shè)備，但從頻譜角度看只進(jìn)行了簡(jiǎn)單的搬移，所以就可以把混頻器所對(duì)應(yīng)的調(diào)制稱為“線性調(diào)劑”。運(yùn)用不同的載波對(duì)不同的信號(hào)進(jìn)行調(diào)制，把載波帶寬劃分為不同頻帶的子信道，每個(gè)并行信號(hào)用一個(gè)子信道傳送。則有 Hs ff 2? （ 3） 當(dāng) Hs ff 2? 是奈奎斯特頻率 [1]。分別為： 8000Hz， 11025Hz， 22050Hz， 44100Hz以保證信號(hào)的質(zhì)量，有利用實(shí)驗(yàn)的進(jìn)行。 Rp 為帶通濾波器最大的衰減系數(shù)， Rs為阻帶最小衰減系數(shù)， Wp 為帶通邊界頻率的歸一化值， Ws 為阻帶邊界頻率的歸一化值。 當(dāng)采 樣頻率為 441KHz 時(shí)， Ws1 和 Wp1 分別為 Ws1=[4100 7600]/22050 Wp1=[4200 7500]/22050 當(dāng)帶通濾波器的 4個(gè)參數(shù)都確定后，就可使用 MATLAB 的 cheb2ord 函數(shù)求第一個(gè)濾波器的截止頻率 Wn 和最小階數(shù) n。 通信系統(tǒng)中的噪聲是存在于信道中疊加在信號(hào)上不需要的電信號(hào)?？芍谡麄€(gè) figure 中 共有 m個(gè)以 n 個(gè)圖排列成的行，而p表示在哪個(gè)圖上話曲線。 在設(shè)計(jì)過(guò)程中通過(guò)使用 MATLAB 的 awgn 函數(shù)對(duì)信號(hào)加入高斯白噪聲，它的調(diào)用形式為： y = awgn(x,SNR) [6] 式中 SNR 為以 dB 為單位的信噪比， x 為加入信號(hào)中的高斯白噪聲。然后巢湖學(xué)院 2020 屆畢業(yè)論文（設(shè)計(jì)） 11 通過(guò) MATLAB 軟件中的 fft 進(jìn)行快速傅里葉變換。 圖 15 信號(hào)恢復(fù)后的頻譜圖 通 過(guò)上面的仿真我們可以恢復(fù) 1 路語(yǔ)音信號(hào)，然后使用 MATLAB 軟件對(duì)信號(hào)進(jìn)行播放，比較與原信號(hào)的區(qū)別。離別在即，站在人生的又一個(gè)轉(zhuǎn)折點(diǎn)上，心中難免思緒萬(wàn)千，一種感恩之情油然而生，感謝 巢湖學(xué)院電子工程與電氣自動(dòng)化 學(xué)院的所有專業(yè)老師，沒(méi)有你們的辛勤勞動(dòng)，就沒(méi)有我們今日的滿載而歸，感謝大學(xué)四年曾經(jīng)幫助過(guò)我的所有同學(xué)。按任意鍵開(kāi)始錄音 1： \n39。 sd1=wavrecord(duration*fs,fs)。 fprintf(39。按任意鍵開(kāi)始錄音 2： \n39。 sd2=wavrecord(duration*fs,fs)。 fprintf(39。按 任意鍵開(kāi)始錄音 3： \n39。 sd3=wavrecord(duration*fs,fs)。 fprintf(39。按任意鍵開(kāi)始聲音樣本的時(shí)域分析： \n39。 %打開(kāi)保存的錄音文件 [sd2,fs]=wavread(39。單位： s39。 subplot(312) plot(t,sd2)。xlabel(39。)。幅度 39。xlabel(39。)。按任意鍵開(kāi)始信號(hào)的調(diào)制和復(fù)用信號(hào)頻域分析： \n39。 s=x1+x2+x3。幅度 39。 %加入高斯白噪聲 pause ys=awgn(s,20)。單位： Hz39。 %（ 5）帶通濾波器的設(shè)計(jì) fprintf(39。 [b1,a1]=cheby2(n1,Rs,Wn1)。 [h2,w2]=freqz(b2,a2)。 mag3=abs(h3)。)。 plot(w2/pi,db2)。 subplot(3,1,3)。)。 y01=y1.*cos(2*pi*4000*t/fs)。)。.39。 subplot(313) stem(t,abs(fft(y03)),39。)。 [n1,Wn1]=cheb2ord(Wp1,Ws1,Rp,Rs)。 xlabel(39。)。xlabel(39。)。 subplot(313) plot(t,yy3)。按任意鍵開(kāi)始 3 路信號(hào)頻域分析： \n39。)。ylabel(39。幅度 39。)。 wavplay(yy2,fs)。)。xlabel(39。 subplot(312) stem(t,abs(fft(yy2)))。單位： Hz39。ylabel(39。)。幅度 39。 pause yy1=filter(b1,a1,y01)。20lg|H(ejw)|39。 db1=20*log10((mag1+eps)/max(mag1))。 pause Rp=。單位 ： Hz39。)。 title(39。.39。 %(6)解調(diào) fprintf(39。w/pi39。)。 title(39。 plot(w1/pi,db1)。 Ws3=[14000 19000]/22050。 Wp2=[9000 13000]/22050。 Rs=40。)。 snr %計(jì)算信噪比 figure(4) stem(t,abs(fft(ys)),39。)。xlabel(39。 x2=4*sd239。ylabel(39。幅度 39。)。 %fft 對(duì)聲音信號(hào)進(jìn)行快速傅里葉變換 xlabel(39。幅度 39。ylabel(39。)。39。 t=0:duration*fs1。 wavwrite(sd3,f