如何对信号进行快速傅里叶变换?

%产生峰值为1的三角波，分析其0~63次谐波的幅值谱和相位谱

clf;

Fs =128;                   %采样频率

T = 1/Fs;                  % 采样周期

N = 128;                   % 采样点数

t = (0:N-1)*T;             % 时间，单位：S

x=zeros(N);

for n=0:N-1

    b=fix((n)/(N/4));

    m=n+1;

    A=1/(N/4);

    if b==0 

        x(m)=A*n;

    elseif b==1||b==2

        x(m)=A*(N/2-n);

    elseif b==3

        x(m)=A*(n-N);

    end;

end;

n=0:N-1;

subplot(3,1,1)

plot(t,x);

xlabel('时间/S');

ylabel('振幅');title('时域波形');grid on;

y=fft(x,N);   %对信号进行快速Fourier变换

mag=abs(y)*2/N;   %求取Fourier变换的振幅;*2/N转变为真实幅值

f=n*Fs/N;

subplot(3,1,2)

plot(f(1:N/2),mag(1:N/2)); %绘出Nyquist频率之前随频率变化的振幅

xlabel('频率/Hz');

ylabel('振幅');title('幅值谱');grid on;

p=mod(angle(y)*180/pi,360);

subplot(3,1,3)

plot(f(1:N/2),p(1:N/2)); %绘出Nyquist频率之前随频率变化的相位

xlabel('频率/Hz');

ylabel('振幅');title('相位谱');grid on;

本文如未解决您的问题请添加抖音号：51dongshi（抖音搜索懂视），直接咨询即可。