%(6) 請用Matlab程序完成一個完整的數(shù)字基帶系統(tǒng)的仿真
%要求：
%1. 畫出輸入信號的碼形，并作出輸入信號的功率譜，眼圖。
%2. 畫出輸出信號的碼形，并作出輸出信號的功率譜，眼圖及誤碼率曲線。
%3.對數(shù)字基帶系統(tǒng)加入噪聲，進(jìn)行分析和探討。
global dt t f df N T %全局變量
close all
clear Eb_N0 Pe%清除以前對于信噪比和誤碼率的定義

N=2^13; 		%采樣點(diǎn)數(shù)
L=32; 			%每碼元的采樣點(diǎn)數(shù)  
M=N/L			%碼元數(shù)   
Rb=2;			%碼速率是2Mb/s		
Ts=1/Rb; 		%碼元間隔      
dt=Ts/L;		%時域采樣間隔					        
df=1/(N*dt)		%頻域采樣間隔     
T=N*dt 			%截短時間        
Bs=N*df/2 		%系統(tǒng)帶寬     
alpha=0.5		%滾降系數(shù)=0.5
t=linspace(-T/2,T/2,N);
%時域橫坐標(biāo)
f=linspace(-Bs,Bs,N)+eps; 
%頻域橫坐標(biāo)


figure(1)
set(1,'Position',[10,50,700,300])%設(shè)定窗口位置及大小

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(T2F(hr));%傅立葉變換取模，為開方
GT=sqrt(HR);%最佳基帶傳輸發(fā)送，接收濾波器的頻域應(yīng)當(dāng)共軛匹配且總體傳輸特性應(yīng)該為升余弦滾降函數(shù)
GR=GT;%發(fā)送接收濾波器系統(tǒng)函數(shù)匹配

for loop1=1:20%;循環(huán)20次
 Eb_N0(loop1)=(loop1-1) %信噪比
eb_n0(loop1)=10^(Eb_N0(loop1)/10)%分貝值變?yōu)檎嬷? Eb=1;	
 n0=Eb/eb_n0(loop1);%信道噪聲單邊功率譜密度
sita=n0*Bs;  %噪聲功率
 n_err=0;     %誤碼計數(shù)，不相同的符號數(shù)
for loop2=1:5  
  a=sign(randn(1,M));%發(fā)送碼元
  imp=zeros(1,N);    
  imp(L/2:L:N)=a/dt;%產(chǎn)生沖擊信號序列
  IMP=T2F(imp);
  n_ch=sqrt(sita)*randn(size(t)); %信道噪聲
  nr=real(F2T(T2F(n_ch).*GR)); %噪聲通過接收濾波器后時域信號，用頻域的方法產(chǎn)生   
  sr=real(f2t(IMP.*HR))+nr; %傳輸后的波形      
  y=sr(L/2:L:N);						%以L/4為起點(diǎn)，步長L,取樣點(diǎn)N  ;符合4/T的取樣偏差
    aa=sign(y);    %接收碼元    
  n_err=n_err+length(find(aa~=a)) ; %誤符號數(shù)為兩者不相同的符號數(shù)
end
Pe(loop1)=n_err/(M*loop2);%誤碼率Pe'=n_err/(M*loop2)
figure(1)

  semilogy(Eb_N0,Pe,'r');
   
  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')
legend('實(shí)際曲線','理論曲線')
title('最佳基帶系統(tǒng)的信噪比/誤碼率曲線(無取樣偏差)')
 
 end