| 查看: 1134 | 回复: 4 | |||
[求助]
数值解非线性方程组时找不到解，希望各位能帮帮我 已有2人参与
|
|
function y=fun(x) alpha=0.125*pi/180;x1=0;z1=5;x2=0;z2=25;L1=5;L2=5;f=30;h0=20;n2=1.49386; y(1) = tan(alpha) - x(1)/x(2); y(2) = n2*sin(atan((x1 - x(1))/(z1 - x(2))) - atan((x(2) - x(4))/(x(1) - x(3)))) - sin(pi/2 - alpha + atan((x1 - x(1))/(z1 - x(2)))) ; y(3) = n2*sin(atan((x2 - x(3))/(z2 - x(4))) - atan((x(2) - x(4))/(x(1) - x(3)))) - sin(pi/2 + atan((x2 - x(3))/(z2 - x(4)))); y(4) = L1 + L2 - f + h0 - h0*n2 + n2*((x(1) - x(3))^2 + (x(2) - x(4))^2)^(1/2); end ---------------------------------------------------------------------- 命令行输入 >> fsolve('fun',[0.1 5.1 0.1 25.1]) ---------------------------------------------------------------------- No solution found. fsolve stopped because the relative size of the current step is less than the default value of the step size tolerance squared, but the vector of function values is not near zero as measured by the default value of the function tolerance. 合理的解应该是在初值附近的，可为什么解不出来呢 |
回复此楼» 猜你喜欢
化工0817调剂
已经有7人回复
-
一志愿211，0860总分286食品工程专业求调剂
已经有8人回复
-
一志愿天津大学材料与化工275求调剂
已经有18人回复
-
标题:捡漏预警|08工科/09农学调剂!英语要求低，过线即有机会!
已经有11人回复
-
0856材料与化工353分求调剂
已经有11人回复
-
308求调剂
已经有7人回复
-
085600 材料与化工 295 求调剂
已经有8人回复
-
求调剂材料专硕293
已经有6人回复
-
欢迎211本科同学，过A区国家线，A区非偏远一本，交叉学科课题组
已经有37人回复
-
收调剂
已经有7人回复
2楼2019-01-30 17:49:49
hzlhm
至尊木虫 (著名写手)
- 应助: 387 (硕士)
- 金币: 17808.5
- 红花: 53
- 帖子: 2879
- 在线: 606.3小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
月只蓝: 金币+10, 感谢帮助！ 2019-01-31 18:04:04
感谢参与，应助指数 +1
月只蓝: 金币+10, 感谢帮助！ 2019-01-31 18:04:04
|
利用最小二乘法（lsqnonlin），可以求得 x1=0.0109083348225495 ;x2=4.99999947195261;x3=0.0384214997768871;4=25.0000321470904 验证： y1 = tan(alpha) - x(1)/x(2)=-2.16860202711944e-09 y2 = n2*sin(atan((x1 - x(1))/(z1 - x(2))) - atan((x(2) - x(4))/(x(1) - x(3)))) - sin(pi/2 - alpha + atan((x1 - x(1))/(z1 - x(2)))) =5.90246506399546e-06 y3 = n2*sin(atan((x2 - x(3))/(z2 - x(4))) - atan((x(2) - x(4))/(x(1) - x(3)))) - sin(pi/2 + atan((x2 - x(3))/(z2 - x(4))))=-3.15643589759254e-05 |
» 本帖已获得的红花（最新10朵）
XiaoGuads
3楼2019-01-31 14:39:33
4楼2019-02-04 16:56:15
dingd
铁杆木虫 (职业作家)
- 应助: 1641 (讲师)
- 金币: 15037.3
- 散金: 101
- 红花: 234
- 帖子: 3410
- 在线: 1223.7小时
- 虫号: 291104
- 注册: 2006-10-28
5楼2019-02-28 18:14:16
