24小时热门版块排行榜    

查看: 1839  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

as23p

木虫 (小有名气)

[求助] 关于Matlab算符号积分

请问大家,用Matlab算积分是不是要有什么技巧才能保证正确性啊?

工作中需要计算这个积分的解析表达:int('(a*cos(x)-b)^28'),其中我对于a和b的选取是可以保证括号里的量小于1的。但用Matlab算了后,代入数值发现结果相当大。然后我干脆用Matlab执行命令:int('(6*sqrt(10)*cos(x)-18)^28',0,acos(3/sqrt(10))),结果高达10^26,明显错误,因为被积函数是小于1的,因此积分值按说是不应大于acos(3/sqrt(10))。

有同学知道是什么原因导致这种结果以及如何纠正么?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

as23p

木虫 (小有名气)

jjdg: 你贴个数学形式的公式出来可能更方便大家帮你看看 2012-01-16 10:11:28
我觉得不是你说的原因。我手动算了一下,结果和Matlab给出的解析表达式是一致的,但代入积分限后就不对 我也自己编程算了,用四精度一样是错误。由于这个积分只是我整个计算的一小环,所以如果需要专门计算它从总计算量上来讲是无法忍受的。
3楼2012-01-16 06:06:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

zhfzh

木虫 (正式写手)

【答案】应助回帖


感谢参与,应助指数 +1
jjdg(金币+1): 感谢参与 2012-01-15 10:04:01
as23p(金币+10): 有帮助 是不老靠谱的:( 2012-02-09 07:13:38
matlab的默认精度是10的-16
你直接算误差是很大的
建议自己编程计算比较靠普
2楼2012-01-14 12:04:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

as23p

木虫 (小有名气)

没人理也无所谓了,我已经自己解决了
4楼2012-01-20 07:33:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

snoopyzhao

至尊木虫 (职业作家)

★ ★
jjdg(金币+1): 感谢参与 2012-01-24 22:16:45
jjdg(金币+1): 春节快乐 2012-01-24 22:16:52
引用回帖:
4楼: Originally posted by as23p at 2012-01-20 07:33:21:
没人理也无所谓了,我已经自己解决了

可否给出你的解决方案?

我用 maxima 算的结果也异常大……
5楼2012-01-24 09:22:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856调剂,是学校就去 +5 sllhht 2026-03-19 6/300 2026-03-19 19:15 by ms629
[考研] 288求调剂 +13 于海海海海 2026-03-19 13/650 2026-03-19 19:11 by 枫桥ZL
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +13 rare12345 2026-03-18 13/650 2026-03-19 17:30 by jean5056
[考研] 梁成伟老师课题组欢迎你的加入 +9 一鸭鸭哟 2026-03-14 11/550 2026-03-19 17:22 by !本暗一次!
[考研] 321求调剂 +8 何润采123 2026-03-18 10/500 2026-03-19 16:46 by 何润采123
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] 317求调剂 +3 申子申申 2026-03-19 6/300 2026-03-19 14:16 by 申子申申
[考研] 一志愿吉林大学材料学硕321求调剂 +6 Ymlll 2026-03-18 9/450 2026-03-19 10:28 by 星空星月
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 085410人工智能专硕317求调剂(0854都可以) +3 xbxudjdn 2026-03-18 3/150 2026-03-18 22:14 by zhq0425
[考研] 302求调剂 +10 呼呼呼。。。。 2026-03-17 10/500 2026-03-18 12:45 by Linda Hu
[考研] 290求调剂 +3 p asserby. 2026-03-15 4/200 2026-03-17 16:35 by wangkm
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 275求调剂 +4 太阳花天天开心 2026-03-16 4/200 2026-03-17 10:53 by 功夫疯狂
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[基金申请] 今年的国基金是打分制吗? 50+3 zhanghaozhu 2026-03-14 3/150 2026-03-16 17:07 by 北京莱茵润色
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 0703一志愿211 285分求调剂 +5 ly3471z 2026-03-13 5/250 2026-03-16 16:16 by 哦哦123
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
信息提示
请填处理意见