shengyuxian.m

利用matlab编的对升余弦系统进行的仿真

%实验二，升余弦眼图比较
%本程序中时间单位是微秒
%频率单位为MHz
%码速率单位是Mb/s

global dt t df N
close all

N=2^14; %采样点数
L=64;   %每码元的采样点数
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('滚降系数1=[0.05]');
if alpha==[],alpha=0.05;end
alpha1=input('滚降系数2=[0.5]');
if alpha1==[],alpha1=0.5;end
alpha2=input('滚降系数3=[1]');
if alpha2==[],alpha2=1;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);        % 升余弦脉冲信号的傅氏变换
g3=cos(alpha1*pi*t/Ts)./(1-(2*alpha1*t/Ts).^2);
g4=g1.*g3;
G4=t2f(g4);
g5=cos(alpha2*pi*t/Ts)./(1-(2*alpha2*t/Ts).^2);
g6=g1.*g5;
G6=t2f(g6);
%--------------------------------------
figure(1)
 set(1,'Position',[10,50,500,300])
                %设定窗口位置及大小 
figure(2)
 set(2,'Position',[10,300,200,200])
                %设定窗口位置及大小 
 hold on
 grid 
 xlabel('t in us')
 ylabel('s(t) in V')
 
 figure(3)
 set(3,'Position',[220,300,200,200])
                %设定窗口位置及大小 
 hold on
 grid 
 xlabel('t in us')
 ylabel('s(t) in V')
 
 figure(4)
 set(4,'Position',[420,300,200,200])
                %设定窗口位置及大小 
 hold on
 grid 
 xlabel('t in us')
 ylabel('s(t) in V')
 %--------------------------------------
EP=zeros(size(f))+eps;
EP1=zeros(size(f))+eps;
EP2=zeros(size(f))+eps;
for ii=1:60
 a=sign(randn(1,M))+1;%产生的信号幅度为0、2

 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;%累计平均
 %--------------------------------------------
 S1=t2f(imp).*G4;       %升余弦信号的傅氏变换
 s1=f2t(t2f(imp).*G4);  %升余弦信号的时域波形
 s1=real(s1);
 P1=S1.*conj(S1)/T;      %升余弦信号的功率谱
 EP1=(EP1*(ii-1)+P1+eps)/ii;%累计平均
 %-------------------------------------------
 S2=t2f(imp).*G6;       %升余弦信号的傅氏变换
 s2=f2t(t2f(imp).*G6);  %升余弦信号的时域波形
 s2=real(s2);
 P2=S2.*conj(S2)/T;      %升余弦信号的功率谱
 EP2=(EP2*(ii-1)+P2+eps)/ii;%累计平均
 %-------------------------------------------
 
 figure(1)
 subplot(1,3,1);plot(f,30+10*log10(EP),'g');
 title('滚降系数一的功率谱')
 grid
 axis([-3,+3,-50,50])
 xlabel('f (MHz)')
 ylabel('Ps(f)  (dBm/MHz)')
 
 figure(1)
 subplot(1,3,2);plot(f,30+10*log10(EP1),'g');
 title('滚降系数二的功率谱')
 grid
 axis([-3,+3,-50,50])
 xlabel('f (MHz)')
 ylabel('Ps(f)  (dBm/MHz)')
 
 figure(1)
 subplot(1,3,3);plot(f,30+10*log10(EP2),'g');
 title('滚降系数三的功率谱')
 grid
 axis([-3,+3,-50,50])
 xlabel('f (MHz)')
 ylabel('Ps(f)  (dBm/MHz)')
 
 figure(2)
 tt=[0:dt:Na*L*dt];
 for l=1:Na*L:N-Na*L
  plot(tt,s(l:l+Na*L));
 
 end
 title('滚降系数一的眼图')
 
 figure(3)
 tt1=[0:dt:Na*L*dt];
 for l1=1:Na*L:N-Na*L
  plot(tt1,s1(l1:l1+Na*L));

 end
 title('滚降系数二的眼图')
 
 figure(4)
 tt2=[0:dt:Na*L*dt];
 for l2=1:Na*L:N-Na*L
  plot(tt2,s2(l2:l2+Na*L));

 end
 title('滚降系数三的眼图')
 
end