为什么整个区域最后变成了无数值（Not a number）？:我用FDTD计算耦合微带线，运行一段时间后，整个区域就变成了Not a number，如图所示，这是什么原因？
还有就是算出来的电场值为什么这么大？如表所示（
E_z(i,j,k=k_s_o_u_r_c_e)）
| -7E+24 | 1.4E+25 | -2.5E+25 | 4.21E+25 | -6.4E+25 | 8.8E+25 | -1.1E+26 | 1.27E+26 | -1.3E+26 | 1.27E+26 | -1.1E+26 | 8.78E+25 | -6.4E+25 | 4.2E+25 | -2.5E+25 | 1.4E+25 | -7.1E+24 | 3.27E+24 | -1.4E+24 | 5.41E+23 | -1.9E+23 | 6.37E+22 | -1.9E+22 | 5.37E+21 | -1.4E+21 | 3.26E+20 | -7.1E+19 | 1.43E+19 | -2.6E+18 | 4.5E+17 |
| 7.26E+24 | -1.5E+25 | 2.65E+25 | -4.4E+25 | 6.64E+25 | -9.2E+25 | 1.16E+26 | -1.3E+26 | 1.39E+26 | -1.3E+26 | 1.15E+26 | -9.1E+25 | 6.61E+25 | -4.4E+25 | 2.63E+25 | -1.4E+25 | 7.28E+24 | -3.4E+24 | 1.42E+24 | -5.5E+23 | 1.96E+23 | -6.4E+22 | 1.92E+22 | -5.3E+21 | 1.35E+21 | -3.2E+20 | 6.85E+19 | -1.4E+19 | 2.49E+18 | -4.2E+17 |
| -6.4E+24 | 1.29E+25 | -2.3E+25 | 3.9E+25 | -5.9E+25 | 8.13E+25 | -1E+26 | 1.17E+26 | -1.2E+26 | 1.17E+26 | -1E+26 | 8.04E+25 | -5.8E+25 | 3.83E+25 | -2.3E+25 | 1.26E+25 | -6.3E+24 | 2.91E+24 | -1.2E+24 | 4.7E+23 | -1.7E+23 | 5.39E+22 | -1.6E+22 | 4.41E+21 | -1.1E+21 | 2.58E+20 | -5.5E+19 | 1.08E+19 | -1.9E+18 | 3.22E+17 |
| 5E+24 | -1E+25 | 1.83E+25 | -3E+25 | 4.56E+25 | -6.2E+25 | 7.76E+25 | -8.9E+25 | 9.26E+25 | -8.8E+25 | 7.66E+25 | -6.1E+25 | 4.37E+25 | -2.9E+25 | 1.74E+25 | -9.5E+24 | 4.75E+24 | -2.2E+24 | 9.01E+23 | -3.4E+23 | 1.21E+23 | -3.9E+22 | 1.15E+22 | -3.1E+21 | 7.82E+20 | -1.8E+20 | 3.79E+19 | -7.3E+18 | 1.3E+18 | -2.1E+17 |
| -3.5E+24 | 7.05E+24 | -1.3E+25 | 2.1E+25 | -3.1E+25 | 4.1E+25 | -4.9E+25 | 5.53E+25 | -5.8E+25 | 5.47E+25 | -4.7E+25 | 3.75E+25 | -2.7E+25 | 1.84E+25 | -1.1E+25 | 6.06E+24 | -3E+24 | 1.33E+24 | -5.5E+23 | 2.09E+23 | -7.3E+22 | 2.32E+22 | -6.8E+21 | 1.86E+21 | -4.7E+20 | 1.07E+20 | -2.3E+19 | 4.35E+18 | -7.7E+17 | 1.24E+17 |
| 2.25E+24 | -4.5E+24 | 8.11E+24 | -1.3E+25 | 1.85E+25 | -2.3E+25 | 2.43E+25 | -2.6E+25 | 2.66E+25 | -2.5E+25 | 2.18E+25 | -1.7E+25 | 1.27E+25 | -9.2E+24 | 6.03E+24 | -3.1E+24 | 1.41E+24 | -6.2E+23 | 2.5E+23 | -9.4E+22 | 3.23E+22 | -1E+22 | 3.07E+21 | -8.9E+20 | 2.32E+20 | -5.5E+19 | 1.17E+19 | -2.3E+18 | 4.04E+17 | -6.5E+16 |
| -1.4E+24 | 2.71E+24 | -4.8E+24 | 7.62E+24 | -1.1E+25 | 1.28E+25 | -1.3E+25 | 1.41E+25 | -1.4E+25 | 1.35E+25 | -1.2E+25 | 9.26E+24 | -6.9E+24 | 5.03E+24 | -3.3E+24 | 1.71E+24 | -8E+23 | 3.45E+23 | -1.4E+23 | 5.12E+22 | -1.7E+22 | 5.51E+21 | -1.6E+21 | 4.6E+20 | -1.2E+20 | 2.81E+19 | -6E+18 | 1.18E+18 | -2.1E+17 | 3.45E+16 |
| 7.85E+23 | -1.6E+24 | 2.76E+24 | -4.3E+24 | 6.06E+24 | -7.4E+24 | 8.16E+24 | -8.7E+24 | 8.78E+24 | -8.3E+24 | 7.12E+24 | -5.7E+24 | 4.17E+24 | -2.9E+24 | 1.87E+24 | -1E+24 | 4.88E+23 | -2.1E+23 | 8.58E+22 | -3.2E+22 | 1.08E+22 | -3.4E+21 | 9.94E+20 | -2.7E+20 | 6.82E+19 | -1.6E+19 | 3.37E+18 | -6.6E+17 | 1.18E+17 | -1.9E+16 |
| -4.4E+23 | 8.77E+23 | -1.6E+24 | 2.47E+24 | -3.5E+24 | 4.41E+24 | -5.1E+24 | 5.52E+24 | -5.6E+24 | 5.28E+24 | -4.6E+24 | 3.6E+24 | -2.6E+24 | 1.79E+24 | -1.1E+24 | 6.11E+23 | -3E+23 | 1.33E+23 | -5.4E+22 | 2.01E+22 | -6.9E+21 | 2.16E+21 | -6.3E+20 | 1.68E+20 | -4.2E+19 | 9.47E+18 | -2E+18 | 3.84E+17 | -6.8E+16 | 1.1E+16 |
| 2.43E+23 | -4.8E+23 | 8.65E+23 | -1.4E+24 | 2E+24 | -2.6E+24 | 3.12E+24 | -3.5E+24 | 3.54E+24 | -3.3E+24 | 2.88E+24 | -2.3E+24 | 1.64E+24 | -1.1E+24 | 6.64E+23 | -3.6E+23 | 1.8E+23 | -8.1E+22 | 3.29E+22 | -1.2E+22 | 4.19E+21 | -1.3E+21 | 3.79E+20 | -1E+20 | 2.45E+19 | -5.5E+18 | 1.14E+18 | -2.2E+17 | 3.75E+16 | -6E+15 |
逐步看看是不是吸收边界问题，我以前也出现这样的情况，后来才发现是对边界网格中的E和H位置关系没弄清楚
逐步看看是不是吸收边界问题，我以前也出现这样的情况，后来才发现是对边界网格中的E和H位置关系没弄清楚
xugangthu 发表于 2009-5-12 12:46
请问你说的是前一个问题还是后一个？我用UPML对整个区域（包括主网格和边界）计算
明显的数值发散
检查公式
明显的数值发散
检查公式
cem-uestc 发表于 2009-5-12 15:24
刚才仔细检查了，没问题，我把我的程序发上来，麻烦你帮我看一下，到底哪儿出错了，我仿真的是耦合微带线
没有时间给你调试了
看学生有空没有，让一位学生跟你调试一下
没有时间给你调试了
看学生有空没有，让一位学生跟你调试一下
没有时间给你调试了
看学生有空没有，让一位学生跟你调试一下
看来学校的网速的确有问题
没有时间给你调试了
看学生有空没有，让一位学生跟你调试一下
cem-uestc 发表于 2009-5-12 22:23
谢谢啦。。。。。。。。。。。。。。。。。。
解决了没？
我计算1550nm银球的时候发散，最大值达到10的38次方
计算650nm银球的时候却不发散，能够得到稳定的值。
很奇怪
解决了没？
我计算1550nm银球的时候发散，最大值达到10的38次方
计算650nm银球的时候却不发散，能够得到稳定的值。
很奇怪
lh_lihong 发表于 2009-5-19 15:54
还没有完全解决，不过发现好像跟微带线模型中衬底的介电常数有关，介电常数不能太大，最大只能为1.35，再大刚开始运行是好的，运行足够长时间后得到的场值就会很大
衬底是什么物质，有电导率么？
还没有完全解决，不过发现好像跟微带线模型中衬底的介电常数有关，介电常数不能太大，最大只能为1.35，再大刚开始运行是好的，运行足够长时间后得到的场值就会很大
lyz4986 发表于 2009-5-19 16:08
很奇怪，算法与材料无关。
应该不是材料参数的影响
我有个教训
检查是不是被零除了
更正:不是零除的问题,看样子是发散
检查公式与代码是否一一吻合
衬底是什么物质，有电导率么？
lh_lihong 发表于 2009-5-19 16:20
我只设置了衬底的介电常数，电导率好像是认为是零
检查公式与代码是否一一吻合
cem-uestc 发表于 2009-5-20 21:35
我代码用的是本论坛下载的UPMLfdtd3D这个程序，在它的基础上加了微带线的设置，请问在主网格中，编程时，上面的空气部分和下面的衬底部分有什么区别？我个人认为仅介电常数不一样，用各自的介电常数，分界面上的介电常数用二者的平均值
更正:不是零除的问题,看样子是发散
yangcong 发表于 2009-5-20 21:35
请问哪些原因会导致发散呢？我的程序如附件
% PML边界最外层的PEC边界
% C1ez(1,:,:)=-1.0;
% C1ez(ih_tot,:,:)=-1.0;
% C2ez(1,:,:)=0.0;
% C2ez(ih_tot,:,:)=0.0;
% C3ey(1,:,:)=-1.0;
% C3ey(ih_tot,:,:)=-1.0;
% C4ey(1,:,:)=0.0;
% C4ey(ih_tot,:,:)=0.0;
这里看样子有问题
PEC边界有问题
介质地板接触PML边界没有？
对不起,我只是推测
但觉得好像在边界处有问题
介质地板接触PML边界没有？
cem-uestc 发表于 2009-5-24 18:21
恩，连在一起，这里会出什么问题吗？
对不起,我只是推测
但觉得好像在边界处有问题
yangcong 发表于 2009-5-24 18:26
可能是吧，我们导师也说可能是边界有问题
介质与PML连在一起的话，如果PML没有处理边界媒质信息，是要出问题的。
介质与PML连在一起的话，如果PML没有处理边界媒质信息，是要出问题的。
cem-uestc 发表于 2009-5-24 20:55
哦，我还想问一下微带线和接地板都看成是很薄的金属片，在编程时具体怎么处理？
就把金属面上切向的两个场分量设为零即可
一碰到边界就振荡,不知为什么
就把金属面上切向的两个场分量设为零即可
cem-uestc 发表于 2009-5-24 22:03
只需要设置面上的切向分量吗？面上的法向分量要不要也设为零？
shiyan1fujian.rar (3.39 KB) 
