信号频谱分析_第1页
信号频谱分析_第2页
信号频谱分析_第3页
已阅读5页，还剩3页未读， 继续免费阅读

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

声音信号的频谱分析姓名：苏小平学号：1305080116专业：智能电网信息工程第一部分：goldwave30图一第二部分：关于声音的Matlab分析Matlab30以利于观察。程序如下：[y,Fs,bits]=wavread('至上励合-繁星.wav');y=y(:,1);sigLength=length(y);Y=fft(y,sigLength);Pyy=Y.*conj(Y)/sigLength;halflength=floor(sigLength/2);f=Fs*(0:halflength)/sigLength;figure;plot(f,Pyy(1:halflength+1));xlabel('Frequency(Hz)');t=(0:sigLength-1)/Fs;figure;plot(t,y);xlabel('Time(s)');图二图二为处理前的声道频域波形程序如下：[y,fs,bits]=wavread('至上励合繁星.wav');L=length(y);noise=0.1*randn(L,2);y_z=y+noise;wavwrite(y_z,8000*6,'加噪声音乐');n=length(y);y_zp=fft(y_z,n);f=fs*(0:n/2-1)/n;figure(1);subplot;plot(y_z);title('加噪语音信号时域波形');xlabel('time')ylabel('幅值A')（2）利用低通滤波器将噪声滤除程序如下：fp=1500fc=1700As=100Ap=1wc=2*pi*fc/fswp=2*pi*fp/fswdel=wcwpbeta=0.112*(As8.7)N=ceil((As8)/2.285/wdel)wn=kaiser(N+1,beta)ws=(wp+wc)/2/pi;b=fir1(N,ws,wn);freqz(b,1);[y,fs,bits]=wavread('至上励合-繁星.wav');L=length(y);noise=0.1*randn(L,2);y_z=y+noise;n=length(y);x=fftfilt(b,y_z);X=fft(x,n);f=fs*(0:n/2-1)/n;figure;gridon;plot(abs(X));title('滤波后信号的频谱');gridonwavwrite(x,fs,bits,'滤波后音乐');滤波器的幅值与相位图：通过低通滤波器后有用的声音的频谱：第三部分：LinkinMatlab的程序作出时域波形的频谱图。40040035030025020015010050000.511.52Frequency(Hz)2.5x1041.81.8.6.4.200.20.40.60.8-1051015Time(s)202530000----加入随机噪声以后的频谱图加噪语音信号时域波形1.510.5A值 0幅-0.5-1

人人文库> 全部分类> 教育资料 > 辅导培训

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论