Matlab频域幅值输出

Y = fft(xt);
Y_M1 = abs(Y/L); % 求幅值

P1 = Y_M1(1:L/2+1);
P1(2:end-1) = 2*P1(2:end-1);  % * 2 恢复幅值。
f = fs * (0:L/2) / L; % 一半，因为是对称的
figure(2);
plot(f,P1);
xlabel('f/Hz');
ylabel('原始幅频特性');