matlab課程設(shè)計---matlab語音信號的采集與處理-閱讀頁
2025-06-23 05:22本頁面
【正文】 。 grid。F(jw)39。 subplot(3,1,3) plot(w,P1*180/pi)。 xlabel(39。)。相位 （ 度 ） 39。 信號 f(t)= 21 ()2 tet?? 傅里葉變換的時頻特性如圖所示： ii f(t1)的頻譜。 t=5:r:5。 W=2*pi。 w=k*W/N。 %f(t) F=r*f1*exp(j*t39。 %f(t)的傅里葉變換 F1=abs(F)。 subplot(3,1,1)。 hold on。t39。 ylabel(39。)。f(t1)39。 subplot(3,1,2)。 hold on。w39。 ylabel(39。)。 subplot(3,1,3)。 hold on。 xlabel(39。)。相位 (度 )39。 傅里葉變換的時移特性如圖所示： V 傅里葉變換的頻移特性 信號 f(t)=g2(t)為門函數(shù)，試?yán)L出信號 f1(t)=f(t) 10jte? 以及信號 f2(t)=f(t) 10jte 的 頻譜圖，并與原信號頻譜圖進(jìn)行比較。 t=2:R:2。 f1=f.*exp(j*10*t)。 W1=2*pi*5。 k=N:N。 F1=f1*exp(j*t39。 %f1(t)傅里葉變換 F2=f2*exp(j*t39。 %f2(t)傅里葉變換 F1=real(F1)。 subplot(2,1,1)。hold on。w39。 ylabel(39。)。頻譜 F1(jw)39。 subplot(2,1,2)。 hold on。w39。 ylabel(39。)。頻譜 F2(jw)39。 傅里葉變換的頻移特性如圖所示 ： ② 程序設(shè)計實驗 I 試確定下列信號的傅里葉變換的數(shù)學(xué)表達(dá)式。 t=2:R:2。 f=fourier(Heaviside(t+1)Heaviside(t1)) 運行結(jié)果： f =2/w*sin(w) (b) f(t)= 3 ()tet?? MATLAB 程序 ： R=。 syms t。 MATLAB程序： r= 。 f1=exp(3*t).*stepfun(t,0) 。 subplot(3,1,1) 。grid 。 plot(t,f2) 。 subplot(3,1,3)。grid 。 w2=40000。 As=35。s39。 [z,p,k]=butter(N,Wn,39。)。 w=linspace(1,4000,1000)。 magH=abs(H)。 figure(11)。k39。 xlabel (39。)。幅度 /B39。 濾波器頻率響應(yīng)圖： I 設(shè)計一個 butterworth 高 通濾波器，滿足以下性能指標(biāo)：通帶截止 頻率 w1=40000rad/s，通帶最大衰減 Ap=3dB，阻帶截止頻率 w2=10000rad/s，阻帶最大衰減 As=35dB. w1=40000。 Ap=3。 Fs=10000。s39。 [Z,P,K]=buttap(N)。 [AT,BT,CT,DT]=lp2hp(A,B,C,D,Wc)。 [num2,den2]=bilinear(num1,den1,Fs)。 plot(W*Fs,abs(H))。頻率（ Hz） 39。幅值 39。 濾波器頻率響應(yīng)圖： （ 2） MATLAB 語音信號的采集與處理 MATLAB對原始 語 音 信號進(jìn)行分析，畫出它的時域波形和頻譜 。 [x,fs,Nbits] =wavread(39。) 。 sound(x,fs) 。原始語音信號 39。 xlabel(39。)。fuzhi n39。 figure(2) freqz(x) %繪制原始語音信號的頻率響應(yīng)圖 title(39。) n=length(x)。 %傅里葉變換 y2=fftshift(y1)。 %得出頻點 figure(3) subplot(2,1,1)。原始語音信號 FFT 頻譜 39。 plot(f,abs(y2))。原始語音信號頻譜 39。 %定義 k 值 ,噪聲與原始語音長度一致 c1=10*sin(2*pi*5000*k)。噪聲信號時域圖 39。 %繪制噪聲信號的時域波形圖 figure(5) freqz(c1) %繪制噪聲信號的頻率響應(yīng)圖 title(39。) s=length(c1)。 %對噪聲信號進(jìn)行頻譜變換 xcc=fftshift(xc)。 figure(6) %繪制噪聲語音的 FFT 頻譜圖 plot(f,abs(xcc)) title(39。)。對 加噪后 的 語音 進(jìn)行分析，并畫出其 信號時域和頻譜圖 。 x1=wavread(39。)。 t=0:1/fs:(length(x1)1)/fs。%噪聲幅值 d=[Au*sin(2*pi*5500*t)]39。 sound(x2,22050)。 y2=fft(x2,1024)。 subplot(2,1,1)。 title(39。)。time n39。 ylabel(39。)。 plot(t,x2) title(39。)。time n39。 ylabel(39。)。 plot(f,abs(y1(1:512)))。原始語音信號頻譜 39。 xlabel(39。)。fuzhi39。 subplot(2,1,2)。 title(39。)。Hz39。 ylabel(39。)。另外縱軸上的數(shù)值范圍也發(fā)生了相應(yīng)的變化。 從 看 書的那一刻，我就很郁悶，《 MATLAB 及在電子信息類（課程中的應(yīng)用） 》到底事學(xué)什么的，有什么用。不過，這些感受都是在對這門課程不了解的情況下產(chǎn)生的。 對我們來說這是 一個很陌生的軟件，還是全英文環(huán)境 ， 通 過讀 書，才對 MATLAB 的概念有所了解 。所以， 我不在把 matlab 作為 VC 的附屬品， 認(rèn)真的了解這個軟件。 Matlab的應(yīng)用領(lǐng)域極為廣泛，除數(shù)學(xué)計算和分析外，還被廣泛用于自動控制、系統(tǒng)仿真、數(shù)字信號處理、圖形圖像分析、數(shù)理統(tǒng)計、人工智能、虛擬現(xiàn)實技術(shù)、通信工程、金融系統(tǒng)等領(lǐng)域。 此次課程設(shè)計讓我對我所學(xué)的書本知識進(jìn)行了鞏固，拓展了專業(yè)知識。我用此次機會，把自己把知識在實踐中的到了運用，把理論的抽象具體化，使 我把知識和實踐得到很好的融合。把死板的課本知識變得生動有趣，激發(fā)了學(xué)習(xí)的積極性。常常是再堅持一兩個小時，就能夠解決你的問題。 因而， 在單一的書本學(xué)習(xí)中，有這么一次課程設(shè)計是 放松中學(xué)習(xí)的一個不可多得的機會。它還讓我學(xué)到了很多課本上學(xué)不到的道理，開拓視野的同時增加知識。對于一些辛苦，只要自己堅持，一定可以看到成功的果實。當(dāng)然主要原因還是自己不夠努力，讓自己的缺點一直保持為缺陷。 參考文獻(xiàn)： 向宏，岳恒立，鄭雪峰編著。電子工業(yè)出版社。電子工業(yè)出版社。 ，周濤，基于 MATLAB 的系統(tǒng)分析
畢業(yè)設(shè)計相關(guān)推薦
鄂ICP備17016276號-1