| 查看: 1834 | 回复: 6 | |||
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人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有3人回复
-
物理学I论文润色/翻译怎么收费?
已经有245人回复
-
基底STO,薄膜SRO，XRD里面的振荡，是laue震荡还是kiessig振荡? 怎么判断？
已经有2人回复
-
B区学生调剂-兰州交通大学材料科学与工程学院
已经有11人回复
-
山西大同大学物理学专业还有调剂名额，欢迎调剂！
已经有12人回复
-
桂林理工大学物理学专业招收调剂
已经有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人回复
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
people891108: 金币+20, ★★★很有帮助, 非常感谢 2015-05-13 09:04:34
感谢参与，应助指数 +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朵）
people8911082楼2015-05-12 20:14:30
people891108
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 4105.6
- 散金: 58
- 红花: 7
- 帖子: 468
- 在线: 193.2小时
- 虫号: 2122748
- 注册: 2012-11-12
- 专业: 化学反应工程
3楼2015-05-13 09:42:26
dingd
铁杆木虫 (职业作家)
- 计算强帖: 4
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
4楼2015-05-13 09:51:11
情弦风
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 907.2
- 散金: 145
- 帖子: 95
- 在线: 59.6小时
- 虫号: 3754702
- 注册: 2015-03-21
- 专业: 磁共振结构成像与疾病诊断
5楼2016-02-25 20:48:55
|
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
7楼2017-04-13 23:00:23
