%?clear?all
clear?all;
close?all;
clc;

%?extracted?sources‘?signals
load?alz.mat;
load?fw.mat;

m=min（length（alz）length（fw））;
N1=1;%20000;
N2=m;%N1+3000;
s1=alz（N1:N2）;??%?Adjust?the?power?of?alz
s2=fw（N1:N2）;
s1=（s1-min（s1））/（max（s1）-min（s1））;
s2=（s2-min（s2））/（max（s2）-min（s2））;
fs=16000;
N=8;
auwrite（s1fsN‘s1‘）;
auwrite（s2fsN‘s2‘）;
s1=1.3*s1;
Enegy1=0;
Enegy2=0;
for?k=1:length（s1）
????Enegy1=Enegy1+s1（k）^2;
????Enegy2=Enegy2+s2（k）^2;
end
Ratio=Enegy1/Enegy2;

S1=fft（s1）;
S2=fft（s2）;

sita1=-pi/3;????%?the?angle?of?the?first?source.??range:?from?-pi/2?to?pi/2
sita2=pi/3;?????%?the?angle?of?the?second?source.?range:?from?-pi/2?to?pi/2
r1=1.0;?????????%?the?radius?of?the?first?sound?source
r2=1.0;?????????%?the?radius?of?the?second?sound?source
Nfft=length（S1）;
D=0.433;????????%?distance?between?two?microphones?suppose?the?diameter?of?circle?is?0.5.
V=344.21;???????%23?degrees????%?velocity?of?sounds
fs=1.6e4;???????%?the?sampling?frequency
Nbeg=0.5;
Nend=1.5;
No=11;

%******************************************
%?Variants?for?Reverberation
%******************************************
T60=0.4;????????????%??time?for?degrading?up?to?60db;?seen?[1]
L=0.5*T60*fs;???????%??length?of?filter?for?reverberation

%?Mixed?sound?signals?in?near--field?equation
d11=sqrt（D*D/4+r1*r1+D*r1*sin（sita1））;
d12=sqrt（3*D*D/4+r1*r1-sqrt（3）*D*r1*cos（sita1））;
d13=sqrt（D*D/4+r1*r1-D*r1*sin（sita1））;
d21=sqrt（D*D/4+r2*r2+D*r2*sin（sita2））;
d22=sqrt（3*D*D/4+r2*r2-sqrt（3）*D*r2*cos（sita2））;
d23=sqrt（D*D/4+r2*r2-D*r2*sin（sita2））;
for?f=1:Nfft
??????X1（f）=exp（-i*2*pi*f/Nfft*d11*fs/V）*d11/d11*S1（f）+exp（-i*2*pi*f/Nfft*d21*fs/V）*d21/d21*S2（f）;
??????X2（f）=exp（-i*2*pi*f/Nfft*d12*fs/V）*d11/d12*S1（f）+exp（-i*2*pi*f/Nfft*d22*fs/V）*d21/d22*S2（f）;
??????X3（f）=exp（-i*2*pi*f/Nfft*d13*fs/V）*d11/d13*S1（f）+exp（-i*2*pi*f/Nfft*d23*fs/V）*d21/d23*S2（f）;
end
mid_mix1=real（ifft（X1））;
mid_mix2=real（ifft（X2））;
mid_mix3=real（ifft（X3））;
%*********Add?reverberation***************

%*********Impulse?Response[2]*************

%?%?addtional?noise:?40dB
%?save（‘X1‘‘X1‘）;
%?save（‘X2‘‘X2‘）;
%?save（‘X3‘‘X3‘）;

%?load?X1;
X1=awgn（X120）;
%?load?X2;
X2=awgn（X220）;
%?load?X3;
X3=awgn（X320）;

mixed1=real（ifft（X1））;
mixed2=real（ifft（X2））;
mixed3=real（ifft（X3））;

mixed1=（mixed1-min（mixed1））/（max（mixed1）-min（mixed1））;
mixed2=（mixed2-min（mixed2））/（max（mixed2）-min（mixed2））;
mixed3=（mixed3-min（mixed3））/（max（mixed3）-min（mixed3））;

