function?[M]?=?correlation（XNFFTFSwinOVERLAPN）
%
%	function?[M]?=?correlation（XFSNFFTOVERLAPN）
%	making?correlation?matrices?for?each?frequency
%	X?must?has?2?column
%	
%	Shiro?Ikeda?14July1998?

SPLEN?=?1000;

[LENGTHdim]?=?size（X）;

FL?=?floor（（LENGTH-NFFT）/（NFFT-OVERLAP））+1;
DL?=?（SPLEN-1）*（NFFT-OVERLAP）+NFFT;
NF?=?NFFT/2+1;

repeat=floor（FL/SPLEN）;

MARGIN=zeros（2*NFN）;
te=OVERLAP;?M=zeros（2*NF2*N+2）;?AVE=zeros（2*NF1）;
for?ITER=1:repeat
??ts?=?te-OVERLAP+1;?te?=?ts-1+DL;
??B1=specgram（X（ts:te1）NFFTFSwinOVERLAP）;
??B2=specgram（X（ts:te2）NFFTFSwinOVERLAP）;
??if?ITER==1
????STKH1=B1（:1:N）;?STKH2=B2（:1:N）;
??end
??for?f=1:NF
????f1=2*f-1;?f2=2*f;?
????AVE（f1）?=?AVE（f1）+B1（f:）*ones（SPLEN1）;
????AVE（f2）?=?AVE（f2）+B2（f:）*ones（SPLEN1）;
????TMP=[MARGIN（f1:）B1（f:）;MARGIN（f2:）B2（f:）];
????MARGIN（f1:f2:）=[B1（fSPLEN-N+1:SPLEN）;B2（fSPLEN-N+1:SPLEN）];
????for?tau=0:N
??????TMP1=TMP（:N-tau+1:SPLEN+N-tau）*TMP（:N+1:SPLEN+N）‘;
??????M（f1:f22*tau+1:2*tau+2）=M（f1:f22*tau+1:2*tau+2）+TMP1;
????end
??end
end

ts?=?te-OVERLAP+1;?te?=?LENGTH;
B1=specgram（X（ts:te1）NFFTFSwinOVERLAP）;
B2=specgram（X（ts:te2）NFFTFSwinOVERLAP）;
RL=size（B12）;
if?repeat==0
??STKH1=B1（:1:N）;?STKH2=B2（:1:N）;
end
STKT1=B1（:RL-N+1:RL）;?STKT2=B2（:RL-N+1:RL）;
for?f=1:NF
??f1=2*f-1;?f2=2*f;?
??AVE（f1）?=?AVE（f1）+B1（f:）*ones（RL1）;
??AVE（f2）?=?AVE（f2）+B2（f:）*ones（RL1）;
??TMP=[MARGIN（f1:）B1（f:）;MARGIN（f2:）B2（f:）];
??for?tau=0:N
????TMP1=TMP（:N-tau+1:RL+N-tau）*TMP（:N+1:RL+N）‘;
????M（f1:f22*tau+1:2*tau+2）=M（f1:f22*tau+1:2*tau+2）+TMP1;
??end
end

AVE=AVE/FL;

for?f=1:NF
??f1=2*f-1;?f2=2*f;?
??TMPA?=?AVE（f1:f2）*AVE（f1:f2）‘;
??M（f1:f21:2）=M（f1:f21:2）-TMPA*FL;
??M（f1:f21:2）=M（f1:f21:2）+M（f1:f21:2）‘;
??for?tau=1:N
????t1=2*tau+1;?t2=2*tau+2;
????T1?=?AVE（f1:f2）*FL-[STKT1（fN-tau+1:N）;STKT2（fN-tau+1:N）]*ones（tau1）;
????T2?=?AVE（f1:f2）*FL-[STKH1（f1:tau）;STKH2（f1:tau）]*ones（tau1）;
????TMP1=T1*AVE（f1:f2）‘;?TMP2=AVE（f1:f2）*T2‘;
????M（f1:f2t1:t2）=M（f1:f2t1:t2）-TMP1-TMP2+TMPA*（FL-tau）;
????M（f1:f2t1:t2）=M（f1:f2t1:t2）+M（f1:f2t1:t2）‘;
??end
end

for?tau=0:N
??M（:2*tau+1:2*tau+2）=M（:2*tau+1:2*tau+2）/（2*（FL-tau））;
end

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

?????文件????????333??1998-07-23?16:02??ica_matlab\ica\decorrelation.m

?????文件???????2179??1998-07-23?16:02??ica_matlab\ica\correlation.m

?????文件???????1151??1998-07-23?16:05??ica_matlab\ica\sepfilter.m

?????文件???????1213??1998-07-24?12:08??ica_matlab\ica\ica_f.m

?????文件???????1364??1998-07-23?16:02??ica_matlab\ica\permutation.m

?????文件??????66840??2003-06-23?21:42??ica_matlab\ica\X_room.wav

?????文件????????960??2003-06-23?21:45??ica_matlab\ica\Readme

?????文件?????123244??2003-06-23?21:42??ica_matlab\ica\X_linear.wav

?????目錄??????????0??2011-10-23?12:57??ica_matlab\ica

?????目錄??????????0??2011-10-23?12:57??ica_matlab

-----------?---------??----------?-----??----

???????????????197284????????????????????10