版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

愛(ài)紅茶的貓

木蟲(chóng) (小有名氣)

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

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

求助一個(gè)matlab編程的優(yōu)化問(wèn)題，類似霍夫變換（Hough Transform）的累加問(wèn)題。
一維向量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

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

回復(fù)此樓

» 猜你喜歡

0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章已經(jīng)有15人回復(fù)
一志愿武漢理工材料工程專碩調(diào)劑已經(jīng)有5人回復(fù)
0856調(diào)劑，是學(xué)校就去已經(jīng)有7人回復(fù)
申博26年已經(jīng)有3人回復(fù)
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑已經(jīng)有10人回復(fù)
0703化學(xué)調(diào)劑，六級(jí)已過(guò)，有科研經(jīng)歷已經(jīng)有12人回復(fù)
294求調(diào)劑材料與化工專碩已經(jīng)有13人回復(fù)
288求調(diào)劑已經(jīng)有14人回復(fù)
307求調(diào)劑已經(jīng)有8人回復(fù)
一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑已經(jīng)有5人回復(fù)

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

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

feixiaolin

榮譽(yù)版主 (文壇精英)

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

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

★
感謝參與，應(yīng)助指數(shù) +1
愛(ài)紅茶的貓: 金幣+1, ★有幫助, 可是這段程序還是基于雙重for循環(huán)的，實(shí)際效率并沒(méi)有提高多少，不知道哪位可以給出類似于將a， b構(gòu)造成2維矩陣，然后整體進(jìn)行運(yù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

銅蟲(chóng) (初入文壇)

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

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

★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
愛(ài)紅茶的貓: 金幣+5, ★★★很有幫助, 經(jīng)過(guò)驗(yàn)證，A矩陣如果為零矩陣的話沒(méi)有誤差，但是方法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)過(guò)試驗(yàn)，D的元素有時(shí)候有數(shù)值誤差，但是都在1.0e-14量級(jí)，稍微處理就行

贊一下

回復(fù)此樓

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

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

愛(ài)紅茶的貓

木蟲(chóng) (小有名氣)

應(yīng)助: 4 (幼兒園)
金幣: 5195.6
散金: 818
紅花: 3
帖子: 162
在線: 67.4小時(shí)
蟲(chóng)號(hào): 2130898
注冊(cè): 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ù)，能不能請(qǐng)您給解釋一下方法二里
counta=diff(ida)
countb=diff(idb)
B(a,b)=B(a,b)+counta'*countb
這幾行的含義？
我驗(yàn)證了一下，如果A是零矩陣的話2中方法是沒(méi)有誤差的。
謝謝！

贊一下

回復(fù)此樓

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

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

yisc1981

銅蟲(chóng) (初入文壇)

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

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

引用回帖:

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

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

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

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

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

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

贊一下(1人)

回復(fù)此樓

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

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

愛(ài)紅茶的貓

木蟲(chóng) (小有名氣)

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

引用回帖:

5樓: Originally posted by yisc1981 at 2014-11-28 00:13:53
理解編程思想就不難理解那幾句話了，至于那些函數(shù)是什么意思可以查詢相關(guān)的說(shuō)明，這里將編程思想說(shuō)說(shuō)吧。
因?yàn)閍,b中存儲(chǔ)的是A矩陣的行標(biāo)和列標(biāo)，如果那個(gè)下標(biāo)出現(xiàn)的次數(shù)多，則A中對(duì)應(yīng)元素被加的量就大。所以 ...

非常感謝您的回復(fù)，非常清楚，學(xué)習(xí)了！

贊一下

回復(fù)此樓

6樓2014-11-28 11:10:08

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主愛(ài)紅茶的貓的主題更新

返回列表

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

[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +7	Ymlll 2026-03-18	10/500	2026-03-19 19:43 by maocaozhuxi
[考研] 0703化學(xué)調(diào)劑 +4	18889395102 2026-03-18	4/200	2026-03-19 16:13 by 30660438
[考研] 0703化學(xué)調(diào)劑 +5	pupcoco 2026-03-17	8/400	2026-03-19 13:58 by houyaoxu
[考研] 346求調(diào)劑[0856] +3	WayneLim327 2026-03-16	6/300	2026-03-19 11:21 by WayneLim327
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4	我不是一根蔥 2026-03-18	4/200	2026-03-19 09:11 by 浮云166
[考研] 0703化學(xué)調(diào)劑 +8	妮妮ninicgb 2026-03-15	12/600	2026-03-19 05:55 by anny19840123
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8	Ncdx123456 2026-03-13	9/450	2026-03-18 14:40 by haxia
[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢(mèng)人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 收復(fù)試調(diào)劑生 +4	雨后秋荷 2026-03-18	4/200	2026-03-18 14:16 by elevennnne
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +7	ZY-05 2026-03-13	8/400	2026-03-18 14:13 by 007_lilei
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 一志愿南京大學(xué)，080500材料科學(xué)與工程，調(diào)劑 +4	Jy? 2026-03-16	4/200	2026-03-17 11:02 by gaoqiong
[考研] 275求調(diào)劑 +4	太陽(yáng)花天天開(kāi)心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 東南大學(xué)364求調(diào)劑 +5	JasonYuiui 2026-03-15	5/250	2026-03-16 21:28 by 木瓜膏
[考研] 289求調(diào)劑 +4	這么名字咋樣 2026-03-14	6/300	2026-03-14 18:58 by userper
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs