ami.m

ami信号的波形图。包括时域波形

clear all
close all
k=input('取样点数=2^k,k=[14]');
if isempty(k),
	k=14;
end
N=2^k;	%采样点数
L=64;	%每码元采样点数
M=N/L;	%码元数
Rb=2;	%码元速率
Ts=1/Rb;	%码元间隔
dt=Ts/L;	%时域取样间隔
df=1/(N*dt);	%频域取样间隔
T=N*dt;	%时域截短时间
Bs=N*df/2;	%频域系统带宽
t=[-T/2+dt/2:dt:T/2];	%时域横坐标
f=[-Bs+df/2:df:Bs];		%频域横坐标  
%定义全零距阵，将来作为求功率的累加器
Ep=zeros(size(f));		
for ii=1:50
	n=1;	%赋初值
	while mod(n,2)~=0,	%保证1的个数为偶数个
		a=round(rand(1,M));
		n=length(find(a==1));
	end
	bb=find(a==1);	%找出a中1的位置
	%找出第奇数个1的位置
	cc=bb(1:2:length(bb));
	%将第奇数个位置的1变为-1	
	a(cc)=-ones(size(cc));	
	ami=zeros(L,M);
    for loop=1:L/2
    	ami(loop,:)=a;	%ami第loop行元素
    end
    ami=reshape(ami,1,N);
    AMI=fftshift(fft(ami)).*dt;
    P=AMI.*conj(AMI)/T;
    Ep=(Ep*(ii-1)+P)/ii;
end
subplot(2,1,2)
aa=30+10*log10(Ep+eps);
plot(f*Ts,aa,'r')
grid on
axis([0,3,-1,30])
xlabel('f*Tb')
ylabel('Ps(f)(dBm/MHz)')
title('AMI码功率谱密度')
subplot(2,1,1)
plot(t,ami,'b')
grid on
axis([-5,5,-1.5,1.5])
xlabel('t(us)')
ylabel('s(t)(v)')
title('AMI码时域图')