线性调频模糊函数代码

clearallclcclftaup=1;b=10;up_down=-1;x=lfm_ambg(taup,b,up_down);taux=-1.1*taup:0.01:1.1*taup;fdy=-b:0.01:b;figure(1)mesh(100*taux,fdy./10,x)xlabel('Delay-\mus')ylable('Doppler-MHz')zlable('|\chi(\tau,fd)|')title('模糊函数')figure(2)contour(100.*taux,fdy./10,x)xlabel('Delay-\mus')ylable('Doppler-MHz')title('模糊函数等高线')gridonN_fd_0=(length(fdy)+1)/2;x_tau=x(N_fd_0,:);function[x]=lfm_ambg[taup,b,up_down]eps=0.000001;i=0;mu=up_down*b/2./taupfortau=-1.1*taup:0.01:1.1*taup;i=i+1;j=0;forfd=-b:0.01:bj=j+1;val1=1-abs(tau)/taup;val2=pi*taup*(1.0-abs(tau)/taup);val3=(fd+mu*tau);val=val2*val3+eps;x(j,i)=abs(val1*sin(val)/val);endendclearallclcclftaup=1;b=10;up_down=-1;x=lfm_ambg(taup,b,up_down);taux=-1.1*taup:0.01:1.1*taup;fdy=-b:0.01:b;figure(1)mesh(100*taux,fdy./10,x)xlabel('Delay-\mus')ylabel('Doppler-MHz')zlable('|\chi(\tau,fd)|')title('模糊函数')figure(2)contour(100.*taux,fdy./10,x)xlabel('Delay-\mus')ylable('Doppler-MHz')title('模糊函数等高线')gridonN_fd_0=(length(fdy)+1)/2;x_tau=x(N_fd_0,:);function[x]=lfm_ambg[taup,b,up_down]eps=0.000001;i=0;mu=up_down*b/2./taupfortau=-1.1*taup:0.01:1.1*taup;i=i+1;j=0;forfd=-b:0.01:bj=j+1;val1=1-abs(tau)/taup;val2=pi*taup*(1.0-abs(tau)/taup);val3=(fd+mu*tau);val=val2*val3+eps;x(j,i)=abs(val1*sin(val)/val);endend舒文23:58:41clearallclcclftaup=1;b=10;up_down=-1;%x=lfm_ambg(taup,b,up_down);eps=0.000001;i=0;mu=up_down*b/2./taup;fortau=-1.1*taup:0.01:1.1*taup;i=i+1;j=0;forfd=-b:0.01:bj=j+1;val1=1-abs(tau)/taup;val2=pi*taup*(1.0-abs(tau)/taup);val3=(fd+mu*tau);val=val2*val3+eps;x(j,i)=abs(val1*sin(val)/val);endendtaux=-1.1*taup:0.01:1.1*taup;fdy=-b:0.01:b;figure(1)mesh(100*taux,fdy./10,x)xlabel('Delay-\mus')ylable('Doppler-MHz')zlable('|\chi(\tau,fd)|')title('模糊函数')figure(2)contour(100.*taux,fdy./10,x)xlabel('Delay-\mus')ylabel('Doppler-MHz')title('模糊函数等高线')gridonN_fd_0=(length(fdy)+1)/2;x_tau=x(N_fd_0,:);