双极性归零码 matlab,双极性不归零码的Matlab实现

%双极性不归零码

%时间单位 us 频率单位 MHz 码元速率 Mb/s

close all %关闭所有的窗口

clear all %清除所有的变量

k=14;

N=2^k;  %总采样点数

L=256;  %每码元采样点数

M=N/L;  %码元数

Rb=2;  %码元速率

Ts=1/Rb; %码元间隔

dt=Ts/L; %时域分辨率dt(时域取样间隔)

df=1.0/(N*dt); %频域分辨率df，满足df*dt=1/N

T=N*dt;  %计算时间宽度T

Bs=N*df/2; %系统带宽

t=linspace(-T/2,T/2,N); %产生时间取样点矩阵

f=linspace(-Bs,Bs,N); %产生频域取样点矩阵

EEp=zeros(size(f)); %定义全零距阵，将来作为求功率的累加器

for ii=1:20

n=1;%赋初值

while rem(n,2)~=0;%保证仿真时有偶数个1

nrz=round(rand(1,M));%随机产生单极性归零码

n=length(find(nrz==1));

end

po0=find(nrz==0);%找出nrz中0的位置

nrz(po0)=-1;%将0赋为-1

copy0=ones(1,L);%定义复制的次数L

temp0=nrz(copy0,:); %将nrz第一行复制L次，生成L*M的矩阵

dnrz=reshape(temp0,1,L*M);%将temp0重排成1*L*M的矩阵

Z=t2f(dnrz,dt);%傅式变换

PP=Z.*conj(Z)/T;%求功率谱

EEp=(EEp*(ii-1)+PP)/ii;%求平均功率谱

end

%开窗口1

figure(1)

set(1,'position',[10,50,300,200])%设置窗口位置和大小

plot(t,dnrz)

title('双极性不零码')

axis([0,10,-1.5,1.5])%设置坐标范围

figure(2)

set(2,'position',[350,50,300,200])%设置窗口位置和大小

bb=30+10*log10(EEp+eps);

plot(f,bb)

grid on%显示网格

axis([-10,+10,-50,50])%设置坐标范围

最后修改于 2006-06-04 18:11

阅读(?)评论(0)