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