close all;clear all;clc;
code=[1,1,1,0,1,0,1];%7位巴克碼
tao=10e-6;%脈沖寬度10us
fc=0.5e6;%調頻信號起始頻率
fs=10e6;%采樣頻率l00MHz
B=1e6;%調頻信號帶寬
t_tao=0:1/fs:tao-1/fs;
N=length(t_tao);
k=B/fs*2*pi/max(t_tao);
n=length(code);
pha=0;
s=zeros(1,n*N);
for i=1:n
        if code(i)==1
        pha=pi;
    else    pha=0;
    end       
       s(1,(i-1)*N+1:i*N)=cos(2*pi*fc*t_tao+k*cumsum(t_tao)+pha);
    end  
t=0:1/fs:7*tao-1/fs;
figure,subplot(2,1,1),plot(t,s),xlabel('t(單位：秒)'),
title('混合調制信號（7位巴克碼＋線性調頻');
s_fft_result=abs(fft(s(1:N)));
subplot(2,1,2),plot((0:fs/N:fs/2-fs/N),abs(s_fft_result(1:N/2))),
xlabel('頻率f(單位:Hz)'),title('碼內信號頻譜');