matlab的語音信號(hào)頻譜分析報(bào)告(編輯修改稿)
2025-06-09 22:31 本頁面
【文章內(nèi)容簡介】 gure(4)subplot(2,2,1)。plot(th1,h1)title(39。矩形窗時(shí)域39。)。subplot(2,2,2)。plot(w,abs(fh1))。title(39。矩形窗頻域39。)。subplot(2,2,3)。plot(th2,h2)。title(39。布萊克曼窗時(shí)域39。)。subplot(2,2,4)。plot(w,abs(fh2))。title(39。布萊克曼窗頻域39。)%解調(diào)后信號(hào)n3=length(y3)t2=0:n31。y4=cos(pi*1/2*t2)y5=y3.*y439。%調(diào)制后的信號(hào)函數(shù)figure %濾波處理y6=conv(h1,y5)。%用矩形窗對調(diào)制后信號(hào)進(jìn)行濾波f6=fft(y6)。n4=length(f6)。ty6=(0:n41)/fs。w3=2/n4*[0:n41]。%sound(y6,fs)。figure(5)subplot(2,2,1)。plot(ty6,y6)。title(39。矩形窗濾波后音頻39。)subplot(2,2,2)。plot(w3,abs(f6))。title(39。矩形窗濾波后頻譜39。)y7=conv(h2,y5)。%用布萊克曼窗進(jìn)行濾波f7=fft(y7)。n5=length(f7)。ty7=(0:n51)/fs。w4=2/n5*[0:n51]。%sound(y7,fs)。subplot(2,2,3)。plot(ty7,y7)。title(39。布萊克曼窗濾波后音頻39。)subplot(2,2,4)。plot(w4,abs(f7))。title(39。布萊克曼窗濾波后頻譜39。)實(shí)驗(yàn)五clc。clear。close。[y,fs,bit]=wavread(39。I do39。)。y0=y(:,1)。l=length(y0)。%加三余弦混合噪聲t0=(0:l1)/fs。d0=[*cos(2*pi*3000*t0)]39。t1=(0:l1)/fs。d1=[*cos(2*pi*5000*t1)]39。t2=(0:l1)/fs。d2=[*cos(2*pi*8000*t2)]39。noise=d2+d1+d0。y1=y0+noise。%sound(y1,fs)。a=length(noise)。%繪制三余弦噪聲音頻圖wa=2/a*[0:a1]。Noise=fft(noise)。figure(1)subplot(2,3,4)。plot(noise(1:150))。 xlabel(39。時(shí)間（s）39。)ylabel(39。幅值39。)title(39。三余弦信號(hào)音譜39。) subplot(2,3,1)。%繪制三余弦噪聲頻譜圖plot(wa,abs(Noise))。grid on。 xlabel(39。W39。)title(39。噪聲頻譜39。)w0=2/l*[0:l1]。%繪制加噪信號(hào)音頻Y1=fft(y1)。subplot(2,3,5)plot(w0,abs(Y1))。grid on。 xlabel(39。W39。)title(39。加噪信號(hào)頻譜39。)ly1=length(y1)。ty1=(0:ly11)/fs。subplot(2,3,2)。plot(ty1,y1)。xlabel(39。時(shí)間（s）39。)ylabel(39。幅值39。)title(39。加噪信號(hào)音譜39。) m=rand(l,1)。 %lm=length(m)。y2=m+y0。%將噪聲信號(hào)與原聲音信號(hào)疊加wm=2/lm*[0:lm1]。M=fft(m)。figure(2)subplot(2,2,3)。plot(m(1:150))xlabel(39。時(shí)間（s）39。)ylabel(39。幅值39。)title(39。白噪信號(hào)音譜39。)subplot(2,2,1)。plot(wm,abs(M))。grid on。 xlabel(39。W39。)title(39。噪聲頻譜39。)l=length(y2)。ty2=(0:l1)/fs。w=2/l*[0:l1]。Y2=fft(y2)。subplot(2,2,4)plot(w,abs(Y2))。grid on。 xlabel(3
環(huán)評公示相關(guān)推薦
鄂ICP備17016276號(hào)-1