24小时热门版块排行榜    

查看: 1833  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 26博士申请 +3 1042136743 2026-03-17 3/150 2026-03-17 23:30 by 轻松不少随
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-16 8/400 2026-03-17 19:03 by Wangjingyue
[考研] 304求调剂 +8 小熊joy 2026-03-14 8/400 2026-03-17 17:29 by ruiyingmiao
[考研] 302求调剂 +9 负心者当诛 2026-03-11 9/450 2026-03-17 17:13 by ruiyingmiao
[考研] 332求调剂 +6 Zz版 2026-03-13 6/300 2026-03-17 17:03 by ruiyingmiao
[考研] 26考研求调剂 +6 丶宏Sir 2026-03-13 6/300 2026-03-17 16:13 by 醉在风里
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] 271求调剂 +12 生如夏花… 2026-03-11 14/700 2026-03-17 10:56 by lovewei0727
[硕博家园] 深圳大学硕士招生(2026秋,传感器方向,仅录取第一志愿) +4 xujiaoszu 2026-03-11 9/450 2026-03-17 10:29 by xujiaoszu
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 085600调剂 +5 漾漾123sun 2026-03-12 6/300 2026-03-16 15:58 by 漾漾123sun
[教师之家] 焦虑 +7 水冰月月野兔 2026-03-13 9/450 2026-03-16 10:00 by Quakerbird
[考研] 材料工程327求调剂 +3 xiaohe12w 2026-03-11 3/150 2026-03-14 20:20 by ms629
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 材料与化工(0856)304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 工科调剂 +4 Jiang191123! 2026-03-11 4/200 2026-03-13 15:15 by Miko19
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 一志愿山大07化学 332分 四六级已过 本科山东双非 求调剂! +3 不想理你 2026-03-12 3/150 2026-03-13 14:18 by JourneyLucky
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
信息提示
请填处理意见