fortran或者其他编程高手来帮我编编这个程序~

fortran或者其他编程高手来帮我编编这个程序~
例1：盲波束形成方向图性能.两信源的Doppler频率分别为40Hz和45Hz,DOA分别为20°和-20°.由于它们的Doppler频率为5Hz,小于fs/K(200/32≈6.3Hz)用传统傅氏变换难以分辨.采用式(9)的盲波束形成方法估计得到Doppler频率后,用式(8)计算盲波束形成的权系数Ws,opt(fd),图1给出了对应于第一个信源(fd1=40Hz,θ1=20°)的盲波束形成方向图.由图1可见,方向图在θ1方向上形成主瓣,而在θ2=-20°方向上形成零点,其方向图性能非常理想,其中阵元幅相误差方差为20%,进一步的实验是将θ2向θ1方向靠近,当两方向位于差不多一个波束宽度之内时,盲波束形成的方向图在θ2方向上仍能形成深零陷,但是θ1方向并不是它的峰值点,这是由于波束形成的分辨率的限制,使得方向图在形成零陷的同时将主瓣挤偏.
Ws,opt(fd)=R-1sγs(fd) (8)
{ dl,l=1,2,…,p}=arg maxfd{?ξ(fd)=γHs(fd)R-1sγs(fd)} (9)
线性约束最小方差(LCMV)波束形成是自适应波束形成中最为常用的方法之一,传统的LCMV波束形成器要求已知目标信号的导向矢量a(θ),在导向约束最小方差准则下的最优权矢量可写作：
Wopt=μR-1sa(θ0) (3)
式中μ为任一非零常数,Rs为阵列空域相关矩阵,实际中用阵列快拍矢量来估计,即 ,当a(θ0)未知时,上述自适应波束形成难以进行.阵列信号处理中,同一传感器在一短时内的特性可认为恒定,通常假设仅空域采样有误差,而时域采样较精确,勿需考虑误差,基本符合实际.下面给出不利用导向矢量而利用Doppler信息的目标函数实现盲ADBF,从而完成DOA-Doppler稳健估计.
设N维列矢量Ws是ADBF器的加权系数,则阵列波束形成的第k次输出为
y(k)=WHsXs(k),k=1,2,…,k (4)
如果适当调整权矢量Ws,使其与除a(θl)外的所有其它方向矢量a(θj),(j≠l),都正交,且与a(θl)的内积近似等于s-1l,则y(k)与exp〔j(k-1)ψ(fdl)〕接近.因此,对如下目标函数寻优
(5)
可以求出P个局部最小点,得到P个Doppler估计 dl,l=1,2,…,P和P个空域权矢量Ws,l(l=1,2,…,P),即P个波束,由此得到P个DOA的估计值 l(l=1,2,…,P).
为了求解式(5),展开目标函数得
Fs(Ws,fd)=WHsRsWs+K-WHsγs(fd)-γHs(fd)Ws (6)
式中Rs与式(3)中的相同,而N维列矢量
(7)
式中Δl(fd)=ψ(fd)-ψ(fdl),(l=1,2,…,p);sincΔ= ,当fd已知时,由式(7)可知γs(fd)就是对阵列各单元同时作Doppler滤波后的阵列矢量.由式(6),先固定fd,则式(5)的最优权矢量
Ws,opt(fd)=R-1sγs(fd) (8)
从式(5)可知,计算自适应波束形成的权矢量Ws时,并不像常规方法(见式(3))那样要求a(θ0)已知,而仅要求目标的Doppler已知(或用后面介绍的估计得到),因此它是一个盲波束形成器(Blind beamforming).
将式(8)代入式(6),得到Doppler的估计如下
{ dl,l=1,2,…,p}=arg maxfd{?ξ(fd)=γHs(fd)R-1sγs(fd)} (9)

呵呵,智能天线啊.
不过这个没给参数别人没法给你编啊,你可以找论文作者要程序,前提是人家愿意给你.
波束成型其实就是些向量相乘,程序不会太难,关键是看懂算法,搞清楚物理意义啊.