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

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%???載波調(diào)制
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
x?=?randsrc（1Nu[1-1]）;????????????????????%信息生成
phase0?=?0;??????????????????????????????????%載波相位，在相干解調(diào)情況下可設(shè)為0
add_phase?=?[];??????????????????????????????%附加相位
phase1?=?0;
time?=?ts:ts:Tb;?????????????????????????????%一個(gè)碼元周期內(nèi)采樣點(diǎn)數(shù)
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的點(diǎn)數(shù)
yca?=?cos（2*pi*fc*cpfsk_num?+?add_phase?+?phase0）;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%??加衰落加噪聲
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

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