%%?調(diào)制
J?=?sqrt（-1）;??????????????????????????????%%虛數(shù)符號(hào)
num?=?1000;
data?=?round（rand（1num））;%%產(chǎn)生一行隨機(jī)數(shù)
nsymbol?=?length（data）/2;?????????????%%計(jì)算發(fā)送符號(hào)數(shù)
qdpsk_mod?=?zeros（1nsymbol）;
mid_data?=?zeros（1num+2）;
qdpsk_data?=?zeros（1num）;
%%?QDPSK的調(diào)制模塊??（功率不作歸一化）
for?loop?=?1:nsymbol
????if?xor（mid_data（2*loop-1）mid_data（2*loop））?==?0
????????mid_data（2*loop+1）?=?xor（mid_data（2*loop-1）data（2*loop-1））;
????????mid_data（2*loop+2）?=?~xor（mid_data（2*loop）data（2*loop））;
????elseif?xor（mid_data（2*loop-1）mid_data（2*loop））?==?1
????????mid_data（2*loop+1）?=?xor（mid_data（2*loop）data（2*loop））;
????????mid_data（2*loop+2）?=?~xor（mid_data（2*loop-1）data（2*loop-1））;
????end
end

for?loop?=?1:nsymbol
????qdpsk_data（2*loop-1）?=?mid_data（2*loop+1）;
????qdpsk_data（2*loop）?=?mid_data（2*loop+2）;
end

for?loop?=?1:nsymbol???????????????%%調(diào)制
????symbol1?=?qdpsk_data（2*loop-1）;
????symbol2?=?qdpsk_data（2*loop）;
????
????if?symbol1?==?0?&&?symbol2?==?0?
????????qdpsk_mod（loop）?=?1;
????elseif?symbol1?==?0?&&?symbol2?==?1?
????????qdpsk_mod（loop）?=?J;?
????elseif?symbol1?==?1?&&?symbol2?==?1?
????????qdpsk_mod（loop）?=?-1;
????elseif?symbol1?==?1?&&?symbol2?==?0?
????????qdpsk_mod（loop）?=?-J;
????end
end


r?=?qdpsk_mod;

%%與乘載波相乘模塊
%%定義參數(shù)
RB?=?50*1e3??;????????????%%基帶碼元速率（Hz）?50K
f?=?1*1e7?;?????????????%%載波頻率（Hz）100M
x?=?20?;??????????????????%%采樣頻率與載波頻率之比??（奈奎斯特定律，至少應(yīng)大于2）
A?=?1?;???????????????????%%載波幅值
time?=?[0:1/（x*f）:nsymbol/RB?-?1/（x*f）];??????????????????%%采樣時(shí)間點(diǎn)（單位為秒）?


C?=?exp（2*pi*f*time*J）;????????????%%e指數(shù)載波

for?loop1?=?1:nsymbol
????wavesignal（??（loop1-1）*（x*f/RB）+1:loop1*（x*f/RB）??）?=?real（???r（loop1