bpsk14.m

本课题为给语音识别系统提供语音信号端点检测的功能

%本程序中时间单位是微秒 
%频率单位为MHz 
%码速率单位是Mb/s 
 
global dt t f df N T 
close all 
clear Eb_N0 Pe 
 
N=2^15; %采样点数 
L=4;    %每码元的采样点数 
M=N/L   %码元数 
Rb=2;		%码速率是2Mb/s 
Ts=1/Rb;        %码元间隔 
dt=Ts/L;        %时域采样间隔 
df=1/(N*dt)     %频域采样间隔 
T=N*dt          %截短时间 
Bs=N*df/2       %系统带宽 
alpha=0.5 	%滚降系数 
 
t=linspace(-T/2,T/2,N); %时域横坐标 
f=linspace(-Bs,Bs,N)+eps;   %频域横坐标 
 
figure(1) 
 set(1,'Position',[60,50,700,200]) 
                %设定窗口位置及大小  
   
%升余弦 
hr1=sin(pi*t/Ts)./(pi*t/Ts); 
hr2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).^2); 
hr=hr1.*hr2; 
HR=abs(fft(hr));%取模是为了忽略时延 
GT=sqrt(HR); 
GR=GT; 
 
for loop1=1:10 
 Eb_N0(loop1)=(loop1-1)  %Eb/N0 in dB 
 eb_n0(loop1)=10^(Eb_N0(loop1)/10);%原始值 
 Eb=1;	 
 n0=Eb/eb_n0(loop1); %信道的噪声谱密度    
 sita=n0*Bs; %信道中噪声功率 
 
 for loop2=1:5   
  a=sign(randn(1,M)); 
  imp=zeros(1,N);    %产生冲激序列 
  imp(L/2:L:N)=a/dt; 
  IMP=fft(imp); 
   
  n_ch=sqrt(sita)*randn(size(t)); %信道噪声 
  nr=real(fft(ifft(n_ch).*GR));    %输出噪声 
 
  sr=real(fft(IMP.*HR))+nr;       %接收信号 
  y=sr(L/4:L:N);		  %取样延时Ts/4 
  aa=sign(y);       %判决                  
  n_err=length(find(aa~=a))   
 end 
 Pe(loop1)=n_err/(M*loop2) 
 figure(1)  
  semilogy(Eb_N0,Pe,'G'); 
  eb_n0=10.^(Eb_N0/10); 
  hold on 
  semilogy(Eb_N0,0.5*erfc(sqrt(eb_n0))); 
  axis([0,9,1e-4,1]) 
  xlabel('Eb/N0') 
  ylabel('Pe') 
 end 
 
%end