| 查看: 3109 | 回复: 3 | |||
[交流]
【求助】非线性方程组求解 已有2人参与
|
用matlab自带的fsolve函数和拟牛顿法等一系列优化算法求解结果都不理想。方程组形式和系数矩阵算法如下。请教如何求解这个复杂的非线性方程组。 系数矩阵算法如下： function g=g_matrix() global rate n topsize=rate*n*32e-9; basesize=n*32e-9; dtop=rate*32e-9; dbase=32e-9; toploc=[-topsize/2+dtop/2:dtop:topsize/2-dtop/2]; baseloc=[-basesize/2+dbase/2:dbase:basesize/2-dbase/2]; E1=2e4; a0=5e-9; xinf1=rate*2000*32e-9; xinf2=2000*32e-9; c1=-(xinf1+a0)*log((xinf1+a0)^2/(a0)^2)+(xinf1-a0)*log((xinf1-a0)^2/(a0)^2); c2=-(xinf2+a0)*log((xinf2+a0)^2/(a0)^2)+(xinf2-a0)*log((xinf2-a0)^2/(a0)^2); g=zeros(n,n); for i=1:n for j=1:n if (j==i) g(i,j)=-(2*a0*log(4)+c1)/(2*pi*E1*a0*dtop)-(2*a0*log(4)+c2)/(2*pi*E1*(1+0*baseloc(i))*a0*dbase); else g(i,j)=2*(log(abs(xinf1))-log(abs(toploc(i)-toploc(j))))/(pi*E1*dtop)+... 2*(log(abs(xinf2))-log(abs(baseloc(i)-baseloc(j))))/(pi*E1*dbase*(1+0*baseloc(i))); end end end function x=x_vector() global n rate topsize=rate*n*32e-9; basesize=n*32e-9; dtop=rate*32e-9; dbase=32e-9; toploc=[-topsize/2+dtop/2:dtop:topsize/2-dtop/2]; baseloc=[-basesize/2+dbase/2:dbase:basesize/2-dbase/2]; x=zeros(n,1); for i=1:n x(i)=toploc(i)-baseloc(i); x(i)=1e9*x(i); end 请教各位大虫帮帮忙。谢谢啦。 [ Last edited by sunarising on 2010-7-5 at 10:42 ] |
回复此楼» 猜你喜欢
自荐读博
已经有8人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
-
自然科学基金委宣布启动申请书“瘦身提质”行动
已经有4人回复
-
求个博导看看
已经有18人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Matlab求解偏微分方程组
已经有22人回复
- matlab求解非线性方程组
已经有16人回复
- 求高人指点用matlab求解非线性方程组，解决了追加100金币；
已经有11人回复
- matlab的fsove 命令求解非线性方程组
已经有6人回复
- 【求助】多元方程组求解及拟合
已经有9人回复
- 【求助】用mathematica 5.0求解一个非线性方程组失败，特发帖求助！
已经有5人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 非线性方程组的迭代法（数值计算高手请进）
已经有7人回复
- 【求助】求教matlab解非线性方程组
已经有9人回复
- 【求助】求解非线性方程
已经有8人回复
- 【分享】PDECOL-经典数值求解非线性偏微分方程组源程序
已经有142人回复
2楼2010-07-05 16:34:19
3楼2010-07-05 18:30:49
4楼2010-07-05 18:47:12