7位巴克码和13位巴克码

close all;clear all; clc; code7=[1,1,1,0,0,1,0]; code13=[1,1,1,1,1,-1,-1,1,1,-1,1,-1,1]; tao=0.5e-6; fc=10e6; fs=40e6; t_tao=0:1/fs:tao-1/fs; N=length(code7);%N=7 NN=length(code13);%NN=13 pha=0; t=0:1/fs:7*tao-1/fs; tt=0:1/fs:13*tao-1/fs; N2=length(t); NN2=length(tt); s=zeros(1,length(t)); for i=1:N if code7(i)==1 pha=1; else pha=0; end s(1,(i-1)*length(t_tao)+1:i*length(t_tao))=sin(2*pi*fc*t_tao+pha*pi); end ss=zeros(1,length(tt)); for i=1:NN if code13(i)==1 pha=1; else pha=0; end ss(1,(i-1)*length(t_tao)+1:i*length(t_tao))=sin(2*pi*fc*t_tao+pha*pi); end figure subplot(2,1,1) plot(t,s),xlabel('t(单位：秒)'),title('二相码（7位巴克码）'); subplot(2,1,2) plot(tt,ss),xlabel('t(单位：秒)'),title('二相码（13位巴克码）'); f=(-N2/2:N2/2-1)/N2*fs; y_fft=fft(s); ff=(-NN2/2:NN2/2-1)/NN2*fs; yy_fft=fft(ss); figure subplot(2,1,1) plot(f,abs(y_fft)) xlabel('频率f(单位:Hz)'),title('7位巴克码的频谱'); subplot(2,1,2) plot(ff,abs(yy_fft)) xlabel('频率f(单位:Hz)'),title('13位巴克码的频谱');