资 源 简 介
做盲源分离的FASTICA算法clc;clear all;close all;function y=ICA(M1,M2)MixedS = [M1;M2];% 将混合声音显示subplot(3,2,3),plot(data3),title("混合声音1"),subplot(3,2,4),plot(data4),title("混合声音2"),% wavwrite(MixedS(1,:),Fs1,"MixedS1.wav");% wavwrite(MixedS(2,:),Fs1,"MixedS2.wav");MixedS_bak=MixedS; % 将混合后的数据备份，以便在恢复时直接调用%%%%%%%%%%%%%%%%%%%%%%%%%% 中心化 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%MixedS_mean=zeros(2,1);MixedS_mean=mean(MixedS,2);MixedS = MixedS-repmat(MixedS_mean,1,size(MixedS,2));%%%%%%%%%%%%%%%%%%%%%%%%%%% 白化 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%MixedS_cov=cov(MixedS"); % cov为求协方差的函数[E,D]=eig(MixedS_cov); % 对图片矩阵的协方差函数进行特征值分解
