readme.txt

TOYFDTD1 is a stripped-down minimalist, 3D FDTD code demonstrating the basic tasks in implementing a

本人要用FDTD计算一个波导缝隙天线的问题。一开始从网上下了一个叫ToyFDTD的开源的c程序做参考。它求解的是一个简单的矩形波导内的场分布，给出了一种最简单情况下求解FDTD问题的基本框架。然后我就根据这个例子编自己的问题，加上电壁和吸收边界条件。
现在讲一下我遇到的问题，很多，慢慢讲
1 激励是TE10，采用正弦波和高斯脉冲有什么区别呢？我看高本庆的书上讲用正弦的周期函数激励的话一般只要计算3个周期就可以稳定。但是正弦波也是一个随时间变化的场， 应该怎么样判断它已经稳定了呢？
2 激励的场强如果是1，那么由它产生的场的幅度是不是应该都不可能大于1呢？我是这么觉得的，但是ToyFDTD计算的结果却是大于1的？有没有人能解释一下呢？
3 关于吸收边界的问题。因为我这个是天线问题，要计算远区场，所以我的网格空间就不能只包括波导，要用吸收边界来截断。我用了2阶Mur条件，但是在棱边上不好处理，不知道有没有谁有好的方法
4 因为要算远区场，觉得近区远区的转换比较麻烦，有没有可以参考的程序？

下面是ToyFDTD的源程序，很简洁，注释也很详细。不过要注意下标，因为作者只是考虑理想电壁的情形。
但是从程序上看，他除了一个激励面之外，其余五个面都是用同样的方法处理的，感觉有点问题，因为在传播方向上不应该是电壁的。

希望大家和我一起讨论讨论。
 
1、正弦波的频谱是一个频点，而高斯脉冲的频谱仍然是满足高斯分布，是一个带谱。简谐波的情况一般比较简单，按照公式或者经验选取即可，或者设一个场观测点，采两个周期的场值比较一下就行了。
2、只要输出能量不大于输入能量幅度瞬时大于输入是可以出现的。
3、可以参看西电葛德彪的《电磁场时域有限差分方法》相关章节。
4、如果以后要在这个领域安身立命，还是自己编一编，对能力的锻炼是很重要的。
并没有看那个程序。不过按一般的算法，由于吸收边界也只是在无穷远才能将入射波衰减为0，而程序只能算有限区域，所以吸收边界的边界也是要截断的，最简单的处理就是设为电壁。不知道是不是和你问的是一个问题，如果说错了，不要怪我，呵呵。