%rho_r=c/(2*Fr)而不是rho_r=c/(2*Bw)%Hsrcfexp函数里忘记乘pi了clcclearcloseall参数设置距离向参数设置R_eta_c=20e3;%景中心斜距Tr=2.5e-6;%发射脉冲时宽Kr=20e12;%距离向调频率alpha_os_r=1.2;%距离过采样率Nrg=320;%距离线采样数距离向参数计算Bw=abs(Kr)*Tr;%距离信号带宽Fr=alp

clcclearcloseall参数设置距离向参数设置R_eta_c=20e3;%景中心斜距Tr=2.5e-6;%发射脉冲时宽Kr=20e12;%距离向调频率alpha_os_r=1.2;%距离过采样率Nrg=320;%距离线采样数距离向参数计算Bw=abs(Kr)*Tr;%距离信号带宽Fr=alpha_os_r*Bw;%距离向采样率Nr=round(Fr*Tr);%距离采样点数（脉冲序列长度）方

clcclearcloseall参数设置距离向参数设置R_eta_c=20e3;%景中心斜距Tr=2.5e-6;%发射脉冲时宽Kr=20e12;%距离向调频率alpha_os_r=1.2;%距离过采样率Nrg=320;%距离线采样数距离向参数计算Bw=abs(Kr)*Tr;%距离信号带宽Fr=alpha_os_r*Bw;%距离向采样率Nr=round(Fr*Tr);%距离采样点数（脉冲序列长度）方

clcclearcloseall距离向参数R_eta_c=20e3;%景中心斜距Tr=25e-6;%发射脉冲时宽Kr=0.25e12;%距离向调频率Fr=7.5e6;%距离向采样率Nrg=256;%距离线采样点数Bw=abs(Kr*Tr);%距离信号带宽方位向参数Delta_f_dop=80;%多普勒带宽Fa=104;%方位向采样率Naz=256;%距离线数；theta_sq_c=[0,22.8]

clcclearcloseall距离向参数R_eta_c=20e3;%景中心斜距Tr=25e-6;%发射脉冲时宽Kr=0.25e12;%距离向调频率Fr=7.5e6;%距离向采样率Nrg=256;%距离线采样点数Bw=abs(Kr*Tr);%距离信号带宽方位向参数Delta_f_dop=80;%多普勒带宽Fa=104;%方位向采样率Naz=256;%距离线数；theta_sq_c=[0,22.8]

2022下半年到现在，一直在学习仿真《合成孔径雷达成像算法与实现》这本书，因为期末考试和其他的一些杂货，到现在应该是这本书的前八章都仿真出来了，但是我感觉有很多细节和知识点掌握的并不是很好，最近师兄让我复现一篇双基成像的论文的时候

clcclearcloseall距离向参数R_eta_c=20e3;%景中心斜距Tr=25e-6;%发射脉冲时宽Kr=0.25e12;%距离向调频率Fr=7.5e6;%距离向采样率Nrg=256;%距离线采样点数Bw=abs(Kr*Tr);%距离信号带宽方位向参数Delta_f_dop=80;%多普勒带宽Fa=104;%方位向采样率Naz=256;%距离线数；theta_sq_c=[0,22.8]

clcclearcloseall距离向参数R_eta_c=20e3;%景中心斜距Tr=25e-6;%发射脉冲时宽Kr=0.25e12;%距离向调频率Fr=7.5e6;%距离向采样率Nrg=256;%距离线采样点数Bw=abs(Kr*Tr);%距离信号带宽方位向参数Delta_f_dop=80;%多普勒带宽Fa=104;%方位向采样率Naz=256;%距离线数；theta_sq_c=[0,22.8]

《合成孔径雷达成像算法与实现》术语定义目标波束覆盖区星下点雷达轨迹速度方位向零多普勒面最短距离最近位置零多普勒时刻波束宽度目标轨迹信号空间和图像空间距离斜距平面地距斜视角距离横向目标是被SAR照射的地球表面上的一个假想点

clcclearcloseall%参数设置TBP=80;%时间带宽积T=10e-6;%脉冲持续时间N_ZD=60;%零频点位于中点右侧的距离，P58%参数计算B=TBP/T;%信号带宽K=B/T;%线性调频频率alpha_os=5;%过采样率F=alpha_os*B;%采样频率N=2*ceil(F*T/2);%采样点数dt=T/N;%采样时间间隔t_c=N_ZD*dt;%脉冲时间偏移%变量设置t=

clcclearcloseall%%参数设置Ta=128;%脉冲持续时间Ka=-0.085;%方位向调频率alpha_OS=[5,0.25];%过采样因子Delta_f_dop=abs(Ka*Ta);%方位信号带宽Fa=alpha_OS*Delta_f_dop;%方位向采样率N=2*ceil(Fa*Ta/2);%方位向采样点数dt=Ta./N;%采样时间间隔df=Fa./N;%采样频率间隔%%变量

clcclearcloseall%%参数设置Ta=128;%脉冲持续时间ka=-0.085;%方位向调频率alpha_OS=[5,0.25];%过采样因子Delta_f_dop=abs(ka*Ta);%方位信号带宽Fa=alpha_OS*Delta_f_dop;%方位向采样率N=2*ceil(Fa*Ta/2);%方位向采样点数dt=Ta./N;%采样时间间隔df=Fa./N;%采样频率间隔%%变量

clcclearcloseall%参数设置TBP=724;%时间带宽积T=42e-6;%脉冲持续时间t_0=1e-6;%脉冲回波时延Nfft=2^11;%fft长度%参数计算B=TBP/T;%信号带宽K=B/T;%线性调频频率alpha_os=1.07;%过采样率F=alpha_os*B;%采样频率N=2*ceil(F*T/2);%采样点数dt=T/N;%采样时间间隔df=F/N;%采样频率间隔%

