| 查看: 4204 | 回复: 15 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
xiangzigood新虫 (初入文坛)
|
[求助]
matlab fsolve得出的解全部为复数解
|
||
|
我用matlab中的FSOLVE解一个六元非线性方程组，为什么得出的解全部为复数解，可是理论上应该解为实数，本人Matlab新手，代码如下，大神们帮忙看一下吧！ function[x,x0,fval,exitflag]=fsolveyalinjieliu d=0.045;D=0.05;k=1.4;ff=0.004;L=100;P1=1500000;T1=300;Pa=101325;M=29;R=8.314; x0=unifrnd(0,1,1,6); options=optimset('MaxFunEvals',20000,'MaxIter',5000); [x,fval,exitflag]=fsolve(@f,x0,options); function F=f(x) F=[((k+1)/2)*log((x(2)^2*x(3))/(x(1)^2*x(4)))-(1/(x(1)^2)-1/(x(2)^2))+4*ff*L/D; (d^2/D^2)*x(5)*sqrt((2*k*M)*((Pa/x(5))^(2/k)-(Pa/x(5))^((k+1)/k))/(R*x(6)*(k-1)))-x(2)*x(5)*sqrt((k*M)/(R*x(6))); (d^2/D^2)*x(5)*sqrt((2*k*M)*((Pa/x(5))^(2/k)-(Pa/x(5))^((k+1)/k))/(R*x(6)*(k-1)))-x(1)*P1*sqrt((k*M)/(R*T1)); 1+((k-1)*x(1)^2)/2-x(3); 1+((k-1)*x(2)^2)/2-x(4); ] end end |
回复此楼» 猜你喜欢
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗
已经有8人回复
-
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有5人回复
-
论文投稿，期刊推荐
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有4人回复
-
孩子确诊有中度注意力缺陷
已经有14人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
求教fsolve如何查看每次迭代的x值？
已经有6人回复- 求助 MATLAB-fsolve
已经有3人回复
- 求助matlab 非线性方程组全部解！
已经有6人回复
- 使用fsolve求解非线性方程问题
已经有9人回复
- 使用matlab 中的fsolve 老提示出错
已经有4人回复
- 各位朋友，如何用matlab编写这个方程，并求解？
已经有15人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- Matlab,解复数方程
已经有8人回复
- 求助！matlab用fsolve函数求解非线性方程组的问题！
已经有19人回复
- 请教下：如何求一个函数在一定范围内的全部解
已经有10人回复
- 求各位前辈帮忙！matlab解复数方程，出现图中的警告，怎来修改呢？谢谢了啊
已经有5人回复
- 【求助】matlab中fsolve解的精度控制【已解决】
已经有4人回复
- 一个方程组的解法
已经有8人回复
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555150
- 散金: 27828
- 红花: 2147
- 沙发: 66666
- 帖子: 1602244
- 在线: 65200.8小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
3楼2013-06-15 17:53:25
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.5小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2013-06-15 17:38:26
xiangzigood
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 321.5
- 红花: 1
- 帖子: 18
- 在线: 11.7小时
- 虫号: 1891458
- 注册: 2012-07-13
- 专业: 制造系统与自动化
4楼2013-06-16 08:25:52
xiangzigood
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 321.5
- 红花: 1
- 帖子: 18
- 在线: 11.7小时
- 虫号: 1891458
- 注册: 2012-07-13
- 专业: 制造系统与自动化
5楼2013-06-16 08:26:27