1 简介
2 部分代码
function F=sinc_interp(sig_ar1,Kv_R,d_Kv)%sinc_interp complete standard sinc interplation% sig_ar1 : the matrix of raw data% Kv_R : the matrix of RCM% d_Kv: the resolution in interplation axis[Na,Nr]=size(sig_ar1);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%// STOLT MATRIX // %%%%%%%%%sinc_matrix=[-0.003 0.010 -0.024 0.062 0.993 -0.054 0.021 -0.009;-0.007 0.021 -0.049 0.131 0.973 -0.098 0.040 -0.017;-0.012 0.032 -0.075 0.207 0.941 -0.134 0.055 -0.023;-0.016 0.043 -0.101 0.287 0.896 -0.160 0.066 -0.027;-0.020 0.054 -0.125 0.371 0.841 -0.176 0.074 -0.030;-0.024 0.063 -0.147 0.457 0.776 -0.185 0.078 -0.031;-0.027 0.071 -0.165 0.542 0.703 -0.185 0.079 -0.031;-0.030 0.076 -0.178 0.625 0.625 -0.178 0.076 -0.030;-0.031 0.079 -0.185 0.703 0.542 -0.165 0.071 -0.027;-0.031 0.078 -0.185 0.776 0.457 -0.147 0.063 -0.024;-0.030 0.074 -0.176 0.841 0.371 -0.125 0.054 -0.020;-0.027 0.066 -0.160 0.896 0.287 -0.101 0.043 -0.016;-0.023 0.055 -0.134 0.941 0.207 -0.075 0.032 -0.012;-0.017 0.040 -0.098 0.973 0.131 -0.049 0.021 -0.007;-0.009 0.021 -0.054 0.993 0.062 -0.024 0.010 -0.003;-0.000 0.000 -0.000 1.000 0.000 -0.000 0.000 -0.000];interp_len=8;F=zeros(Na,Nr,'single'); % initialize F(kx,ky) array for interpolationF1=zeros(Na,Nr+interp_len,'single');down_r=interp_len/2+1;up_r=down_r+Nr-1;sig_ar2=zeros(Na,Nr+interp_len,'single');sig_ar2(:,down_r:up_r)=sig_ar1;Kv_R1=zeros(Na,Nr+interp_len,'single');Kv_R1(:,down_r:up_r)=Kv_R;%% % % % % % % % % % % % % % % % new interploation % % % % % % % % % % % % % % % % % % % % % % %for m=1:Na % for each k loopfor n=down_r:up_ricKX=ceil(Kv_R1(m,n)/d_Kv)+n; % closest grid point in KX domainif(~((icKX>=down_r)&&(icKX<=up_r)))continue;endinterp_n_stolt=(icKX-4):1:(icKX+3);cKX=(icKX-n)*d_Kv; % and its KX valueres_Kz=-(Kv_R1(m,n)-cKX);row_idx=round(res_Kz/d_Kv*16);if(row_idx == 0)F1(m,n)= sig_ar2(m,icKX);continue;endF1(m,n)=sum(sig_ar2(m,interp_n_stolt).*sinc_matrix(row_idx,:));endend% % % % % % % % % % % % % % % new interploation % % % % % % % % % % % % % % % % % % % % % % %F=F1(:,down_r:up_r);clear F1 Kv_R1 sig_ar2;
3 仿真结果
4 参考文献
[1]黄峣, 尹春雷, 程潇. 合成孔径雷达成像RD算法仿真实验[J]. 兵器装备工程学报, , 35(5):109-111.
博主简介：擅长智能优化算法、神经网络预测、信号处理、元胞自动机、图像处理、路径规划、无人机等多种领域的Matlab仿真，相关matlab代码问题可私信交流。
部分理论引用网络文献，若有侵权联系博主删除。
