% 用于给阵列天线中的线阵的阵元激励进行泰勒加权
function [Itaylor,Imax]=taylor(SLL,n1,Nelement)
%--------------------------------泰勒线源----------------------------------%
% SLL---副瓣电平
% n1---用于控制副瓣数量
MM=Nelement; %阵元数
R0=10^(-SLL./20); %主副瓣电平比
A=(1/pi)*log(R0+sqrt(R0^2-1)); %中间因子 反双曲
minn=2*A.^2+0.5; %n1的最小值
sigma=n1/sqrt(A.^2+(n1-0.5).^2); %波瓣展宽因子
Z=zeros(1,n1-1);ss=zeros(1,n1-1);F=zeros(1,n1-1);
for n=1:1:n1-1
Z(n)=sigma*sqrt(A*A+(n-0.5)*(n-0.5));%确定零点位置
end
for n=1:1:n1-1
ss(n)=(factorial(n1-1))^2/(factorial(n1-1+n)*factorial(n1-1-n));%factorial阶乘
for m=1:1:n1-1
l(m)=1-(n/Z(m))^2;
end
zz=cumprod(l); %cumprod对数组计算累积连乘 得到与原数组维度相同的数组
F(n)=ss(n)*zz(n1-1);%F(m,A,n1)
end
if(mod(MM,2)==1)%阵元数为奇数
ii=(0:1:(MM-1)/2);p=2*ii/(MM-1);%在[0,L/2]的范围内均匀采样 pmax=x/(L/2) p=[0,1]
CC=2*(1:1:(MM+1)/2)/(MM+1); %归一化
for m=1:1:(MM+1)/2;
gx=0;
for n=1:1:n1-1;
U(n)=cos(pi*p(m)*n);
gx=2*F(n)*U(n)+gx;
Itaylor(m)=gx;
end
end
else %阵元数为偶数
ii=(1:1:MM/2);p=(2*ii-1)/(MM-1);%在[0,L/2]的范围内均匀采样 pmax=x/(L/2) p=[0,1]
CC=2*(1:1:(MM)/2)/MM; %归一化
for m=1:1:(MM)/2;
gx=0;
for n=1:1:n1-1;
U(n)=cos(pi*p(m)*n);
gx=2*F(n)*U(n)+gx;
Itaylor(m)=gx;
end
end
end
Itaylor=Itaylor+1; %F(0,A,n1)=1;此处的1为一个单位矩阵 到此I求解完毕
Imax=max(Itaylor);
end
- 1
- 2
- 3
前往页