【正文】 AIC23 的內(nèi)部寄存器。簡(jiǎn)單的講，可以在原聲音流中疊加延遲一段時(shí)間后的聲流，實(shí)現(xiàn)回聲效果。聲音遇到較遠(yuǎn)物體產(chǎn)生的反射會(huì)比遇到較近的反射波晚些到達(dá)聲源位置，所以回聲和原聲的延遲隨反射物體的距離大小改變。另一方面，為了語(yǔ)音識(shí)別實(shí)用化的需要，講者自適應(yīng)、聽覺模型、快速搜索識(shí)別算法以及進(jìn)一步的語(yǔ)言模型的研究等課題倍受關(guān)注。它的強(qiáng)大數(shù)據(jù)處理能力和高運(yùn)行速度，是最值得稱道的兩大特色。 DSP系統(tǒng)特點(diǎn)和設(shè)計(jì)基本原則DSP（digital signal processor）是一種獨(dú)特的微處理器，是以數(shù)字信號(hào)來處理大量信息的器件。語(yǔ)音通信的研究涉及到人工智能、數(shù)字信號(hào)處理、微型計(jì)算機(jī)技術(shù)、語(yǔ)言聲學(xué)、語(yǔ)言學(xué)等許多領(lǐng)域，所以說語(yǔ)音的通信是一個(gè)多學(xué)科的綜合研究領(lǐng)域，其研究成果具有重要的學(xué)術(shù)價(jià)值。通常這些信號(hào)處理的過程要滿足實(shí)時(shí)且快速高效的要求，隨著DSP技術(shù)的發(fā)展，以DSP為內(nèi)核的設(shè)備越來越多，為語(yǔ)音信號(hào)的處理提供了良好的平臺(tái)。FIR數(shù)字濾波器的設(shè)計(jì)出來的是兩種窗的圖形，通過兩種窗的比較，我了解了他們各自的特點(diǎn)，幅頻和相頻特性。對(duì)這種語(yǔ)言環(huán)境也有了新的了解。矩形窗設(shè)計(jì)振幅譜39。H振幅39。xlabel(39。h39。xlabel(39。title(39。%矩形窗函數(shù)設(shè)計(jì)的頻響 magH=abs(H)。h=(1).^n.*sin((i).*Wc)./((i).*pi)。Wc=*pi。p0163。用矩形窗設(shè)計(jì)一個(gè)線性相位高通濾波器。title(39。w2=。grid on。freqz(b,1,512)。w1=。（3）根據(jù)過渡帶寬度和阻帶最小衰減，確定窗函數(shù)類型和窗口長(zhǎng)度N；（4）求FIR濾波器單位脈沖響應(yīng)；（5）分析幅頻特性，若不滿足要求，可適當(dāng)改變窗函數(shù)形式或長(zhǎng)度N，重復(fù)上述設(shè)計(jì)過程，以得到滿意的結(jié)果。對(duì)此，采取兩項(xiàng)措施：一是將h（n）截短；二是將其右移。實(shí)驗(yàn)七 FIR數(shù)字濾波器的設(shè)計(jì)一、實(shí)驗(yàn)?zāi)康恼莆沼么昂瘮?shù)設(shè)計(jì)FIR數(shù)字濾波的原理及其設(shè)計(jì)步驟；熟悉線性相位數(shù)字濾波器的特性。)。plot(W/pi,20*log10(abs(H)))。頻率39。subplot(2,1,1)。[Bap,Aap]=zp2tf(Z,P,K)。[N,Wn]=buttord(wp2,ws2,rp,rs,39。wp1=wp*ts。rp=1。)。頻率/Hz39。figure(2)。[H,W]=zp2tf(Z,P,K)。)。179。209。201。%通帶最大衰減 rs=35。基于chebyshev I型模擬濾波器原型使用沖激不變轉(zhuǎn)換方法設(shè)計(jì)數(shù)字濾波器，要求參數(shù)為通帶截止頻率Wp=；通帶最大衰減Ap=1dB；阻帶截止頻率Ws=；阻帶最小衰減As=35dB。MATLAB提供了一組標(biāo)準(zhǔn)的數(shù)字濾波器設(shè)計(jì)函數(shù)，大大簡(jiǎn)化了濾波器的設(shè)計(jì)工程。一般來說，在要求時(shí)域沖激響應(yīng)能模仿模擬濾波器的場(chǎng)合，一般使用沖激響應(yīng)不變法。二、實(shí)驗(yàn)原理在IIR濾波器的設(shè)計(jì)中，常用的方法是：先根據(jù)設(shè)計(jì)要求尋找一個(gè)合適的模擬原型濾波器，然后根據(jù)一定的準(zhǔn)則將此模擬原型濾波器轉(zhuǎn)換為數(shù)字濾波器。)。相位部分39。xlabel(39。振幅39。title(39。零極點(diǎn)圖39。angx=angle(x).*180/pi。x2=*exp(1*j*w)+*exp(2*j*w)+*exp(3*j*w)+*exp(4*j*w)+*exp(5*j*w)。（一）程序b=[,]。MATLAB提供了求離散時(shí)間系統(tǒng)頻響特性的函數(shù)freqz，調(diào)用freqz的格式主要有兩種。離散系統(tǒng)中z變化建立了時(shí)域函數(shù)h（n）與z域函數(shù)H（z）之間的對(duì)應(yīng)關(guān)系。在MATLAB中系統(tǒng)函數(shù)的零極點(diǎn)可通過函數(shù)roots得到，也可借助函數(shù)tf2zp得到，tf2zp的語(yǔ)句格式為[Z,P,K]=tf2zp(B,A)其中，B與A分別表示H（z）的分子與分母多項(xiàng)式的系數(shù)向量。)。subplot(2,1,2)。k=0:1:N1。)。subplot(2,1,2)。k=0:1:N1。axis([0 11 ])。)。frequency in pi unit39。subplot(224)。title(39。)。plot(w/pi,realX)。相位 part39。ylabel(39。grid。)。magnitude39。xlabel(39。imagX=imag(H)。w=[0:500]*pi/500 H=freqz(x,y,w)。X(e)=+(n)=cosjwjw++++（２）計(jì)算32點(diǎn)序列5pn16，0≦n≦31的32點(diǎn)和64點(diǎn)DFT，分別繪出幅度譜圖形，并繪出該序列的DTFT圖形。在本實(shí)驗(yàn)中，將學(xué)習(xí)利用MATLAB計(jì)算離散時(shí)間信號(hào)的DTFT和DFT，并加深對(duì)其相互關(guān)系的理解。grid on。stem(n,x)。(2)N=10。title(39。n=0:N1。grid on。grid on。%學(xué)號(hào)01 x=exp((2+j*w)*n)。分析:該序列具有周期性,、編制程序產(chǎn)生復(fù)正弦序列x(n)=e(2+j學(xué)號(hào)后兩位）n，并繪出其圖形。title(39。phi=0。title(39。%學(xué)號(hào)01 phi=0。grid on。stem(n,x)。w=2*pi。學(xué)號(hào)后兩位248。246。單位階躍序列39。y=[zeros(1,NM+1)ones(1,NM)]。分析:在上圖的基礎(chǔ)上平移了1個(gè)單位.(3)N=6。title(39。M=1。單位階躍序列39。n=(N1):N1。0編制程序產(chǎn)生單位階躍序列u(n)=237。grid on。stem(n,y)。(2)N=6。title(39。程序:(1)N=4。1n=0編制程序產(chǎn)生單位沖激序列d(n)=237。n=0：N1 x=a.^n stem(n,x)單邊指數(shù)序列n的取值范圍為n=0。與連續(xù)的正弦信號(hào)不同，正弦序列的自變量n必須為整數(shù)。（3）矩形序列矩形序列有一個(gè)重要的參數(shù)，就是序列的寬度N。離散序列通常用x（n）來表示，自變量必須是整數(shù)。而要研究離散時(shí)間信號(hào)，首先需要產(chǎn)生出各種離散時(shí)間信號(hào)。在以后的學(xué)習(xí)中還需要深入了解這方面的內(nèi)容。在設(shè)計(jì)的過程中遇到了很多問題，剛剛開始時(shí)曾天真的認(rèn)為只要把以前的程序改了參數(shù)就可以用了，可是問題沒有我想象中的那么簡(jiǎn)單，單純的搬程序是不能解決問題的。頻率圖中，女生高 27頻成分較多。，把高頻濾除，聲音變得比較低沉。d:39。%顯示標(biāo)表示的意義 grid。)。低通濾波器濾波后信號(hào)的頻譜圖39。%顯示標(biāo)表示的意義 grid。)。subplot(3,2,5)。%butter函數(shù)（num為分子系數(shù)den為分母系數(shù)）[h,w]=freqz(num,den)。%阻帶截至頻率 Rp=1。grid。)。疊加信號(hào)的頻譜圖39。y0=fft(yy)。ylabel(39。)。LineWidth39。f=[(n11)*(2*pi/N)]/pi。%將兩個(gè)聲音疊加：wavplay(yy,fs)。%對(duì)錄入信號(hào)進(jìn)行采樣[x1,fs,bits]=wavread(39。fs=8000。)。%快速傅立葉變換subplot(3,2,2)。ylabel(39。%信號(hào)的時(shí)域波形 title(39。f=[(n11)*(2*pi/N)]/pi。%對(duì)語(yǔ)音信號(hào)進(jìn)行采樣wavplay(x1,fs)。)。pause。)。double39。%暫停命令 fprintf(39。%錄音時(shí)間 fprintf(39。)。頻率w/pi39。%原始信號(hào)的頻譜圖 title(39。%顯示縱坐標(biāo)的表示意義 grid。)。%頻帶寬度 subplot(2,2,3)。%播放原始信號(hào) N=length(y)。%加網(wǎng)格[y,fs,bits]=wavread(39。%顯示橫坐標(biāo)表示的意義 ylabel(39。)。%快速傅立葉變換 subplot(2,2,2)。幅值/A39。%錄制信號(hào)的時(shí)域波形title(39。f=[(n11)*(2*pi/N)]/pi。%播放原始信號(hào),解決落后半拍wavplay(x,fs)。)。pause。)。double39。%暫停命令 fprintf(39。%錄音時(shí)間 fprintf(39。%sound(x, Fs, bits)。%網(wǎng)格wavplay(ys1,8000)。%顯示標(biāo)表示的意義 title(39。)。%在四行兩列的第八個(gè)窗口顯示圖形plot(f, abs(fftwave(n1)))。)。頻率w/pi39。%繪制后信號(hào)的時(shí)域的圖形 title(39。ys1= fftfilt(h,x)。h=fir1(n,wn/pi,39。wn1=(ws1+wc1)/2。wc2=2*pi*3200/fs。%加網(wǎng)格wavplay(ys,8000)。%顯示標(biāo)表示的意義title(39。)。%在四行兩列的第七個(gè)窗口顯示圖形 plot(f,abs(fftwave(n1)))。)。高通濾波后信號(hào)的時(shí)域波形39。ylabel(39。%在四行兩列的第五個(gè)窗口顯示圖形 plot(ys)。,kaiser(n+2,beta))。wn=(ws+wc)/2。beta=。)。幅值/A39。%繪制模值 xlabel(39。)。頻率w/pi39。%信號(hào)的時(shí)域波形title(39。%信號(hào)送入濾波器濾波，ys為輸出 fftwave=fft(ys)。band39。width=wswc。%加網(wǎng)格3．1低通濾波器濾波 fs=8000。%顯示縱坐標(biāo)表示的意義 title(39。)。原始信號(hào)的頻譜圖39。%加網(wǎng)格y0=fft(x)。%加標(biāo)題 ylabel(39。plot(x)。%采樣間隔 n1=1:N/2。)。)。頻率（HZ）39。title(39。[Bz,Az]=bilinear(B,A,Fs1)。)。wap=2*tan(wp/2)。wp=2*pi*fp/Fs。fs=[1000,3200]。)。頻率（HZ）39。title(39。[Bz,Az]=bilinear(B,A,Fs1)。high39。10[N,wc]=buttord(wap,was,rp,rs,39。T=1。rp=1。9 巴特沃思低通濾波器MATLAB程序如下。ylabel(39。)。[h,w]=freqz(Bz,Az,512,Fs1*22050)。s39。[N,wc]=buttord(wap,was,rp,rs,39。ws=2*pi*fs/Fs。Fs=22050。grid。)。截取后的聲音時(shí)域圖形39。C:39。%加標(biāo)題grid。)。頻率w/pi39。%原始信號(hào)的頻譜圖 title(39。%加網(wǎng)格y0=fft(x)。%加標(biāo)題 ylabel(39。plot(x)。%采樣間隔 n1=1:N/2。C:39。)。錄音結(jié)束39。x=wavrecord(time*Fs,Fs,39。%文字提示 pause。%抽樣頻率 time=3。奈奎斯特采樣定理：要使實(shí)信號(hào)采樣后能夠不失真還原，采樣頻率必須大于信號(hào)最高頻率的兩倍：Ws179。a239。?(jW)與X(jW)的關(guān)系：理想抽樣信號(hào)的頻譜是連續(xù)信號(hào)頻譜的Xaa周期延拓,重復(fù)周期為Ws(采樣角頻率)。?Xa(jW)=229。xXa(jW)171。dT(t)=頻域描述：利用傅氏變換的性質(zhì)，時(shí)域相乘頻域卷積，若n=165。xa(nT)d(tnT)xn=165。 DF線性相位的概念，即線性相位對(duì)h(n)、H(w)及零點(diǎn)的約束，了解四種FIR DF的頻響特點(diǎn)。希望以后每一次實(shí)習(xí)都可以學(xué)到更多更有用的知識(shí)。做完這個(gè)數(shù)字圖像處理的實(shí)習(xí)，不僅僅讓我對(duì)本學(xué)期學(xué)過的圖像處理方法的原理了解得更透徹了，而且對(duì)MATLAB的操作也變得更加熟練。在老師的指導(dǎo)和幫助下，我們也都將之逐一地完善了。這些基本的圖像處理方法都可以用我們學(xué)過的MATLAB語(yǔ)句來實(shí)現(xiàn)，只需要設(shè)計(jì)好版面，并處理好軟件的細(xì)節(jié)功能就可以了。第一篇：數(shù)字信號(hào)處理實(shí)習(xí)感想數(shù)字圖像處理個(gè)人體會(huì)這次實(shí)習(xí)讓我學(xué)會(huì)了如何用MATLAB語(yǔ)句來實(shí)現(xiàn)一些功能，雖然用起來還不是很熟練，但收獲不少。這些功能包括數(shù)字圖像處理中所學(xué)到的各種處理圖像的方法，我們將它們劃分為四個(gè)模塊：1）圖像增強(qiáng)與恢復(fù)；2）彩色圖像處理；3）形態(tài)學(xué)圖像處理；4）圖像提取與分割。此外，有些軟件總是會(huì)存在很多的漏洞，像是圖像處理完后下一次打開原圖像部分還會(huì)出現(xiàn)上一次處理的圖像之類的種種問題。而我比較擅長(zhǎng)的就是版面設(shè)計(jì)，做一二界面需要的時(shí)間很短，我們?cè)谧鲞@個(gè)軟件的過程中沒有因?yàn)榉止?zhēng)吵，一切都進(jìn)行得很順利！實(shí)習(xí)初期我們遇到了很多困難，便會(huì)有種不知所措的感覺，不過在接下去的實(shí)習(xí)過程中，困難總是會(huì)被一一克服，原本頭痛的問題最后都會(huì)給我們新的收獲?；蛟S這就是學(xué)習(xí)的精神吧！最后，我們要感謝院里給我們提供了這次實(shí)習(xí)的機(jī)會(huì)，感謝黃忠朝老師和玉罡老師指導(dǎo)我們完成了這次實(shí)習(xí)。=fft(x)對(duì)采集的信號(hào)做快速傅立葉變換，并用[h1,w]=freqz(h)進(jìn)行DTFT變換。xa(t)d(tnT)=229。165。?a(t)Xa(jW)171。2p1165。TTk=165。Xa(jW)=237。Ws/2即連續(xù)信號(hào)是帶限的，且信號(hào)最高頻率不超過抽樣頻率的二分之一，則可不失真恢復(fù)。程序如下：Fs=8000。,time)。)。%錄制語(yǔ)音信號(hào) fprintf(39。按Enter鍵回放錄音39。%按任意鍵播放語(yǔ)音信號(hào)wavwrite(x,Fs,39。%返回采樣點(diǎn)數(shù) df=fs/N。subplot(2,1,1)。)。%顯示縱坐標(biāo)的表示意義 grid。plot(f,abs(y0(n1)))。%加標(biāo)題 xlabel(39。幅值 39。)。x=wavread(39。%畫出截取后的時(shí)域圖形 title(39。頻率39。)。fs=1200。wp=2*pi*fp/Fs。was=2*tan(ws/2)。[B,A]=butter(N,wc,39。figure(1)。巴特沃斯低通濾波器39。)。gridon。Fs=22050。ws=2*pi*fs/Fs。was=2*tan(ws/2)。[B,A]=butter(N,wc,39。)。plot(w,abs(h))。xlabel(39。耗損（dB）39。%帶通濾波fp=[1200,3000]。rs=100。Fs1=1。s39。)。plot(w,abs(h))。xlabel(39。耗損（dB）39。C:39。%返回采樣點(diǎn)數(shù) df=fs/N。subplot(4,2,1)。)。%顯示縱坐標(biāo)的表示意義 grid。%原始信號(hào)的頻譜圖 title(39。頻率w/pi39。)。%加標(biāo)題 grid。ws=2*pi*1200/fs。h=fir1(n,wn/pi,39。ys=fftfilt(h,x)。%在四