毕业设计 Matlab仿真……

毕业设计，很急很急，用Matlab仿真的
辐射场的公式是

 
公式里，k=16，a分别是1.24m和0.114m，V0=30伏特，r是一个任意值，是相对空间内任意一点P(r,theta,fai)如果不好弄就随便定个值，定成inf也行。theta1=theta2=pi/12和pi/4，如果可以最好能做theta1与theta2不等的，和例图类似，只是辐射方向最大值在仰角大约5°左右，theta1=pi/12,theta2=pi/4.
然后theta取值从0到2*pi,对式子求theta'的（theta'是点Q(a,theta',fai')处的方位角，圆锥单臂长度）积分从sita1到pi-sita2
我自己编的干脆出不来，我也不懂这个theta'和theta怎么区分……
然后是贝塞尔函数，

最后要仿真的就是这个有限长双锥天线方向性函数

仿真出来类似这样

这个图里theta a就是之前的theta1,theta b就是之前的theta2
拜托哪位会编程，能不能把程序用文本给我发到邮箱，必有重谢……555……
lava_1987@163.com