ip_06_09.m

来自「这是通信系统仿真的第六章的例题的仿真程序代码！挺好的」· M 代码 · 共 29 行

% MATLAB script for Illustrative Problem 6.9.
echo on
d=[1 1 1 0 1 0 0 1 0 0 0 1];
p(1)=0;
for i=1:length(d)
  p(i+1)=rem(p(i)+d(i),2);
  echo off ;
end
echo on ;
a=2.*p-1;
b(1)=0;
dd(1)=0;
for i=1:length(d)
  b(i+1)=a(i+1)+a(i);
  d_out(i+1)=rem(b(i+1)/2+1,2);
  echo off ;
end
echo on ;
d_out=d_out(2:length(d)+1);
subplot(2,1,1);
n=1:length(d)
stairs(n,d)

axis([-0.1 length(d)+1 -0.1 1.1])
title('d')
subplot(2,1,2);
stairs(n,d_out)
axis([-0.1 length(d)+1 -0.1 1.1])
title('d-out')