单极性归零不归零

%单极性归零码srz
x=[1 1 1 0 0 0 0 0 0 0 0 1 1 0 0];
t0=200;
t=0:1/t0:length(x); %给出相应的时间序列
for i=1:length(x) %计算码元的值
if x(i)==1 %如果输入信息为1
for j=1:t0/2
y((2*i-2)*t0/2+j)=1; %定义前半段时间值为1
y((2*i-1)*t0/2+j)=0; %定义后半段时间值为0
end
else
for j=1:t0 %如果输入信息为0
y((i-1)*t0+j)=0; %定义所有时间值为0
end
end
end
y=[y,x(i)];
plot(t,y);
title('1 1 1 0 0 0 0 0 0 0 0 1 1 0 0'); grid on;
axis([0,i, -0.1,1.1]);
legend('单极性归零码')
xlabel('时间t/s')
ylabel('幅度A')
% 单极性不归零码snrz
function y=snrz(x)
x=[1 1 1 0 0 0 0 0 0 0 0 1 1 0 0];
grid=300;
t=0:1/grid:length(x);
for i=1:length(x) %计算码元的值
if(x(i)==1) %如果信息为1
for j=1:grid
y((i-1)*grid+j)=1;
end
else
for j=1:grid %反之，信息元为0，码元对应点值取0
y((i-1)*grid+j)=0;
end;end;
end;
y=[y,x(i)]; %为了画图，将y序列加上最后一位
M=max(y) ;
m=min(y);
plot(t,y);
axis([0,i,m-0.1,M+0.1]);
title('1 1 1 0 0 0 0 0 0 0 0 1 1 0 0');
%序列111000000001100对应的单极性不归零码波形
legend('单极性不归零码')
xlabel('时间t/s')
ylabel('幅度A')。