| 查看: 3168 | 回复: 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 ] |
回复此楼» 猜你喜欢
一志愿985，本科211，0817化学工程与技术319求调剂
已经有4人回复
-
311求调剂
已经有6人回复
-
0703化学调剂 290分有科研经历，论文在投
已经有4人回复
-
化学工程321分求调剂
已经有6人回复
-
本人考085602 化学工程 专硕
已经有9人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有7人回复
-
326求调剂
已经有3人回复
-
焦虑
已经有8人回复
-
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
