barkcode.m

bark码产生程序!!!!!!!!很有用

function BARK=barkcode(t,fc,fp,fs) %t调制时间；fc载频；fp重复频率；fs采样频率；
% close all;clc;
%t=0.00065,fc=10E5,fp=10E3,fs=10E8;
n = 0:1/fs:1/fp;
  w=1/fp*1/10;
  m=t/(1/fp);
code=[1,1,1,0,0,1,0];         %7位巴克码,脉冲宽度:7*码宽=7us.
tao=1E-6;       %码宽1us
t_tao=0:1/fs:tao-1/fs;
%t_tao=50E-6:1/fs:50E-6+tao-1/fs;
N=length(code);
%N=length(n);
pha=0;
T=0:1/fs:7*tao-1/fs;
%T=50E-6:1/fs:7*tao-1/fs+50E-6;
s=zeros(1,length(T));
for i=1:N
    if code(i)==1
        pha=pi;
    else pha=0;
    end
 s(1,(i-1)*length(t_tao)+1:i*length(t_tao))=cos(2*pi*fc*t_tao+pha);
 %s(1,(i-1)*length(n)+1:i*length(n))=cos(2*pi*fc*n+pha);
end
slp=(square(2*pi*fp*T,7)+1)/2;
sp2=slp.*s;
n1=[0:7000];
n2=[0:10E4];
x1n=sp2;
x2n=0;
nn=[min(min(n1),min(n2)):max(max(n1),max(n2))];
y1=zeros(1,length(nn));
y2=y1;
x1a=zeros(1,min(n1)-min(nn));
x1b=zeros(1,max(nn)-max(n1));
x2a=zeros(1,min(n2)-min(nn));
x2b=zeros(1,max(nn)-max(n2));
y1=[x1a x1n x1b];
y2=[x2a x2n x2b];
y=y1+y2;
BARK=repmat(y,1,m);
%figure;
%plot(T,s);
%subplot(311),plot(slp);%axis([0 0.00065 -1 1]);
%subplot(312),plot(y);
%subplot(313),plot(BARK);