clc;clear?all;close?all;
%?讀入混合前的原始圖片并顯示
[S1Fs1]=wavread?（‘music.wav‘）;?
S1?=?S1‘;
[S2Fs2]=wavread?（‘man.wav‘）;
S2?=?S2‘;
subplot（321）plot（S1）title（‘輸入聲音1‘）
subplot（322）plot（S2）title（‘輸入聲音2‘）
%?將其組成矩陣
S=[S1;S2];?????
%?weight=rand（size（S1））;???????????????%?取一隨機矩陣，作為信號混合的矩陣
%?MixedS=weight*S;??????????????????????%?模擬麥克風獲取信號
%?讀取麥克風獲取的信號
M1?=?wavread?（‘MixedS1.wav‘）‘;
M2?=?wavread?（‘MixedS2.wav‘）‘;
MixedS?=?[M1;M2];
%?將混合聲音顯示
subplot（323）plot（MixedS（1:））title（‘混合聲音1‘）
subplot（324）plot（MixedS（2:））title（‘混合聲音2‘）
%?wavwrite（MixedS（1:）Fs1‘MixedS1.wav‘）;
%?wavwrite（MixedS（2:）Fs1‘MixedS2.wav‘）;
MixedS_bak=MixedS;?????????????????????????%?將混合后的數據備份，以便在恢復時直接調用
%%%%%%%%%%%%%%%%%%%%%%%%%%??中心化??%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
MixedS_mean=zeros（21）;
MixedS_mean=mean（MixedS2）;
MixedS?=?MixedS-repmat（MixedS_mean1size（MixedS2））;
%%%%%%%%%%%%%%%%%%%%%%%%%%%??白化??%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

MixedS_cov=cov（MixedS‘）;????????????????????%?cov為求協方差的函數
[ED]=eig（MixedS_cov）;??????????????????????%?對圖片矩陣的協方差函數進行特征值分解
Q=inv（sqrt（D））*（E）‘;????????????????????????%?Q為白化矩陣
MixedS_white=Q*MixedS;??????????????????????%?MixedS_white為白化后的圖片矩陣
IsI=cov（MixedS_white‘）;?????????????????????%?IsI應為單位陣????????????

%%%%%%%%%%%%%%%%%%%%%%%%　FASTICA算法??%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
X=MixedS_white;????????????????????????????%?以下算法將對X進行操作
[VariableNumSampleNum]=size（X）;
numofIC=VariableNum;???????????????????????%?在此應用中，獨立元個數等于變量個數
B=zeros（numofICVariableNum）;??????????????%?初始化列向量w的寄存矩陣B=[b1??b2??...???bd]
for?r=1:numofIC
????i=1;maxIterationsNum=100;???????????????%?設置最大迭代次數（即對于每個獨立分量而言迭代均不超過此次數）
????b=rand（numofIC1）-.5;??????????????????%?隨機設置b初值
????b=b/norm（b）;???????????????????????????%?對b標準化?norm（b）:向量元素平方和開根號
????while?i<=maxIterationsNum+1
????????bOld=b;??????????????????????????
????????t=X‘*b;
????????g=t.*exp（-t.^2/2）;
????????dg=（1-t.^2）.*exp（-t.^2/2）;
????????b=X*g/SampleNum-mean（dg）*b;
????????b=b-B*B‘*b;????????????????????????%?對b正交化
????????b=b/norm（b）;?
????????if?abs（abs（b‘*bOld）-1）<1e-9????????%?如果收斂，則
?????????????B（:r）=b;?????????????????????%?保存所得向量b
?????????????break;
?????????end
????????i=i+1;????????
????end
end
if?i?==?maxIterationsNum+1??????????%?循環結束處理
??????fprintf（‘\n第%d分量在%d次迭代內并不收斂?！?rmaxIterationsNum）;
??????break;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%??ICA計算的數據復原并構圖??%%%%%%%%%%%%%%%%%%%%%%%%%
ICAedS=B‘*Q*MixedS_bak;?????????????????????%?計算ICA后的矩陣

%?將混合矩陣重新排列并輸出
subplot（325）plot（ICAedS（1:））title（‘ICA解混聲音1‘）
subplot（326）plot（ICAedS（2:））title（‘ICA解混聲音2‘）
wavwrite（ICAedS（1:）Fs1‘SS1.wav‘）;
wavwrite（ICAedS（2:）Fs1‘SS2.wav‘）;

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

?????文件???????3027??2011-05-29?00:28??BBS\ICADemo.m

?????文件?????882044??2006-05-28?15:31??BBS\man.wav

?????文件?????882044??2011-05-29?00:20??BBS\MixedS1.wav

?????文件?????882044??2011-05-29?00:20??BBS\MixedS2.wav

?????文件?????882142??2006-05-28?15:32??BBS\music.wav

?????文件?????882044??2011-05-29?00:32??BBS\SS1.wav

?????文件?????882044??2011-05-29?00:32??BBS\SS2.wav

?????目錄??????????0??2011-05-29?00:32??BBS

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

??????????????5295389????????????????????8