| 查看: 3438 | 回复: 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 ] |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有4人回复
-
怎么看青基中了没有啊
已经有6人回复
-
国社科又开始会评了，不知道这次命运如何
已经有11人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
基金未中，这种答复是模板吗？
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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人回复
4楼2010-07-05 18:47:12
3楼2010-07-05 18:30:49
投票:
10