clc;
close?all;
clear?all;

[hcfs]=wavread（‘clean.wav‘）;%調(diào)出語音信號(hào)
s=hc;

%?參數(shù)設(shè)置
M=5;?%?陣元數(shù)目
N=length（s）;?%?采樣快拍數(shù)?18476
c=340;?%?信號(hào)的傳播速度，采用語音信號(hào)的傳播速度
f0=fs;?%?信號(hào)的中心頻率?8000
fj=1000;?%聚焦頻率
lamda=c/f0;?%?信號(hào)的波長(zhǎng)
d=0.02;?%?陣元間距采取信號(hào)波長(zhǎng)的一半
snr_dB=-5;?%?信噪比
snr=10^（snr_dB/10）;?%?線性信噪比
sir_dB=-10;?%?信干比均為-10dB
sir=10^（sir_dB/10）;?%?線性信干比??0.1000
theta_s=0*pi/180;?%?信號(hào)到達(dá)方向
theta_i1=45*pi/180;?%?干擾1到達(dá)方向
angle=[theta_s?theta_i1];
degrad=pi/180;

%?生成信號(hào)源、干擾源及噪聲信號(hào)
power_s=0;
for?t=1:N
????power_s=power_s+（s（t））^2;
end????????????????%?151.7425
power_s=power_s/N;?%?信號(hào)源能量?0.0082
power_i=power_s/sir;?%?干擾源能量??0.0821
power_n=power_s/snr;?%?噪聲信號(hào)能量??0.0260
%?i1=normrnd（0sqrt（power_i/2）1N）+j*normrnd（0sqrt（power_i/2）1N）;?%
%?干擾1的隨機(jī)包絡(luò)
noise=0.08*wgn（MNpower_n）;?%?噪聲信號(hào)??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=s