关于信号发生器的相位误差

%关于相位误差
clear
close all
clc

fs=100; %采样率
fc=10;  %载波
fpe=1;  %相位噪声

t=0:1/fs:100;
c_pe=cos(2*pi*fpe*t);
c_c=cos(2*pi*fc*t);
c_cpe=cos(2*pi*fc*t+c_pe);

C_C=fft(c_c);
C_C=abs(C_C)/length(C_C);
C_CPE=fft(c_cpe);
C_CPE=abs(C_CPE)/length(C_CPE);

F=linspace(0,fs,length(t));
N=floor(length(F)/2);

plot(F(1:N),C_C(1:N),'r');
hold on
plot(F(1:N),C_CPE(1:N));
set(gca,'xtick',[0:1:F(N)]);
grid on;