FASTica算法matlab代码

clear
clc
 
%观察信号矩阵 X 为 K*N 矩阵
%源信号矩阵 S 为 K*N 矩阵
%混合矩阵 A 为 K*K 矩阵
 
K=2;
%生成 K 个源信号, 存储到 S
f1=30;f2=15;f3=120;
Fs=3000;
N=1000;
t=0:1/Fs:(N-1)/Fs;
 
figure(1);
%S(1,:)=sin(2*pi*f3*t).*[1+1.5*sin(2*pi*f2*t)]; 
S(1,:)=0.3*cos(2*pi*f1*t+0.5*sin(2*pi*f2*t))+0.2*sin(2*pi*f3*t);
subplot(2,1,1);
plot(t,S(1,:),'k');
set(gca,'fontname','Times New Roman','fontsize',9);
title('信号1');
xlabel('time/ms','fontname','Times New Roman','fontsize',9);
ylabel('amplitude','fontname','Times New Roman','fontsize',9);
hold on;
 
S(2,:)=sin(2*pi*f1*t)+sin(2*pi*f2*t); 
subplot(2,1,2);
plot(t,S(2,:),'k');
set(gca,'fontname','Times New Roman','fontsize',9);
title('信号2');xlabel('time/ms','fontname','Times New Roman','fontsize',9);