clc;clear;close all;
Adc=2;%直流分量幅值
A1=3;%频率F1信号的幅值
A2=1.5;%频率F2信号的幅值
F1=50;%信号1频率（Hz）
F2=75;%信号2频率
Fs=256;%采样频率
P1=-30;%信号1相位
P2=90;%信号2相位
N=256;%采样点数
t=[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);
%显示原始信号
plot(S);
title('原始信号');
figure;
Y=fft(S,N);%做FFT变换
Ayy=(abs(Y));%取模
plot(Ayy(1:N));%显示原始的FFT模值结果
title('FFT模值');
figure;
Ayy=Ayy/(N/2);%换算成实际的幅值
Ayy(1)=Ayy(1)/2;
F=([1:N]-1)*Fs/N;%换算成实际的频率值
plot(F(1:N/2),Ayy(1:N/2));%显示换算后的模值结果
title('幅度-频率曲线图');
figure;
Pyy=[1:N/2];
for i=1:N/2
Pyy(i)=phase(Y(i)); %计算单位
Pyy(i)=Pyy(i)*180/pi; %换算为角度
end
plot(F(1:N/2),Pyy(1:N/2));%显示相位图
title('相位-频率曲线图');
