GSC 采用广义旁瓣抵消（GSC）自适应波束形成方法实现时域和频域滤波

采用广义旁瓣抵消（GSC）自适应波束形成方法实现时域和频域滤波，采用LMS自适应算法，最终实现语音增强。（文件中包含纯净语音及不同信噪比的带噪语音）

% 
clear;
clc

M=4;   % 麦克风数目

clean = wavread（‘clean.wav‘）;
x1 = wavread（‘noisy0.wav‘）;
x2 = wavread（‘noisy10.wav‘）;
x3 = wavread（‘noisy15.wav‘）;
x4 = wavread（‘noisy20.wav‘）;
% x5 = wavread（‘noisy25.wav‘）;
Len = length（clean）;
mic1 = x1（1:Len）;
mic2 = x2（1:Len）;
mic3 = x3（1:Len）;
mic4 = x4（1:Len）;
% mic5 = x5（1:Len）;
x = [mic1‘;mic2‘;mic3‘;mic4‘];%;mic5‘
FBFout=sum（x）/M;

B = [1 -1 0 0 
     0 1 -1 0 
     0 0 1 1 ];
Bout = B*x;
x1 = Bout（1:）‘;
x2 = Bout（2:）‘;
x3 = Bout（3:）‘;
% x4 = Bout（4:）‘;

y1 = zeros（1Len）;
y2 = y1;
y3 = y1;
y4 = y1;
MCout = y1;
N = 64;
h1 = zeros（1N）;
h2 = h1;
h3 = h1; 
h4 = h1;
GSCout = zeros（1Len）;
u = 0.0002;
weight = [];
for i=1:fix（Len/N）-1
    X1 = fft（x1（（i-1）*N+1:（i+1）*N））;
    H1 = fft（[h1zeros（1N）]）;
 

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件       36996  2010-04-06 14:39  clean.wav
     文件        1919  2010-06-27 10:52  FGSC.m
     文件       36996  2010-04-06 14:44  noisy0.wav
     文件       36996  2010-04-06 14:44  noisy10.wav
     文件       36996  2010-04-06 14:41  noisy15.wav
     文件       36996  2010-04-06 14:40  noisy20.wav
     文件       36996  2010-06-23 09:17  noisy5.wav
     文件         755  2010-06-11 08:43  shiyuchengxu.m