| 查看: 438 | 回复: 2 | ||
[求助]
求 crazy climber的完整程序~有哪个大神指导一下，不胜感激！！
|
| 要crazy climber的完整程序，不是只有out=climber（c）的那个程序~包括后面的脊线链接 什么的都要~求助~~ |
回复此楼» 猜你喜欢
319求调剂
已经有4人回复
-
081700化学工程与技术 一志愿中海洋 323 求调剂学校
已经有15人回复
-
085600，320分求调剂
已经有14人回复
-
085600调剂
已经有3人回复
-
求调剂：085600材料与化工，考材科基，总分319
已经有26人回复
-
0703化学
已经有5人回复
-
085600，专业课化工原理，320分求调剂
已经有8人回复
-
一志愿南农090401，268，求调剂
已经有3人回复
-
材料334求调剂
已经有11人回复
-
323求调剂
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 哪位大神来指导下数字信号处理怎么学，感激不尽
已经有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
