音频信号去噪报告.docxVIP

实验目的：1、熟悉MATLAB语言的基本用法；2、掌握MATLAB语言中音频数据与信息的读取与播放方法；3、掌握在MATLAB中设计滤波器的方法；4、掌握噪声产生的方法；5、掌握MATLAB语言中信号频谱的绘制方法。实验原理：选取一段.wav格式的音乐或歌曲，用MATLAB将选取的音频文件读出来，在读出的音频信号中加入随机高斯噪声，再将带有噪声的音乐信号通过低通滤波器滤除噪声，还原音乐。实验步骤：仿真程序audio_signal.m:clear all%% 读取原始音频文件[x,Fs]=audioread(C:\Users\yhtcp\Desktop\音频信号去噪\music.wav);%默认采样频率Fs为44100Hzx1=x(:,1);%1声道数据x2=x(:,2);%2声道数据t=0:1/Fs:(length(x1)-1)/Fs;%取时域横轴tN = 1024; %fft点数f=10^3;%取频率轴单位为KHzX0=fft(x1,N);X1=10*log(abs(fftshift(X0)));%用对数表示原始音频信号的频谱figure(1);%原始音频信号时域图、频谱图subplot(2,1,1);plot(t,x1,r);axis([0 95 -1 1]);%设置坐标轴范围set(gca,XTick,0:5:95),grid on;%设置x坐标轴刻度，绘制方格线set(gca,YTick,-0.8:0.2:0.8),grid on;%设置y坐标轴刻度，绘制方格线title(原始声音信号波形);xlabel(时间轴，单位s);subplot(2,1,2);plot((-N/2+1:N/2)*Fs/N/f,X1,b);axis([-25 25 -60 40]);set(gca,XTick,-25:5:25),grid on;set(gca,YTick,-60:20:40),grid on;title(原始声音信号频谱图);xlabel(频率轴，单位kHz);ylabel(单位:dB);%% 原始音频加噪处理SNR=16;y1=awgn(x1,SNR,measured);%对音频信号叠加高斯白噪声，信噪比为16Ym=max(max(max(y1)),max(abs(min(y1))));%找出极值Y=y1/Ym;%归一化处理audiowrite(C:\Users\yhtcp\Desktop\音频信号去噪\music_noise.wav,Y,Fs);%生成原始信号叠加噪声后的wav格式音频文件，试听叠加噪声效果Y0=fft(y1,N);Y1=10*log(abs(fftshift(Y0)));%用对数表示噪声信号的频谱figure(2);%原始音频信号叠加噪声后时域图、频谱图subplot(2,1,1);plot(t,y1,r);axis([0 95 -1 1]);set(gca,XTick,0:5:95),grid on;set(gca,YTick,-0.8:0.2:0.8),grid on;title(叠加噪声后信号的波形);xlabel(时间轴，单位s);subplot(2,1,2);plot((-N/2+1:N/2)*Fs/N/f,Y1,b);axis([-25 25 -60 40]);set(gca,XTick,-25:5:25),grid on;set(gca,YTick,-60:20:40),grid on;title(叠加噪声后信号的频谱);xlabel(频率轴，单位kHz);ylabel(单位:dB);%% 设计低通滤波器wn=0.15;%滤波器归一化截止频率为0.15n=128;%滤波器阶数为128hh=fir1(n,wn,low);figure(3);freqz(hh);%绘制滤波器的频率响应图%% 对加噪音频信号进行滤波处理z1=filter(hh,1,y1);%使含有噪声的信号通过一个已设计的低通滤波器hhZ0=fft(z1,N);Z1=10*log(abs(fftshift(Z0)));、z=[z1 x2];、Zm=max(max(max(z)),max(abs(min(z))));、Z=z/Zm;figure(4);subplot(2,1,1);plot(t,z1,r);axis([0 95 -1 1]);set(gca,XTick,0:5:95),grid on;set(gca,YTick,-0.8:0.2:0.8),grid on;title(滤波后信号的波形);xlabel(时间轴，单位s);subplot(2,1,2);plot((-N/2+1:N/2)*Fs/N/f,Z1,b);axis([-25 25 -60 40]);set(gca,XTick,-25:5:25),grid on;s