clcclearcloseall%参数设置TBP=724;%时间带宽积T=42e-6;%脉冲持续时间t_0=1e-6;%脉冲回波时延Nfft=2^11;%fft长度%参数计算B=TBP/T;%信号带宽K=B/T;%线性调频频率alpha_os=1.07;%过采样率F=alpha_os*B;%采样频率N=2*ceil(F*T/2);%采样点数dt=T/N;%采样时间间隔df=F/N;%采样频率间隔%

代码复现如下：clcclearcloseall%参数设置TBP=100;%时间带宽积T=7.2e-6;%脉冲持续时间t_0=1e-6;%脉冲回波时延%参数计算B=TBP/T;%信号带宽K=B/T;%线性调频频率alpha_os=1.25;%过采样率F=alpha_os*B;%采样频率N=2*ceil(F*T/2);%采样点数dt=T/N;%采样时间间隔df=F/N;%采样频率间隔%变量设置t=-T

clcclearallcloseall%参数设置TBP=100;%时间带宽积T=10e-6;%脉冲持续时间%参数计算B=TBP/T;%信号带宽K=B/T;%信号调频频率alpha_os=8;%过采样率F=alpha_os*B;%采样率N=2*ceil(F*T/2);%采样点数dt=T/N;%采样时间间隔df=F/N;%采样频率间隔%变量设置t=-T/2:dt:T/2-dt;%时间变量f=-F/2:

与图3.7的代码区别只在于原始信号的表达式对了一个时间偏移代码复现如下：clcclearallcloseall%参数设置TBP=100;%时间带宽积T=10e-6;%脉冲持续时间tc=1e-6;%时间偏移量%参数计算B=TBP/T;%信号带宽K=B/T;%信号调频频率alpha_os=8;%过采样率F=alpha_os*B;%采样率N=2*ceil(F*T/2);%采样点数dt=T/N;%采样时间

代码复现如下：clcclearallcloseall%参数设置TBP=100;%时间带宽积T=10e-6;%脉冲持续时间%参数计算B=TBP/T;%信号带宽K=B/T;%信号调频频率alpha_os=8;%过采样率F=alpha_os*B;%采样率N=2*ceil(F*T/2);%采样点数dt=T/N;%采样时间间隔df=F/N;%采样频率间隔%变量设置t=-T/2:dt:T/2-dt;%时间变量

代码复现如下：clcclearallcloseall%参数设置TBP=100;%时间带宽积T=10e-6;%脉冲持续时间%参数计算B=TBP/T;%信号带宽K=B/T;%信号调频频率alpha_os=8;%过采样率F=alpha_os*B;%采样率N=2*ceil(F*T/2);%采样点数dt=T/N;%采样时间间隔df=F/N;%采样频率间隔%变量设置t=-T/2:dt:T/2-dt;%时间变量

代码复现如下：clcclearcloseall%参数设置TBP=100;%时间带宽积T=7.2e-6;%脉冲持续时间t_0=1e-6;%脉冲回波时延%参数计算B=TBP/T;%信号带宽K=B/T;%线性调频频率alpha_os=1.25;%过采样率，使用较高的过采样率是为了提高采样频率F=alpha_os*B;%采样频率N=2*ceil(F*T/2);%采样点数dt=T/N;%采样时间间隔df=F

代码复现如下：clccloseallclearallN=[4,8,16,8];n=256;n1=2^10;fori=1:4%t=-N(i)/2:1/n:N(i)/2-1/n;t=linspace(-N(i)/2,N(i)/2,n);y=sinc(t);title_1=[num2str(N(i)),'点sinc插值核'];ifi==4y=sinc(t).*kaiser(length(t),2.5)'

代码复现如下：clccloseallclearallN=500;t=-8:1/N:8-1/N;st=sin(pi*t)./(pi*t);figurep1=plot(t,st,'k');holdonp2=plot(t,st.*kaiser(length(t),2.5)','r-.');title('Kaiser窗加权后的Sinc函数，\beta=2.5'),xlabel('时间（采样点）'),yla

代码复现如下：clccloseallclearall%参数设置B=5.80e+6;%信号带宽T=7.26e-6;%脉冲持续时间K=B/T;%线性调频频率alpha=5;%过采样率F=alpha*B;%采样频率N=F*T;%采样点数dt=T/N;%采样时间间隔%变量设置t=-T/2:dt:T/2-dt;%时间变量%另一种变量设置方式%t=linspace(-T/2,T/2-dt,N);%信号表达st

%%插值核clccloseallclearallN=40;x1=-4:1/N:4-1/N;hx=sin(pi*x1)./(pi*x1);ifx1==0hx=1;endfigureplot(x1,hx,'k')axis([-44,-0.41.2])gridonarrow([-4,0],[4,0],'Color','k','Linewidth',1);arrow([0,-0.4],[0,1.2],'C

代码复现如下：x=linspace(-20,20,32);y=linspace(-20,20,32);SINC_1=zeros(length(x),length(y));fori=1:length(x)forj=1:length(y)SINC_1(i,j)=sinc(x(i))*sinc(y(j));endend%SINC_1=imrotate(SINC_1,8,'bilinear','crop'

本人也是初学《合成孔径雷达成像算法与实现》，在GitHub上下载的关于CS算法的代码，仔细读了读，发现中间有一些不合理之处，动手改了一改，重新写了一个代码供初学者学习，有什么问题留言即可。