24小时热门版块排行榜    

查看: 2132  |  回复: 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的回帖

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

木虫 (小有名气)

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

as23p

木虫 (小有名气)


jjdg(金币+1): 感谢参与 2012-02-09 09:57:15
引用回帖:
: Originally posted by snoopyzhao at 2012-01-24 09:22:44:
可否给出你的解决方案?

我用 maxima 算的结果也异常大……

Matlab给出的解析结果是正确的,但数值结果不对,应该是因为各项的系数过大,代入数值的话变成了很大的数乘以很小的数,造成了很大的数值误差?最后在我的程序中改用高斯积分公式了。
6楼2012-02-09 06:53:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

LYR90

新虫 (初入文坛)

请问,符号矩阵运算时,由于四到五个符号矩阵进行乘法以及求逆运算,所以一直求解不出结果来,有没有什么好的简化方法?
7楼2014-05-12 10:34:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 as23p 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 10/500 2026-08-17 23:43 by sytp
[基金申请] 93BebMhtakh前后11位开头都是大写 +3 且听虎啸 2026-08-17 4/200 2026-08-17 22:31 by ccec101
[基金申请] filecode=后面第一个是大写字母 +8 wangze12014 2026-08-14 10/500 2026-08-17 17:05 by xter9665
[基金申请] 时间戳变了,能看出什么问题? +4 基诺咪客 2026-08-17 4/200 2026-08-17 16:10 by Vivilian
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +12 Ldrop2023 2026-08-13 15/750 2026-08-17 15:02 by 小豌豆_发芽
[基金申请] 时间戳又变了8-15 +13 archvillain 2026-08-15 25/1250 2026-08-16 20:13 by zhaosm1982
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见,大家看看几个A,几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
信息提示
请填处理意见