| 查看: 508 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求 crazy climber的完整程序~有哪个大神指导一下，不胜感激！！
|
||
| 要crazy climber的完整程序，不是只有out=climber（c）的那个程序~包括后面的脊线链接 什么的都要~求助~~ |
回复此楼» 猜你喜欢
filecode
已经有12人回复
-
不应该看fileCode
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有33人回复
-
FileCode能看出啥？
已经有19人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有5人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有13人回复
-
是这周出结果还是下周出结果？
已经有3人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 哪位大神来指导下数字信号处理怎么学，感激不尽
已经有30人回复
- HKU申请提交RP问题，求各位大神指教，不胜感激！
已经有10人回复
3楼2012-10-29 15:39:16
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
