| 查看: 417 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
cw277木虫 (正式写手)
|
[求助]
求助，请看一个matlab程序~~
|
|
|
看有没有什么地方可以优化？ clear all; fp=fopen('c:\\1.txt','w'); tic min=6; num=0; for i=1:500 l=GETBIT(i,0) + GETBIT(i,1)*sym('x1') + GETBIT(i,2)*sym('x2') + GETBIT(i,3)*sym('x3') + GETBIT(i,4)*sym('x4') + GETBIT(i,5)*sym('x5') + GETBIT(i,6)*sym('x6') + GETBIT(i,7)*sym('x7') + GETBIT(i,8)*sym('x1*x2') + GETBIT(i,9)*sym('x1*x3') + GETBIT(i,10)*sym('x1*x4') + GETBIT(i,11)*sym('x1*x5') + GETBIT(i,12)*sym('x1*x6') + GETBIT(i,13)*sym('x1*x7') + GETBIT(i,14)*sym('x2*x3') + GETBIT(i,15)*sym('x2*x4') + GETBIT(i,16)*sym('x2*x5') + GETBIT(i,17)*sym('x2*x6') + GETBIT(i,18)*sym('x2*x7') + GETBIT(i,19)*sym('x3*x4') + GETBIT(i,20)*sym('x3*x5') + GETBIT(i,21)*sym('x3*x6') + GETBIT(i,22)*sym('x3*x7') + GETBIT(i,23)*sym('x4*x5') + GETBIT(i,24)*sym('x4*x6') + GETBIT(i,25)*sym('x4*x7') + GETBIT(i,26)*sym('x5*x6') + GETBIT(i,27)*sym('x5*x7') + GETBIT(i,28)*sym('x6*x7') ; f=sym(' x1 + x1*x3+ x2*x3+x4 + x1*x4 + x2*x4 + x1*x2*x4 + x3*x4 + x1*x3*x4 + x1*x2*x3*x4 + x1*x3*x5 + x4*x5 + x1*x2*x4*x5 + x1*x2*x3*x4*x5 + x6 + x2*x6 + x1*x2*x6 + x1*x2*x3*x6 + x1*x2*x4*x6 + x1*x2*x3*x4*x6 + x1*x5*x6 + x3*x5*x6 + x1*x4*x5*x6 + x3*x4*x5*x6 + x7 + x2*x7 + x1*x2*x7 + x2*x3*x7 + x1*x4*x7+ x1*x2*x4*x7+ x1*x2*x3*x4*x7 +x5*x7 +x1*x5*x7 +x1*x3*x5*x7 +x1*x2*x3*x5*x7 +x2*x4*x5*x7 +x2*x3*x4*x5*x7 +x6*x7 +x1*x2*x6*x7 +x1*x3*x6*x7 +x1*x2*x3*x6*x7 +x2*x4*x6*x7 +x1*x3*x4*x6*x7 +x2*x3*x4*x6*x7 +x5*x6*x7 +x2*x5*x6*x7 +x1*x2*x5*x6*x7 +x2*x3*x5*x6*x7 +x1*x4*x5*x6*x7 +x3*x4*x5*x6*x7'); s=char(expand(l*f)); s=strrep(s,'^2',''); s=strrep(s,' 2*',' 0*'); s=strrep(s,' 4*',' 0*'); s=strrep(s,' 6*',' 0*'); s=strrep(s,' 8*',' 0*'); s=char(expand(sym(s))); s=strrep(s,' 2*',' 0*'); s=strrep(s,' 4*',' 0*'); s=strrep(s,' 6*',' 0*'); s=strrep(s,' 8*',' 0*'); s=strrep(s,' 10*',' 0*'); s=char(expand(sym(s))); s=strrep(s,'x1','x'); s=strrep(s,'x2','x'); s=strrep(s,'x3','x'); s=strrep(s,'x4','x'); s=strrep(s,'x5','x'); s=strrep(s,'x6','x'); s1=strrep(s,'x7','x'); fprintf(fp,'%d\r\n',i); if(poly_degree(s1)<=min&poly_degree(s1)>0) min=poly_degree(s1); num=i; end end fprintf(fp,'%d,%d\r\n',min,num); fclose(fp); toc; |
回复此楼» 猜你喜欢
确定了，国自然21号放榜
已经有4人回复
-
静等基金结果
已经有16人回复
-
filecode
已经有6人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有11人回复
-
关于Filecode分析方法
已经有3人回复
-
国自然结果
已经有8人回复
-
FileCode能看出啥？
已经有8人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有22人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有4人回复
-
8月时间戳变的，举个手。玩一下，释放压力
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助啊！求一程序，用matlab程序做，用蒙特卡洛方法模拟
已经有10人回复
- 求助matlab一程序。
已经有4人回复
- 求助matlab编程
已经有13人回复
- 【求助】请教一个matlab程序，谢谢帮忙！！
已经有15人回复
- 【求助】大家看看，matlab安装过程的截图，找出错
已经有14人回复
- 【求助】请教把java程序转化成matlab程序
已经有9人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
3楼2012-10-08 13:39:47
Toauto
至尊木虫 (著名写手)
- 应助: 23 (小学生)
- 金币: 13024.8
- 散金: 1904
- 红花: 25
- 帖子: 1925
- 在线: 417.2小时
- 虫号: 339744
- 注册: 2007-04-07
- 专业: 文献学
2楼2012-10-08 12:04:58
100