Untitled.m

clear;
clf;
close all
Fs=800e3;
fz=80e3;
fj=1e3;
i=4096;
A=1;
Kf=0.2;             
t=（1:i）/Fs;         
tf=（1:i-1）/Fs;      
f=（（1:i）-1）*Fs/i;   
tz=1;               

Mod_Sig_c=A*cos（2*pi*fj*t）;
%画图基带信号
subplot（311）;plot（tMod_Sig_c）;title（‘基带信号‘）;xlabel（‘时间（S）‘）;ylabel（‘幅度（V）‘）;grid on;legend（‘基带‘）;

Carrier_c=cos（2*pi*fz/Fs*（1:i））;%载波cos
%画载波
subplot（312）;plot（t（1:300）Carrier_c（1:300））;title（‘载波信号‘）;xlabel（‘时间（S）‘）;ylabel（‘幅度（V）‘）;grid on;legend（‘载波‘）;

int_m（1）=0;      %对m_fun积分
for p=1:length（1:i）-1
    int_m（p+1）=int_m（p）+Mod_Sig_c（p）*tz;
end

SFm_c=A*cos（2*pi*fz/Fs*（1:i）+int_m*Kf）;
SFm   =   SFm_c;
% SFm=awgn（SFm_c0.00000022）;
%画调制信号
subplot（313）;plot（tSFm）;title（‘调制信号‘）;xlabel（‘时间（S）‘）;ylabel（‘幅度（V）‘）;grid on;legend（‘调制信号‘）;
ff_SFm=fft（SFm）;
ff_S