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变化的图形。
但计算组合数时可能会出现溢出问题，导致计算结果错误，该怎么处理呢？

回复此楼

» 猜你喜欢

参与限项已经有3人回复
假如你的研究生提出不合理要求已经有7人回复
实验室接单子已经有4人回复
全日制（定向）博士已经有4人回复
对氯苯硼酸纯化已经有3人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有12人回复
不自信的我已经有12人回复
所感已经有4人回复
要不要辞职读博？已经有7人回复
北核录用已经有3人回复

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-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);
可以出来图形，如下，但感觉不对，是不是程序出错了？