>> x = randn(1, 10000); >> A = sum(A.^2) - sum(fft(x).^2)/length(x); % 根据帕斯瓦尔定理,傅里叶变换系数必须除以 N,才会等于真正的频谱值; % 本例中才会接近于 0;