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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 新手重金求高手解一积分问题
41/1返回列表
查看: 528  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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的值。
求高手解答~~~~
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
低调做人，高调做事
1楼 2012-10-27 22:31:26
已阅   回复此楼   关注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
赞一下(1人) 回复此楼
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 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