rz.m

matlab 对通信原理的一个仿真

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];   %频域横坐标

EP=zeros(size(f))+eps;
for ii1=1:30
b=round(rand(1,M)); %产生随机码
s=zeros(1,N);
for ii=1:L/2;
s(ii+[0:M-1]*L)=b;
end
S=t2f(s); %S是s的傅氏变换
a=f2t(S); %a是S的傅氏反变换
a=abs(a);
P=S.*conj(S)/T; %功率谱
EP=(EP*(ii1-1)+P+eps)/ii1;

end
figure(1) %输出的功率画图
plot(f,10*log10(EP),'b');
grid
axis([-15,15,-100,20]);
xlabel('f (MHz)');ylabel('P (V/MHz)') ;title('功率谱图')
figure(2) %输出信号波形画图
plot(t,a,'b')
grid
axis([-3,3,-0.3,1.3]);
xlabel('t (us)');ylabel('s(t) (V)');title('单极性RZ的时域图')