%?keystone變換，根據理想的頻域距離模型進行計算驗證。未產生信號收發流程。
clear?all;
close?all;

N?=?106;?%?快時間采樣點數
M?=?256;??%?慢時間采樣點數
L?=?N;
K?=?M;

B?=?100e6;?%?信號帶寬
c?=?3e8;?%?光速
fc?=?1e9;?%?中心頻率

eta?=?B/（fc*L）;

S_lm?=?zeros（LM）;?%?原始距離向頻譜
S_lk?=?zeros（LM）;?%?中間變量
S_lmp?=?zeros（LM）;?%?keystone變換后的距離向頻譜

l?=?（-L/2:L/2-1）‘;
m?=?（-M/2:M/2-1）‘;
%?m?=?（0:M-1）‘;

T?=?25e-6;?%?脈沖重復周期


vi?=?[1000];?%?目標速度
R0?=?[20];?%?初始位置
fdci?=?2*fc*vi/c;
bi?=?2*R0/c?*?B/L;

%?產生目標回波距離向頻譜
for?ll=-L/2:L/2-1
???for?mm=-M/2:M/2-1
???????S_lm（ll+1+L/2mm+1+M/2）?=?sum（exp（-1i*2*pi.*bi*（ll））?.*?exp（1i*2*pi*（1+eta*（ll））.*fdci*（mm）*T））+??randn;
???????S_lm（ll+1+L/2mm+1+M/2）?=?sum（exp（-1i*2*pi.*bi*（ll））?.*?exp（1i*2*pi*（1+eta*（ll））.*fdci*（mm）*T））;
???end
end



%?求?S（lk‘）
for?ll=-L/2:L/2-1
????for?kp=-M/2:M/2-1
????????S_lk（ll+1+L/2kp+1+M/2）?=?S_lm（ll+1+L/2:）*exp（-1i*2*pi*（1+eta*ll）/M*（kp）*m）;?%?
????end
end

%?求?S（lm‘）
%?S_lmp?=?fft（S_lk[]2）;

for?ll=-L/2:L/2-1
????for?mp=-M/2:M/2-1
????????kp?=?m;
????????S_lmp（ll+1+L/2mp+1+M/2）?=?S_lk（ll+1+L/2:）*exp（1i*2*pi/M*kp*（mp））;
????end
end

varef?=?M;
%?IFFT求距離像
s_lm?=?ifft（S_lm）;??%?fftshift
s_lmp?=?ifft（S_lmp）;
fig_position?=?[1050450310];????
set（gcf‘position‘fig_position）;
fig_fontsize?=?12;??%?字號

figure;
contour（abs（s_lm））??%?原始信號
hold?on;
figure;
contour（abs（s_lmp））?%?keystone變換后的信號

set（gcf‘Color‘‘none‘）??%?背景色
set（gca‘LineWidth‘0.5）;?%?坐標軸線寬?該寬度影響gird線寬
set（gca‘FontSize‘fig_fontsize‘FontName‘‘Times?New?Roman‘）;
%?h=legend（‘\itK\rm?=?1‘‘\itK\rm?=?2‘‘\itK\rm?=?4‘‘\itK\rm?=?8‘‘Location‘‘Best‘‘Orientation‘‘vertical‘）;?%?方向：vertical?或?horizontal
%?set（h?‘Box‘?‘off‘）????%?或者?legend（gca‘boxoff‘）??或者??legend（‘boxoff‘）
%?set（h?‘Color‘?‘none‘）?%?圖注背景透明
xlabel（‘\fontname{Times?New?Roman}\fontsize{12}雷達發射脈沖數?\itM‘）;
ylabel（‘\fontname{Times?New?Roman}\fontsize{12}距離分辨單元‘）;
set（gca‘MinorGridLinestyle‘‘--‘）;
set（gca‘GridLinestyle‘‘--‘）;??%?網格線型
%?set（gca‘XGrid‘‘off‘）;?%?x軸網格
%?set（gca‘YGrid‘‘off‘）;?%?y軸網格
set（gca‘XMinorGrid‘‘off‘）;?%?x軸子網格
set（gca‘YMinorGrid‘‘off‘）;?%?y軸子網格
set（gca‘XMinorTick‘‘on‘）;?%?x軸子網格標識
set（gca‘YMinorTick‘‘on‘）;?%?y軸子網格標識


s_lmp_FFT?=?fft（s_lmp[]2）;
figure;
contour（abs（s_lmp_FFT））
[vaktvidkt]=max（max（abs（s_lmp_FFT/M）））;

s_lm_FFT?=?fft（s_lm[]2）;
figure;
contour（abs（s_lm_FFT））
[vafftvidfft]=max（max（abs（s_lm_FFT）））;
figure;?
imagesc（abs（S_lk））
figure;?
imagesc（angle（S_lmp））

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????632??2018-07-11?07:33??gold_gener.m

?????文件???????3437??2018-09-17?15:03??keystoneSinc.m

?????文件???????3040??2018-07-27?20:21??Kstone10MCXFultSpectrum.m

?????文件???????2673??2018-08-06?16:02??ch4_keystone_contour.m

?????文件????????363??2018-07-26?15:02??gold_gen.m

-----------?---------??----------?-----??----

????????????????10145????????????????????5