有没有人用有限元算过空波导的本征值问题?
12-08
想问一下,均匀的空波导,
用顶点元方法算出来的本征值,
会不会有伪解?
我觉得应该不会有,但是算出来的值相差很大。
查程序查了很久,得出的矩阵应该不会错的。
用顶点元方法算出来的本征值,
会不会有伪解?
我觉得应该不会有,但是算出来的值相差很大。
查程序查了很久,得出的矩阵应该不会错的。
会不会分割的有问题
或是你的程序有小的致命的问题,比如某个系数错了
现在算的结果是,大部分都能和精确解对上,
但是多出了一些多余的解(大概每10个就有2,3个是多余的),
还有一些解是重复的,但却不是简并模式的,比如最低次模式
得到了两个一模一样的值,很奇怪。
为什么会多出那么许多东东呢?
而且前十个模式解中,丢掉了一个没有算出来。
你用什么算法求解的,有的算法,比如lanz,你选的参数的不正确,容易丢模
如果其他的模式求解的是对的,单单多出来几个莫名其妙的,应该是假模了八
lanz是什么算法?
你说的算法,是指求解矩阵本征值的算法吗?
我直接用了matlab里面的eig函数,
用Cholesky法。
不过好象用不用得到的结果都差不多。
哦,我用的是fortran,
在netlib上面找的一个剧烦的程序
matlab真好,呵呵
lanz 就是lanczos
你说的那个eeig我试过,好像用的是anordi(是不是这么拼呀?)
我感觉lanz是它的一个特列,有朋友用过,求介质谐振器的谐振频率
据说还行呀
估计还是你的程序问题
但是就是这个朋友提醒我假模问题,好像他是由这个问题的,过两天我闲下来,将它推荐的一个paper帖一下,你自己先找找
恩,好的,多谢了,呵呵。
不过我的程序从道理上说应该不会有假模啊,
因为是求解纵向场,在顶点上是连续的,
所以用顶点元是符合物理意义的。
难道还可能有别的原因导致假模?
heeh, 这个问题俺还不好回答,因为我也是球的纵向场
在欧的程序中从来没有出现过假模
等欧交了论文,偶也要研究一下这个问题,欢迎到时候继续讨论相关问题 ^_^
相关文章:
- 求助:本征值问题的有限元计算 (05-08)
- 请大侠们推荐一本关于有限元理论入门的中文书籍 (05-08)
- 请问有限元或者有限差分中这样一个问题!(05-08)
- 谁有有限元求解任意截面波导问题的代码?(05-08)
- 波导中的模式与反射的关系?(05-08)
- 请教如何用FDTD分析波导膜片 (05-08)
射频专业培训教程推荐