| 查看: 416 | 回复: 2 | |||
[求助]
求 crazy climber的完整程序~有哪个大神指导一下，不胜感激！！
|
| 要crazy climber的完整程序，不是只有out=climber（c）的那个程序~包括后面的脊线链接 什么的都要~求助~~ |
回复此楼» 猜你喜欢
拟解决的关键科学问题还要不要写
已经有11人回复
-
救命帖
已经有5人回复
-
限项规定
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
招博士
已经有3人回复
-
存款400万可以在学校里躺平吗
已经有35人回复
-
最失望的一年
已经有18人回复
-
求推荐博导
已经有4人回复
-
求推荐英文EI期刊
已经有5人回复
-
疑惑？
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 哪位大神来指导下数字信号处理怎么学，感激不尽
已经有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