单通道盲源分离程序

基于ICA的单通道盲源分离算法，程序代码方便读者理解，里面包含了相关数据源和程序包

function [ykalman] = ex1kalman（）
    load（‘Data/simdata.mat‘‘t‘‘y‘‘ynf‘‘ysin‘）;
    % State = （y ydot）. We only observe （y）.
    % X（t+1） = F X（t） + noise（Q）
    % Y（t） = H X（t） + noise（R）
    fprintf（1 ‘-----------------Kalman-----------------\n‘）;
    ss = 2; % state size
    os = 1; % observation size
    F = [1 0; 0 1]; 
    H = [1 0];
    Q = 0.1*eye（ss）;
    R = 1*eye（os）;
    initx = [0 0]‘;
    inity = 1*eye（ss）;

    %[xfilt Vfilt VVfilt loglik] = kalman_filter（y F H Q R initx inity）;
    ykalman = kalman_smoother（y F H Q R initx inity）;

    % dfilt = x（1:） - xfilt（1:）;
    % mse_filt = sqrt（sum（sum（dfilt.^2）））

%     dsmooth = x（1:） - xsmooth（1:）;
%     mse_smooth = sqrt（sum（sum（dsmooth.^2）））
    ykalman = ykalman（1:）;
    
    % 带阻滤波，消除50Hz工频干扰
     % 经频谱分析可知，有用信号的频率成分可能落在工频带，因此不能使用此方法
%     Hd = bandstop;
%     ykalman = filter（Hdykalman）;
    
    %计算信噪比
    eta = ynf - ykalman; % Residual Signal
    snrpkalman = 10*log（ （ynf*ynf‘）/（eta*eta‘） ）;   
%     figure
%     plot（tynf‘r‘tynf+ysin‘m‘tykalman‘b‘）;
%     legend（‘ynf‘ ‘ynf+ysin‘ ‘ykalman‘）
%     title（‘Kalman‘）;
    
    save（‘Data/kalman.mat‘‘ykalman‘‘snrpkalman‘）;
end

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

     文件      69460  2013-07-11 22:15  ircbp2.1_wkp\Data\dictCA.mat

     文件      15987  2013-07-11 22:15  ircbp2.1_wkp\Data\realdata.mat

     文件      16063  2013-07-11 22:15  ircbp2.1_wkp\Data\simdata.mat

     文件       1337  2013-07-11 22:15  ircbp2.1_wkp\ex1kalman.m

     文件        748  2013-07-11 22:15  ircbp2.1_wkp\ex2sgolay.m

     文件       1065  2013-07-11 22:15  ircbp2.1_wkp\ex3bp.m

     文件       9925  1980-01-01 14:04  ircbp2.1_wkp\ex4ircbp.m

     文件       9955  1980-01-01 14:04  ircbp2.1_wkp\ex4ircbp.m~

     文件       1453  2013-07-11 22:15  ircbp2.1_wkp\getdata.m

     文件       2145  2013-07-11 22:15  ircbp2.1_wkp\getdictCA.m

     文件   16393591  2013-12-07 23:32  ircbp2.1_wkp\toolbox\cvx-a64.tar.gz

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

             16521729                    11