matlab如何实现单极性归零码型设计?

matlab如何实现单极性归零码型设计?

归零码
　　信号电平在一个码元之内都要恢复到零的编码方式，它包括曼彻斯特编码和差分曼彻斯特编码两种编码方式。　
　　return to zero code (RZ)
　　一种二进制信息的编码，用极性不同的脉冲分别表示二进制的“1”和“0”，在脉冲结束之后要维持一段时间的零电平。能够自同步，但信息密度低。
　　它是码元中间的信号回归到0电平,判断的方法是"从正电平到零电平的转换边表示的码元为0,而从负电平到零电平转换边表示码元1",使得这种编码为自定时的编码.
单极性归零码
　　单极性归零码（RZ）即是以高电平和零电平分别表示二进制码1 和0，而且在发送码1 时高电平在整个码元期间T 只持续一段时间τ，其余时间返回零电平．在单极性归零码中，τ／T 称为占空比．单极性归零码的主要优点是可以直接提取同步信号，因此单极性归零码常常用作其他码型提取同步信号时的过渡码型．也就是说其他适合信道传输但不能直接提取同步信号的码型，可先变换为单极性归零码，然后再提取同步信号．
(一)用到的函数: f2t函数和t2f函数。
function x=f2t(X)
global dt df t f T N
%x=f2t(X)
%x为时域的取样值矢量
%X为x的傅氏变换
%X与x长度相同并为2的整幂
%本函数需要一个全局变量dt(时域取样间隔)
X=[X(N/2+1:N),X(1:N/2)];
x=ifft(X)/dt;
function X=t2f(x)
global dt df N t f T
%X=t2f(x)
%x为时域的取样值矢量
%X为x的傅氏变换
%X与x长度相同,并为2的整幂。　
%本函数需要一个全局变量dt(时域取样间隔)
H=fft(x);
X=[H(N/2+1:N),H(1:N/2)]*dt;
(二),主程序
global dt t f df N T
close all
k=input('取样点数＝2^k, k取10左右');
if isempty(k), k=13; end
z=input('每个信号取样点数＝2^z,zif isempty(z), z=5; end
f0=input('f0=取1(kz)左右');
if isempty(f0), f0=1; end
N=2^k
L=2^z;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 %系统带宽
f=[-Bs+df/2:df:Bs]; %频域横坐标
t=[-T/2+dt/2:dt:T/2]; %时域横坐标
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;
figure(1) %输出的功率画图
plot(f,10*log10(EP),'b');
grid
axis([-15,15,-100,20]);
xlabel('f (KHz)');ylabel('P (V/KHz)') ;title('功率谱图')
figure(2) %输出信号波形画图
plot(t,a,'b')
grid
axis([-3,3,0,1.3]);
xlabel('t (ms)');ylabel('s(t) (V)');title('单极性RZ的时域图')
end