function [p]=syssim(snr_in_dB)
E=1;
SNR=exp(snr_in_dB*log(10)/10); %信噪比
sgma=E/sqrt(2*SNR);	 %噪聲的標(biāo)準(zhǔn)偏差
M=10000;  %碼元數(shù) 
L=32;  %每個(gè)碼元采樣點(diǎn)數(shù)
N=M*L;  %采樣點(diǎn)數(shù)
Rb=2;   %碼速率
Ts=1/Rb;  %碼元間隔
dt=Ts/L;  %時(shí)域采樣間隔
df=1/(N*dt);  %頻域采樣間隔
d=16;  %半個(gè)碼元采樣點(diǎn)數(shù)
T=N*dt;  %截短時(shí)寬
Bs=N*df/2;  %頻帶寬度
t=[-T/2+dt/2:dt:T/2];   %時(shí)域坐標(biāo)
t1=[-T/2+dt/2:32*dt:T/2];   %時(shí)域坐標(biāo)
f=[-Bs+df/2:df:Bs];   %頻域坐標(biāo)

alpha=0.5;	%滾降系數(shù)
g1=sin(pi*t/Ts)./(pi*t/Ts);
g2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).^2);
g=g1.*g2;
G=sqrt(t2f(g,dt));  %發(fā)送和接收濾波器頻譜
Cf=1-0.5*exp(-2*i*pi*f*d);  %	多徑信道
Gf=G.^2.*Cf;  %系統(tǒng)總的頻譜
%%  產(chǎn)生序列
a = sign(randn(1,M));
b = sign(randn(1,M));
%%  a路經(jīng)過發(fā)送端的成型濾波器和信道
imp=zeros(1,N);
imp(L/2:L:N)=a/dt;
S=t2f(imp,dt).*Gf;
s=f2t(S,dt);

s1=zeros(1,N);       %此段程序?qū)⒔?jīng)傅立葉變換得到的數(shù)字序列
s1(N/2+1:N)=s(1:N/2);%在中軸反轉(zhuǎn),如111-1-11,經(jīng)變換變成-1-11111
s1(1:N/2)=s(N/2+1:N);%因此應(yīng)把它反過來,這主要是由
s1=real(s1);         %f2t和t2f函數(shù)引起的.
gr=s1+sgma*randn(1,N);  %信號(hào)加噪聲，實(shí)現(xiàn)SNR
gr1=gr(16:32:N);  %取每個(gè)碼元中間的一個(gè)樣本
%%  b路經(jīng)過發(fā)送端的成型濾波器和信道
imp=zeros(1,N);
imp(L/2:L:N)=b/dt;
S=t2f(imp,dt).*Gf;
s=f2t(S,dt);

s1=zeros(1,N);       %此段程序?qū)⒔?jīng)傅立葉變換得到的數(shù)字序列
s1(N/2+1:N)=s(1:N/2);%在中軸反轉(zhuǎn),如111-1-11,經(jīng)變換變成-1-11111
s1(1:N/2)=s(N/2+1:N);%因此應(yīng)把它反過來,這主要是由
s1=real(s1);         %f2t和t2f函數(shù)引起的.
gr=s1+sgma*randn(1,N);  %信號(hào)加噪聲，實(shí)現(xiàn)SNR
gr2=gr(16:32:N);  %取每個(gè)碼元中間的一個(gè)樣本

%%  判決和誤碼概率的計(jì)算
numoferr=0;
for ii=1:M,
	%隨后為判決
    if (gr1(ii)<0),
		decis=-1;	 %判決為“0”
	else 
		decis=1;		%判決為“1”
	end;
	if(decis~=a(ii)),
		numoferr=numoferr+1;  %如果出錯(cuò)，計(jì)數(shù)器加1
	end;
	if (gr2(ii)<0),
		decis=-1;	 %判決為“0”
	else 
		decis=1;		%判決為“1”
	end;
	if(decis~=b(ii)),
		numoferr=numoferr+1;  %如果出錯(cuò)，計(jì)數(shù)器加1
	end;
end;
p=numoferr/2/M;		  %誤差概率估算