24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lhfx_313

至尊木虫 (文坛精英)

应助: 1 (幼儿园)
贵宾: 0.155
金币: 29747.1
散金: 2852
红花: 20
沙发: 102
帖子: 11130
在线: 1568.7小时
虫号: 550265
注册: 2008-04-23
专业: 通信理论与系统

[交流] 【求助】matlab级数求和问题已有2人参与

只有X为未知数，x取值范围为[50:1:700]，目的是画出pr随X变化的图形。
但计算组合数时可能会出现溢出问题，导致计算结果错误，该怎么处理呢？

回复此楼

» 猜你喜欢

售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有7人回复
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急已经有8人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有6人回复
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急已经有7人回复
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急已经有6人回复
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急已经有6人回复
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急已经有7人回复
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急已经有5人回复
售SCI文章，我:8O.5.5.1O.54,科目全,可十急已经有5人回复
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急已经有6人回复

1楼 2010-04-09 11:38:22

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

lhfx_313

至尊木虫 (文坛精英)

应助: 1 (幼儿园)
贵宾: 0.155
金币: 29747.1
散金: 2852
红花: 20
沙发: 102
帖子: 11130
在线: 1568.7小时
虫号: 550265
注册: 2008-04-23
专业: 通信理论与系统

引用回帖:

Originally posted by hopingzou at 2010-04-11 10:02:38:
看了一下你的程序，觉得没有编错，只是gamma函数求的值太大了，数值大小超过了计算机所容忍的最大值，如果取精度的话，这本身就是一种错误。所以，能否自己简化一下表达式，再编程。仅供参考。

我也觉得应该是错在这里了，用斯特灵公式简化了下感觉还是不对谢谢你噢

赞一下

回复此楼

7楼2010-04-11 10:34:39

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

查看全部 8 个回答

hopingzou

金虫 (著名写手)

应助: 0 (幼儿园)
金币: 2332.3
散金: 256
红花: 1
帖子: 1130
在线: 432.4小时
虫号: 832201
注册: 2009-08-22
性别: GG
专业: 机械动力学

★
小木虫(金币+0.5):给个红包，谢谢回帖交流

觉得数组不大，理论上应该不会有溢出，不妨把错误贴上！

赞一下(1人)

回复此楼

2楼2010-04-09 12:38:36

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

lhfx_313

至尊木虫 (文坛精英)

应助: 1 (幼儿园)
贵宾: 0.155
金币: 29747.1
散金: 2852
红花: 20
沙发: 102
帖子: 11130
在线: 1568.7小时
虫号: 550265
注册: 2008-04-23
专业: 通信理论与系统

引用回帖:

Originally posted by hopingzou at 2010-04-09 12:38:36:
觉得数组不大，理论上应该不会有溢出，不妨把错误贴上！

t=4;
S=50;
L=t/S;
p=zeros(1,66);
Ye=0;
i=1;
for x=50:10:700

for k=50:1:x
Ye=(gamma(x+1)/(gamma(k+1)*(gamma(x-k+1))))*((t/S)^k)*
((1-t/S)^(x-k))+Ye;
end
p(1,i)=S*Ye;
i=i+1;
Y=0;
end
A=[50:10:700];
plot(A,p);
可以出来图形，如下，但感觉不对，是不是程序出错了？