www.pudn.com > targer_point.zip > hamming_vs_rectangle.m, change:2013-10-09,size:2076b

```%%汉明窗加权以后的AASR与未加权函数比较

clc;clear all;
%%恒值参数
cj=sqrt(-1);
C=3e8;  %光速
%%典型星载雷达参数
La=10;
V=7100;
B1=3500;

%%未加权AASR
for n=110:1:150
m=n/100;
for i=-10:0.01:10
j=int32(i.*100+1001);
g(n-109,j)=sinc((La*B1)/(4*V).*(i)).^2;
end
end

s1=zeros(41,1);
s2=zeros(41,1);
s3=zeros(41,1);
AASR1=zeros(41,1);
%%旁瓣求积分
for n=110:1:150
m=n/100;
for i=-1:0.01:1
j=int32(i.*100+1001);
s1(n-109,1)=0.01.*g(n-109,j+200.*m)+s1(n-109,1);
s2(n-109,1)=0.01.*g(n-109,j-200.*m)+s2(n-109,1);
s3(n-109,1)=0.01.*g(n-109,j)+s3(n-109,1);
AASR1(n-109,1)=20.*log10((s1(n-109,1)+s2(n-109,1))./s3(n-109,1));
end
end

%%未加权AASR另一种方法
% M=@(f)sinc((La*B1)/(4*V).*(f)).^4;
% q1=zeros(1,31);
% q2=zeros(1,31);
%  for n=110:1:140
%   m=n/100;
%   A1=@(f)(sinc((La*B1)/(4*V).*(f+2.*m))).^4;
%   A2=@(f)(sinc((La*B1)/(4*V).*(f-2.*m))).^4;
%   AASR1=20.*log10((q1+q2)./q3);
%  end

%%加权AASR
%%方向图加窗
for n=110:1:150
m=n/100;
for i=-10:0.01:10
j=int32(i.*100+1001);
y1(n-109,j)=sinc((La*B1)/(4*V).*(i)).^2;
end
y(n-109,:)=conv(y1(n-109,:),hamming(64)','same');
y(n-109,:)=y(n-109,:)./max(y(n-109,:));
end

sum1=zeros(41,1);
sum2=zeros(41,1);
sum3=zeros(41,1);
AASR2=zeros(41,1);
%%加窗后的旁瓣求积分
for n=110:1:150
m=n/100;
for i=-1:0.01:1
j=int32(i.*100+1001);
sum1(n-109,1)=0.01.*y(n-109,j+200.*m)+sum1(n-109,1);
sum2(n-109,1)=0.01.*y(n-109,j-200.*m)+sum2(n-109,1);
sum3(n-109,1)=0.01.*y(n-109,j)+sum3(n-109,1);
AASR2(n-109,1)=20.*log10((sum1(n-109,1)+sum2(n-109,1))./sum3(n-109,1));
end
end

t=1.1:0.01:1.5;
figure
plot(t,AASR1,t,AASR2)
grid on;
legend('未加权','加权');
xlabel('\fontsize{9}过采样率');
ylabel('\fontsize{9}AASR（db）');
title('\fontsize{9}旁瓣抑制');

figure
plot(-10:0.01:10,g(1,:),-10:0.01:10,y(1,:))
legend('未加权','加权');

```