13位巴克码二相相位编码模糊图

close all;clear all;clc;
t=10e-6;          %线性调频信号的时间
fs=100e6;         %采样频率
fc=10e6;          %线性调频信号起始频率
B=2e6;            %调频带宽
ft=0:1/fs:t-1/fs;  %调制信号，也就是时间信号
N=length(ft);
k=B/fs*2*pi/max(ft);  %计算调制灵敏度，也就是理解为线性调频的步进系数
y=modulate(ft,fc,fs,'fm',k);
y_fft_result=fft(y);
figure,subplot(2,1,1),plot(ft,y);
subplot(2,1,2),plot((0:fs/N:fs/2-fs/N),abs(y_fft_result(1:N/2)));
 

D-83