clc;clear;Adc=2; %直流分量幅度A1=3; %频率F1信号的幅度A2=1.5; %频率F2信号的幅度A3=2.2;%频率F3信号的幅度F1=50; %信号1频率(Hz)F2=75; %信号2频率(Hz)F3=75.5; %信号2频率(Hz)Fs=256; %采样频率(Hz)P1=-30; %信号1相位(度)P2=90; %信号相位(度)P3=60; %信号相位(度)N=1024; %采样点数%采样总时长4s ; 频率分辨率0.25Hz ; 相应的谐频点在201，301，303t=0:1/Fs:N/Fs; %采样时刻%信号S=Adc+A1*cos(2*pi*F1*t+pi*P1/180)+A2*cos(2*pi*F2*t+pi*P2/180)+A3*cos(2*pi*F3*t+pi*P3/180);figure;plot(t,S);%显示原始信号title(原始信号);Y = fft(S,N); %做FFT变换figure;Ayy = abs(Y); %取模n = 1:N;%plot(n,Ayy);stem(Ayy(1:N)); %显示原始的FFT模值结果title(FFT 模值);Ayy=Ayy/(N/2); %换算成实际的幅度Ayy(1)=Ayy(1)/2;F=([1:N]-1)*Fs/N; %换算成实际的频率值，Fn=(n-1)*Fs/Nfigure;plot(F,Ayy);%stem(Ayy(1:N/2)); %显示换算后的FFT模值结果title(幅度-频率曲线图);Pyy=[1:N/2];for i=1:N/2Pyy(i)=angle(Y(i)); %计算相位Pyy(i)=Pyy(i)*180/pi; %换算为角度end;figure;plot(F(1:N/2),Pyy(1:N/2));%stem(F(1:N/2),Pyy(1:N/2)); %显示相位图title(相位-频率曲线图);& 转自：http://ass.3daixie.com/2018061815675074.html

原著是一个有趣的人,若有侵权,请通知删除