| 查看: 734 | 回复: 11 | |||
| 当前主题已经存档。 | |||
ltqly木虫 (小有名气)
|
[交流]
【求助】方程求解
|
||
|
Sample Text tan（x*d/sqrt(a))=x/sqrt(a)*(H1+H2)/(x*x-H1*H2) 求解这个方程的解x， a=5.1e-6 H1=60.0/30.3 H2=30.0/30.3 本人已经初步试着，用牛顿法解过，求出的x*d/sqrt(a)接近于k*pi，pi=3.1415926，x有无穷解。 谢谢 |
回复此楼» 猜你喜欢
0856求调剂
已经有13人回复
-
0703化学
已经有11人回复
-
329求调剂
已经有7人回复
-
321求调剂
已经有8人回复
-
0856材料化工调剂 总分330
已经有12人回复
-
309求调剂
已经有3人回复
-
070300化学354求调剂
已经有4人回复
-
330分求调剂
已经有3人回复
-
299求调剂
已经有8人回复
-
一志愿北京理工大学本科211材料工程294求调剂
已经有6人回复
★ ★ ★ ★ ★ ★
sunxiao(金币+3,VIP+0):谢谢参与，欢迎常来仿真编程版 4-8 23:39
ltqly(金币+3,VIP+0):谢谢 4-9 12:07
sunxiao(金币+3,VIP+0):谢谢参与，欢迎常来仿真编程版 4-8 23:39
ltqly(金币+3,VIP+0):谢谢 4-9 12:07
|
d是多少啊，我用MATLAB搞了一下，设d=a算的,把方程的左边减右边设成一个函数，画图 function y=fun(x) a=5.1e-6; H1=60.0/30.3; H2=30.0/30.3; y=tan(x.*a./sqrt(a))-x./sqrt(a).*(H1+H2)./(x.*x-H1.*H2); 可以看出，原方程只有零解，正负1.4左右有两个无穷间断点，x有无穷解 [ Last edited by fspdlh on 2009-4-8 at 23:17 ] |
2楼2009-04-08 22:55:37
ltqly
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2863.3
- 散金: 716
- 帖子: 177
- 在线: 250.4小时
- 虫号: 686537
- 注册: 2009-01-02
- 性别: GG
- 专业: 零件成形制造
3楼2009-04-09 12:08:05
ltqly
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2863.3
- 散金: 716
- 帖子: 177
- 在线: 250.4小时
- 虫号: 686537
- 注册: 2009-01-02
- 性别: GG
- 专业: 零件成形制造
4楼2009-04-09 12:13:25
★ ★ ★ ★ ★
ltqly(金币+3,VIP+0):谢谢! 4-9 17:52
sunxiao(金币+2,VIP+0):谢谢参与，欢迎常来仿真编程版 4-9 23:10
ltqly(金币+3,VIP+0):谢谢! 4-9 17:52
sunxiao(金币+2,VIP+0):谢谢参与，欢迎常来仿真编程版 4-9 23:10
|
看了一下，确实应该是在每个区间都有解，上由于计算精度的原因画图得不到想要的曲线。 由于函数存在无穷间断点，牛顿法应该是不行的，可以仿照这个帖子中的方法，对每个有解区间分别用二分法来做。 帖子地址：http://muchong.com/bbs/viewthread.php?tid=1251344 [ Last edited by fspdlh on 2009-4-9 at 13:02 ] |
5楼2009-04-09 12:23:09
ltqly
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2863.3
- 散金: 716
- 帖子: 177
- 在线: 250.4小时
- 虫号: 686537
- 注册: 2009-01-02
- 性别: GG
- 专业: 零件成形制造
6楼2009-04-09 17:54:23
7楼2009-04-09 18:11:05
ltqly
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2863.3
- 散金: 716
- 帖子: 177
- 在线: 250.4小时
- 虫号: 686537
- 注册: 2009-01-02
- 性别: GG
- 专业: 零件成形制造
8楼2009-04-09 22:20:32
9楼2009-04-09 22:46:13
ltqly
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2863.3
- 散金: 716
- 帖子: 177
- 在线: 250.4小时
- 虫号: 686537
- 注册: 2009-01-02
- 性别: GG
- 专业: 零件成形制造
10楼2009-04-10 22:19:13
20