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('原始幅频特性');
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('原始幅频特性');