| 查看: 602 | 回复: 4 | ||
[求助]
matlab计算 已有1人参与
|
|
b=Integrate(f(x),0,a). 当a 为一n维向量时，b也同样为一n维向量。 在matlab中是如何实现运算的呢，求高手帮帮忙啊！ |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : matlab计算.docx
2014-04-10 09:16:08, 17.69 K
» 猜你喜欢
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有7人回复
-
2026年国自然面上资助率
已经有12人回复
-
小木虫看见有人已经知道结果了
已经有20人回复
-
微信指数没变化，科研之友没阅读
已经有18人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有8人回复
-
HXDI做水性聚氨酯乳液，是不是特别容易出渣
已经有3人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有15人回复
-
申请2027年材料/化学类博士
已经有5人回复
-
你们的时间戳变了吗
已经有4人回复
-
今年的WR进展到哪一步了？
已经有9人回复
mathstudy
金虫 (正式写手)
- 数学EPI: 2
- 应助: 161 (高中生)
- 金币: 4320.1
- 散金: 38
- 红花: 16
- 帖子: 446
- 在线: 141.1小时
- 虫号: 2515489
- 注册: 2013-06-20
- 性别: GG
- 专业: 计算数学与科学工程计算
kexiaolong2楼2014-04-10 09:38:34
3楼2014-04-10 09:51:23
mathstudy
金虫 (正式写手)
- 数学EPI: 2
- 应助: 161 (高中生)
- 金币: 4320.1
- 散金: 38
- 红花: 16
- 帖子: 446
- 在线: 141.1小时
- 虫号: 2515489
- 注册: 2013-06-20
- 性别: GG
- 专业: 计算数学与科学工程计算
|
大体上是这样 备注：运行中会有提示说Warning: Maximum function count exceeded; singularity likely. 这是因为你的被积函数值很大>exp(38)=3.1856e+016；(下面的程序放在M文件里，在命令窗口中输入a 调用Integral(a)就可以了) function b=Integral(a) % main program nn=length(a); b=zeros(nn,1); for n=1:nn b(n)=quad(@myfun,0,a(n)); end end function f=myfun(x) % integral function f=exp(38+32.633./(2.586+0.0138*x+0.834*x.^0.372)); end |
4楼2014-04-10 22:29:34
mathstudy
金虫 (正式写手)
- 数学EPI: 2
- 应助: 161 (高中生)
- 金币: 4320.1
- 散金: 38
- 红花: 16
- 帖子: 446
- 在线: 141.1小时
- 虫号: 2515489
- 注册: 2013-06-20
- 性别: GG
- 专业: 计算数学与科学工程计算
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kexiaolong: 金币+30, ★★★★★最佳答案, 谢谢帮助！ 2014-04-12 12:49:45
kexiaolong: 金币+30, ★★★★★最佳答案, 谢谢帮助！ 2014-04-12 12:49:45
|
大体上是这样 备注：运行中会有提示说Warning: Maximum function count exceeded; singularity likely. 这是因为你的被积函数值很大>exp(38)=3.1856e+016；(下面的程序放在M文件里，在命令窗口中输入a 调用Integral(a)就可以了) function b=Integral(a) % main program nn=length(a); b=zeros(nn,1); for n=1:nn b(n)=quad(@myfun,0,a(n)); end end function f=myfun(x) % integral function f=exp(38+32.633./(2.586+0.0138*x+0.834*x.^0.372)); end |
5楼2014-04-10 22:30:09