24小时热门版块排行榜    

查看: 2067  |  回复: 19
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

yung1347

银虫 (初入文坛)

[求助] 求解非线性方程组 已有3人参与

论文需要求解下面方程组,希望各位能给小弟提供点解决方法和思路,最近卡在这里有点烦了。

求解非线性方程组


@月只蓝 发自小木虫IOS客户端
回复此楼
不畏将来,不念过往
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hmmwjs

木虫 (著名写手)

你猜

引用回帖:
10楼: Originally posted by hmmwjs at 2016-08-20 09:01:12
突然发现 fitness太大,于是重新修正了一下条件,得到如下结果,建议你试试,优化这个东东本身就不是精确解,估计你想找精确解也找不到。这个还是比较不错的结果【L1 = 0.01718092158585438        L2 = 0.01718248636476 ...

我用的是遗传算法进行最优化的,Function函数代码如下:
function sum = Fqiu( L )
%FQIU Summary of this function goes here
%   Detailed explanation goes here
aa = 0.7;
bb = 0.7;
r = 25000;
b = 3;
h = 2;
rr = 5;
sgm1 = 4e6;
sgm2 = 4e4;
%%
% syms L(1) L(2)
G1 = -aa^(-1/bb)*(r/sgm1)^((1-bb)/bb)*(L(2) - L(1))^(1/bb)+...
    h+rr+b-sqrt(rr^2-L(1)^2)-sqrt(rr^2-L(2)^2)-...
    1/(bb+1)*aa^(-1/bb)*sgm1^((bb-1)/bb)*r^(1/bb)*(L(2)-L(1))^((1+bb)/bb)+...
    (r*h+r*rr+r*b-r*sqrt(b^2-L(1)^2)-r*sqrt(rr^2-L(2)^2)-sgm2)*(L(2)-L(1));
G2 = -r*L(1)/2*sqrt(rr^2-L(1)^2)-r*rr^2/2*asin(L(1)/rr)+r*L(1)*(h+b+rr-sqrt(rr^2-L(2)^2))+...
    r*L(2)/2*sqrt(rr^2-L(2)^2)-r*rr^2/2*asin(L(2)/rr);
sum = G1.^2+G2.^2;
end

» 本帖已获得的红花(最新10朵)

亲身体验通常是成长所必须的
11楼2016-08-20 09:02:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答

yung1347

银虫 (初入文坛)

刚刚上传的图片有点问题,这张才是正确.
求解非线性方程组-1



发自小木虫IOS客户端
不畏将来,不念过往
2楼2016-08-19 23:32:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hmmwjs

木虫 (著名写手)

你猜

引用回帖:
2楼: Originally posted by yung1347 at 2016-08-19 23:32:15
刚刚上传的图片有点问题,这张才是正确.

这个用solve解不了吗,我没试过,提醒一下而已

发自小木虫Android客户端
亲身体验通常是成长所必须的
3楼2016-08-20 06:48:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hmmwjs

木虫 (著名写手)

你猜

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
2楼: Originally posted by yung1347 at 2016-08-19 23:32:15
刚刚上传的图片有点问题,这张才是正确.

刚在等车的过程中掏出笔记本计算了一下,发现不论是solve还是写Function后用fsolve(@Qiujie,[0,0])求解的过程,结果都是 L1=0;L2=0;
试验了一下,改变里面的值,求解的结果会变,由此基本可以肯定,L1=0,L2=0是该方程的解
求解非线性方程组-2
001.jpg


求解非线性方程组-3
002.jpg

» 本帖已获得的红花(最新10朵)

亲身体验通常是成长所必须的
4楼2016-08-20 07:37:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见