用matlab处理CST结果的问题（关于metamaterial): untitled.JPG

不清楚啊！来看看啊！

一下是我仿的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方面的问题。所以建议重发吧

软件结合是个不错的方式，貌似也是发展只需

这个问题一直都有，参考的资料可以看DRSmith组还有JAKong组的提取方法，分别做个试试。我没做过这个，但以前看到过类似的问题，有争论没结果。

回复 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应该没有影响的呀

有一个稳健性算法的，smith他们都指出这个问题了，设置合适的阻抗容差带基本可以解决这个问题@@
Robust method to retrieve the constitutive effective parameters of metamaterials @Xudong Chen.PHYSICAL REVIEW E.2004

这篇文章的结构参数我还没搞懂 是说正方形边长2.5mm 介质厚度0.25mm么？请教各位大大

我按照smith05年文章建模如图但是结果却跟文章结果大不一样

回复 16# liuguanglg 的帖子
您的偏振方向设置对了么，我用的CST算的，边界条件设置是E-wall and H-wall