%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%?Simulation?program?to?realize?BPSK??system
%
%?pro?by?pengsuyu????stu.id=1015231037
%?參考文獻(xiàn)：<>?yamamura?&&?harada
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
close?all;
clear?all;
clc;
%%%%%%%%%開(kāi)始的參數(shù)設(shè)置%%%%%%%%%%%%%%%%
sr=256000.0;?%?碼率
ml=1;????????%?調(diào)制電平數(shù)
br=sr.*ml;???%?比特率，由于是BPSK，所以在這里比特率等于碼率
nd?=?1000;???%?碼數(shù)，即仿真一次的符號(hào)數(shù)目，考慮到數(shù)據(jù)量大浪費(fèi)時(shí)間，數(shù)據(jù)量小又不精確，所以選擇2000
snr=3;?????%?信噪比?單位DB。按照誤差曲線，測(cè)試時(shí)不要太大，否則，不易看到誤碼。要想容易看到誤碼，就增加碼數(shù)，比如去nd=10000，或者增加仿真次數(shù)
IPOINT=8;????%?上采樣速率，按照要求取8
irfn=21;?????%?濾波器抽頭數(shù)帶一定的時(shí)域均衡效果。???????
alfs=0.1;????%?滾降因子
%%%%%%%%%%%滾降濾波器的參數(shù)設(shè)置并產(chǎn)生該濾波器的時(shí)域沖擊響應(yīng)序列%%%%%%%%%%%%%%%
[xh]?=?hrollfcoef（irfnIPOINTsralfs1）;???%其他人寫(xiě)的滾降濾波器函數(shù)-這個(gè)濾波器函數(shù)在工程目錄下供調(diào)用，該函數(shù)最后一個(gè)參數(shù)設(shè)置是發(fā)送濾波，還是接收濾波。
[xh2]?=?hrollfcoef（irfnIPOINTsralfs0）;??%最后一個(gè)參數(shù)是1，則是發(fā)送濾波，0則是接收濾波
????
%%%%%%%%%%%計(jì)算部分%%%%%%%%%%%%%%
nloop=100;????%?仿真次數(shù)
noe?=?0;????%?發(fā)生錯(cuò)誤的比特?cái)?shù)
nod?=?0;????%?一共發(fā)送的比特?cái)?shù)

for?ii=1:nloop????
%%%%%%%%%%%產(chǎn)生數(shù)據(jù)?%%%%%%%%%%%%%%%%

	data=rand（1nd）>0.5;??%隨機(jī)產(chǎn)生0-1的比特?cái)?shù)據(jù)，并記錄在data中
????data1=data.*（-2）+1;??????%將0-1的比特轉(zhuǎn)換成-1?和?1的電平??其中，0對(duì)應(yīng)1??1對(duì)應(yīng)-1
	data2=?oversamp（?data1?nd??IPOINT）?;??????%其他人寫(xiě)的上采樣模塊，調(diào)用起來(lái)很方便???????????????????????????????????????????
???
	data3?=?conv（data2xh）;??%通過(guò)卷積實(shí)現(xiàn)滾降低通濾波
%%%%%%%%%%%%%加高斯噪聲%%%%%%%%%%%%
	
????data4=awgn（data3snr）;??????%用AWGN函數(shù)制作高斯噪聲。
	data5=conv（data4xh2）;??????%接收卷積濾波

	sampl=irfn*IPOINT+1;????????%下采樣
	data6?=?data5（sampl:IPOINT:IPOINT*nd+sampl-1）;????
????
%%%%%%%%%%%%解調(diào)%%%%%%%%%%%%%%
????demodata=data6?%%%%%%%%%%%誤碼率計(jì)算%%%%%%%%%
????noe2=sum（abs（data-demodata））;??%?計(jì)算不一樣的碼的個(gè)數(shù)，即誤碼個(gè)數(shù)
	nod2=length（data）;??????????????%計(jì)算發(fā)送的碼的個(gè)數(shù)?
	noe=noe+noe2;??????????????????%若仿真多次，則累加起來(lái)
	nod=nod+nod2;???????????????????%若仿真多次，則累加起來(lái)

