简介：面向信号处理、语音分离与盲源分离入门者及MATLAB实践者，这份资料围绕JADE盲源分离算法的原理与程序实现展开。内容从盲信号分离的基本问题切入，说明观测混合信号、源信号统计独立性以及“盲”的双重含义，并重点介绍JADE利用信号四阶累积量、借助源信号非高斯性完成分离的思路，同时提及三阶累积量的用法及最多允许一个高斯源等约束，便于读者理解HOS类算法的适用边界。压缩包为zip格式，约426KB，文件总数与类型明细上游暂未提供，需以实际下载内容为准；资料兼具MATLAB实现与算法讲解属性，适合结合代码对照理论进行实验复现和课程学习。已有953人学习，可作为语音信号处理与盲源分离方向的入门参考。
1. 盲源分离为什么绕不开“四阶累积量”：JADE 的切入点
先给一个反直觉的结论：JADE 能把混在一起的信号拆开，靠的核心依据，并不是“信号本来应该长什么样”，真正起作用的判据是“信号有多不像高斯分布”。二阶统计量只能说明两路信号之间有没有相关性，而任何正交旋转都不改变相关性；四阶累积量对高斯变量恒等于零，只有非高斯源才会在上面留下非零投影。JADE 的全部结构，包括白化、四阶张量、特征矩阵和联合对角化，本质上都是在把这一点残留的投影稳定地提取出来。一份完整的 matlab 算法实现里，真正需要做判断的其实只有累积量构造和联合对角化这两处，其余步骤主要是在为它们准备一个更干净的坐标系。把这条主线理顺之后，理解后面每一步的动机就会顺畅很多。
2. 白化：为什么把混合矩阵掰成正交之后问题少了一半
盲源分离的线性模型写成 x = A s，其中 x 是 n 维观测，s 是 m 维源，A 是 n×m 的混合矩阵。直接估计 A 要从 m×n 个自由参数里搜索，规模一大基本就失去方向。JADE 第一步做的是白化，把“估计任意矩阵 A”变成“估计一个正交阵 U”。参数从 m² 个降到 m(m-1)/2 个，m=8 时就是从 64 降到 28，而且正交性天然带来一系列可以代数求解的性质。
2.1 协方差特征分解与 W = Λ^{-1/2}Vᵀ 的来历
若源信号已经归一化为零均值、单位协方差，那么观测协方差就是
R_x = E[xxᵀ] = A E[ssᵀ] Aᵀ = A Aᵀ
对 R_x 做特征分解 R_x = V Λ Vᵀ，把特征值从大到小排好，取前 m 个。定义
W = Λ_m^{-1/2} V_mᵀ
W 是 m×n 的矩阵。此时
E[(Wx)(Wx)ᵀ] = W R_x Wᵀ = Λ_m^{-1/2}V_mᵀ (VΛVᵀ) V_m Λ_m^{-1/2} = I
也就是白化后的数据 z = Wx 各通道之间已经完全去相关、并且功率归一。更关键的是 z = W A s，记 U = W A，则 U Uᵀ = I，所以 U 是一个 m×m 的正交阵。也就是说，白化之后剩下的唯一未知量只是一个“旋转”。
这一步之所以成立，前提是源的二阶统计量是单位阵。源功率各不相同也不影响结论，只是最终的 U 会多吸收一个对角缩放，分离出来的信号变成单位功率版本。这个缩放后面要靠一个参考通道或者主观判断拉回来。
2.2 MATLAB 里白化的三行核心代码与特征值排序陷阱
function [W, Z] = jade_whiten(X, m) % X: n x T 观测矩阵，每行一个通道 % m: 估计的源数目（必须 <= n） % W: m x n 白化矩阵；Z: m x T 白化数据 [n, T] = size(X); % 1) 去均值。ICA 模型里源是零均值的，观测必须先对齐 X = X - mean(X, 2); % 2) 样本协方差 Rx = (X * X') / T; % 3) 特征分解。eig 不保证特征值有序，必须手动排 [V, D] = eig(Rx); [d, idx] = sort(diag(D), 'descend'); V = V(:, idx); % 4) 只保留前 m 个主方向 d = d(1:m); V = V(:, 1:m); % 5) 白化矩阵：先乘 Λ^{-1/2}，再乘 V' W = diag(1 ./ sqrt(d)) * V'; <p> <a href="https://download.csdn.net/download/wzz110011/87607398" style="color:#ec7500;font-size:14px;"> 本文还有配套的精品资源，点击获取 </a> <img alt="menu-r.4af5f7ec.gif" src="https://csdnimg.cn/release/wenkucmsfe/public/img/menu-r.4af5f7ec.gif" style="width:16px;margin-left:4px;vertical-align:text-bottom;cursor:text;"> </p>)
