24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

people891108

木虫 (正式写手)

应助: 0 (幼儿园)
金币: 4105.6
散金: 58
红花: 7
帖子: 468
在线: 193.2小时
虫号: 2122748
注册: 2012-11-12
专业: 化学反应工程

[求助] 用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人回复
0702一志愿吉大B区求调剂本科期间发表一篇Sci 已经有3人回复
物理学I论文润色/翻译怎么收费? 已经有203人回复
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？已经有2人回复
B区学生调剂-兰州交通大学材料科学与工程学院已经有11人回复
山西大同大学物理学专业还有调剂名额，欢迎调剂！已经有13人回复
桂林理工大学物理学专业招收调剂已经有13人回复
VASP 的一组 GPU / CPU 基准测试记录已经有0人回复
美国顶级物理期刊《应用物理快报》（APL）的编辑欺骗和歧视作者及AIP的官僚主义傲慢已经有6人回复
津理工大学晶体材料全国重点实验室刘红军教授课题组招收博士生一名已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

用1stopt自定义函数来拟合非线性曲线，函数中有复数，为什么得不到结果啊？？？已经有5人回复
用1stOpt解出非线性方程组的解作为初值用Matlab求解，解不出来已经有4人回复
看了下1stOpt的介绍，感觉略显浮夸已经有14人回复
MATLAB解非线性方程组已经有9人回复
请教高手，如何解非线性方程组！！！已经有13人回复
非线性方程的线性化问题，大神们进来指教下！！！！感谢了！！！！！！！！！已经有9人回复
1stopt解方程组已经有12人回复
牛顿迭代法解非线性矩阵方程已经有4人回复
用matlab求解一个非线性方程组的解已经有4人回复
用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？已经有11人回复
sos!!1stopt解非线性方程组，一直出现常数定义错误的问题！求大神帮助！已经有12人回复
matlab解非线性方程组解不开已经有3人回复
用1stopt拟合非线性方程结果与其他软件拟合结果差异大已经有5人回复
1stopt运行程序，总是出现字符函数错误的提示，请高手帮助已经有5人回复
用1stopt或者origin非线性拟合磁导率已经有5人回复
1stopt4阶非线性常微分方程，帮小弟看一下，跪求。已经有7人回复
1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？已经有13人回复
【求助】求教matlab解非线性方程组已经有9人回复

1楼 2015-05-12 17:21:57

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

仰望星空7

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1237.7
帖子: 42
在线: 21.3小时
虫号: 4127268
注册: 2015-10-08
专业: 工程热物理相关交叉领域

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的回帖

查看全部 7 个回答

dingd

铁杆木虫 (职业作家)

计算强帖: 4
应助: 1641 (讲师)
金币: 15037.3
散金: 101
红花: 234
帖子: 3410
在线: 1223.7小时
虫号: 291104
注册: 2006-10-28

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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

赞一下

回复此楼