| 查看: 4738 | 回复: 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 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有59人回复
-
filecode=后面第一个是大写字母
已经有5人回复
-
时间戳又变了8-15
已经有21人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有10人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求教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人回复
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
2楼2013-06-15 17:38:26
nono2009
超级版主 (文学泰斗)
No gains, no pains.
-
专家经验: +21105 - 应助: 28684 (院士)
- 贵宾: 513.911
- 金币: 2555230
- 散金: 27828
- 红花: 2148
- 沙发: 66666
- 帖子: 1602255
- 在线: 65200.9小时
- 虫号: 827383
- 注册: 2009-08-13
- 性别: GG
- 专业: 工程热物理与能源利用
- 管辖: 科研家筹备委员会
3楼2013-06-15 17:53:25
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
6楼2013-06-16 11:01:05
xiangzigood
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 321.5
- 红花: 1
- 帖子: 18
- 在线: 11.7小时
- 虫号: 1891458
- 注册: 2012-07-13
- 专业: 制造系统与自动化
7楼2013-06-17 08:26:26
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
8楼2013-06-17 08:31:02
9楼2013-06-17 10:40:40
xiangzigood
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 321.5
- 红花: 1
- 帖子: 18
- 在线: 11.7小时
- 虫号: 1891458
- 注册: 2012-07-13
- 专业: 制造系统与自动化
|
function[x,x0,fval,exitflag]=fsolveyalinjieliu d=0.04;D=0.05;k=1.4;ff=0.004;L=1000;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)*sqrt(2*((Pa/x(5))^(2/k)-(Pa/x(5))^((k+1)/k))/(k-1))-x(2); (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); x(6)-(x(3)*T1)/x(4); ] end end 这个应该比较完整，六个方程六个未知数 |
10楼2013-06-17 11:21:58
