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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 用1stopt解非线性方程组问题
71/1返回列表
查看: 1742  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

people891108

木虫 (正式写手)

[求助] 用1stopt解非线性方程组问题 已有1人参与

本人刚开始接触1stopt,最开始用matlab求解结果不好，打算用1stopt求解非线性方程组，下面是我的程序，解的方法选用包维尔法，每次解一次方程组需要迭代很长时间，求大神帮我优化这个程序和算法，第二个方程0.7491=x3+x4+x5中0.7491是个变量，每次我都手动修改数值，希望大神能帮我将其他数值[0.5549; 0.6009; 0.7181; 0.7491; 1.0403]也能编进去，谢谢

程序
Title "concentration";
Parameters x(1:8)[0,];
//Variable ;
Function 1.678=x1+x2;
0.749151=x3+x4+x5;
50=x6+x7+x8+x4+x5;
0=x2+x7-x4-2*x5-x8;
3.82024*10^(-13)= x7*x8;
0.43996*10^(-6)=(x7*x4/(x3));
0.74128*10^(-10)=(x7*x5/(x4));
0.039833*10^(-6)=(x7*x1/(x2));
回复此楼

» 猜你喜欢

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

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
people891108: 金币+20, ★★★很有帮助, 非常感谢 2015-05-13 09:04:34
循环常数 a        x1        x2        x3        x4        x5        x6        x7        x8
0.5549        1.14173198863636        0.536268011363644        0.0224757584546091        0.528580867483308        0.00384337406208307        49.4675753252834        1.87075303955625E-8        4.14463700340161E-7
0.6009        1.10180642808357        0.576193571916431        0.027075916130922        0.571866314491743        0.00195776937733538        49.4257641557996        2.08306045808516E-8        0.000411739500621744
0.7181        0.999546948870126        0.678453051129874        0.0414689554283631        0.674809322011355        0.00182172256028892        49.3233686175059        2.70367977046525E-8        3.11034738591439E-7
0.7491        0.973269657335869        0.704730342664131        0.0459880433975497        0.701493763743082        0.00161819285935604        49.2968877925099        2.88425929463931E-8        2.22044929482337E-7
1.0403        0.743218579157401        0.934781420842599        0.106317357324403        0.933709132758929        0.000273509916853229        49.0654919888832        5.00963406693268E-8        0.000525318346303248
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

people891108
2楼2015-05-12 20:14:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

people891108

木虫 (正式写手)
送红花一朵
引用回帖:
2楼: Originally posted by dingd at 2015-05-12 20:14:30
循环常数 a        x1        x2        x3        x4        x5        x6        x7        x8
0.5549        1.14173198863636        0.536268011363644        0.0224757584546091        0.528580867483308        0.00384337406208307        49.4675753252834        1.87075303955625E-8        4.14463700340161E-7
0.6009 ...

请问能把程序代码发给我么，我用循环常数试了一下，结果只出一个解
赞一下 回复此楼
3楼2015-05-13 09:42:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)
引用回帖:
3楼: Originally posted by people891108 at 2015-05-13 09:42:26
请问能把程序代码发给我么，我用循环常数试了一下，结果只出一个解...

如果你用的是老版本就别试了，无法进行这种循环计算。
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

情弦风
4楼2015-05-13 09:51:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

情弦风

新虫 (小有名气)
送红花一朵
引用回帖:
4楼: Originally posted by dingd at 2015-05-13 09:51:11
如果你用的是老版本就别试了，无法进行这种循环计算。...

大神 我需要求解5600个循环的非线性方程组，可以把代码贴出来吗
赞一下 回复此楼
5楼2016-02-25 20:48:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

仰望星空7

新虫 (初入文坛)
function F=myfun(X)
T1=X(1);T2=X(2);T3=X(3);T4=X(4);Tf=X(5);q=X(6);
c1=0.003;
c2=0.003;
h=0.002;
w1=0.002;
w2=0.002;
v=5;
T0=3287;
%燃气对流
syms sigma At A Pr hc1
dt=0.013;d=0.03;Pc=3.5*10^6;c=2000;C0=9.18*10^5;Cp=1500;u=1.56*10^(-5);
m=0.022;k=1.26;Ma=3;
sigma=(T1*(1+(k-1)*Ma*Ma/2)/(2*T0)+0.5)^(-0.68)*(1+(k-1)*Ma*Ma/2)^(-0.12);
At=pi*dt*dt/4;A=pi*d*d/4;
Pr=u*Cp/m;
hc1=(0.026/dt^0.2)*(u^0.2*Cp/Pr^0.6)*(Pc/c)^0.8*(At/A)^0.9*sigma;
F(1)=hc1*(T0-T1)+C0-q;
%内壁导热
syms k2
m1=387;
k2=m1/c1;
F(2)=k2*(T1-T2)-q;
%冷却剂对流
syms Cpl ml lf lt pf pt uf ut Ref Prf hc20 hc2 yita
Cpl=4.1*Tf+1892.2;
ml=-0.0001*Tf+0.1383;lf=10^(-8)*(-5.3*Tf+683.7);lt=10^(-8)*(-5.3*T2+683.7);
pf=-0.8102*Tf+831.4661;pt=-0.8102*T2+831.4661;uf=pf*lf;ut=pt*lt;
de=2*h*w1/(h+w1);Ref=v*w1/lf;Prf=uf*Cpl/ml;
hc20=0.023*Ref^0.8*Prf^0.34*(uf/ut)*ml/de;
m2=202.4;m=(2*hc20/(m2*w2))^0.5;yita=tanh(m*h)/(m*h);
hc2=hc20*yita;
F(3)=hc2*(T2-Tf)-q;
F(4)=T2+T3-2*Tf;
%外壁导热
m3=387;
F(5)=(T3-T4)/(c2/m3)-q;
%空气对流
syms av Gra hc3
g=9.8;Ta=293;av=1/((T4+Ta)/2);d=0.2;pa=1.06;Pra=0.7;ma=0.04;la=4*10^(-5);
Gra=g*av*(T4-Ta)*d^3/la^2;hc3=0.48*(Gra*Pra)^0.25*ma/d;
F(6)=hc3*(T4-Ta)-q;
end
赞一下 回复此楼
6楼2016-03-10 15:22:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小双子蓝胖子

新虫 (初入文坛)
引用回帖:
4楼: Originally posted by dingd at 2015-05-13 09:51:11
如果你用的是老版本就别试了，无法进行这种循环计算。...

你好 有没有新的版本下载地址呀 求助谢谢 需要用到循环
赞一下 回复此楼
7楼2017-04-13 23:00:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 people891108 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