% gaussian_sinc.m - Gauusian windowed sinc waveform% written by Nadav Levanon on 4 June 2003% single bit shapenn=201;  nn2=(nn-1)/2; nn22=nn2/2;small=.00000001;nnn=-nn2:nn2;arg_bit=small+4*pi/nn*nnn;amp_bit=sin(arg_bit)./arg_bit; s_gauss=.7; gauss_weight=exp(-0.5*(nnn./(nn/4*s_gauss)).^2);ab_a=amp_bit.*gauss_weight;sig_type=input(' Barker13 =1, P4 =2, single bit =3,  = ? ');if sig_type==1    phase_vec=pi*[0 0 0 0 0 1 1 0 0 1 0 1 0];  % Barker 13elseif sig_type==2    mm=input(' No. of elements of P4 signal = ? ');    m=1:mm;    phase_vec=pi*(1/mm*(m-1).^2-(m-1)); % P4elseif sig_type==3    phase_vec=[0 0 0 0 0 0 0 0 0 0];endlb=length(phase_vec);vec_length=(lb+3)*nn22+1;ab=zeros(lb,vec_length);bpv=ones(1,nn); % bit phase vectorfor k=1:lb  ab(k,:)=[zeros(1,(k-1)*nn22), ab_a.*exp(j*phase_vec(k)*bpv) , zeros(1,(lb-k)*nn22)];endu_amp_complex=sum(ab);u_amp1=abs(u_amp_complex);u_phase1=1/pi*angle(u_amp_complex);t_axis=nnn/nn22;figure(10)plot( u_amp1,'k','linewidth',1.5)gridfigure(11), clf, hold offplot( t_axis,ab_a,'k','linewidth',2.5)hold onplot(t_axis,gauss_weight,'k--','linewidth',1.5)plot(t_axis,amp_bit,'k:','linewidth',1.5)gridaxis([ -inf inf -.3 1.1])xlabel('{\it\tau} / {\itt_b}')disp(' ')disp('  With ambfn7 ')disp('    use    "u_amp = u_amp1;"     and     "u_phase = u_phase1;" ')disp(' ')