【正文】 1)=sum((xn1.*xm2)39。)ylabel(39。)for N=1:11xm1=[fliplr(xn1(1:N)) fliplr(xn1(N+1:12))]。x1(Nm)39。)。y(n)39。stem(n,xm1)xlabel(39。)title(39。figure()stem(n,yn)xlabel(39。)title(39。(3) 聲音文件線性卷積序列x1(n)=[2 0 1 2 5 7 0 5],讀取一段聲音數(shù)據(jù)（Windows XP 菜單命令），當(dāng)循環(huán)卷積長度大于或等于兩序列長度之和時(shí)，循環(huán)卷積等于線性卷積。[xn1,fs,nbits]=wavread(39。 zeros(1,7)]。)。)endxm3=[zeros(1,686) fliplr(xn2)]。m39。n=N39。G:\Windows XP 菜單命令39。X=fft(xn2,1024)。 n=12。 a=sqrt(2)*n*pi。 subplot(5,1,1) plot(x1)。 y=fft(yn,fs)。 subplot(5,1,5) plot(yn)； 采樣定理演示1② 要求輸入采樣頻率fs（根據(jù)程序處理需要指定范圍）后，在時(shí)域演示信號(hào)波形、采樣脈沖及采樣后信號(hào)；在頻域演示不同采樣頻率下對(duì)應(yīng)信號(hào)的頻譜。 a=sqrt(2)*n*pi。 subplot(5,1,1) plot(x1)。 y=fft(yn,fs)。 subplot(5,1,5) plot(yn)；. 采樣定理演示2 課程設(shè)計(jì)原理 頻譜分析原理時(shí)域分析只能反映信號(hào)的幅值隨時(shí)間的變化情況，除單頻率分量的簡單波形外，很難明確提示信號(hào)的頻率組成和各頻率分量大小，而頻譜分析能很好的解決此問題。但上式的計(jì)算效率很低，因?yàn)橛写罅康闹笖?shù)(等價(jià)于三角函數(shù))運(yùn)算，故實(shí)際中多采用快速Fourier變換(FFT)。以頻率f為橫坐標(biāo)，為縱坐標(biāo)，可以得到幅值譜；以頻率f為橫坐標(biāo)，為縱坐標(biāo)，可以得到相位譜；以頻率f為橫坐標(biāo)，為縱坐標(biāo)，可以得到實(shí)頻譜；以頻率f為橫坐標(biāo)，為縱坐標(biāo)，可以得到虛頻譜。它有兩種形式：振幅特性在通帶內(nèi)是等波紋的，在阻帶內(nèi)是單調(diào)下降的切比雪夫I型；振幅特性在阻帶內(nèi)是等波紋的，在通帶內(nèi)是單調(diào)下降的切比雪夫II型。脈沖響應(yīng)不變法是從S平面映射到z平面，這種映射不是簡單的代數(shù)映射，而是S平面的每一條寬為2π/T的橫帶重復(fù)地映射到整個(gè)z平面。有限長單位沖激響應(yīng)(FIR)數(shù)字濾波器具有以下優(yōu)良的特點(diǎn)：(1) 可在設(shè)計(jì)任意幅度頻率特性濾波器的同時(shí)，保證精確、嚴(yán)格的線性相位特性。常用的窗函數(shù)[4]有以下幾種：矩形窗(Rectangular window)、三角窗(Triangular window)、漢寧窗(Hanning window)、海明窗(Hamming window)、布拉克曼窗(Blackman window)、切比雪夫窗(Chebyshev window)、巴特里特窗(Bartlett window)及凱塞窗(Kaiser window)。當(dāng)sflag＝symmetric時(shí)，為對(duì)稱采樣；當(dāng)sflag＝periodic時(shí)，為周期采樣，此時(shí)hann函數(shù)計(jì)算N+1個(gè)點(diǎn)的窗，但是僅返回前N個(gè)點(diǎn)。在MATLAB中，實(shí)現(xiàn)布拉克曼窗的函數(shù)為blackman，調(diào)用格式如下： w=blackman (N)w=blackman (N,39。各種窗函數(shù)的性能比較可見表21，在設(shè)計(jì)FIR濾波器的過程中可以根據(jù)要求選擇合適的窗函數(shù)：表21 瓣若所設(shè)計(jì)的數(shù)字濾波器是低通的，那么上述設(shè)計(jì)工作可以結(jié)束，若所設(shè)計(jì)是高通、帶通或帶阻濾波器，那么還需進(jìn)行以下步驟。雙線性Z變換滿足以上4個(gè)條件的映射關(guān)系，其變換公式為 FIR課程設(shè)計(jì)步驟和過程(1)根據(jù)對(duì)阻帶衰減及過渡帶的指標(biāo)要求，選擇窗函數(shù)的類型，并估計(jì)窗口長度N；(2)構(gòu)造希望逼近的頻率響應(yīng)函數(shù)Hd(ejw)，即Hd(ejw)=Hdg(w)ejw（N1）/2；(3)計(jì)算hd(n)；(4)加窗得到設(shè)計(jì)結(jié)果：h(n)=hd(n)w(n)。當(dāng)相移為177。當(dāng)無相移時(shí)，N為奇數(shù)時(shí)，所設(shè)計(jì)的FIR數(shù)字帶通濾波器為Ⅰ型濾波器；當(dāng)N為偶數(shù)時(shí)，為Ⅱ型濾波器。在MATLAB中通過wavplay(‘filename’)讀取語音信號(hào)數(shù)據(jù)。 原語音信號(hào)時(shí)域圖和頻譜圖。設(shè)計(jì)指標(biāo)：wpz=4500Hz。其設(shè)計(jì)程序如下：wpz=2*4500/fs。rp=2。)。[bz,az]=bilinear(b,a,fs)。 濾波后語音信號(hào)的時(shí)域圖和頻譜圖 切比雪夫高通程序的調(diào)試和運(yùn)行結(jié)果 。)。x=[rot90(x1,3),rot90(x1,3)]。 wsz=500Hz。wsz=2*500/fs。rs=100。[b,a]=cheby2(N,rs,wc,39。)。其濾波程序如下：Yk=Z.*Hk。G:\Windows 39。x1=*sin(2*pi*100*t)+*sin(2*pi*220*t)+ *sin(2*pi*6000*t)。設(shè)計(jì)指標(biāo)：wpz=4500Hz。 rp=20。wp=2*fs*tan(wpz*pi/2)。[N,wc]=cheb2ord(wp,ws,rp,rs,39。s39。 帶通濾波器的幅頻圖。其程序如下：[y,fs,nbits]=wavread(39。通過老師的指導(dǎo)，和自己的思考，總結(jié)如下：問題1：Subscript indices must either be real positive integers or logicals.解決辦法：在訪問矩陣（包括向量、二維矩陣、多維數(shù)組，下同）的過程中，下標(biāo)索引要么從 0 開始，要么出現(xiàn)了負(fù)數(shù)。2. 如果 a 是函數(shù)，則這是因?yàn)閙atlab在所有已添加的路徑中都無法找到該函數(shù)對(duì)應(yīng)的m文件而導(dǎo)致的，對(duì)此，把該m文件移動(dòng)到當(dāng)前路徑下再運(yùn)行便可[7]。解決方法：檢查打開文件目錄，一般打開文件需要在MATBLE下的work內(nèi)。(7) 數(shù)組元素求和：sum直方圖：hist(10) 矩陣分析：左右翻轉(zhuǎn)：fliplr線性變化數(shù)組：linspace，聚合矩陣：cat/horzcat/vertcat[8]二項(xiàng)分布：binornd剛開始自己編的程序運(yùn)行就錯(cuò)誤很多，當(dāng)時(shí)都沒信心能把這個(gè)課程設(shè)計(jì)完成，但什么都要自己面對(duì)的。G:\Windows 39。xlabel(39。y39。)Y=fft(y,32768)。xlabel(39。Y39。)t=0::。figure(2)subplot(2,1,1)。)。title(39。subplot(2,1,2)。)。title(39。subplot(2,1,2)。)。title(39。)。wpz=2*4500/fs。rp=2。)。[bz,az]=bilinear(b,a,fs)。f39。)。Hk=[H fliplr(H)]。f=22050*(0:16384)/32768。f39。)。subplot(2,1,1)。n39。)。 (2) 切比雪夫高通：[y,fs,nbits]=wavread(39。plot(abs(y))。ylabel(39。原語音信號(hào)時(shí)域波形39。plot(f,abs(Y(1:16385)))。ylabel(39。原語音信號(hào)頻譜圖39。y=y(1:20191,:)z=x+y。n39。)。f=22050*(0:16384)/32768。f39。)。sync39。)。ws=2*fs*tan(wsz*pi/2)。s39。,39。figure(3)Hk=freqz(bz,az,16384,fs)。)。title(39。H=Hk39。subplot(2,1,2)。)。title(39。plot(abs(yk))。)。title(39。G:\Windows 39。xlabel(39。y39。)Y=fft(y,32768)。xlabel(39。Y39。)t=0::。figure(2)subplot(2,1,1)。)。title(39。subplot(2,1,2)。)。title(39。)。wpz=[2*700/fs,2*4500/fs]。rp=20。)。[bz,az]=bilinear(b,a,fs)。f39。)。Hk=[H fliplr(H)]。f=22050*(0:16382)/32768。f39。)。subplot(2,1,1)。n39。)。 (4) hamming濾波器完整程序：[x1,fs,bits]=wavread(39。plot(x1)。X1=fft(x1)。title(39。x2=x1+*cos(2*pi*10000*t/22050000)。title(39。mag2=abs(X2)。)。Btg=wpgwsg。hng=fir1(Ng1,wcg,39。subplot(2,2,1)。)。title(39。Yg=fft(yg)。濾波后信號(hào)頻譜39。)。subplot(2,2,3)。)。加噪聲聲音濾波后波形39。wpd=pi*。Nd=N0d+mod(N0d+1,2)。subplot(2,2,1)。)。title(39。Yd=fft(yd)。濾波后信號(hào)頻譜39。)。subplot(2,2,3)。)。加噪聲聲音濾波后波形39。wpdai2=pi*。Mdai=ceil(*pi/Btdai)。subplot(2,2,1)。)。title(39。Ydai=fft(ydai)。濾波后信號(hào)頻譜39。)。subplot(2,2,3)。)。加噪聲聲音濾波后波形39。盡我所知，除文中特別加以標(biāo)注和致謝的地方外，不包含其他人或組織已經(jīng)發(fā)表或公布過的研究成果，也不包含我為獲得 及其它教育機(jī)構(gòu)的學(xué)位或?qū)W歷而使用過的材料。除了文中特別加以標(biāo)注引用的內(nèi)容外，本論文不包含任何其他個(gè)人或集體已經(jīng)發(fā)表或撰寫的成果作品。本人授權(quán)　　 　 　大學(xué)可以將本學(xué)位論文的全部或部分內(nèi)容編入有關(guān)數(shù)據(jù)庫進(jìn)行檢索，可以采用影印、縮印或掃描等復(fù)制手段保存和匯編本學(xué)位論文。、圖表要求：1）文字通順，語言流暢，書寫字跡工整，打印字體及大小符合要求，無錯(cuò)別字，不準(zhǔn)請(qǐng)他人代寫2）工程設(shè)計(jì)類題目的圖紙，要求部分用尺規(guī)繪制，部分用計(jì)算機(jī)繪制，所有圖紙應(yīng)符合國家技術(shù)標(biāo)準(zhǔn)規(guī)范。作者簽名： 日期： 年 月 日導(dǎo)師簽名： 日期： 年 月 日注 意 事 項(xiàng)（論文）的內(nèi)容包括：1）封面（按教務(wù)處制定的標(biāo)準(zhǔn)封面格式制作）2）原創(chuàng)性聲明3）中文摘要（300字左右）、關(guān)鍵詞4）外文摘要、關(guān)鍵詞 5）目次頁（附件不統(tǒng)一編入）6）論文主體部分：引言（或緒論）、正文、結(jié)論7）參考文獻(xiàn)8）致謝9）附錄（對(duì)論文支持必要時(shí)）：理工類設(shè)計(jì)（論文）正文字?jǐn)?shù)不少于1萬字（不包括圖紙、程序清單等）。本人完全意識(shí)到本聲明的法律后果由本人承擔(dān)。作 者 簽 名：　　 　 　　日　 期：　　 　 　　 指導(dǎo)教師簽名：　　 　 　　 日　　期：　　 　 　　 使用授權(quán)說明本人完全了解 大學(xué)關(guān)于收集、保存、使用畢業(yè)設(shè)計(jì)（論文）的規(guī)定，即：按照學(xué)校要求提交畢業(yè)設(shè)計(jì)（論文）的印刷本和電子版本；學(xué)校有權(quán)保存