chbe.m

来自「线列阵加权程序」· M 代码 · 共 28 行

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)]%归一化