不同波形的模糊函数

%abs：不确定图；
%abs.^2：模糊图；
%% 常规脉冲信号
if 0
	tp = 100e-6;
	fc = 10e9;
	N = 90;
	[t fd] = meshgrid（linspace（-1 1 N） linspace（-10 10 2*N））;
	A = abs（sin（pi*fd/tp*tp.*（1-abs（t*tp）/tp））/pi./fd/tp/tp./（1-abs（t*tp）/tp））...
		.*（1-abs（t*tp）/tp）;
	A_norm = A/max（max（A））;
	if 1
% 		figure;surf（t fd A_norm ‘EdgeColor‘ ‘none‘）;
		figure;mesh（t fd A_norm）;colormap jet;
		title（‘载频10GHz、脉宽100\mus的单脉冲‘）;
		xlabel（{‘归一化时延‘;‘（t_p的倍数）‘}）;
		ylabel（{‘归一化多普勒频移‘;‘（1/t_p的倍数）‘}）;
		zlabel（‘接收机幅度响应A（tF_D）‘）;
	end
	if 0
		figure;imagesc（A_norm）;colormap jet;
		set（gca ‘XTick‘ [1 N/2 N]）;
		set（gca ‘XTickLabel‘ {-1 0 1}）;
		set（gca ‘YTick‘ [1 N N*2]）;
		set（gca ‘YTickLabel‘ {-10 0 10}）;
		title（‘载频10GHz、脉宽100\mus的单脉冲‘）;
		xlabel（‘归一化时延（t_p的倍数）‘）;
		ylabel（‘归一化多普勒频移（1/t_p的倍数）‘）;
	end
	if 1
		figure;contour（t fd A_norm）;colormap jet;
		title（‘载频10GHz、脉宽100\mus的单脉冲‘）;
		xlabel（‘归一化时延（t_p的倍数）‘）;
		ylabel（‘归一化多普勒频移（1/t_p的倍数）‘）;
	end
	if 1
		[row col] = size（A_norm）;
		range_ambiguity = A_norm（row/2 :）;
		velocity_ambiguity = A_norm（: col/2）;
		figure;plot（t range_ambiguity ‘b‘）;grid on;
		title（‘载频10GHz、脉宽100\mus单脉冲的距离模糊图‘）;
		xlabel（‘归一化时延（t_p的倍数）‘）;
		ylabel（‘|\chi（\tau0）|‘）;
		figure;plot（fd velocity_ambiguity ‘b‘）;grid on;
		title（‘载频10GHz、脉宽100\mus单脉冲的速度模糊图‘）;
		xlabel（‘归一化多普勒频移（1/t_p的倍数）‘）;
		ylabel（‘|\chi（0f_d）|‘）;
	end
end
%% 线性调频信号
if 0
	tp = 10e-6;
	B = 1e6;
	fc = 10e9;
	k = B/tp;
	N = 90;
	[t fd] = meshgrid（linspace（-1 1 N） linspace（-20 20 2*N））;
	A = （1-abs（t*tp）/tp）.*abs（sin（pi*（fd/tp*tp+B*t*tp）.*（1-abs（t*tp）/tp））...
		/pi./（fd/tp*tp+B*t*tp）./（1-abs（t*tp）/tp））;
	A_norm = A/max（max（A））;
	if 1
		figure;mesh（t fd A_norm）;colormap jet;
		title（‘载频10GHz、脉宽10\mus、调制带宽1MHz的线性调频信号‘）;
		xlabel（{‘归一化时延‘;‘（t_p的倍数）‘}）;
		ylabel（{‘归一化多普勒频移‘;‘（1/t_p的倍数）‘}）;
		zlabel（‘接收机幅度响应A（tF_D）‘）;
	end
	if 0
		figure;imagesc（A_norm）;colormap jet;
		set（gca ‘XTick‘ [1 N/2 N]）;
		set（gca ‘XTickLabel‘ {-1 0 1}）;
		set（gca ‘YTick‘ [1 N N*2]）;
		set（gca ‘YTickLabel‘ {-20 0 20}）;
		title（‘载频10GHz、脉宽10\mus、调制带宽1MHz的线性调频信号‘）;
		xlabel（‘归一化时延（t_p的倍数）‘）;
		ylabel（‘归一化多普勒频移（1/t_p的倍数）‘）;
	end
	if 1
		figure;contour（t fd A_norm）;colormap jet;
		title（‘载频10GHz、脉宽10\mus、调制带宽1MHz的线性调频信号‘）;
		xlabel（‘归一化时延（t_p的倍数）‘）;
		ylabel（‘归一化多普勒频移（1/t_p的倍数）‘）;
	end
	if 1
		[row col] = size（A_norm）;
		range_ambiguity = A_norm（row/2 :）;
		velocity_ambiguity = A_norm（: col/2）;
		figure;plot（t range_ambiguity ‘b‘）;grid on;
		title（‘载频10GHz、脉宽10\mus、调制带宽1MHz的线性调频信号的距离模糊图‘）;
		xlabel（‘归一化时延（t_p的倍数）‘）;
		ylabel（‘|\chi（\tau0）|‘）;
		figure;plot（fd velocity_ambiguity ‘b‘）;grid on;
		title（‘载频10GHz、脉宽10\mus、调制带宽1MHz的线性调频信号的速度模糊图‘）;
		xlabel（‘归一化多普勒频移（1/t_p的倍数）‘）;
		ylabel（‘|\chi（0f_d）|‘）;
	end
end
%% 相干脉冲串
if 0
	clear xxxxxx;	%不clear的话结果会一直叠加，维度一直增大
	taup = 0.2;
	pri = 1;
	n = 5;
	if taup > pri/2
		disp ‘ERROR. Pulsewidth nust be less than the PRI/2.‘;
		return;
	end
	g