rcos.m

通信系统仿真：单极性归零（RZ）波形及其功率谱

global dt t df N
close all
N=2^13; %采样点数
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=0.5;% 滚降系数为0.5
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',[350,50,300,200])%设定窗口位置及大小 
hold on
grid 
xlabel('t(us)')
ylabel('s(t)(V)')
 
EP=zeros(size(f))+eps;
for ii=1:50
a=sign(randn(1,M));
a=a+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;%累积平均

figure(2)          %作升余弦信号的眼图
tt=[0:dt:Na*L*dt];
for jj=1:Na*L:N-Na*L
plot(tt,s(jj:jj+Na*L));
title('升余弦滚降波形的眼图')
end
end                  
figure(1)           %作升余弦信号的功率谱
plot(f,30+10*log10(EP),'g');
grid
axis([-3,+3,-50,50])
xlabel('f (MHz)')
ylabel('Ps(f)  (dBm/MHz)')
title('升余弦滚降波形的功率谱')