matlab fastica工具箱安装

尝试使用独立成分分析和主成分分析，将融合到一起的两组正弦函数分开。但独立成分分析无法运行，需下载工具箱。因经典的ICA算法已经基本不再使用，而是被一种名为FastICA的改进算法替代。顾名思义，该算法的优点在与Fast，即运算速度快。
FastICA工具箱下载地址：https://research.ics.aalto.fi/ica/fastica/code/dlcode.shtml

%FastICA程序
 % 清空环境
clear;
clc;

% 参数设置
fs = 1000;          % 采样频率 (Hz)
t = 0:1/fs:1-1/fs;  % 时间向量 (1秒的时间)

% 第一个正弦波的参数
A1 = 1;             % 幅度
f1 = 5;             % 频率 (Hz)
phi1 = 0;           % 相位 (弧度)

% 第二个正弦波的参数
A2 = 0.5;           % 幅度
f2 = 50;            % 频率 (Hz)
phi2 = pi/4;        % 相位 (弧度)

% 生成正弦波
y1 = A1 * sin(2 * pi * f1 * t + phi1);
y2 = A2 * sin(2 * pi * f2 * t + phi2);

% 叠加两个正弦波
y = y1 + y2;

% 生成混合信号矩阵
X = [y1; y2]';

% 独立成分分析
[icasig, ~, ~] = fastica(X');

% 绘制结果
figure;
subplot(3, 1, 1);
plot(t, y);
title('叠加后的信号');
xlabel('时间 (s)');
ylabel('幅度');
grid on;

subplot(3, 1, 2);
plot(t, icasig(1, :));
title('ICA 第一个独立成分');
xlabel('时间 (s)');
ylabel('幅度');
grid on;

subplot(3, 1, 3);
plot(t, icasig(2, :));
title('ICA 第二个独立成分');
xlabel('时间 (s)');
ylabel('幅度');
grid on;



%主成分分析
% 清空环境
clear;
clc;

% 参数设置
fs = 1000;          % 采样频率 (Hz)
t = 0:1/fs:1-1/fs;  % 时间向量 (1秒的时间)

% 第一个正弦波的参数
A1 = 1;             % 幅度
f1 = 5;             % 频率 (Hz)
phi1 = 0;           % 相位 (弧度)

% 第二个正弦波的参数
A2 = 0.5;           % 幅度
f2 = 50;            % 频率 (Hz)
phi2 = pi/4;        % 相位 (弧度)

% 生成正弦波
y1 = A1 * sin(2 * pi * f1 * t + phi1);
y2 = A2 * sin(2 * pi * f2 * t + phi2);

% 叠加两个正弦波
y = y1 + y2;

% 生成特征矩阵
X = [y1; y2]';

% 主成分分析
[coeff, score, ~] = pca(X);

% 绘制结果
figure;
subplot(3, 1, 1);
plot(t, y);
title('叠加后的信号');
xlabel('时间 (s)');
ylabel('幅度');
grid on;

subplot(3, 1, 2);
plot(t, score(:, 1));
title('PCA 第一个主成分');
xlabel('时间 (s)');
ylabel('幅度');
grid on;

subplot(3, 1, 3);
plot(t, score(:, 2));
title('PCA 第二个主成分');
xlabel('时间 (s)');
ylabel('幅度');
grid on;