分享自己编写的一维有限元计算实例“金属衬底介质片对平面波的反射”代码

cem-uestc · 发表于 2010-6-9 23:09:00

你最后一层的反射系数是如何得到？

kerbcurb · 发表于 2010-6-9 23:13:19

本帖最后由 kerbcurb 于 2010-6-9 23:14 编辑

回复 31# cem-uestc 的帖子
最后的公式，就是计算介质区域以外的发射率，前面给出的Doc_1那个描述，最后一个公式，需要修改，不过，我自己的方法好像不需要改变那个，我的方法也是解析法

kerbcurb · 发表于 2010-6-13 11:20:46

本帖最后由 kerbcurb 于 2010-6-13 11:23 编辑

下面的是Hz模式所谓的计算结果，没有采取金的书上迭代求解R的方法，是求得全部的Am,Bm后，计算R。实际上，分为N段后，每一段有A,B两个未知数，总共2N个未知数，分为N段提供N + 1个匹配点，两个端点【边界点】提供两个方程，其余N - 1个点【连接点】提供2N - 2个方程，因此总共2N个方程，因此可以解出所有的Am,Bm。蓝色圆点是金的书上的数据，红色+号是用这种方法计算的结果，迭代方法速度要快一些，因为没有关于矩阵的计算，不过本方法速度也不慢