愛紅茶的貓

木蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 5195.6
散金: 818
紅花: 3
帖子: 162
在線: 67.4小時
蟲號: 2130898
注冊: 2012-11-16
專業(yè): 信息處理方法與技術(shù)

[求助] MATLAB編程優(yōu)化問題已有2人參與

求助一個matlab編程的優(yōu)化問題，類似霍夫變換（Hough Transform）的累加問題。
一維向量a里面為矩陣A的行下標(biāo)信息，一維向量b里面是A的列下標(biāo)信息。
使用for循環(huán)的賦值程序如下：

for i = 1 ：m
   for j = 1 : n
      count = A(a(i), b(j)) + 1;
      A(a(i), b(j)) = count;
   end
end

問題是matlab的for循環(huán)效率低下，如果m和n的值非常大的話（>10000），計算時間很長。
在此求助matlab高手，能否使用矩陣實現(xiàn)的方法代替for循環(huán)來提高編程效率？
期待確實有幫助的應(yīng)助回帖

回復(fù)此樓

» 猜你喜歡

346求調(diào)劑[0856] 已經(jīng)有7人回復(fù)
一志愿山大07化學(xué) 332分四六級已過本科山東雙非求調(diào)劑！已經(jīng)有3人回復(fù)
310求調(diào)劑已經(jīng)有3人回復(fù)
機(jī)械專碩299求調(diào)劑至材料已經(jīng)有4人回復(fù)
070300化學(xué)319求調(diào)劑已經(jīng)有7人回復(fù)
08工科 320總分求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有12人回復(fù)
307求調(diào)劑已經(jīng)有3人回復(fù)
265求調(diào)劑已經(jīng)有3人回復(fù)
301求調(diào)劑已經(jīng)有10人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

求matlab編程代碼，求幫助，謝謝！已經(jīng)有11人回復(fù)
關(guān)于有約束非線性函數(shù)的優(yōu)化問題已經(jīng)有8人回復(fù)
有償求助會用matlab對粒子群算法編程的同學(xué) 已經(jīng)有5人回復(fù)
【課件】最優(yōu)化方法及其Matlab程序設(shè)計已經(jīng)有390人回復(fù)
LARS算法問題求助已經(jīng)有6人回復(fù)
matlab編程約束條件問題已經(jīng)有10人回復(fù)
我用matlab編寫了一個風(fēng)力機(jī)葉片的優(yōu)化設(shè)計程序，有點問題已經(jīng)有13人回復(fù)
請教下優(yōu)化類的算法處理帶約束問題已經(jīng)有7人回復(fù)
matlab程序中調(diào)用m文件的問題已經(jīng)有13人回復(fù)
求助：遺傳算法的matlab編程求解資源優(yōu)化問題已經(jīng)有7人回復(fù)
求助：資源優(yōu)化問題的遺傳算法matlab編程已經(jīng)有5人回復(fù)
Matlab最小二乘參數(shù)優(yōu)化已經(jīng)有3人回復(fù)
matlab編程中遇到的問題已經(jīng)有4人回復(fù)
求助 matlab交點問題已經(jīng)有19人回復(fù)
matlab 問題求助，是否計算結(jié)果錯誤，如何改正？已經(jīng)有3人回復(fù)
求助有關(guān)Matlab有約束非線性最優(yōu)化問題已經(jīng)有9人回復(fù)
求助使用MATLAB計算綠色葉片的覆蓋度的問題已經(jīng)有24人回復(fù)
求助matlab編程已經(jīng)有13人回復(fù)
【求助】用matlab最優(yōu)化方法進(jìn)行參數(shù)擬合已經(jīng)有17人回復(fù)

1樓 2014-11-24 11:31:33

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yisc1981

銅蟲 (初入文壇)

應(yīng)助: 2 (幼兒園)
金幣: 731.7
帖子: 42
在線: 16.6小時
蟲號: 2194115
注冊: 2012-12-18
性別: GG
專業(yè): 人工智能與知識工程

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
愛紅茶的貓(nono2009代發(fā)): 金幣+14 2014-11-28 17:03:05

