非线性调频matlab仿真

close all;clear;clc;
B=2e6;
tau=（200）/B; %时宽带宽积
betaL=50/tau;
betaC=20/tau;
pi = 3.141592654;
j = sqrt（-1）;
N=1024;
t=linspace（-tau/2tau/2N）;
m=0;
%雷达参数

ft=（t./tau）.*（betaL+betaC*1./sqrt（（1-3.9868*t.^2./tau^2）））;
%非线性调频频率函数
ft_lfm=（t./tau）*B;
%线性调频频率函数
figure（1）
plot（t/tauft*tau‘b‘t/tauft_lfm*tau‘r‘）;

thetat=pi*betaL*t.^2/tau-pi*betaC*tau*sqrt（1-3.9868*t.^2/（tau^2））/2;
%非线性调频相位函数
thetat_lfm=pi*B*t.^2/tau;
%线性调频相位函数
figure（2）
plot（t/tauthetat‘b‘t/tauthetat_lfm‘r‘）;

xt=cos（thetat）+j*sin（thetat）;
xtfft=fft（xtN）;
%非线性调频频谱
xt_lfm=cos（thetat_lfm）+j*sin（thetat_lfm）;
xtfft_lfm=fft（xt_lfmN）;
%非线性调频频谱
figure（3）
plot（t/taufftshift（abs（xtfft））‘b‘t/taufftshift（abs（xtfft_lfm））‘r‘）;

htfft=conj（xtfft）;
Sot=fftshift（ifft（htfft.*xtfft））;
Z=abs（Sot）;
Z=Z/max（Z）;
Z=20*log10（Z+10e-6）;
%非线性匹配滤波器输出幅度
htfft_lfm=conj（xtfft_lfm）;
Sot_lfm=fftshift（ifft（htfft_lfm.*xtfft_lfm））;
Z_lfm=abs（Sot_lfm）;
Z_lfm=Z_lfm/max（Z_lfm）;
% Z_lfm=20*log10（Z_lfm+10e-6）;
Z_lfm=20*log10（Z_lfm）;
%线性匹配滤波器输出幅度
figure（4）
plot（t/tauZ‘b‘t/tauZ_lfm‘r‘）;

Dot=xt.*（cos（2*pi*7/tau*t）+j*sin（2*pi*7/tau*t））;
Dotfft=fft（Dot）;
Dot_nlfm=fftshift（ifft（htfft.*Dotfft））;
Dot_nlfm_am=abs（Dot_nlfm）;
Dot_nlfm_am=Dot_nlfm_am/max（Dot_nlfm_am）;
Dot_nlfm_am=20*log10（Dot_nlfm_am+10e-6）;
%非线性调频信号加频移
Dot_lfm=xt_lfm.*exp（j*2*pi*7/tau*t）;
Dot_lfmfft=fft（Dot_lfm）;
Dot_lfm=fftshift（ifft（htfft_lfm.*Dot_lfmfft））;
Dot_lfm_am=abs（Dot_lfm）;
Dot_lfm_am=Dot_lfm_am/max（Dot_lfm_am）;
Dot_lfm_am=20*log10（Dot_lfm_am+10e-6）;
figure（9）
plot（t