taylor_Taylor加权_线阵泰勒加权_泰勒线阵_泰勒加权_

% 用于给阵列天线中的线阵的阵元激励进行泰勒加权 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