clear?all;?
close?all;
clc;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%???變量定義、賦值
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
SNR=102;??????????????????????????????????????%定義信噪比，單位DB
fc=10e6;?????????????????????????????????????%最小載頻??最大*2
h?=?0.6;?????????????????????????????????????%調制指數
fz=1000;?????????????????????????????????????%衰落頻率
fs=50e6;?????????????????????????????????????%采樣頻率
Rb=1e6;??????????????????????????????????????%碼元速率
Nu=1000;??????????????????????????????????????%碼元個數
Nn=500;??????????????????????????????????????%RMS取值個數?
Vre=1.5;?????????????????????????????????????%比較標準值
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
N=fs/fc;?????????????????????????????????????%采樣點數??
Tc=1/fc;?????????????????????????????????????%載波周期
Tb=1/Rb;?????????????????????????????????????%碼元周期
ts?=?1/fs;

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%???載波調制
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x?=?randsrc（1Nu[1-1]）;????????????????????%信息生成
phase0?=?0;??????????????????????????????????%載波相位，在相干解調情況下可設為0
add_phase?=?[];??????????????????????????????%附加相位
phase1?=?0;
time?=?ts:ts:Tb;?????????????????????????????%一個碼元周期內采樣點數
for?i?=?1:length（x）
????q?=?time/（2*Tb）;
????if?i?~=?1
????????theta_k?=?pi*h*sum（x（1:i?-?1））;
????else
????????theta_k?=?0;
????end
????phase1?=?theta_k?+?2*pi*h*x（i）*q;
????add_phase?=?[add_phasephase1];
end
cpfsk_num?=?0:ts:（length（add_phase）?-?1）*ts;?%cpfsk的點數
yca?=?cos（2*pi*fc*cpfsk_num?+?add_phase?+?phase0）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%??加衰落加噪聲
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%?ynoise=awgn（sqrt（2）*ycaSNR‘measured‘）;???????%加