| 查看: 268 | 回复: 3 | |||
| 当前主题已经存档。 | |||
lshuihui62铜虫 (小有名气)
|
[交流]
【求助】方程组求解
|
||
|
非线性方程组（见图）：其中Dζ、Dγ、Dβ、Yζ、Yγ、Yβ为未知数，对应到程序里面依次为x(1)...x(6)，其余均为常数，我想问一下在MATLAB中，可不可以用解一般非线性方程组的方法解这个非线性方程组，下面是我的程序代码，但是没有解出来，我想请教一下我的问题在哪里？未知数的最小值是零，但是为确保分母不为零，初值就不能赋为零。 function F=myfun(x) F=[((7.22-3.56)*sqrt(x(1))*exp(-(x(4))^2/(4*x(1))))/(erf((x(4))/(2*sqrt(x(1)))))-(sqrt(pi)*x(4)*(3.56-3.05)/2)-(((3.05-0.924)*sqrt(x(2))*exp(-(x(4))^2/(4*x(2))))/(erf((x(5))/(2*sqrt(x(2))))-erf((x(4))/(2*sqrt(x(2)))))); ((0.263-0.016)*sqrt(x(3))*exp(-(x(5))^2/(4*x(3))))/(erf((x(5))/(2*sqrt(x(3)))-erf((x(6))/(2*sqrt(x(3))))))-(sqrt(pi)*x(5)*(0.924-0.263)/2)-(((3.05-0.924)*sqrt(x(2))*exp(-(x(5))^2/(4*x(2))))/(erf((x(4))/(2*sqrt(x(2))))-erf((x(5))/(2*sqrt(x(2))))))]; >> x0=[]; >> options=optimset('Display','iter'); >> [x,fval]=fsolve(@myfun,x0,options) |
回复此楼» 猜你喜欢
重庆交大26年硕士生招生拟调剂通知已出！欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
UJN物理学专业调剂
已经有12人回复
-
物理学I论文润色/翻译怎么收费?
已经有236人回复
-
0702一志愿吉大B区求调剂 本科期间发表一篇Sci
已经有0人回复
-
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组招收2026硕士研究生1名
已经有0人回复
-
山西大同大学物理学硕士研究生招收调剂生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
［调剂信息］211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
[调剂信息]211智能人工感知方向国家青年特聘专家课题组招收调剂研究生
已经有0人回复
-
一志愿物理所凝聚态物理352求调剂
已经有0人回复
q68
铁杆木虫 (著名写手)
- 应助: 3 (幼儿园)
- 贵宾: 1.615
- 金币: 14800.2
- 散金: 10
- 红花: 3
- 帖子: 2323
- 在线: 149.8小时
- 虫号: 325340
- 注册: 2007-03-17
- 性别: GG
- 专业: 金属材料的合金相、相变及
2楼2009-06-01 18:34:52
xjxj1234xjxj
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1154.7
- 散金: 100
- 帖子: 304
- 在线: 52.9小时
- 虫号: 534354
- 注册: 2008-03-27
- 专业: 人力资源开发与管理
3楼2009-06-04 22:54:30
lshuihui62
铜虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 42.2
- 散金: 10
- 帖子: 102
- 在线: 19.1小时
- 虫号: 762831
- 注册: 2009-05-04
- 专业: 金属功能材料
4楼2009-06-07 10:55:30