引用回帖:

4樓: Originally posted by 愛紅茶的貓 at 2014-11-26 13:57:16
十分感謝回復(fù)，能不能請您給解釋一下方法二里
counta=diff(ida)
countb=diff(idb)
B(a,b)=B(a,b)+counta'*countb
這幾行的含義？
我驗證了一下，如果A是零矩陣的話2中方法是沒有誤差的。
謝謝！...

理解編程思想就不難理解那幾句話了，至于那些函數(shù)是什么意思可以查詢相關(guān)的說明，這里將編程思想說說吧。
   因為a,b中存儲的是A矩陣的行標(biāo)和列標(biāo)，如果那個下標(biāo)出現(xiàn)的次數(shù)多，則A中對應(yīng)元素被加的量就大。所以首先是找出a，b中所含不同的下標(biāo)到底有多少，只對這些不同下標(biāo)進(jìn)行一次尋址。找到A中對應(yīng)的元素后，加上其因該被訪問的次數(shù)就行了。從而不必要反復(fù)用for循環(huán)。

下面幾行得到a,b中到底還有多少個不同元素，
   a=sort(a) %為什么要先排序是為了后面通過diff能求出不同元素的計數(shù)，也就是a中不同元素應(yīng)該被訪問的次數(shù)。
   b=sort(b)
[a,ida]=unique(a)  %找出不同的元素，經(jīng)過這一步后a中就只含有不同元素了，不會有重復(fù)的。
[b,idb]=unique(b)

下面幾行是找出a,b不同元素出現(xiàn)的次數(shù)，也就是訪問的計數(shù)
ida=[0,ida]
idb=[0,idb]
counta=diff(ida)
countb=diff(idb)
為什么能得到計數(shù)，是因為ida記錄的是a中不同元素在a中最后一次出現(xiàn)的位置。

B(a,b)得到的是A的一個子矩陣（“上面講了，B=A)，也就是這些元素才能被你給的那些兩個for循環(huán)訪問到。

D=counta'*countb  得到一個與B(a,b)同樣大小的矩陣，他們的對應(yīng)元（行被訪問次數(shù)乘以列被訪問次數(shù)，就是其交叉點處總的訪問次數(shù)）
D(i,j)就是A(a(i),b(i))應(yīng)該被最終加上的數(shù)字，也就是你那兩個for循環(huán)中A(a(i),b(i))應(yīng)該被訪問到的次數(shù)。
我也覺得很難用語言表達(dá)清楚，如果還是不明白你一步步用紙驗算，就看出其中的奧秘了。

贊一下(1人)

回復(fù)此樓

5樓2014-11-28 00:13:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

查看全部 6 個回答

feixiaolin

榮譽版主 (文壇精英)

專家經(jīng)驗: +518
信息EPI: 3
應(yīng)助: 942 (博后)
貴賓: 1.275
金幣: 3430
散金: 58785
紅花: 532
沙發(fā): 11
帖子: 24215
在線: 2601.8小時
蟲號: 2139575
注冊: 2012-11-21
專業(yè): 光學(xué)信息獲取與處理
管轄: 數(shù)學(xué)

【答案】應(yīng)助回帖

★
感謝參與，應(yīng)助指數(shù) +1
愛紅茶的貓: 金幣+1, ★有幫助, 可是這段程序還是基于雙重for循環(huán)的，實際效率并沒有提高多少，不知道哪位可以給出類似于將a， b構(gòu)造成2維矩陣，然后整體進(jìn)行運算累加的結(jié)果。 2014-11-24 14:59:52

for i = 1 ：m
   m1=a(i);
   for j = 1 : n
      m2= b(j);
      A(m1, m2) += 1;
   end
end

贊一下(1人)

回復(fù)此樓

2樓2014-11-24 14:45:06

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yisc1981

銅蟲 (初入文壇)

