咨询：MATLAB中使用ifft还原时域信号时的缩放问题修复方法

问题核心原因 #

1. 负频率构造错误：原代码将包含0频率的完整X翻转共轭，导致频域信号重复了0频率分量，破坏了频域对称性，同时使信号长度错误。
2. 缺少幅值缩放：MATLAB的ifft默认会对结果除以点数N，但连续傅里叶逆变换的数值近似需要结合频率步长domega进行缩放，才能匹配解析解的幅值。
3. 时间轴计算偏差：原时间轴基于Fn=2*max(f)的经验值，未结合频域信号的点数和频率步长，导致时间轴与频域信号不匹配。

修正后的完整代码 #

%system parameters
m=1;
k=1;
omega_n=sqrt(k/m);
c=0.2;
xi=c/(2*sqrt(k*m));
omega_d=omega_n*sqrt(1-xi^2);

%frequency range
domega = 0.1;
omega_pos = 0:domega:10; % 正频率分量（包含0频率）
X_pos = 1./(-omega_pos.^2*m + 1i*c.*omega_pos + k); 

% 构造完整频域信号：负频率为正频率（去掉0点）的共轭翻转
X_neg = conj(flip(X_pos(2:end)));
X_full = [X_pos, X_neg];
N = length(X_full);

% 计算匹配的时间轴：由DFT频率步长推导时间步长dt=2π/(N*domega)
dt = 2*pi/(N*domega);
t = 0:dt:(N-1)*dt;

% IFFT重构并施加幅值缩放：补偿ifft的默认除法，匹配连续傅里叶逆变换
time_domain_signal = ifft(X_full, 'symmetric') * N * domega/(2*pi);

%analytical solution
x_impl= exp(-xi*omega_n.*t).*sin(omega_d*t)/(m*omega_d); 

figure
plot(t,time_domain_signal)
hold on
plot(t,x_impl,'--')
legend('IFFT重构信号','解析解')
xlabel('时间 t')
ylabel('幅值')

关键修正点说明 #

• 负频率处理：仅对正频率中除0点外的部分取共轭翻转，作为负频率分量，保证频域信号的对称性和长度正确。
• 幅值缩放：通过N*domega/(2*pi)因子，将ifft的离散结果转换为符合连续傅里叶逆变换的幅值，抵消ifft默认的除以N操作，并补偿积分近似的频率步长。
• 时间轴匹配：基于DFT的频率-时间关系推导时间步长，确保时间轴与频域信号的采样特性完全匹配。

内容的提问来源于stack exchange，提问作者Eni