matlab的語音信號頻譜分析報告-在線瀏覽
2025-06-30 22:31本頁面
【正文】 el(39。)。調(diào)制后信號的頻譜(低頻率調(diào)制)39。grid on。 實驗四clear all。close all[a,fs,bit]=wavread(39。)。%去單列數(shù)據(jù)進(jìn)行分析f1=fft(y1)。tn=(0:n1)/fs。%sound(y1,fs)。plot(tn,y1)。title(39。)。時間39。ylabel(39。)。plot(w,abs(f1))。title(39。)。頻率/pi39。ylabel(39。)。y2=cos(pi*1/2*t)。%信號調(diào)制ty3=(0:length(y3)1)/fs。plot(ty3,y3)。title(39。)。時間39。ylabel(39。)。n2=length(f3)。subplot(2,3,5)。%繪制調(diào)制后信號頻譜圖grid on。AM調(diào)制頻譜39。xlabel(39。)。幅度39。 %解調(diào)后信號n3=length(y3)t2=0:n31。%解調(diào)后的信號函數(shù)subplot(2,3,3)。grid on。時間39。ylabel(39。)。AM解調(diào)音頻信號39。f5=fft(y5)。subplot(2,3,6)。%繪制解調(diào)后信號頻譜圖grid on。AM解調(diào)頻譜39。xlabel(39。)。幅度39。 %設(shè)計巴特沃斯濾波器進(jìn)行濾波去噪[N1,wc1]=buttord(,1,15)。%確定低通濾波器分子分母系數(shù)[H,W]=freqz(b,a)。%低通濾波器波形xlabel(39。)。H39。title(39。)。wavplay(m,fs)。%濾波后信號長度tm=(0:lm1)/fs。M=fft(m)。plot(tm,m)。title(39。)。t39。ylabel(39。)。plot(wm,abs(M))。濾波后波形頻譜圖39。xlabel(39。)。幅度39。%矩形窗和布萊克曼窗N=33。%基于經(jīng)驗的指標(biāo)，其中N為理想低通濾波器階數(shù)，wc為截止頻率hd=ideal(N,wc)。%產(chǎn)生各種窗函數(shù)w2=blackman(N)。%加窗設(shè)計各種FIR濾波器h2=hd.*w239。th2=(0:32)/fs。fh1=fft(h1,M)。fh2=fft(h2,M)。plot(th1,h1)title(39。)。plot(w,abs(fh1))。矩形窗頻域39。subplot(2,2,3)。title(39。)。plot(w,abs(fh2))。布萊克曼窗頻域39。y4=cos(pi*1/2*t2)y5=y3.*y439。%用矩形窗對調(diào)制后信號進(jìn)行濾波f6=fft(y6)。ty6=(0:n41)/fs。%sound(y6,fs)。plot(ty6,y6)。矩形窗濾波后音頻39。plot(w3,abs(f6))。矩形窗濾波后頻譜39。%用布萊克曼窗進(jìn)行濾波f7=fft(y7)。ty7=(0:n51)/fs。%sound(y7,fs)。plot(ty7,y7)。布萊克曼窗濾波后音頻39。plot(w4,abs(f7))。布萊克曼窗濾波后頻譜39。clear。[y,fs,bit]=wavread(39。)。l=length(y0)。d0=[*cos(2*pi*3000*t0)]39。d1=[*cos(2*pi*5000*t1)]39。d2=[*cos(2*pi*8000*t2)]39。y1=y0+noise。a=length(noise)。Noise=fft(noise)。plot(noise(1:150))。時間（s）39。幅值39。三余弦信號音譜39。%繪制三余弦噪聲頻譜圖plot(wa,abs(Noise))。 xlabel(39。)title(39。)w0=2/l*[0:l1]。subplot(2,3,5)plot(w0,abs(Y1))。 xlabel(39。)title(39。)ly1=length(y1)。subplot(2,3,2)。xlabel(39。)ylabel(39。)title(39。) m=rand(l,1)。y2=m+y0。M=fft(m)。plot(m(1:150))xlabel(39。)ylabel(39。)title(39。)subplot(2,2,1)。grid on。W39。噪聲頻譜39。ty
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號-1