%%E7_1_DPSKModem.m程序清單
%仿真DPSK信號(hào)的產(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;??%中頻信號(hào)處理帶寬
Fs=8*10^6;??%采樣速率為8MHz
fc=2*10^6;???%載波頻率為2MHz
N=20000;??????%仿真數(shù)據(jù)的長(zhǎng)度

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

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

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

%與相干載波相乘