%產生原碼
M=10;
t1=linspace（00.001*M60*M）;%時間取0-0.01s，共有600個點
A=rand（1M）;????????????????%取一行十個數，即a[10]
for?i=1:M
????if?A（i）>=0.5
????????A（i）=1;
????else?A（i）=0;????????????%大于0.5取1小于取0
????end
end
source=A;???????????????????%源碼取A的值



for?i=1:M
????if?source（i）==1
????????for?j=1:60
????????????sourcecode（j+60*（i-1））=1;%60個點一格代表一個0/1信號
????????end
????else?for?j=1:60
????????????sourcecode（j+60*（i-1））=0;%60個點一格代表一個0/1信號
????????end?
????end
end??????????????????????????????????%形成輸入波形



aftershift=[];??????????????????????%創建一個空矩陣
aftershift（1）=1;????????????????????%初始碼元為1
for?i=1:M???????????????????????????
????if?source（i）==1?????????????????
????????if?aftershift（i）==1
????????????aftershift（i+1）=0;?????
????????else?aftershift（i+1）=1;????
????????end
????else?aftershift（i+1）=aftershift（i）;?%兩碼元相同則是0%兩碼元不相同則是1
????end
end

t=linspace（00.001*（M+1）（M+1）*60）;
for?i=1:（M+1）
????if?aftershift（i）==1
????????for?j=1:60
????????????aftershiftcode（j+60*（i-1））=1;
????????end
????else?for?j=1:60
????????????aftershiftcode（j+60*（i-1））=-1;
????????end
????end
end??????????????????????????????????????%轉換后的碼形狀?0變-1

t2=linspace（00.00160）;
carrier=cos（2000*pi*t2+0.5*pi）;%載波
for?i=1:（M+1）
????for?j=1:60
????????encode（j+60*（i-1））=aftershiftcode（j+60*（i-1））*carrier（j）;????????%加載波后的波形
???????nencode（j+60*（i-1））=awgn（aftershiftcode（j+60*（i-1））*carrier（j）25）;%加高斯噪聲的波形
???????
??????%?h=rayleighchan（0.0010-1-3）;
??????%?nray（j+60*（i-1））=filter（hencode（j+60*（i-1）））;
????end
end???????????????????????????????????



for?i=1:M*60????????????????????
????D（i）=nencode（i+60）;%取600個點編碼加載波后波形形
end
for?i=1:60*M
???E（i）=nencode（i）;
end
for?i=1:M??
????F（i）=encode（i*60-45）;
end
for?i=1:M*60
???decode（i）=D（i）*E（i）;%差分相干解調
end

load?mypass?lowpasser;
LEN=fix（length（lowpasser）/2）;
%測出濾波器延時的點數?fix?向零靠攏取整數?length（）即lowapsser取點個數
decoded=[];
for?i=1:M*60
????decoded（i）=decode（i）;
end
for?i=（M*60+1）:（M*60+LEN）
????decoded（i）=0;
end
afterfilter=filter（lowpasser1decoded）;
%lowpasser/1為濾波器系數?
%比如說你的采樣頻率Fs=1000Hz，
%??選用8階的濾波器，期望的低通截止頻率是200Hz
%[ba]=butter（80.4‘low‘）;????%0.4=200/（Fs/2）
%y=filter（bax）???????????????%x是原始數據，y是濾波后的數據;

for?i=1:60*M
????if?afterfilter（i+LEN）<=0
????????recieve（i）=1;
????else?recieve（i）=0;
????end
end???????????????????????????????????%接收碼


for?i=1:M
????recieved（i）=recieve（（i-1）*60+30）;
end
source
aftershift;
recieved;
t3=linspace（00.001*（60*M+LEN/60）（60*M+LEN））;

n=0;
for?j=1:M
????for?i=1:60
????????if?recieve（60*j+i-60）~=sourcecode（60*j+i-60）
????????n=n+1;
????????end
????end
end
rate=n/（60*M）;%誤碼率

figure（1）
plot（t1sourcecode）axis（[0?0.001*M?-0.2?1.2]）title（‘原碼‘）xlabel（‘t（sec）‘）;
figure（2）
plot（taftershiftcode）axis（[0?0.001*（M+1）?-1.2?1.2]）title（‘碼形轉換‘）xlabel（‘t（sec）‘）;
figure（3）
plot（tnencode）axis（[0?0.001*（M+1）?-1.2?1.2]）title（‘DBPSK‘）xlabel（‘t（sec）‘）;
figure（4）
plot（t1decode）axis（[0?0

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????3561??2015-01-06?21:47??DBPSKn.m

?????文件????????873??2010-01-03?22:52??mypass.mat

-----------?---------??----------?-----??----

?????????????????4434????????????????????2