应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】非线性方程组求解
51/1返回列表
查看: 3477  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

sunarising

铜虫 (小有名气)

[交流] 【求助】非线性方程组求解 已有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 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2010-07-05 10:38:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunarising

铜虫 (小有名气)
结果应该是递减的，好像还是有些问题。有时间我加你QQ向你请教吧。谢谢啦。
赞一下 回复此楼
3楼2010-07-05 18:30:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

anyuezhiji

银虫 (正式写手)

星空行者
sunarising(金币+30): 2010-07-06 08:41:52
引用回帖:
Originally posted by sunarising at 2010-07-05 18:30:49:
结果应该是递减的，好像还是有些问题。有时间我加你QQ向你请教吧。谢谢啦。

增加遗传代数试试

或许是精度还不够

继续算很多代以后看看结果是否合适


另外

f=f+((li(ii)^2+(x(ii)+u)^2)*(li(ii)-Fi(ii)/k)^2-li(ii)^2*l^2)^2*100;

f=f+sum((x+u)./li.*Fi)^2*100;

后面都乘了100

可以把现在的程序下求出的解代入这两组表达式
如果对这俩组约束方程的重视程度相当的话
应该在后面乘上合适的权重 使他们的值大致相当
当然也可以根据实际情况调整权重，侧重某组约束
赞一下 回复此楼
暗月下没有留下风的痕迹，但它已经寂然飘逝。。By<暗月之寂>:tiger38:
4楼2010-07-05 18:47:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:34 by ZvyPK8n6Nfki
[找工作] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:18 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-28 23:19 by b8fCuTHqckEt
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[公派出国] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:45 by ez6fGg9abYaj
[有机交流] 同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？ +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[找工作] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:28 by ez6fGg9abYaj
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-28 22:21 by ez6fGg9abYaj
[考博] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
[找工作] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 17:30 by ZYXdhzDAy9ZX
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意！ +7 zju2000 2026-09-26 8/400 2026-09-28 16:42 by beefly
[教师之家] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 15:32 by ZYXdhzDAy9ZX
信息提示
关闭
请填处理意见
关闭