请教MATLAB相关问题 如图
程序如下:
%下程序为ICA的调用函数,输入为观察的信号,输出为解混后的信号
function Z=ICA(X)
%以下为主程序,主要为原始信号的产生,观察信号和解混信号的作图
clear all;
clc;
N=10000;
n=1:N;%N为采样点数
s1=load(‘lvbo1.mat’);
s2=load(‘lvbo2.mat’);
s3=load(‘lvbo3.mat’);
s4=load(‘lvbo6.mat’);
S=[s1;s2;s3;s4];%信号组成4*N
A=rand(4,4);
X=A*S;%观察信号
%-----------去均值---------
[M,T] = size(X); %获取输入矩阵的行/列数,行数为观测数据的数目,列数为采样点数
average=mean(X’)’; %均值
for i=1:M
X(i,:)=X(i,:)-average(i)*ones(1,T);
end
%---------白化/球化------
Cx =cov(X’,1); %计算协方差矩阵Cx
[eigvector,eigvalue]= eig(Cx); %计算Cx的特征值和特征向量
W=eigvalue^(-1/2)eigvector’; %白化矩阵
Z=WX; %正交矩阵
%----------迭代-------
Maxcount=10000; %最大迭代次数
Critical=0.00001; %判断是否收敛
m=M; %需要估计的分量的个数