figure;
subplot（311）;
plot（mixed1‘-r‘）;
xlabel（‘signal?points‘）
ylabel（‘Amplitude‘）
title（‘Mixed?signal1‘）
subplot（312）;
plot（mixed2‘-r‘）;
xlabel（‘signal?points‘）
ylabel（‘Amplitude‘）
title（‘Mixed?signal2‘）
subplot（313）;
plot（mixed3‘-r‘）;
xlabel（‘signal?points‘）
ylabel（‘Amplitude‘）
title（‘Mixed?signal3‘）

%?%DSBF
%?B1=0;
%?L=zeros（19No）;
%?for?belta=-pi/2:pi/18:pi/2
%?????B1=B1+1;
%?????R1=0;
%?????z（B1）=0;
%?????for?focus=Nbeg:0.1:Nend
%?????????R1=R1+1;
%?????????Z（B1

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????7380??2007-02-07?17:31??超聲算法\1\DSBF.m

?????文件???????2287??2015-10-06?08:36??超聲算法\10\createSignal.m

?????文件???????7017??2018-12-28?15:50??超聲算法\10\examples.m

?????文件???????1515??2015-10-06?08:36??超聲算法\10\license.txt

?????文件???????4176??2015-10-06?08:36??超聲算法\10\plotSteeredResponse.m

?????文件???????1668??2015-10-06?08:36??超聲算法\10\sphericalCoordinateSystem.m

?????文件???????2183??2015-10-06?08:36??超聲算法\10\steeredResponseDelayAndSum.m

?????文件???????1656??2015-10-06?08:36??超聲算法\10\steeringVector.m

?????文件?????472907??2015-10-06?08:36??超聲算法\10\TN?beamformers.pdf

?????文件???????1112??2006-07-19?09:39??超聲算法\11\線陣DBF\linearityarray.m

?????文件???????1056??2006-07-19?09:40??超聲算法\11\線陣DBF\linearityarray2.m

?????文件???????1219??2018-12-28?15:45??超聲算法\12\beamforming.m

?????文件????????993??2012-02-03?20:18??超聲算法\13\diffusenoisefield.m

?????文件???????1938??2008-11-24?23:10??超聲算法\13\DSB.m

?????文件????????311??2015-05-12?00:48??超聲算法\13\DSB_BeamPattern.m

?????文件????????664??2010-03-27?16:00??超聲算法\14\example21.m

?????文件????????790??2013-12-24?20:19??超聲算法\15\circular_array.m

?????文件???????1060??2013-03-05?13:54??超聲算法\16\beamform\circle_beampattern.m

?????文件???????1125??2012-10-22?11:20??超聲算法\16\beamform\cylinder_beampattern.m

?????文件??????10408??2013-03-03?15:27??超聲算法\16\beamform\Fr_bin2dec.m

?????文件??????22789??2013-03-03?15:24??超聲算法\16\beamform\fr_dec2bin.m

?????文件???????1053??2013-03-05?16:27??超聲算法\16\beamform\ULAbeamformzuoye.m

?????文件???????1324??2012-10-22?11:20??超聲算法\16\beamform\URAbeampattern.m

?????文件???????1659??2018-12-28?18:14??超聲算法\17\Beamforming.m

?????文件???????9079??2013-12-25?08:57??超聲算法\18\概率數據互聯（pdaf）\PDAF-demo\pdaf.m

?????文件??????26806??2013-12-22?17:39??超聲算法\18\概率數據互聯（pdaf）\UF_EF\blackscholes.m

?????文件????????386??2013-12-22?17:39??超聲算法\18\概率數據互聯（pdaf）\UF_EF\bsffun.m

?????文件????????790??2013-12-22?17:39??超聲算法\18\概率數據互聯（pdaf）\UF_EF\bshfun.m

?????文件???????5600??2013-12-22?17:39??超聲算法\18\概率數據互聯（pdaf）\UF_EF\data\c2925.prn

?????文件???????3560??2013-12-22?17:39??超聲算法\18\概率數據互聯（pdaf）\UF_EF\data\c3025.prn

............此處省略411個文件信息