exp1.m

一些比较基本的程序

%实验一:单极性归零(RZ)波形的功率谱
%本程序中时间单位是微秒
%频率单位为MHz
%码元速率单位是Mb/s
clear all;
close all
global dt df t f N
N=2^13;%采样点数
L=64;%每码元的采样点数
M=N/L;%码元数
Rb=2;
Ts=1/Rb;%元长度
dt=Ts/L;%频域采样间隔
df=1/(N*dt);%时域采样间隔 
T=N*dt;%截短时间
Bs=N*df/2 ;%系统带宽 
Rt=0.5;% 占空比
t=[-T/2+dt/2:dt:T/2];%时域横坐标
f=[-Bs+df/2:df:Bs];  %频域横坐标
figure(1)
 set(1,'position',[10,50,300,200])%设定窗口位置及大小 
figure(2)
 set(2,'position',[400,50,300,200])%设定窗口位置及大小 
 EP=zeros(size(f));
 for loop=1:50
    rz=zeros(1,N);
    a=round(rand(1,M));
    for ii=1:Rt*Ts/dt
         rz(ii+[0:M-1]*L)=a;
    end                           %产生单极性归零码
    RZ=t2f(rz);  
    P=RZ.*conj(RZ)./T;            %单极性归零码的功率谱    
    EP=(EP*(loop-1)+P)/loop;      %计算累计平均
 end
 
 figure(1)
 aa=30+10*log10(EP+eps);%加eps以避免除以零，取对数运算以对应单位
 plot(f,aa,'r');
 grid
 axis([-15,+15,-50,50])
 xlabel('f (MHz)')
 ylabel('Ps(f) (dBm/MHz)')
 
 figure(2)
 plot(t,rz,'r')
 grid
 axis([-2,+2,-1.5,1.5])
 xlabel('t (us)')
 ylabel('s(t) (V)')
 clear all;