一下是我仿的smith 2005年的文章,见附件,结果如图,也会有这种跳跃的情况,当时感觉是计算误差的问题,也很好判断,但算自己的结构发现很复杂啊。。。
本帖最后由 feymanloveyou 于 2010-1-23 18:48 编辑
回复 3# austin555
在用hfss提取s参数,用matlab处理s参数得到metamaterial有效介电常数和磁导率的过程中,我也出现过相位正反跳荡情况。问题好像出在两个复数根,选取哪个进行运算。
另外建议此贴发在metamaterial板块,谢谢
的S-parameter导入matlab算结构的参数, Can you send you code to me :
pkuwuli@163.com
回复 5# pkuwuli
程序不是我一个人写的,不好意思。建议您参考smith2005年的文章
回复 4# feymanloveyou
谢谢您的回复,问题应该是在求Z的时候开了根号,因为强行规定了Z的实部为正,当根号内的数值在二象限贴近x轴时候,它的根我们取在一象限贴近y轴,而如果根号内数值在三象限贴近x轴时候(注意和前一种情况比较,这个根号内的数值时连续的),而开出的根却因为要把Z实部定为正,这个根要在四象限贴近y轴,于是这两个开出的根在实部会连续,但虚部差很多。我的问题是,怎样解决这个情况呢?
Ps:怎么把帖子移动过去啊?还是重新发?
回复 7# austin555
这个问题显然是metamaterial方面的问题。所以建议重发吧
回复 3# austin555
阻抗实部跟折射率虚部都要大于零,另外关于仿真中相位的问题希望可以讨论下。QQ:401783556
回复 7# austin555
在1楼途中,我发现那个频段很高。所以在高频段请注意有效介质理论是否还成立。所以我们从3楼的图讨论解决方法
我的观点是物理意义+函数连续条件。
Smith那样表示z的开方是没有任何问题的。是我们在编程的时候出了问题。MATLAB在复数开根号时候,自动在复平面取了一个解,而真正的是有无穷多组解(2*pi*n*i)。所以我的观点是多取一些n,带到程序里面,多算几条曲线。再加以连续性条件和物理意义方面的综合考虑,几组曲线的综合就是真正提取的曲线了。
不知道我说清楚了没。
本帖最后由 austin555 于 2010-1-27 11:19 编辑
回复 12# feymanloveyou
我大概懂您的意思,但是在对z开方的时候并不是三角函数,为什么会有无穷多个解呢?二次方根不就只是两个解吗还有,我看一些资料时,注意到scale law在THz的低段的时候都还是有效的,也就是说在几十个GHz应该没有影响的呀
回复 16# liuguanglg 的帖子
您的偏振方向设置对了么,我用的CST算的,边界条件设置是E-wall and H-wall