24小时热门版块排行榜    

查看: 632  |  回复: 4

kexiaolong

金虫 (正式写手)

[求助] 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

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mathstudy

金虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
循环一下....b(i)=Integrate(f(x),0,a(i))).

» 本帖已获得的红花(最新10朵)

2楼2014-04-10 09:38:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kexiaolong

金虫 (正式写手)

送红花一朵
引用回帖:
2楼: Originally posted by mathstudy at 2014-04-10 09:38:34
循环一下....b(i)=Integrate(f(x),0,a(i))).

我对matlab 编程不在行啊,能否给出求解过程呢?谢谢!
3楼2014-04-10 09:51:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mathstudy

金虫 (正式写手)

引用回帖:
3楼: Originally posted by kexiaolong at 2014-04-10 09:51:23
我对matlab 编程不在行啊,能否给出求解过程呢?谢谢!...

大体上是这样 备注:运行中会有提示说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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mathstudy

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
kexiaolong: 金币+30, ★★★★★最佳答案, 谢谢帮助! 2014-04-12 12:49:45
引用回帖:
2楼: Originally posted by mathstudy at 2014-04-10 09:38:34
循环一下....b(i)=Integrate(f(x),0,a(i))).

大体上是这样 备注:运行中会有提示说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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kexiaolong 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 07:26 by C79jjtjAKjEn
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:38 by 5BDX0d0WFp7t
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 05:50 by 5BDX0d0WFp7t
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-15 05:49 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 05:22 by 5BDX0d0WFp7t
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 05:13 by 5BDX0d0WFp7t
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-15 05:00 by mibUvS8DDCwf
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 01:21 by mibUvS8DDCwf
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-14 23:24 by vZfe6xYu34yj
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:13 by vZfe6xYu34yj
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:36 by vZfe6xYu34yj
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +4 两三岁ss 2026-09-14 4/200 2026-09-14 20:29 by 班主任的同桌
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 QUjhNVAcOSff 2026-09-13 3/150 2026-09-14 15:09 by Aj1rhIDL5ixY
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见