clear?all;
SNR=0:2:30;
snr=10.^（SNR/10）;
sigma=sqrt（1./snr）;%噪聲與信號(hào)的幅度比值
K=200000;%未調(diào)制信號(hào)長(zhǎng)度，信號(hào)越長(zhǎng)，曲線越精確
data=randint（1K）;
%%%%%%%%%%%%%%%%%%%%設(shè)置參數(shù)%%%%%%%%%%%%%%%%%%%%%%%%
nt=2;%發(fā)射天線數(shù)
nr=2;%接收天線數(shù)
n=4;?%調(diào)制方式nPSK

%%%%%%%%%%%%%%%%%%%%NPSK信號(hào)調(diào)制%%%%%%%%%%%%%%%%%%%%
npskdata=npsk（ndata）;%調(diào)用調(diào)制方式子函數(shù)
x=npskdata;???????????%調(diào)制后信號(hào)

%%%%%%%%%%%%%%%%%%%Alamouti%%%%%%%%%%%%%%%%%%%%%%%%
for?k=1:length（SNR）
????error_symbol=0;%對(duì)每一個(gè)信噪比下的誤符號(hào)率清零
	for?i=1:length（npskdata）/2
????	x1=x（2*i-1）;
????????x2=x（2*i）;?????????%對(duì)信源進(jìn)行空時(shí)操作，每?jī)蓚€(gè)數(shù)據(jù)符號(hào)一起
????????X=[x1?-x2‘;x2?x1‘];%空時(shí)發(fā)射信號(hào)矩陣
		H=rey（ntnr）;??????%調(diào)用信道子函數(shù)
		noise=sigma（k）.*rey（ntnr）;
????????R=sqrt（1/2）.*H*X+noise;%接收信號(hào)矩陣
????????y1=0;y2=0;
		for?i=1:nr
		????y1=H（i1）‘*R（i1）+H（i2）*R（i2）‘+y1;
		????y2=H（i2）‘*R（i1）-H（i1）*R（i2）‘+y2;
		end??

	????%進(jìn)行最大似然檢測(cè)
	????DEC_X=judge（ny1y2）;

????????%統(tǒng)計(jì)符號(hào)錯(cuò)誤的總數(shù)
	????if?DEC_X（1）~=x1
	????error_symbol=error_symbol+1;
	????end
	????if?DEC_X（2）~=x2
	????error_symbol=error_symbol+1;
	????end
	end
?????SER（k）=error_symbol/length（npskdata）;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%畫圖%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
semilogy（SNRSER‘b*-‘）;
grid?on
legend（‘Alamouti?QPSK（2TX?2RX）‘）
title（‘QPSK調(diào)制下的Alamouti碼方案在慢瑞利衰減信道上的誤符號(hào)率性能‘）
xlabel（‘信噪比/dB‘）;
ylabel（‘誤符號(hào)率‘）;
hold?on

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

?????文件??????71168??2009-03-30?10:47??運(yùn)行結(jié)果.doc

?????文件????????449??2009-03-29?21:07??judge.m

?????文件????????304??2009-03-29?21:02??npsk.m

?????文件????????241??2009-03-29?21:00??rey.m

?????文件???????1461??2009-03-29?21:34??alamouti.m

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

????????????????73623????????????????????5