未标题-2.jpg

无意中在网上得到的，因为没有权限进入讨论区，在这里发一下，供大家讨论讨论

下载下来看看

不错，奖励
用C语言计算MOM的解析公式，实部、虚步分开写的，用复数数据类型要方便一些

怎么用啊~~~:love :'(
好东西不会用~~

:31bb:31bb:31bb

本帖最后由 kerbcurb 于 2010-2-18 13:49 编辑

这个程序是算金属圆柱散射的吧，我感觉稍微有点问题，
角频率：w = 2*pi*f;
频率：f = 1.0e9;
光速：c = 3.0e8;
波长：lambda = c / f = 0.1m;
半径：a = 30 * lambda / ( 2 * pi ) = 3 / (2 * pi );
周长：L = 2 * pi * a = 3m;
分为Nu=100段，每一段的长度：l = 3 / 100 = 0.03 = 0.3 * lambda =  3 * lambda /10;也就是说每个波长取不到4个点，这样的话精度要降低不少。

另外，很多人喜欢以频率作为基本输入参量，这样不好，较好的办法是以波长作为输入基本参量，这样便于分析问题，实际上电磁学的很多公式里面，描述的几何尺寸前面都有一个系数k【=2pi / lambda】，实际上在计算之前可以对自变量做一个代换：kx -- > x，这样就会省去很多麻烦。

再者，看里面代码的注释是C++的，C++支持复数，所以代码还可以写的更精简一些，比如cgas的高斯消去法里面的矩阵可以定义为复数，汉克尔函数[贝塞尔函数：bs1、bs2]的返回值也可以是复数

看看别人做的，半径分别是1波长，2波长

对于时谐的散射问题，通常以波长为控制参数来计算

每个波长取样4个点，很少啊，可以先设置每波长的取样参数数，来控制离散数，这样避免改动方程尺寸时，引起较大的计算误差

MOM算射问题，每个波长20点，算RCS能看出一些误差，如果取到25～30个点，RCS基本上看不出来，有一位同学在群里问他的FDTD算得问什么差的那么多，我说是精度问题，他改变采样点数，他取到50～60个点，图形才与MOM算得差不多，但是误差还是明显能看出来的，我没算过FDTD，不知道FDTD算RCS每个波长需要取多少个点。

{:7_1234:}{:7_1234:}{:7_1234:}