一、简介
基于matlab不同信噪比下的麦克风接收信号
二、源代码
clear all;
close all;
clc
%****************************************************%
[s,fs]=wavread('C9_1_y.wav');
s=s/max(abs(s));
sound_length = 6400;
%装载房间冲激响应
load h.mat;
%麦克风1的信号
s1=conv(s,h1);
s1=s1(1:sound_length);
%麦克风2的信号
s2=conv(s,h2);
s2=s2(1:sound_length);
%****************************************************%
%添加白噪声
%SNR=20dB
s1_20db = awgn(s1,20,'measured','db');
s2_20db = awgn(s2,20,'measured','db');
%SNR=10dB
s1_10db = awgn(s1,10,'measured','db');
s2_10db = awgn(s2,10,'measured','db');
%SNR=0dB
s1_0db = awgn(s1,0,'measured','db');
s2_0db = awgn(s2,0,'measured','db');
%SNR=-5dB
s1_m5db = awgn(s1,-5,'measured','db');
s2_m5db = awgn(s2,-5,'measured','db');
figure(1);
%显示原始输入语音信号波形
subplot(5,2,1),plot(s1,'k'),title('麦克风1（原始）'),xlabel('采样点'),ylabel('幅度');
subplot(5,2,2),plot(s2,'k'),title('麦克风2 （原始）'),xlabel('采样点'),ylabel('幅度');
subplot(5,2,3),plot(s1_20db,'k'),title('麦克风1（SNR=20dB）'),xlabel('采样点'),ylabel('幅度');
subplot(5,2,4),plot(s2_20db,'k'),title('麦克风2（SNR=20dB）'),xlabel('采样点'),ylabel('幅度');
subplot(5,2,5),plot(s1_10db,'k'),title('麦克风1（SNR=10dB）'),xlabel('采样点'),ylabel('幅度');
subplot(5,2,6),plot(s2_10db,'k'),title('麦克风2（SNR=10dB）'),xlabel('采样点'),ylabel('幅度');
subplot(5,2,7),plot(s1_0db,'k'),title('麦克风1（SNR=0dB）'),xlabel('采样点'),ylabel('幅度');
function [h]=rir(fs, mic, n, r, rm, src);
% 房间脉冲响应
% fs 采样频率
% mic 麦克风坐标（行向量）
% n 虚拟声源个数 (2*n+1)^3
% r 墙壁反射系数（-1<R<1）
% rm 房间尺寸（行向量）
% src 声源坐标（行向量）
%
% h 房间脉冲响应
nn=[-n:1:n];
rms=nn+0.5-0.5*(-1).^nn;
srcs=(-1).^(nn);
xi=[srcs*src(1)+rms*rm(1)-mic(1)]; % 式（9-2）
yj=[srcs*src(2)+rms*rm(2)-mic(2)]; % 式（9-3）
zk=[srcs*src(3)+rms*rm(3)-mic(3)]; %式（9-4）
三、运行结果
四、备注
版本：2014a
