求教叉状全息图的画法！！急！！！

kilsyou

                               
登录/注册后可看大图

一般参考光选择传播方向与z轴成一定角度的
平面波u,其数学表达式如下:
u =exp(-ikxx-ikzz),                                                 (2)
在z =0平面我们可以得到两束光的相干图样:
I =1+E20+2E0cos(kxx-qθ),                                   (3)
由此干涉图样决定的全息光栅周期为Λ=2π/kx。
这时周期透过率函数T(x)可以用傅里叶级数表
示:
T(x) =∑A0exp(i2πnx/Λ),                                      (4)
式中A0=1Λ∫T(x)exp(-i2nπx/Λ)dx（积分上下限为Λ/2  -Λ/2）,是n级衍
射的振幅。
可把T(x)简化为闪耀型光栅的透过函数
T(x) =exp[imMod(x,Λ)],                                        (5)
m是相位调制的振幅,
Mod(a,b) = a-bint(a/b),
最后(4)式的强度分布简化为
I =Mod(qθ-kx,2π),                                                （6）
(6)式的强度分布图

                               
登录/注册后可看大图

            
（5）式透射函数图
各位大侠帮忙把（5）全息图程序做出来！！！！

[ 本帖最后由 kilsyou 于 2009-5-1 15:24 编辑 ]

kilsyou

谁有计算全息图的程序，发过来参考一下