| 查看: 495 | 回复: 2 | ||
[求助]
求 crazy climber的完整程序~有哪个大神指导一下，不胜感激！！
|
| 要crazy climber的完整程序，不是只有out=climber（c）的那个程序~包括后面的脊线链接 什么的都要~求助~~ |
回复此楼» 猜你喜欢
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
-
听说今天filecode变了
已经有45人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有7人回复
-
关于豆爷回答的JTJC与%2F数量
已经有6人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有31人回复
-
filecode
已经有19人回复
-
基金中了
已经有11人回复
-
固定端突然变了，今天
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 哪位大神来指导下数字信号处理怎么学，感激不尽
已经有30人回复
- HKU申请提交RP问题，求各位大神指教，不胜感激！
已经有10人回复
18662515085
银虫 (正式写手)
- 应助: 142 (高中生)
- 金币: 419.8
- 散金: 24
- 红花: 11
- 帖子: 487
- 在线: 220.4小时
- 虫号: 1955839
- 注册: 2012-08-25
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
【答案】应助回帖
★ ★
感谢参与，应助指数 +1
xiegangmai: 金币+2, 谢谢应助 2012-10-29 20:08:21
感谢参与，应助指数 +1
xiegangmai: 金币+2, 谢谢应助 2012-10-29 20:08:21
|
%crazy climber算法 function out=ridge ( c ) %c是要提取脊线的时频矩阵,c元素为非负 [B,A]=size(c); N=floor(A*B/4); %生成N个Climber temp=zeros(B,A); for i=1:1:A*B if mod(i,4)==0 temp(i)=1; end end T=max(max(c))-min(min(c));%系统初始化温度 Tt=T; %系统当前温度 t=2; %系统当前时间 while Tt>=T/1000 %对时间t做循环 for i=4:4:A*B %对每个climber做移动 if mod(i,B)==0 heng=mod(i,B)+B; else heng=mod(i,B); end %计算climber的横坐标 zong=ceil(i/B) %计算climber的纵坐标 p=sign(2*rand-1); %横坐标以0.5的等概率分别向左和右移动 if heng==1 p=1; elseif heng==A p=-1; end %排除边界条件 heng_new=heng+p; %纵坐标按规则移动 p=sign(2*rand-1); if zong==1 p=1; elseif zong==A p=-1; end %排除边界条件 zong_new=zong+p; if c(heng_new,zong_new)>c(heng_new,zong) temp(heng,zong)=0; temp(heng_new,zong_new)=1; else pt=exp((c(heng_new,zong_new)-c(heng_new,zong))/Tt); if(rand<=pt) zong_new=zong+p; temp(heng,zong)=0; temp(heng_new,zong_new)=1; else zong_new=zong; temp(heng,zong)=0; temp(heng_new,zong_new)=1; end end end Tt=T/log2(t); t=t+1; end |
2楼2012-10-29 13:23:16
3楼2012-10-29 15:39:16
100