【正文】 lot(6,3,1)。plot(I1)。title(39。輸入信號(hào) 139。)。 set(gca,39。ylim39。,[4,4])。 subplot(6,3,4)。a=I1。 y1=fft(a, 200)。 %對(duì)信號(hào)做 200 點(diǎn) FFT 變換 y=abs(y1)。%取模，整合幅頻和相頻 plot(y,39。r39。) title(39。輸入信號(hào) 1 頻譜 39。)。 axis([0 300 0 200])。 subplot(6,3,2)。plot(I2)。title(39。輸入信號(hào) 239。)。 set(gca,39。ylim39。,[4,4])。 subplot(6,3,5)。b=I2。 y1=fft(b, 200)。 %對(duì)信號(hào)做 200 點(diǎn) FFT 變換 y=abs(y1)。%取模，整合幅頻和相頻 plot(y,39。r39。) title(39。輸入信號(hào) 2 頻譜 39。)。 axis([0 300 0 200])。 subplot(6,3,3)。plot(I3)。title(39。輸入信號(hào) 339。)。 set(gca,39。ylim39。,[4,4])。 subplot(6,3,6)。c=I3。 y1=fft(c, 200)。 %對(duì)信號(hào)做 200 點(diǎn) FFT 變換 y=abs(y1)。%取模，整合幅頻和相頻 武漢理工大學(xué)《信息處理課群綜合訓(xùn)練與設(shè)計(jì)》報(bào)告 26 plot(y,39。r39。) title(39。輸入信 號(hào) 3 頻譜 39。)。 axis([0 300 0 200]) % 將其組成矩陣 S=[I1。I2。I3]。 % 因此 S_all 是一個(gè)變量個(gè)數(shù)＊采樣個(gè)數(shù)的矩陣 Sweight=rand(size(S,1)) % 取一隨機(jī)矩陣，作為信號(hào)混合的權(quán)矩陣 MixedS=Sweight*S。 % 得到三個(gè)信號(hào)的混合信號(hào)矩陣 % 將混合矩陣重新排列并輸出 subplot(6,3,7),plot(MixedS(1,:)),title(39。混合信號(hào) 139。), set(gca,39。ylim39。,[4,4])。 subplot(6,3,10)。d=MixedS(1,:)。 y1=fft(d, 200)。 %對(duì)信號(hào)做 200 點(diǎn) FFT 變換 y=abs(y1)。%取模，整合幅頻和相頻 plot(y,39。r39。) title(39?；旌闲盘?hào) 1 頻譜 39。)。 axis([0 300 0 200])。 subplot(6,3,8),plot(MixedS(2,:)),title(39?；旌闲盘?hào) 239。), set(gca,39。ylim39。,[4,4])。 subplot(6,3,11)。e=MixedS(2,:)。 y1=fft(e, 200)。 %對(duì)信號(hào)做 200 點(diǎn) FFT 變換 y=abs(y1)。%取模，整合幅頻和相頻 plot(y,39。r39。) title(39?；旌闲盘?hào) 2 頻譜 39。)。 axis([0 300 0 200])。 subplot(6,3,9),plot(MixedS(3,:)),title(39。混合信號(hào) 339。), set(gca,39。ylim39。,[4,4])。 subplot(6,3,12)。f1=MixedS(3,:)。 y1=fft(f1, 200)。 %對(duì)信號(hào)做 200 點(diǎn) FFT 變換 y=abs(y1)。%取模，整合幅頻和相頻 plot(y,39。r39。) title(39?；旌闲盘?hào) 3 頻譜 39。)。 axis([0 300 0 200])。 武漢理工大學(xué)《信息處理課群綜合訓(xùn)練與設(shè)計(jì)》報(bào)告 27 MixedS_bak=MixedS。 % 將混合后的數(shù)據(jù)備份，以便在恢復(fù)時(shí)直接調(diào)用 %%%%%%%%%%%%%%%% 標(biāo)準(zhǔn)化 %%%%%%%%%%%%% MixedS_mean=zeros(3,1)。 for i=1:3 MixedS_mean(i)=mean(MixedS(i,:))。 end % 計(jì)算 MixedS 的均值 for i=1:3 for j=1:size(MixedS,2) MixedS(i,j)=MixedS(i,j)MixedS_mean(i)。 end end %%%%%%%%% 白化 %%%%%%%% MixedS_cov=cov(MixedS39。)。 % cov 為求協(xié)方差的函數(shù) [E,D]=eig(MixedS_cov)。 % 對(duì)矩陣的協(xié)方差函數(shù)進(jìn)行特征值分解 Q=inv(sqrt(D))*(E)39。 % Q 為白化矩陣 MixedS_white=Q*MixedS。 % MixedS_white 為白化后的矩陣 IsI=cov(MixedS_white39。)。 % IsI 應(yīng)為單位陣 %%%%%% FASTICA 算法 %%%%%%% X=MixedS_white。 % 以下算法將對(duì) X 進(jìn)行操作 [VariableNum,SampleNum]=size(X)。 numofIC=VariableNum。 % 在此應(yīng)用中，獨(dú)立元個(gè)數(shù)等于變量個(gè)數(shù) B=zeros(numofIC,VariableNum)。 % 初始化列向量 w 的寄存矩陣 ,B=[b1 ...] for r=1:numofIC i=1。maxIterationsNum=100。 % 設(shè)置最大迭代次數(shù) IterationsNum=0。 b=rand(numofIC,1)。 % 隨機(jī)設(shè)置 b 初值 b=b/norm(b)。 % 對(duì) b 標(biāo)準(zhǔn)化 norm(b):向量元素平方和開根號(hào) while i=maxIterationsNum+1 if i == maxIterationsNum % 循環(huán)結(jié)束處理 fprintf(39。\n 第 %d 分量在 %d 次迭代內(nèi)并不收斂。 39。, r,maxIterationsNum)。 武漢理工大學(xué)《信息處理課群綜合訓(xùn)練與設(shè)計(jì)》報(bào)告 28 break。 end bOld=b。 a2=1。 u=1。 t=X39。*b。 g=t.*exp(a2*t.^2/2)。 dg=(1a2*t.^2).*exp(a2*t.^2/2)。 b=((1u)*t39。*g*b+u*X*g)/SampleNummean(dg)*b。 % 核心公式，參見理論部分公式 b=bB*B39。*b。 % 對(duì) b 正交化 b=b/norm(b)。 if abs(abs(b39。*bOld)1)1e9 % 如果收斂，則 B(:,r)=b。 % 保存所得向量 b break。 end i=i+1。 end B(:,r)=b % 保存所得向量 b end %%%%%%%%% ICA 計(jì)算的數(shù)據(jù)復(fù)原并構(gòu)圖 %%%%%%%% ICAedS=B39。*Q*MixedS_bak。 % 計(jì)算 ICA 后的矩陣 % 將混合矩陣重新排列并輸出 subplot(6,3,13),plot(ICAedS(1,:)),title(39。ICA 解混信號(hào) 139。), set(gca,39。ylim39。,[4,4])。 subplot(6,3,16)。g=ICAedS(1,:)。 y1=fft(g,200)。 %對(duì)信號(hào)做 100 點(diǎn) FFT 變換 y=abs(y1)。%取模，整合幅頻和相頻 plot(y,39。r39。) title(39。分解信號(hào) 1 頻譜 39。)。 axis([0 300 0 200])。 subplot(6,3,14),plot(ICAedS(2,:)),title(39。ICA 解混信號(hào) 239。), set(gca,39。ylim39。,[4,4])。 武漢理工大學(xué)《信息處理課群綜合訓(xùn)練與設(shè)計(jì)》報(bào)告 29 subplot(6,3,17)。h=ICAedS(2,:)。 y1=fft(h,200)。 %對(duì)信號(hào)做 100 點(diǎn) FFT 變換 y=abs(y1)。%取模，整合幅頻和相頻 plot(y,39。r39。) title(39。分解信號(hào) 2 頻譜 39。)。 axis([0 300 0 200])。 subplot(6,3,15),plot(ICAedS(3,:)),title(39。ICA 解混信號(hào) 339。), set(gca,39。ylim39。,[4,4])。 subplot(6,3,18)。i=ICAedS(3,:)。 y1=fft(i,200)。 %對(duì)信號(hào)做 100 點(diǎn) FFT 變換 y=abs(y1)。%取模，整合幅頻和相頻 plot(y,39。r39。) title(39。分解信號(hào) 3 頻譜 39。)。 axis([0 300 0 200])。 武漢理工大學(xué)《信息處理課群綜合訓(xùn)練與設(shè)計(jì)》報(bào)告 本科生綜合訓(xùn)練成績評(píng)定表 姓 名 性 別 專業(yè)、班級(jí) 綜合訓(xùn)練題目： 綜合訓(xùn)練答辯或質(zhì)疑記錄： 成績評(píng)定依據(jù)： 最終評(píng)定成績（以優(yōu)、良、中、及格、不及格評(píng)定） 指導(dǎo)教師 簽字： 年 月 日