clc;
close all;
clear all;
[hc,fs]=audioread('clean1.wav');%调出语音信号
s=hc;
% 参数设置
M=5; % 阵元数目
N=length(s); % 采样快拍数 18476
c=340; % 信号的传播速度,采用语音信号的传播速度
f0=fs; % 信号的中心频率 8000
fj=1000; %聚焦频率
lamda=c/f0; % 信号的波长
d=0.02; % 阵元间距采取信号波长的一半
snr_dB=-5; % 信噪比
snr=10^(snr_dB/10); % 线性信噪比
sir_dB=-10; % 信干比均为-10dB
sir=10^(sir_dB/10); % 线性信干比 0.1000
theta_s=0*pi/180; % 信号到达方向
theta_i1=45*pi/180; % 干扰1到达方向
angle=[theta_s theta_i1];
degrad=pi/180;
% 生成信号源、干扰源及噪声信号
power_s=0;
for t=1:N
power_s=power_s+(s(t))^2;
end % 151.7425
power_s=power_s/N; % 信号源能量 0.0082
power_i=power_s/sir; % 干扰源能量 0.0821
power_n=power_s/snr; % 噪声信号能量 0.0260
% i1=normrnd(0,sqrt(power_i/2),1,N)+j*normrnd(0,sqrt(power_i/2),1,N); %
% 干扰1的随机包络
noise=0.08*wgn(M,N,power_n); % 噪声信号 5*18476
% 固定波束形成
tao1=d*sin(theta_s)/c;
tao2=d*sin(theta_i1)/c;
Ts=1.399/N;
L1=ceil(tao1/Ts);
L2=ceil(tao2/Ts);
s1=s';
i1=0.5*s1; %干扰源
x1=s1+i1+noise(1,:);
x2=[zeros(1,L1),s1(1:N-L1)]+[zeros(1,L2),i1(1:N-L2)]+noise(2,:); %各麦克风接收到的信号
x3=[zeros(1,2*L1),s1(1:N-2*L1)]+[zeros(1,2*L2),i1(1:N-2*L2)]+noise(3,:);
x4=[zeros(1,3*L1),s1(1:N-3*L1)]+[zeros(1,3*L2),i1(1:N-3*L2)]+noise(4,:);
x5=[zeros(1,4*L1),s1(1:N-4*L1)]+[zeros(1,4*L2),i1(1:N-4*L2)]+noise(5,:);
X1=[x1;x2;x3;x4;x5];
X2=1/8*(x1+x2+x3+x4+x5);
% 输出信噪比
e=X2-s1;
% ps1=10*log10(sum((s1).^2)/N);
% pnout=10*log10(sum((e).^2)/N);
% snr1=ps1-pnout;
ps1=sum((s1).^2)/N;
pnout=sum((e).^2)/N;
snr1=10*log10(ps1/pnout) %增强后的语音
pnout1=sum(i1).^2/N+power_n;
snr2=10*log10(ps1/pnout1) %未经增强的语音
snr3=snr1-snr2;
subplot(3,1,1);plot(s1);title('(a)原始语音');xlabel('样点数');ylabel('归一化幅度');axis([0 2*10^4 -1 1]);
subplot(3,1,2);plot(s1+i1+noise(1,:));title('(b)带噪语音');xlabel('样点数');ylabel('归一化幅度');axis([0 2*10^4 -1 1]);
subplot(3,1,3);plot(real(X2));title('(c)固定波束法-增强语音');xlabel('样点数');ylabel('归一化幅度');axis([0 2*10^4 -1 1]);
BeamF.rar_固定波束_固定波束形成_语音增强_麦克风阵列_麦克风阵列语音增强
版权申诉
5星 · 超过95%的资源 88 浏览量
2022-07-14
11:50:05
上传
评论 3
收藏 1KB RAR 举报
我虽横行却不霸道
- 粉丝: 75
- 资源: 1万+
评论2