whiten.m

来自「这是一个基于峰度的ICA算法」· M 代码 · 共 13 行

function [V,z]=whiten(x)
%这是一个对信号进行白化的函数
%x要求是观察信号的列向量x=[x1;x2;x3;...;xn],返回值中V是白化矩阵，z是白化后的
%结果。
[vectorSize, numSamples] = size(x);
for i=1:numSamples
    x(:,i)=x(:,i)-(mean(x'))';
end
B=cov(x');
[E,D]=eig(B);
V=E*D^(-1/2)*E';
z=V*x;