| 查看: 512 | 回复: 2 | |||
[求助]
求 crazy climber的完整程序~有哪个大神指导一下，不胜感激！！
|
| 要crazy climber的完整程序，不是只有out=climber（c）的那个程序~包括后面的脊线链接 什么的都要~求助~~ |
回复此楼» 猜你喜欢
哪位老哥知道今年的国自然具体哪一天放榜？
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有72人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有26人回复
-
我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人
已经有16人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有12人回复
-
招聘 金属平磨液，抛光液研发工程师
已经有3人回复
-
时间戳又变了8-15
已经有23人回复
-
filecode=后面第一个是大写字母
已经有7人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 哪位大神来指导下数字信号处理怎么学，感激不尽
已经有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
