nolinear_ica 非线性卷积盲源分离程序

非线性卷积盲源分离程序，适用于初学盲源分离内容的学生，很有帮助

clear all
clc
%产生非线性卷积盲信号分离

%1.产生源信号
T=0.1;        %取时间区间为[0s  0.1s]
Num=1000;     %10k采样率下，采了1000个点
M=3;          %信源数=传感器数=3
Lum=50;       %卷积的移位值
Kum=100;      %STFT的窗长度
a=-0.1;b=0.1;
t=linspace（0TNum）;
S1=sign（cos（2*pi*155*t））;   %符号信号
S2=sin（2*pi*800*t）;         %高频正弦信号
S3=sin（2*pi*9*t）.*sin（2*pi*300*t）;      %幅度调制信号
S=[S1;S2;S3];

%2.产生混合矩阵A，且A为元胞数组
A=cell（MM）;  %分别取不同时刻的混合矩阵完成卷积
for i=1:M
    for j=1:M
        A{ij}=a+（b-a）*rand（1Lum）;
    end
end

X=zeros（size（S））;
ans1=0;ans2=0;ans3=0;
for i=1:M
    for k=1:Num
        for j=1:M
            for p=0:Lum-1
                if （k-p>0）&（k-p<=Num）
                    ans1=A{ij}（p+1）*S（jk-p）;
                else ans1=0;
                end
            ans2=ans2+ans1;
            end
         

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

     文件       1667  2008-10-08 10:54  feixianxinghunhe.m

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

                 1689                    2