【正文】 020.[6] 劉海軍,彭紹雄,高傳斌,[J]. 兵工自動(dòng)化，2012年4期.[7] 苗京,黃紅星,程衛(wèi)生,[J]. 武漢大學(xué)數(shù)學(xué)學(xué)院，2005年第05期.[8] 于勇,[J]. 電子與信息學(xué)報(bào)，2008,28(5):1236~1239.[9] 朱玲,施心陵,劉亞杰,[J]. 計(jì)算機(jī)工程，2006,32(24):178~179,239.[10] [D]. 成都：電子科技大學(xué)，2010.[11] 高德威,陳天煌,[J]. 武漢理工大學(xué)計(jì)算機(jī)科學(xué)與技術(shù)學(xué)院, 2009年第01期.[12] 陳昌敏,謝維成,[J]. 西華大學(xué)學(xué)報(bào)（自然科學(xué)版）, 2011年3期.[13] 劉濤,[J]. 科技創(chuàng)新導(dǎo)報(bào), 2010年14期.[14] , . A Regularized Solution To Edge of plexity, 1988(4):106123.[15] [J]. 現(xiàn)代計(jì)算機(jī)（專業(yè)版）, 2013年17期.[16] [D]. 武漢理工大學(xué)，2013.[17] [D]. 長(zhǎng)安大學(xué)，2007.[18] [D]. 天津大學(xué)，2005.[19] 郭金泉,伍宏濤,[J]. 通信市場(chǎng), 2008(3):74~78.[20] 吳小峰,[J]. 福建電腦, 2009(4):39~40.[21] 徐紅梅,陳義保,劉加光,[J]. 山東理工大學(xué)學(xué)報(bào), 2008年1月第22卷第一期.function mainclose all。 clear all。 clc。img = double(imread(39。H:\baby\39。))./255。[nrow, ncol] = size(img)。 /讀取行和列/fprintf(39。Wele to demo program.\nPlease wait......\n39。)。v = zeros(size(img))。v_norm = 0。 /置零/for rr =1:nrow for cc=1:ncol %defination of clique temp1 = [rr2 cc1。 rr2 cc+1。 rr1 cc2。 rr1 cc1。 rr1 cc。 rr1 cc+1。 rr1 cc+2。 rr cc1]。 /取8個(gè)像素/ temp2 = [rr+2 cc+1。 rr+2 cc1。 rr+1 cc+2。 rr+1 cc+1。 rr+1 cc。 rr+1 cc1。 rr+1 cc2。 rr cc+1]。 temp0 = find(temp1(:,1)=1 amp。 temp1(:,1)=nrow amp。 temp1(:,2)=1 amp。 temp1(:,2)=ncol amp。 temp2(:,1)=1 amp。 temp2(:,1)=nrow amp。 temp2(:,2)=1 amp。 temp2(:,2)=ncol)。 temp11 = temp1(temp0, :)。 temp22 = temp2(temp0, :)。 temp00 = zeros(size(temp11,1))。 for kk = 1:size(temp11,1) temp00(kk) = abs(img(temp11(kk,1), temp11(kk,2))img(temp22(kk,1), temp22(kk,2)))。 end if size(temp11,1) == 0 v(rr, cc) = 0。 v_norm = v_norm + v(rr, cc)。 else lambda = 10。 temp00 = sin(pi .* temp00./2./lambda)。 v(rr, cc) = sum(sum(temp00.^2))。 v_norm = v_norm + v(rr, cc)。 end endendv = v./v_norm。 v = v.*100。% pheromone function initializationp = .* ones(size(img))。 %paramete setting /圖像預(yù)處理/alpha = 10。 beta = 。 rho = 。 phi = 。 ant_total_num = round(sqrt(nrow*ncol))。ant_pos_idx = zeros(ant_total_num, 2)。 % record the location of ant% initialize the positions of antsrand(39。state39。, sum(clock))。temp = rand(ant_total_num, 2)。ant_pos_idx(:,1) = round(1 + (nrow1) * temp(:,1))。 %row indexant_pos_idx(:,2) = round(1 + (ncol1) * temp(:,2))。 %column indexsearch_clique_mode = 39。839。memory_length = 40。% record the positions in ant39。s memory, convert 2D positionindex (row, col) into /二維轉(zhuǎn)一維行列/% 1D positionindexant_memory = zeros(ant_total_num, memory_length)。total_step_num = 300。total_iteration_num = 4。for iteration_idx = 1: total_iteration_num %record the positions where ant have reached in the last 39。memory_length39。 iterations delta_p = zeros(nrow, ncol)。 for step_idx = 1: total_step_num delta_p_current = zeros(nrow, ncol)。 for ant_idx = 1:ant_total_num ant_current_row_idx = ant_pos_idx(ant_idx,1)。 ant_current_col_idx = ant_pos_idx(ant_idx,2)。 % find the neighborhood of current position if search_clique_mode == 39。439。 rr = ant_current_row_idx。 cc = ant_current_col_idx。 ant_search_range_temp = [rr1 cc。 rr cc+1。 rr+1 cc。 rr cc1]。 elseif search_clique_mode == 39。839。 rr = ant_current_row_idx。 cc = ant_current_col_idx。 ant_search_range_temp = [rr1 cc1。 rr1 cc。 rr1 cc+1。 rr cc1。 rr cc+1。 rr+1 cc1。 rr+1 cc。 rr+1 cc+1]。 end %remove the positions our of the image39。s range temp = find(ant_search_range_temp(:,1)=1 amp。 ant_search_range_temp(:,1)=nrow amp。 ant_search_range_temp(:,2)=1 amp。 ant_search_range_temp(:,2)=ncol)。 ant_search_range = ant_search_range_temp(temp, :)。 %calculate the transit prob. to the neighborhood of current %position ant_transit_prob_v = zeros(size(ant_search_range,1),1)。 ant_transit_prob_p = zeros(size(ant_search_range,1),1)。 for kk = 1:size(ant_search_range,1) temp = (ant_search_range(kk,1)1)*ncol + ant_search_range(kk,2)。 if length(find(ant_memory(ant_idx,:)==temp))==0 %not in ant39。s memory ant_transit_prob_v(kk) = v(ant_search_range(kk,1), ant_search_range(kk,2))。 ant_transit_prob_p(kk) = p(ant_search_range(kk,1), ant_search_range(kk,2))。 else %in ant39。s memory ant_transit_prob_v(kk) = 0。 ant_transit_prob_p(kk) = 0。 end end % if all neighborhood are in memory, then the permissible search range is REcalculated. if (sum(sum(ant_transit_prob_v))==0) || (sum(sum(ant_transit_prob_p))==0) for kk = 1:size(ant_search_range,1) temp = (ant_search_range(kk,1)1)*ncol + ant_search_range(kk,2)。 ant_transit_prob_v(kk) = v(ant_search_range(kk,1), ant_search_range(kk,2))。 ant_transit_prob_p(kk) = p(ant_search_range(kk,1), ant_search_range(kk,2))。 end end ant_transit_prob = (ant_transit_prob_v.^alpha) .* (ant_transit_prob_p.^beta) ./ (sum(sum(((ant_transit_prob_v+eps).^alpha) .* (ant_transit_prob_p.^beta))))。 % generate a random number to determine the next position. rand(39。state39。, sum(100*clock))。 temp = find(cumsum(ant_transit_prob)=rand(1), 1)。 ant_next_row_idx = ant_search_range(temp,1)。 ant_next_col_idx = ant_search_range(temp,2)。 if length(ant_next_row_idx) == 0 ant_next_row_idx = ant_cu