QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 9743|回复: 8
打印 上一主题 下一主题

[代码资源] matlab排列组合问题

[复制链接]
字体大小: 正常 放大
tansuozhe 实名认证       

1

主题

0

听众

7

积分

升级  2.11%

该用户从未签到

自我介绍
地方第三方士大夫多少
跳转到指定楼层
1#
发表于 2009-12-26 14:18 |只看该作者 |倒序浏览
|招呼Ta 关注Ta
各位老师:8 v5 Z1 J+ k" }2 S! o5 h9 M- M3 T4 F
     我是初学matlab。现在遇到一个难题,用matlab的函数combntns()* @( |* J( T& l. Z: L$ Y4 ~
求1-210共210个自然数的组合,我算从210个自然数中取3个进行组合(不重复),
. X. G' ], R8 q5 H' a在matlab命令窗口输入如下函数:
/ R8 G/ m* j" Y》combntns(1:210,3)
( s: A/ M- @$ k) Z& e. B  O/ t运行计算都正常的列出了所有组合数,可是想算从210个自然数中取20个进行组合(不重复)" F% S1 C9 G( _& U0 \" v
我知道这样的组合数非常大为:4.49333E+27,我在在matlab命令窗口输入如下函数:
& B5 i  {! W2 v4 y- @+ f; e% s" Z》combntns(1:210,20)
! l1 y6 b8 }1 x* u& J0 c( @提示出错信息如下:! h+ y2 j8 `( c, Y" o
Error in ==> combntns at 71' ^  a6 W" W. q0 A
  out=[out; choicevec(i)*ones(size(tempout,1),1) tempout], z. d6 m$ ]5 k: R2 `0 p1 y
意思是内存不够吧,现在想请教一下几个问题:
- y; A( x' E" n8 k4 |% _1  提示的内存不够是指内存条不够大,还是硬盘空间不够,这样的大的求组合数,matlab是9 t3 y" L4 Y! o# \( D
   否有函数能胜任,如果有,请问是什么函数?3 W0 i# i+ ]7 a7 H) x' ?( v
2,如果matlab的函数不能胜任,请问我如何用M文件编程计算从210个自然数中取20个的组合数并列出,' K1 i) Q! c" R
   (我想可能一次性算出不大可能),如果分段计算,应该如何编程,烦请各位老师提供一下代码,
2 P& i5 o- w+ `2 D$ g( k$ w    期盼各位老师指教......
zan
转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信

1341

主题

738

听众

2万

积分

数学中国总编辑

  • TA的每日心情

    2016-11-18 10:46
  • 签到天数: 206 天

    [LV.7]常住居民III

    社区QQ达人 邮箱绑定达人 元老勋章 发帖功臣 新人进步奖 原创写作奖 最具活力勋章 风雨历程奖

    群组2011年第一期数学建模

    群组第一期sas基础实训课堂

    群组第二届数模基础实训

    群组2012第二期MCM/ICM优秀

    群组MCM优秀论文解析专题

    对于你的问题一:是内存不够大,非外存不够大;不管是任何函数都不行,因为一次算出所有组合数的可能性为0
    ! i. Y) @9 b, x                        所以任何函数都不可能在一次执行这么多,这是硬件资源所**的,而非软件所能左右" E) w, y" O9 w1 m+ S' ^
    对于你的问题二:你可以先将210分成7段,即:0-25,26-50,.....,201-210;! M  G9 @3 F* I2 r$ m) E1 I
                    然后按下列几步进行:
    , R; U3 s- W5 _. Y+ O(1)在第一段内搜索排列,然后在第二段内搜索排列,然后在第三段内搜索排列。。。。。。。。。这样依次搜索各段的排列,搜索完后,保存数据,进入下一步;
    5 `" e, a- K: S. J4 s(2)然后两段组合搜索排列,注意这次是在两段内搜索,所以得到的排列需要包含两端内的数据
    4 i/ k% N- t6 l) G0 P5 }8 @(3)然后三段组合
    ; c: J! M9 [1 o( }1 x% H" o$ X依次下去& a2 _' ?* U0 L8 o; p% \) Q2 J
    (。。。)最后到10段组合;; T1 F2 W2 d) V1 C$ m4 O* y1 y1 W2 a! G
    这样才能处理过来。虽然比较麻烦,但是只能这样啦,除非在大型机上就不存在这种问题啦;
    已有 1 人评分金币 收起 理由
    大笨象 + 6 写的不错--

    总评分: 金币 + 6   查看全部评分

    回复

    使用道具 举报

    1341

    主题

    738

    听众

    2万

    积分

    数学中国总编辑

  • TA的每日心情

    2016-11-18 10:46
  • 签到天数: 206 天

    [LV.7]常住居民III

    社区QQ达人 邮箱绑定达人 元老勋章 发帖功臣 新人进步奖 原创写作奖 最具活力勋章 风雨历程奖

    群组2011年第一期数学建模

    群组第一期sas基础实训课堂

    群组第二届数模基础实训

    群组2012第二期MCM/ICM优秀

    群组MCM优秀论文解析专题

    我之所以选择每段的长度是25,是因为我的机子在长度为25时,反应快,而且效果比较好
    回复

    使用道具 举报

    1341

    主题

    738

    听众

    2万

    积分

    数学中国总编辑

  • TA的每日心情

    2016-11-18 10:46
  • 签到天数: 206 天

    [LV.7]常住居民III

    社区QQ达人 邮箱绑定达人 元老勋章 发帖功臣 新人进步奖 原创写作奖 最具活力勋章 风雨历程奖

    群组2011年第一期数学建模

    群组第一期sas基础实训课堂

    群组第二届数模基础实训

    群组2012第二期MCM/ICM优秀

    群组MCM优秀论文解析专题

    这么大量的计算和数据,用小型机计算有点太面前啦,应该用大型机计算比较合适: W2 U$ U3 ^: k' n4 y
    : z) |5 f# \$ ?: x6 h! X$ n
    我试了个每段长为27,我的机子都要计算将近两分钟(我的内存2G),更不用说你把所有的都计算完,估计至少得一天吧,而且还没有休息
    回复

    使用道具 举报

    0

    主题

    0

    听众

    9

    积分

    升级  4.21%

    该用户从未签到

    自我介绍
    大家好!我是酱油王!
    很好 顶!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    回复

    使用道具 举报

    兄弟 实名认证       

    0

    主题

    5

    听众

    556

    积分

    升级  85.33%

    该用户从未签到

    新人进步奖

    群组西南大学建模组

    顶!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    回复

    使用道具 举报

    hupanfeng 实名认证       

    0

    主题

    3

    听众

    149

    积分

    升级  24.5%

    该用户从未签到

    自我介绍
    hello~大家好~
    回复

    使用道具 举报

    0

    主题

    4

    听众

    304

    积分

    升级  1.33%

  • TA的每日心情
    开心
    2014-1-1 12:58
  • 签到天数: 16 天

    [LV.4]偶尔看看III

    回复

    使用道具 举报

    0

    主题

    4

    听众

    304

    积分

    升级  1.33%

  • TA的每日心情
    开心
    2014-1-1 12:58
  • 签到天数: 16 天

    [LV.4]偶尔看看III

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-11 01:04 , Processed in 0.494378 second(s), 96 queries .

    回顶部