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