clearclose?all
bit=1000;
n=16;
p=0.6;
signal=rand（1n）<=p;?%產(chǎn)生n位隨機(jī)二進(jìn)制信號(hào)
receive=0;
j=1;
while?j<=n???????????%差分編碼
????if?signal（j）~=receive（j）
????????bi=1;
????????
????else
????????bi=0;
????end
????receive=[receive?bi];
????j=j+1;
end
difference=receive（2:n+1）;?%差分碼，去掉參考位
m=0:1/bit:（n-1）/bit;
figure（1）
subplot（311）?????????????%畫出基帶原碼
stairs（msignal）axis（[0n/bit-0.51.5]）;
title（‘基帶原碼‘）xlabel（‘時(shí)間/秒‘）ylabel（‘幅值‘）;
grid?on
subplot（312）?????????????%畫出差分碼
stairs（mdifference）axis（[0n/bit-0.51.5]）;
title（‘差分碼‘）xlabel（‘時(shí)間/秒‘）ylabel（‘幅值‘）;?
grid?on
for?i=1:n;?????????????????%將差分碼變成雙極性差分碼biploar?????????????
????if?difference（i）==0
????????biploar（i）=-1;
????else
????????biploar（i）=1;
????end
end
t=linspace（0n/10003200）;?%載波時(shí)間切片
carrier=sin（2*pi*1500*t）;?%載波
biploar=repmat（biploar2001）;?
biploar=reshape（biploar1numel（biploar））;
modulate=biploar.*carrier;??%模擬調(diào)制，產(chǎn)生2DPSK調(diào)制波
subplot（313）plot（tmodulate）axis（[0n/bit-1.51.5]）;
title（‘2DPSK調(diào)制波‘）xlabel（‘時(shí)間/秒‘）ylabel（‘幅值‘）grid?on
figure（2）
%（（（（（（信道傳輸））））））））））））
channel=awgn（modulate10）;
subplot（311）plot（tchannel）
axis（[0n/bit-1.51.5]）;
title（‘信道輸出‘）xlabel（‘時(shí)間/秒‘）;ylabel（‘幅值‘）;grid?on
%?（（（（（（相干調(diào)制））））））））））））））））））））
coherent=channel.*carrier;
subplot（312）
plot（tcoherent）;
axis（[0n/bit-1.51.5]）;
title（‘相乘輸出‘）xlabel（‘時(shí)間/秒‘）;ylabel（‘幅值‘）;grid?on
%（（（（（（（低通橢圓濾波器）））））））））））））））））
fp=60;fs=150;Fs=4000;
rp=3;rs=60;?
wp=fp/Fs;
ws=fs/Fs;
[NWn]