[求助]电磁场数值计算中点匹配法和muller迭代运算

如题 Hcwfe=K&/  
列出了相关的电磁场的场量方程，没办法求出解析解，想利用点匹配法和迭代运算求解，但又对这些一窍不通，尤其是要怎么利用软件求解，具体怎么入手 pm,xGo2  
拜 .. "GQ Q8rQ  
|5tZ*$nGa  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册

怎么没朋友出手相救啊？？？ E+@Q u "W  
急~~~~~

兄弟啊，你的问题题名太大，不好回答。我试试回答，看是不是你要的: E!9WZY  
:<J7g`f  
（1）学习计算机编程（如Matlab或C/C++）； {=Zy;Er  
（2）学习矩阵计算（如迭代松弛法求解）； A& iv  
（3）开始计算机软件编程。

天哪，这样的回答也太大了吧 _Sr}3  
$[|(&8+7  
具体的问题是要求半圆拱形波导内的传播常数。列出了波导内传播的场量分布，以及波导外介质的场量分布，根据波导的边界条件，即边界切线场量连续列出n个等式，现在不知道该怎么解这些等式？

这个好像是本征值问题？关于传播常数的方程是线性的，还是非线性的？

线性的 ~Jf(M^E  
现在关键是我不知道该怎么转换成矩阵的形式，方程里面还包括有贝塞尔函数什么的 =.s0"[%   
类似于这样的方程：

我在一些电磁场数值计算的书里看到过解类似的问题要用到各种迭代运算，例如muller迭代等等，但关键的问题是我现在不知道该怎么将列出来的方程转换成矩阵或是要求解的形式，编程的问题就更不知道了，愁~~~

贝塞尔和汉克尔函数有许多书上可以找到计算程序,如Numerical Recipes。 ~gQYgv<7  
4MzPm~Ct  
我的理解如下： a3A3mB
w  
（1）你的未知量是An，Bn，Cn，Dn（n=0，。。。，N-1），和gamma。 sk:B;.z  
（2）左手边第二项的无限求和必须截断成N-1。 /<-=1XJI  
（3）以A1，。。。，AN-1，B1，。。。，BN-1，C1，。。。，CN-1，D1，。。。，DN-1为未知量形成未知矢量，它们的系数就是矩阵元素。矩阵大小是4N x 4N。 0`v-pL0|  
（4）矩阵元素中包含有未知量gamma。 / $s(OFbi#  
（5）为了得到非零An，Bn，Cn，Dn （n=0，。。。，N-1），矩阵的行列式必须等于零，这就是你需要求解的gamma的方程，它是gamma的高阶非线性方程。

谢谢楼上的解答，好像明白了很多 0L,!o[L*  
R7!v=X]i  
既然要矩阵的行列式要等于零，应该就是高阶线性方程吧？就是说这个方程是以行列式等于零来列的？未知数是gamma？ M`@ASL:u  
那这样的方程能具体写出来么？具体在软件计算里应该怎么表示？

既然要矩阵的行列式要等于零，应该就是高阶线性方程吧？ $/5\Hg1  
VlSM/y5  
对，是一元高阶方程。 ^6F, lS_t  
S#HeOPRL  
`-R&4%t%  
就是说这个方程是以行列式等于零来列的？ aJF/y3  
:F#^Q%-IS  
是的。 ?9!9lSH6%  
fkImX:|q  
b!Nr  
未知数是gamma？ >j$a
Y  
E?bv<L,"  
是的。 5JK'2J&  
@p}_"BHYWt  
o:"(\$  
那这样的方程能具体写出来么？ Ex|Z@~T12  
-1NR]#P'  
不行，这样的方程无法具体写出来。 BafNFPc  
zxD=q5in  
6A|XB3  
具体在软件计算里应该怎么表示？ DHuvHK0#  
Li0+%ijM  
编一个计算行列式子程序。 "L1LL iS  
方程可以用求根算法（如New-Raphson迭代法）求解。