調制解調程序清單
%
仿真
DPSK
信號的產生及相幹解調過程
%
將成形濾波器系數量化為
10
比特後，寫入
tra_lpf.coe
文件中
%
將接收濾波器系數量化為
10
比特後，寫入
rec_lpf.coe
文件中
%
將
DPSK
已調數據量化為
8
比特後，寫入
Dpsk.txt
文件中
ps=1*10^6;
%
碼速率為
1MHz
a=0.8;
%
成形濾波器系數為
0.8
B=(1+a)*ps;
%
中頻信號處理帶寬
Fs=8*10^6;
%
采樣速率為
8MHz
fc=2*10^6;
%
載波頻率為
2MHz
N=20000;
%
仿真數據的長度
t=0:1
/Fs:(N*Fs/ps-1)/Fs;%
產生長度為
N,
頻率為
fs
的時間序列
s=randint(N,1,2);
%
產生隨機數據作為原始數據
,
%
並將絕對碼變換為相對碼
ds=ones(1,N);
for i=2:N
if s(i)==1
ds(i)=-ds(i-1);
else
ds(i)=ds(i-1);
end
end
%
對相對碼數據以
Fs
頻率采樣
Ads=upsample(ds,Fs/ps);
%
設計平方升余弦濾波器
n_T=[-2 2];
rate=Fs/ps;
T=1;
Shape_b = rcosfir(a,n_T,rate,T);%figure(4);freqz(Shape_b)
%
對采樣後的數據進行升余弦濾波
;
rcos_Ads=filter(Shape_b,1,Ads);
%
產生載頻信號
f0=sin(2*pi*fc*t);
%
產生
DPSK
已調信號
dpsk=rcos_Ads.*f0;
%
與相幹載波相乘，實現相幹解調
demod_mult=dpsk.*f0;
%
設計接收端低通濾波器
fc=[ps 3.1*10^6];
%
過渡帶
mag=[1 0];
%
窗函數的理想濾波器幅度