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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » 【求助】matlab级数求和问题
51/1返回列表
查看: 639  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

lhfx_313

至尊木虫 (文坛精英)

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



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

» 猜你喜欢
1楼 2010-04-09 11:38:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hopingzou

金虫 (著名写手)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
adu886886(金币+1):谢谢交流，欢迎常来仿真模拟板块 2010-04-11 12:00
看了一下你的程序，觉得没有编错，只是gamma函数求的值太大了，数值大小超过了计算机所容忍的最大值，如果取精度的话，这本身就是一种错误。所以，能否自己简化一下表达式，再编程。仅供参考。
赞一下(2人) 回复此楼
6楼2010-04-11 10:02:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

hopingzou

金虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
觉得数组不大，理论上应该不会有溢出，不妨把错误贴上！
赞一下(1人) 回复此楼
2楼2010-04-09 12:38:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lhfx_313

至尊木虫 (文坛精英)
引用回帖:
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);
可以出来图形，如下 ，但感觉不对，是不是程序出错了？
赞一下 回复此楼
3楼2010-04-09 14:48:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hopingzou

金虫 (著名写手)
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
nono2009(金币+2):鼓励应助！ 2010-04-10 07:59
大致看了下，只是提个建议，或许不是知识层面上的。
赞一下(2人) 回复此楼
4楼2010-04-10 00:40:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 国社科又开始会评了，不知道这次命运如何 +3 雨打竹帘 2026-08-30 4/200 2026-08-30 12:35 by 白大多亮
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 9/450 2026-08-30 11:41 by l0VvVHGBGRLv
[公派出国] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +6 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 11:32 by l0VvVHGBGRLv
[论文投稿] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 6/300 2026-08-30 11:07 by l0VvVHGBGRLv
[基金申请] 面上意见出来了 +9 黄鸟于飞Chao 2026-08-29 17/850 2026-08-30 10:28 by 孤独的英雄6
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:28 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 4/200 2026-08-30 07:17 by ZPa0EcMwuECS
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 5/250 2026-08-30 07:05 by ZPa0EcMwuECS
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 06:22 by ZPa0EcMwuECS
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-30 05:48 by ZPa0EcMwuECS
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +7 ASdOkHsho7FD 2026-08-28 11/550 2026-08-30 05:40 by ZPa0EcMwuECS
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 基金不中，共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
信息提示
关闭
请填处理意见
关闭