线列阵加权程序

				function q=chbe(n,zpb)
				q=zeros(n);
				zp=10^(zpb/20);
				z0=[(2*zp)^(1/(n-1))+(1/(2*zp))^(1/(n-1))]/2;
				[M,k,y,x,r]=jiou(n);%调用子函数
				A=[];
				
				for m=k:-1:1
				    a=n-1;
				    b=fact(m-1);%调用子函数
				    c=fact(n-m-1);
				    s=[];
				    for p=0:(m-1)
				    d=fact(n-p-2);
				    e=fact(m-p);
				    f=fact(p);
				    g=fact(m-p-1);
				    h=a*b*d/(e*f*g*c)*z0^(n-(2*m+1))*(z0^2-1)^(m-p);
				    s=[s,h];
				end
				A=[A,sum(s)];%某一阵元权值
				end
				A(y)=z0^(2*M-1);%公式中需单独计算的值
				B=A(y:-1:r);%对称取值
				q=[[B,A]/A(x)]%归一化