CW HFM LFM宽带模糊函数仿真代码.m

close all; clear all; clc; 

T = 0.26;
fs = 20e3;
fc = 5e3;
B = fc/3.2;
t = 0:1/fs:T-1/fs;

type_index = 3; % Choosing signal: 1 for pcw 2 for lfm 3 for hfm 4 for hfm+pcw
type_name = {‘PCW‘‘LFM‘‘HFM‘‘HFM+PCW‘};

switch type_index
    case 1
    % 1. pcw
    x = sin（2*pi*fc*t）;
    
    case 2
    % 2. lfm
    k = B/T;
    f0 = fc-B/2;
    x = sin（2*pi*（f0*t+k/2*t.^2））;
    
    case 3
    % 3. hfm
    f0 = fc+B/2;
    beta = B/f0/（fc-B/2）/T;
%     f0=f0*（fc-B/2）/fc;
    x = sin（2*pi/beta*log（1+beta*fc*t））;
    
    case 4
    % 4. hfm + pcw
    f0 = fc+B/2;
    beta = B/f0/（fc-B/2）/T;
    x = sin（2*pi/beta*log（1+beta*f0*t）） + sin（2*pi*fc*t）;
end

re_fs = 0.9*fs:10:1.1*fs;
alpha = re_fs/fs;       % Doppler ratio alpha = 1-2*v/c
doppler = （1-alpha）*fc;	% Doppler = 2v/c*fc = （1-alpha）*fc

N_a = length（ resample（xfsmin（re_fs）） ）;
N = N_a + length（x）-1;
afmag = zeros（length（alpha）N）;

hh = waitbar（0 ‘Calculating the Wide-band Ambiguity Function Waiting please! ‘）;
for i = 1:length（alpha）
    x_alpha = resample（xfsre_fs（i））;
    x_temp = zeros（1length（x_alpha））;
    for j = 1:lengt