[求助]请教前辈：问什么用FDTD方法计算出来的二维光子晶体的透射率会大于1？

我参照资料用FDTD方法改编的计算二维光子晶体的透射率程序，为什么出来的结果在有些频谱上透射率居然大于1，这是不可能的啊！！！ dH|^\IQ  
大家帮我分析分析这是什么情况 .. I&8m5F?$`  
Nl`ry2"<  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

光在光子晶体中传播时，沿传播方向（即非周期方向）的能流是守恒的。也就是坡印廷矢量是守恒的，排除数值误差，透射率一定小于1。看你的说法，应该是宽频脉冲激励，然后傅里叶变换求解频谱吧，要注意是否时间足够长，保证透射衰减的差不多了。而且不知道你的入射波如何引入的，如果是线源的话，因为会向两个方向辐射，所以要用两步法计算。总之理论上说，求解坡印廷矢量（注意是垂直于周期方向的）应该是守恒的。

非常感谢m4a1ken的回答，谢谢您。 g1ZV&X=2  
因为我的滤波器只是结构上的设置，并非有源器件，所以透射率一定小于1，这点非常赞同您的观点。 uR#'lb`3  
我加入的是点源（调制高斯脉冲激励源），是宽频频脉冲激励源。 IQ3n
@  
可是我认为不管时间是否够长，透射波都不应该会比入射波的值大啊，除非是真空，否则一定会有衰减吧。 .OmQ'  
希望可以得到您的指点，谢谢！

我把图附上。 n#BvW,6J  
86%%n?"}  
第一幅图入射波是正弦波，左图是入射波波形，右图是透射波波形。 G}VDEC  
第二幅图入射波是调制高斯脉冲，左图是入射波波形，右图是透射波波形。  `?|Rc  
GW3>&j_!d  
可以看到，如果入射波是正弦波，则透射波的（时域）幅值一直小于入射波。 ;1 |x  
可是当入射波是高斯调制脉冲时，则透射波的（时域）幅值就会有大于入射波的时候。 $2KK:{VX  
D+sQPymI  
还请好心人多多指教，谢谢！

晕死，图怎么传不上去啊！

SN7"7joP<  
N,8/Y  
图的顺序调换一下。

你的图看不到啊，是否上传的地址不许外链。 gt@m?w(  
我认为时域有透射大于入射的时候不能说明有问题，你要做傅里叶变换到频域，然后看频域的透射，就是能流的透射，每个频率下的透射应该是小于1的。能流守恒是针对每一个频率的，即每一个平面波。 '<"s \,  
而且用点源激励的话，透射的计算可能要在包围点源和光子晶体的一个闭合曲线上积分吧。 C{U?0!^  
如果可能，把你的计算方法和结构再说具体点。

还是看不到图啊，郁闷了。 ~~D{spMVO  
我昨晚又想了很久，对您的指点有些参透，您看对不对哈。因为脉冲信号持续的时间很短，之后的值几乎为0，可是经过滤波以后，出来的信号待稳定后基本上幅度变化的很小，但是值不等于0。而我又把所有运行时间的点都采了样，然后都做了FFT，所以出来透射率会有大于1的部分。 8{^kQ/]'|  
我觉得您说的对，首先应该看频域的透射率，其次我觉得是只看脉冲信号频率范围内的透射率。因为一旦超过脉冲信号频率范围，还是会有大于1的值，但是不能说明程序错误。可能是我的采样时间或者FFT的点数没有选好。关于这个问题我还想继续请教，呵呵，希望不会太打扰。 - YEZ]:"  
582660910这是我的QQ号，如果不介意，请加我吧！ q+yQwX{  
我想在QQ上应该说的清楚一些，至少可以看见图，呵呵！ ;AG8C#_  
谢谢分享您的宝贵经验，感激！

请问你能告诉我如何计算二维光子晶体的折射率嘛？我完全不懂，但是又急需哎。

我也遇到这种情况，是声子晶体的透射率大于1的情况。大致求解模型是这样的，二维声子晶体在中间，左边给一个高斯脉冲线激励源f1，在声子晶体的右边收集透射过来的信号f2，四个边界用PML吸收法，将fft(f2)/fft(f1),得透射系数。但有些频率信号的透射系数就是大于1，很难理解。求高手指点。