應(yīng)助: 2 (幼兒園)
金幣: 731.7
帖子: 42
在線: 16.6小時
蟲號: 2194115
注冊: 2012-12-18
性別: GG
專業(yè): 人工智能與知識工程

【答案】應(yīng)助回帖

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
愛紅茶的貓: 金幣+5, ★★★很有幫助, 經(jīng)過驗證，A矩陣如果為零矩陣的話沒有誤差，但是方法2中最后5行的意思不太明白，如果能解釋一下的話一定采納為最佳答案。 2014-11-27 11:01:52

clear all
clc

A=rand(10);
B=A; %留著比較
a=[2 2 6 8 5 5 5 5];
m=length(a)
b=[1 1 1 3 3 9 9 9 4 4];
n=length(b);
for i =1:m
   for j = 1 : n
      count = A(a(i), b(j)) + 1;
      A(a(i), b(j)) = count;
   end
end
% 第二種方法
a=sort(a)
b=sort(b)
[a,ida]=unique(a)
[b,idb]=unique(b)
ida=[0,ida]
idb=[0,idb]
counta=diff(ida)
countb=diff(idb)
B(a,b)=B(a,b)+counta'*countb

%檢查兩端程序是否相等
D=A-B
%經(jīng)過試驗，D的元素有時候有數(shù)值誤差，但是都在1.0e-14量級，稍微處理就行

贊一下

回復(fù)此樓

3樓2014-11-26 09:50:41

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

愛紅茶的貓

木蟲 (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 5195.6
散金: 818
紅花: 3
帖子: 162
在線: 67.4小時
蟲號: 2130898
注冊: 2012-11-16
專業(yè): 信息處理方法與技術(shù)

引用回帖:

3樓: Originally posted by yisc1981 at 2014-11-26 09:50:41
clear all
clc

A=rand(10);
B=A; %留著比較
a=;
m=length(a)
b=;
n=length(b);
for i =1:m
   for j = 1 : n
      count = A(a(i), b(j)) + 1;
      A(a(i), b(j)) = count;
   end
...

十分感謝回復(fù)，能不能請您給解釋一下方法二里
counta=diff(ida)
countb=diff(idb)
B(a,b)=B(a,b)+counta'*countb
這幾行的含義？
我驗證了一下，如果A是零矩陣的話2中方法是沒有誤差的。
謝謝！

贊一下

回復(fù)此樓

4樓2014-11-26 13:57:16

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑 +12	yangfz 2026-03-17	12/600	2026-03-21 03:30 by JourneyLucky
[考研] 材料專碩英一數(shù)二306 +7	z1z2z3879 2026-03-18	7/350	2026-03-20 23:48 by JourneyLucky
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9	Ncdx123456 2026-03-19	9/450	2026-03-20 23:41 by lovewei0727
[考研] 考研調(diào)劑求學(xué)校推薦 +3	伯樂29 2026-03-18	5/250	2026-03-20 22:59 by JourneyLucky
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4	@taotao 2026-03-19	4/200	2026-03-20 22:14 by JourneyLucky
[考研] 086500 325 求調(diào)劑 +3	領(lǐng)帶小熊 2026-03-19	3/150	2026-03-20 18:38 by 盡舜堯1
[考博] 招收博士1-2人 +3	QGZDSYS 2026-03-18	3/150	2026-03-20 11:58 by 呱呱呱呱叫
[考研] 085600材料與化工求調(diào)劑 +6	緒幸與子 2026-03-17	6/300	2026-03-19 13:27 by houyaoxu
[考研] 311求調(diào)劑 +11	冬十三 2026-03-15	12/600	2026-03-18 14:36 by 星空星月
[考研] 生物學(xué)071000 329分求調(diào)劑 +3	我愛生物生物愛?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6	薛云鵬 2026-03-15	6/300	2026-03-17 11:05 by 學(xué)員h26Tkc
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3	木子momo青爭 2026-03-15	3/150	2026-03-17 07:52 by laoshidan
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +3	自由煎餅果子 2026-03-16	3/150	2026-03-16 14:10 by 運氣yunqi
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper