frames.m

来自「对数字信号进行分析」· M 代码 · 共 30 行

function [frame1]=frames(system1,tb)
g1=[0,0,0,0,0,0,0,0,0,1];
N1=595;
frame2=zeros(1,N1);
for i=1:1:N1
    p1=xor(g1(1,3),g1(1,10));
    frame2(1,i)=g1(1,10);
     for k=10:-1:2
         g1(1,k)=g1(1,k-1);
     end
     g1(1,1)=p1;
 end
 runum1=length(frame2);
 for i=1:1:runum1
      if frame2(1,i)==0
           ww=-4.5;
      else
          ww=4.5;
      end
      frame2(1,i)=ww;
  end

runum2=length(system1);
frame1=zeros(1,runum1+runum2);
frame1(1,1:runum1)=frame2;
frame1(1,runum1+1:runum1+runum2)=system1;
t1=tb:tb:tb*(runum1+runum2);
figure
stairs(t1,frame1)