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; 
% q3=quadl(M,-1,1); 
% 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; 
%   q1(n-109)=quad(A1,-1,1); 
%   A2=@(f)(sinc((La*B1)/(4*V).*(f-2.*m))).^4; 
%   q2(n-109)=quad(A2,-1,1); 
%   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('未加权','加权');