%例三 升余弦滾降系統

%本程序中時間單位是微秒
%頻率單位為MHz
%碼速率單位是Mb/s

global dt t df N
close all

N=2^13; %采樣點數
L=32;   %每碼元的采樣點數
M=N/L   %碼元數
Rb=2;		%碼速率是2Mb/s
Ts=1/Rb;        %碼元間隔
dt=Ts/L;        %時域采樣間隔
df=1/(N*dt)     %頻域采樣間隔
T=N*dt          %截短時間
Bs=N*df/2       %系統帶寬
Na=4;           %示波器掃描寬度為4個碼元
alpha=input('滾降系數=[0.5]');
if isempty(alpha),alpha=0.5;end

t=[-T/2+dt/2:dt:T/2]; %時域橫坐標
f=[-Bs+df/2:df:Bs];   %頻域橫坐標

g1=sin(pi*t/Ts)./(pi*t/Ts);
g2=cos(alpha*pi*t/Ts)./(1-(2*alpha*t/Ts).^2);
g=g1.*g2;		 %升余弦脈沖波形
G=t2f(g);

figure(1)
 set(1,'Position',[10,50,300,200])
                %設定窗口位置及大小 
figure(2)
 set(2,'Position',[400,50,300,200])
                %設定窗口位置及大小 
 %hold on
 grid 
 xlabel('t in us')
 ylabel('s(t) in V')
 
EP=zeros(size(f))+eps;
for ii=1:100
 a=sign(randn(1,M));

 imp=zeros(1,N);    %產生沖激序列
 imp(L/2:L:N)=a/dt;

 S=t2f(imp).*G;       %升余弦信號的傅氏變換
 s=f2t(t2f(imp).*G);  %升余弦信號的時域波形
 s=real(s);
 P=S.*conj(S)/T;      %升余弦信號的功率譜
 EP=(EP*(ii-1)+P+eps)/ii;

 figure(1)
 plot(f,30+10*log10(EP),'g');
 grid
 axis([-3,+3,-50,50])
 xlabel('f (MHz)')
 ylabel('Ps(f)  (dBm/MHz)')

figure(2)    %畫眼圖
hold on
 grid 
 tt=[0:dt:Na*L*dt];
 for jj=1:Na*L:N-Na*L
  plot(tt,s(jj:jj+Na*L));
 end

end