24小时热门版块排行榜    

查看: 639  |  回复: 3
本帖产生 1 个 博学EPI ,点击这里进行查看

sqking

金虫 (正式写手)

[求助] 求解方程

各位虫友,我在读文献的过程中想验证一个模型的正确性,需要解一个方程,具体如下:tan(X)=3/(3/(X)+0.2092*X) (变量为X)。请大家帮忙解答,给出其中的10个大于0的解,谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

灰头灰雀

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
sqking: 金币+10, 博学EPI+1, ★有帮助 2013-07-22 19:26:44
这个问题可以用matlab来解决。具体的解法是:1、首先,令f=3./(3./(X)+0.2092*X) -tan(X),则问题转化为求解f函数的零点问题。2、然后确定区间:取x的取值区间为[0,20],在matlab中输入为x=0:0.1:20。然后做出图形:plot(x,f)。在图形中可以看出区间正好包含了10个所需零点位置。3、求零点值:在图形中可以大致看出零点的位置则用fzero函数可以求解具体位置。比如第一个零点在2附近,则输入fzero('3./(3./(X)+0.2092*X) -tan(X)',2),则可以得到第一个零点为1.5708,依此类推可以得出前十个零点值为1.5708、4.2204、4.7124、7.2811、7.8540、10.3111、10.9956、13.3489、14.1372、16.4008。
坚持就是胜利
2楼2013-07-22 14:24:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lam987

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
sqking: 金币+10, ★有帮助 2013-07-22 19:26:55
楼上的方法有问题,第一个数验证就不对。
通过matlab求解
syms x
y=solve('tan(x)*(3/x+0.2029*x)=3','x')
y=
   [-10.31114532]
   [                   ]
   [-13.34888952]
   [                   ]
   [-7.281075701]
   [                   ]
   [       0.          ]
有四个结果,分别是-10.31114532,-13.34888952,-7.281075701和0,而“0”有没有意义还需另说。
3楼2013-07-22 19:23:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sqking

金虫 (正式写手)

谢谢,实验室有高人,已经把这个问题搞定了。
4楼2013-07-22 19:27:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 sqking 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 BMTdnp4OJO9f 2026-11-06 4/200 2026-11-07 19:36 by Xc2FxHsHkBNc
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 2rKfVQg2jQww 2026-11-06 5/250 2026-11-07 17:23 by Xc2FxHsHkBNc
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 17:01 by Xc2FxHsHkBNc
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 5sJFDKp2kwuI 2026-11-06 5/250 2026-11-07 16:06 by Xc2FxHsHkBNc
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 5sJFDKp2kwuI 2026-11-06 6/300 2026-11-07 15:44 by Xc2FxHsHkBNc
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +8 ncTVBM6AOJHF 2026-11-04 23/1150 2026-11-07 06:49 by Equinoxhua
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-05 7/350 2026-11-06 16:43 by 5sJFDKp2kwuI
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +6 ELBZ6WwEaRj6 2026-11-05 8/400 2026-11-06 16:34 by 5sJFDKp2kwuI
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 ELBZ6WwEaRj6 2026-11-04 14/700 2026-11-06 16:32 by 5sJFDKp2kwuI
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 ELBZ6WwEaRj6 2026-11-04 10/500 2026-11-06 16:23 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 ncTVBM6AOJHF 2026-11-04 15/750 2026-11-06 13:49 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ncTVBM6AOJHF 2026-11-04 11/550 2026-11-06 13:47 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +8 ncTVBM6AOJHF 2026-11-04 14/700 2026-11-06 13:36 by BkA9Bi5cMhkD
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +8 ncTVBM6AOJHF 2026-11-04 18/900 2026-11-06 13:36 by BkA9Bi5cMhkD
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 kvewwGtKktqB 2026-11-04 16/800 2026-11-06 05:00 by qL1cU46YIdK0
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ELBZ6WwEaRj6 2026-11-04 8/400 2026-11-06 02:24 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 hNpp4Y2CGvSC 2026-11-04 12/600 2026-11-06 01:53 by qL1cU46YIdK0
[基金申请] 各位基础老师们,看完中央《从严治党》的决定后,对学术界更加有信心了吧 +5 iaeyuan 2026-11-01 5/250 2026-11-04 14:43 by jinchengzi
信息提示
请填处理意见