tvamb.m

黑白全电视信号仿真

clear
Fst=35000000;                 %Generate "thumbtack-shape" ambiguity function of "sync-plus-white" signal.
Fsf=0.00025;                  %using 1-D convolution method.
t=0:1/Fst:0.000064;
t1=ones(size(t));
i=0;
for k=-64:64;
     x1=0.325*square(2*pi*(15625)*t,20);
     y1=(0.425*t1)+x1;
     x2=0.125*square(2*pi*(15625)*(t-(1.55e-006)),7.34);
     y2=(0.125*t1)+x2;
     y3=y1+y2;                                      %"sync-plus-white" signal.
     x31=0.10075*square(2*pi*(31200)*(t-(0.5e-006)),70); 
     y4=0.10075*t1+x31;
     x32=0.0612*square(2*pi*(62500)*(t-(2.3e-006)),50); 
     y5=0.0612*t1+x32;
     x33=0.0455*square(2*pi*(156250)*(t-(3e-006)),60); 
     y6=0.0455*t1+x33;
     x34=0.0455*square(2*pi*(390000)*(t-(0.9e-006)),30); 
     y7=0.0455*t1+x34;
     x35=0.052*square(2*pi*(1562500)*(t-(0.4e-006)),69); 
     y8=0.052*t1+x35;
     x36=0.0195*square(2*pi*(3125000)*t,10); 
     y9=0.0195*t1+x36;
     y10=y4+y5+y6+y7+y8+y9;
     [b,a]=cheby1(5,0.5,0.343);
     y11=filter(b,a,y10);
     x4=0.5*square(2*pi*(15625)*(t-(13e-006)),79.3);
     y12=0.5*t1+x4;                                                             
     y13=y12.*y11;
     y=y3+y13;     
     i=i+1;
     y14=y.*exp(j*2*pi*(k/Fsf)*t);       %y14=y3.*exp(j*2*pi*(k/Fsf)*t);
     y15=conj(fliplr(y));                %y15=conj(fliplr(y3));       
     y16(i,:)=conv(y14,y15);           
  end
  mesh((abs(y16)))
  %subplot(2,1,1),plot((abs(y(128,:))))
  %subplot(2,1,2),plot((abs(y(:,192))))