查看: 570  |  回复: 3

shiwen5006

铁杆木虫 (正式写手)

奔跑的蜗牛

[求助] 新手重金求高手解一积分问题

我想求表达式f=3400*0.71*10^(-5)/(1-exp(-(7000-x)/x))的积分
程序如下:
z=zeros(70,1);
t=100:100:7000;
for m=1:70
f=@(x)3400*0.71*10^(-5)/(1-exp(-(7000-x)/x));
z(m,1)=vpa(int('f',0,t(m)));
end
程序能运行,但是搞笑的是我无论怎么改变f的表达式,z的值都不变;
而且即使没有f的表达式,仍然能得到z的值。
求高手解答~~~~
回复此楼
低调做人,高调做事
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zsy312

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
shiwen5006: 金币+6, ★★★很有帮助, 非常感谢! 2012-10-28 17:55:41
xiegangmai: 金币+2, 谢谢应助 2012-10-29 19:44:25
使用自带的maple求解, 我用的是7.0:
CODE:
str=(maple('seq(evalf(int[numerical](3400*0.71*10^(-5)/(1-exp(-(7000-x)/x)),x=0...100*k)),k=1..69)'));
num=str2num(str)'

输出:
CODE:
num =

    2.4140
    4.8280
    7.2420
    9.6560
   12.0700
   14.4840
   16.8982
   19.3128
   21.7286
   24.1468
   26.5693
   28.9984
   31.4371
   33.8885
   36.3564
   38.8445
   41.3568
   43.8974
   46.4706
   49.0807
   51.7319
   54.4286
   57.1754
   59.9769
   62.8377
   65.7626
   68.7567
   71.8251
   74.9732
   78.2067
   81.5315
   84.9538
   88.4803
   92.1182
   95.8750
   99.7588
  103.7785
  107.9435
  112.2642
  116.7518
  121.4188
  126.2786
  131.3463
  136.6385
  142.1736
  147.9726
  154.0585
  160.4579
  167.2008
  174.3214
  181.8595
  189.8609
  198.3794
  207.4787
  217.2346
  227.7385
  239.1023
  251.4646
  265.0003
  279.9343
  296.5627
  315.2865
  336.6678
  361.5298
  391.1486
  427.6622
  475.0783
  542.3943
  658.3197

2楼2012-10-28 09:48:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zsy312

银虫 (小有名气)

注意2008以后的版本就不支持maple了, 除非你有maple too box for matlab.
如果在matlab中输入:
CODE:
z=zeros(70,1);
t=100:100:7000;
syms x
f=@(x)3400*0.71*10^(-5)./(1-exp(-(7000-x)./x));
for m=1:70
    z(m,1)=double(int(f(x),0,t(m)));
end

输出也没有问题:
CODE:
z =

    2.4140
    4.8280
    7.2420
    9.6560
   12.0700
   14.4840
   16.8982
   19.3128
   21.7286
   24.1468
   26.5693
   28.9984
   31.4371
   33.8885
   36.3564
   38.8445
   41.3568
   43.8974
   46.4706
   49.0807
   51.7319
   54.4286
   57.1754
   59.9769
   62.8377
   65.7626
   68.7567
   71.8251
   74.9732
   78.2067
   81.5315
   84.9538
   88.4803
   92.1182
   95.8750
   99.7588
  103.7785
  107.9435
  112.2642
  116.7518
  121.4188
  126.2786
  131.3463
  136.6385
  142.1736
  147.9726
  154.0585
  160.4579
  167.2008
  174.3214
  181.8595
  189.8609
  198.3794
  207.4787
  217.2346
  227.7385
  239.1023
  251.4646
  265.0003
  279.9343
  296.5627
  315.2865
  336.6678
  361.5298
  391.1486
  427.6622
  475.0783
  542.3943
  658.3197
       Inf

3楼2012-10-28 09:57:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgdxbing121

新虫 (初入文坛)

求!!!急求!!!
f(x,t)=积分(下限0,上限0.00021t,积分函数为(exp(-k-55*x^2/k)kdk
积分内部有函数变量和积分变量 且积分上限有变量  想得到f(x,t) 在t=3600*24*2时的图像怎么弄啊!!
找到组织了
4楼2012-10-30 19:22:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 shiwen5006 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 听说今天filecode变了 +19 布布和一二 2026-08-06 33/1650 2026-08-07 01:25 by 云上清扬
[基金申请] 基金中了 +6 laoda193707 2026-08-06 6/300 2026-08-06 23:34 by dragonxp
[基金申请] filecode +10 等待解的谜 2026-08-06 14/700 2026-08-06 21:16 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +3 archvillain 2026-08-06 7/350 2026-08-06 21:05 by archvillain
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 6/300 2026-08-06 20:19 by 苏知砚
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 12/600 2026-08-06 15:11 by eddyin
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[教师之家] 咨询面上基金 +4 李长云 2026-07-31 7/350 2026-08-06 11:20 by 李长云
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +22 低垂的野花 2026-07-31 30/1500 2026-08-05 18:03 by 低垂的野花
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[论文投稿] 十年后又回来了,论文投稿求助 +3 哈哈114477 2026-08-01 3/150 2026-08-04 15:40 by tegsgjy20
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见