exampl_2.m

AMI码的功率谱,通过仿真产生AMI码并且绘制其功率谱图。

%例二 AMI码的功率谱

%本程序中时间单位是微秒
%频率单位为MHz
%码元速率单位是Mb/s

global dt df t f N
close all

N=2^13;   %采样点数
L=64;    %每码元的采样点数
M=N/L     %码元数
Rb=2;
Ts=0.5		%码元宽度是0.5us
dt=Ts/L;        
df=1/(N*dt)     %MHz
T=N*dt          %截短时间
Bs=N*df/2       %系统带宽

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 ii=1:50
 n=1;  %每次有偶数个1
 while rem(n,2)~=0,
   a=round(rand(1,M));
   n=length(find(a==1));
 end
 bb=find(a==1);    %bb是a中1的位置
 cc=bb(1:2:length(bb)); 
    %cc是a中奇序数1的位置
 b=a;              
 b(cc)=-ones(size(cc)); %所有位于cc的1变成-1

 ami=zeros(L,M);
 for loop=1:L,ami(loop,:)=b; end
 ami=reshape(ami,1,N);
 
 AMI=t2f(ami);       %AMI信号的傅氏变换
 P=AMI.*conj(AMI)/T;        %AMI信号的功率谱
 EP=(EP*(ii-1)+P)/ii;       %累计平均

 figure(1)
 aa=30+10*log10(EP+eps);%加eps以避免除以零
 plot(f,aa,'r');
 grid
 axis([-5,+5,-50,50])
 xlabel('f  (MHz)')
 ylabel('Ps(f)   (dBm/MHz)')

 figure(2)
 plot(t,ami,'r')
 grid
 axis([-2,+2,-1.5,1.5])
 xlabel('t (us)')
 ylabel('s(t) (V)')

end