ＩＣＡ算法实现代码

非常实用，独立成分分析（ｉｃａ）ｍａｔｌａｂ实现，希望能够对你有用。

%%%%%%%%%%%%%%%%%%%%%%%%%%  初始化  %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

clc;clear all;close all;

%%%%%%%%%%%%%%  读入原始图像，混合，并输出混合图像  %%%%%%%%%%%%%%%%%%

% 读入混合前的原始图片并显示
t=0:1/100:9;
I1=sin（t）; 
I2=randn（1901）;
I3=square（4*t）;
subplot（431）plot（I1）title（‘输入信号1‘）
subplot（432）plot（I2）title（‘输入信号2‘）
subplot（433）plot（I3）title（‘输入信号3‘）

% 将其组成矩阵
S=[I1;I2;I3];                          % 图片个数即为变量数，图片的像素数即为采样数
                                       % 因此S_all是一个变量个数＊采样个数的矩阵

Sweight=randn（size（S1））;              % 取一随机矩阵，作为信号混合的权矩阵
MixedS=Sweight*S;                      % 得到三个混合信号矩阵

% 将混合矩阵重新排列并输出
subplot（434）plot（MixedS（1:））title（‘混合信号1‘）
subplot（435）plot（MixedS（2:））title（‘混合信号2‘）
subplot（436）plot（MixedS（3:））title（‘混合信号3‘）

MixedS_bak=Mixed

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       4108  2006-06-08 11:12  ica示例\signalOK.m

     目录          0  2008-09-03 13:03  ica示例

----------- ---------  ---------- -----  ----

                 4108                    2