应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 解方程，希望大家帮帮忙！谢谢！
162/2返回列表上一页12
查看: 1888  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

chyanog

金虫 (小有名气)
引用回帖:
7楼: Originally posted by walk1997 at 2013-12-24 09:03:42
又试了下
不用分数代替也可以 不过是复数解(精度也只能是机器精度之内左边==0)
{t1 -> 1174.72 - 207.014 I, t2 -> 772.469 + 196.068 I}
也许还有其他复数解...

1stOpt是全局求解吗，可以算复数解吗
赞一下 回复此楼
11楼2013-12-24 17:33:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

walk1997

金虫 (著名写手)
引用回帖:
8楼: Originally posted by dingd at 2013-12-24 10:03:03
1stOpt试试，很简单，数值解的话比mathematica强的多：

Function
1*(1273 - t1) + 1.76*10^(-4)*5.67*10^(-8)*(1273^4 - (t1)^4) - ((1.5 + 0.6 + (
2*1.38*10^(-23)*t1)/(1.6*10^(-19)))*1.2*10^6*(t1)^2*
E ...

谢谢提醒 LZ的方程是有实数解
NSolve应该还是不够强大 对于无穷多解方程 在给定解参数范围下
没给出解 这点是不够的(对有些方程能给出给定范围的全部解)
不过FindRoot能给出上面这解，另外单独这个的t1 t2带入表达式中，
给出的精度是10^-5 相对于t1,t2也够了，不过总感觉不太满意
mathematica给出的解是
{t1->1721.1348184716746049, t2 -> 279.71114194601986316}
精度在10^-16
另外也有其他实数解 比如
{t1 -> 1937.1675296975197254, t2 -> 355.62015391388776230}

不知道有没办法给出对全部实数解的讨论。。。。
赞一下 回复此楼
12楼2013-12-24 22:17:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chyanog

金虫 (小有名气)
引用回帖:
12楼: Originally posted by walk1997 at 2013-12-24 22:17:28
谢谢提醒 LZ的方程是有实数解
NSolve应该还是不够强大 对于无穷多解方程 在给定解参数范围下
没给出解 这点是不够的(对有些方程能给出给定范围的全部解)
不过FindRoot能给出上面这解，另外单独这个的t1 t2带入表 ...

超越方程的全局求解本身就不容易，Matlab里的fsolve、1stOpt也都不是全局求解的
赞一下 回复此楼
13楼2013-12-24 23:28:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Simewe

新虫 (初入文坛)

【答案】应助回帖

引用回帖:
8楼: Originally posted by dingd at 2013-12-24 10:03:03
1stOpt试试，很简单，数值解的话比mathematica强的多：

Function
1*(1273 - t1) + 1.76*10^(-4)*5.67*10^(-8)*(1273^4 - (t1)^4) - ((1.5 + 0.6 + (
2*1.38*10^(-23)*t1)/(1.6*10^(-19)))*1.2*10^6*(t1)^2*
E ...

1stOpt算数值解强在哪里，是求解速度还是精度？
赞一下(1人) 回复此楼
14楼2013-12-24 23:34:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

walk1997

金虫 (著名写手)
引用回帖:
13楼: Originally posted by chyanog at 2013-12-24 23:28:09
超越方程的全局求解本身就不容易，Matlab里的fsolve、1stOpt也都不是全局求解的...

嗯 是的 不过我想的是指定解的范围 然后给出全部解
这点NSolve对有些超越方程应该是可以做到的，比如
f = Exp[x/100] - Sin[x] == 0
NSolve[{f, -20 < x < 20}, x]
能给出解，
但是对LZ的方程组，失效了，和自己对比NSolve还是不够强
原则上 这样指定范围内的解应该是可以全部给出....
赞一下 回复此楼
15楼2013-12-25 08:00:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingd

铁杆木虫 (职业作家)

【答案】应助回帖

引用回帖:
14楼: Originally posted by Simewe at 2013-12-24 23:34:04
1stOpt算数值解强在哪里，是求解速度还是精度？...

求解能力方面吧！采用的是全局优化算法，数值解方程能力比Mathematica的NSolve和Matlab的fsolve强不少。
赞一下 回复此楼
16楼2013-12-25 09:38:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yangyang1231 的主题更新
162/2返回列表上一页12
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 4/200 2026-02-23 06:46 by jsjzfl
[硕博家园] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 8rmuugja8q 2026-02-22 6/300 2026-02-23 06:39 by w4l55oybr1
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 w89i99eaeh 2026-02-22 4/200 2026-02-23 06:36 by w4l55oybr1
[博后之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 khieu8v8m0 2026-02-22 5/250 2026-02-23 06:34 by w4l55oybr1
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 5/250 2026-02-23 06:29 by w4l55oybr1
[博后之家] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +6 3dfhjxgsh7 2026-02-22 8/400 2026-02-23 06:21 by w4l55oybr1
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 usprnugpzw 2026-02-21 10/500 2026-02-23 04:58 by 5jlh3qtdvx
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 khieu8v8m0 2026-02-22 6/300 2026-02-23 02:08 by 5jlh3qtdvx
[考博] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +5 3dfhjxgsh7 2026-02-22 6/300 2026-02-23 02:04 by 5jlh3qtdvx
[教师之家] 版面费该交吗 +7 苹果在哪里 2026-02-22 8/400 2026-02-22 22:37 by otani
[基金申请] 面上可以超过30页吧？ +4 阿拉贡aragon 2026-02-22 4/200 2026-02-22 21:22 by 山西悬空寺空悬�
[教师之家] 为什么中国大学教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？ +5 欢乐颂叶蓁 2026-02-21 5/250 2026-02-22 21:15 by 山西悬空寺空悬�
[论文投稿] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +4 usprnugpzw 2026-02-21 6/300 2026-02-22 19:48 by w89i99eaeh
[考研] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 3dfhjxgsh7 2026-02-22 4/200 2026-02-22 16:52 by khieu8v8m0
[找工作] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-22 3/150 2026-02-22 16:37 by khieu8v8m0
[公派出国] 售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急 +3 usprnugpzw 2026-02-21 4/200 2026-02-22 16:27 by khieu8v8m0
[基金申请] “人文社科而论，许多学术研究还没有达到民国时期的水平” +4 苏东坡二世 2026-02-18 5/250 2026-02-22 16:07 by liangep1573
[基金申请] 什么是人一生最重要的？ +4 瞬息宇宙 2026-02-21 4/200 2026-02-22 11:44 by huagongfeihu
[基金申请] 今年春晚有几个节目很不错，点赞！ +11 瞬息宇宙 2026-02-16 12/600 2026-02-21 21:14 by lq493392203
[基金申请] 体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低 +9 瞬息宇宙 2026-02-20 12/600 2026-02-21 10:39 by 欢乐颂叶蓁
信息提示
关闭
请填处理意见
关闭