%%E7_1_DPSKModem.m程序清單
%仿真DPSK信號的產(chǎn)生及相干解調(diào)過程
%將成形濾波器系數(shù)量化為10比特后，寫入tra_lpf.coe文件中
%將接收濾波器系數(shù)量化為10比特后，寫入rec_lpf.coe文件中
%將DPSK已調(diào)數(shù)據(jù)量化為8比特后，寫入Dpsk.txt文件中
ps=1*10^6;???%碼速率為1MHz
a=0.8;???????%成形濾波器系數(shù)為0.8
B=（1+a）*ps;??%中頻信號處理帶寬
Fs=8*10^6;??%采樣速率為8MHz
fc=2*10^6;???%載波頻率為2MHz
N=20000;??????%仿真數(shù)據(jù)的長度

t=0:1/Fs:（N*Fs/ps-1）/Fs;%產(chǎn)生長度為N頻率為fs的時間序列
s=randint（N12）;???????%產(chǎn)生隨機數(shù)據(jù)作為原始數(shù)據(jù)
%并將絕對碼變換為相對碼
ds=ones（1N）;
for?i=2:N
????if?s（i）==1?
????????ds（i）=-ds（i-1）;
????else
????????ds（i）=ds（i-1）;
????end
end
%對相對碼數(shù)據(jù)以Fs頻率采樣
Ads=upsample（dsFs/ps）;

%設(shè)計平方升余弦濾波器
n_T=[-2?2];
rate=Fs/ps;
T=1;
Shape_b?=?rcosfir（an_TrateT）;%figure（4）;freqz（Shape_b）
%對采樣后的數(shù)據(jù)進行升余弦濾波;
rcos_Ads=filter（Shape_b1Ads）;

%產(chǎn)生載頻信號
f0=sin（2*pi*fc*t）;???????
%產(chǎn)生DPSK已調(diào)信號
dpsk=rcos_Ads.*f0;???????

%與相干載波相乘