%?	fprintf（‘%d\t%e\n‘iinoe2/nod2）;???%一次仿真中的誤碼率。
end???

%%%%%%%%%%%輸出誤碼率%%%%%%%%%%%%%%%%%%

ber?=?noe/nod;
fprintf（‘仿真總次數(shù):%d\nSNR:%d?dB\n誤碼數(shù)：%d\n發(fā)送碼的總個(gè)數(shù)：%d\n誤碼率：%e\n‘iisnrnoenodnoe/nod）;

%%%%%%%%%%%%畫(huà)出圖像%一共8幅圖%%%%%%%%%%%%
????figure（1）???????????????%圖1是調(diào)制前后的序列圖
????subplot（211）;
????stem（data‘.‘）;
????subplot（212）;
????stem（data1‘.‘）;?
????title（‘調(diào)制前后序列圖‘）;
????figure（2）;??????????????%圖2是調(diào)制后的信號(hào)的星座圖
????ss=data1+j*eps;?
????plot（ss‘black*‘）;?
????grid?on;?
????axis（[-2?2?-2?2]）;
????xlabel（‘I-channel‘）;?
????ylabel（‘Q-channel‘）;?
????legend（‘Signal?under?White?Gaussian?Noise?channel‘‘Signal?Constellation‘）;??
????title（‘調(diào)制后信號(hào)星座圖‘）;
????figure（3）;
????stem（data2‘.‘）;????????%圖3是上采樣后的序列
????figure（4）;
????plot（data3）;????????%圖4是滾降濾波之后的波形趨于平滑
????eyediagram（data32*IPOINT）;????%圖5是眼圖1--高斯噪聲之前
????title（‘高斯噪聲之前‘）;
????eyediagram（data416）;????%圖6是眼圖2--高斯噪聲之后
????title（‘高斯噪聲之后‘）;
????figure（7）;
????plot（data5）;????????%圖7是匹配濾波之后的波形
????figure（8）;
????ss1=data6+j*eps;?
????plot（ss1‘black*‘）;?
????grid?on;?
????xlabel（‘I-channel‘）;?
????ylabel（‘Q-channel‘）;?
????legend（‘Signal?under?White?Gaussian?Noise?channel‘‘Signal?Constellation‘）;??
????title（‘濾波抽取信號(hào)星座圖‘）;
%%%%%%%%

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----

?????文件????????562??2001-11-30?07:56??qpsk\compconv.m

?????文件????????684??2001-11-30?07:56??qpsk\compoversamp.m

?????文件???????1652??2001-11-27?01:37??qpsk\hrollfcoef.m

?????文件???????4548??2013-06-10?13:07??qpsk\qpsk.asv

?????文件???????3970??2013-06-11?10:19??qpsk\qpsk.m

?????文件????????721??2013-06-10?11:20??qpsk\qpskdemod.m

?????文件????????913??2013-06-10?11:19??qpsk\qpskmod.m

?????文件???????4684??2001-11-30?06:31??qpsk\qpsk_fading.m

?????文件????????131??2013-06-11?10:17??版本2改進(jìn).txt

?????文件?????765952??2013-06-11?10:15??仿真圖說(shuō)明.doc

?????文件?????311808??2013-06-11?10:11??結(jié)果分析.doc

?????文件???????3598??2013-06-10?13:02??bpsk\bpsk.asv

?????文件???????3636??2013-06-11?10:05??bpsk\bpsk.m

?????文件???????1652??2001-11-27?01:37??bpsk\hrollfcoef.m

?????文件????????497??2001-11-27?09:39??bpsk\oversamp.m

?????目錄??????????0??2013-06-10?13:28??qpsk

?????目錄??????????0??2013-06-10?13:25??bpsk

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

??????????????1105008????????????????????17