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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » 关于Matlab算符号积分
71/1返回列表
查看: 1830  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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))。

有同学知道是什么原因导致这种结果以及如何纠正么？
回复此楼

» 收录本帖的淘帖专辑推荐

文献

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-01-13 13:06:58
已阅   回复此楼   关注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
你直接算误差是很大的
建议自己编程计算比较靠普
赞一下(1人) 回复此楼
2楼2012-01-14 12:04:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

as23p

木虫 (小有名气)
jjdg: 你贴个数学形式的公式出来可能更方便大家帮你看看 2012-01-16 10:11:28
我觉得不是你说的原因。我手动算了一下，结果和Matlab给出的解析表达式是一致的，但代入积分限后就不对 我也自己编程算了，用四精度一样是错误。由于这个积分只是我整个计算的一小环，所以如果需要专门计算它从总计算量上来讲是无法忍受的。
赞一下(1人) 回复此楼
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 算的结果也异常大……
赞一下(2人) 回复此楼
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给出的解析结果是正确的，但数值结果不对，应该是因为各项的系数过大，代入数值的话变成了很大的数乘以很小的数，造成了很大的数值误差？最后在我的程序中改用高斯积分公式了。
赞一下(1人) 回复此楼
6楼2012-02-09 06:53:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

LYR90

新虫 (初入文坛)
请问，符号矩阵运算时，由于四到五个符号矩阵进行乘法以及求逆运算，所以一直求解不出结果来，有没有什么好的简化方法？
赞一下 回复此楼
7楼2014-05-12 10:34:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 as23p 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 326求调剂 +4 上岸的小葡 2026-03-15 5/250 2026-03-16 08:39 by Linda Hu
[考研] 化学调剂0703 +7 啊我我的 2026-03-11 7/350 2026-03-15 23:03 by 凌千颂111
[考研] 调剂 +8 调剂的考研学生 2026-03-09 8/400 2026-03-15 22:14 by Winj1e
[考研] 中科院材料273求调剂 +3 yzydy 2026-03-15 3/150 2026-03-15 21:15 by ms629
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 297求调剂 +4 学海漂泊 2026-03-13 4/200 2026-03-14 11:51 by 热情沙漠
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 0703求调剂 +7 jtyq001 2026-03-10 7/350 2026-03-14 01:06 by JourneyLucky
[考研] 2026考研调剂+本科延边大学+山东大学+生物化学与分子生物学+有项目经验 +3 ccdsscjy 2026-03-10 3/150 2026-03-14 00:12 by JourneyLucky
[考研] 材料与化工（0856）304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 341求调剂 +4 番茄头--- 2026-03-10 4/200 2026-03-13 23:12 by JourneyLucky
[考研] 290求调剂 +9 ADT 2026-03-11 9/450 2026-03-13 21:55 by JourneyLucky
[考研] 285化工学硕求调剂（081700） +6 柴郡猫_ 2026-03-12 6/300 2026-03-13 20:46 by hmn_wj
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
[考研] 收调剂 +7 调剂的考研学生 2026-03-10 7/350 2026-03-10 17:57 by 麦茶汤圆
信息提示
关闭
请填处理意见
关闭