| 1 | 1/1 | 返回列表 |
| 查看: 2075 | 回復(fù): 0 | ||
[求助]
MATLAB 中的kron()函數(shù)怎么精簡(jiǎn)???新手上路多多關(guān)照。!
|
|
h4=-z*k*(kron(kron(kron(a1x,a2),I3),I4)... +kron(kron(kron(a1,a2x),I3),I4)... +kron(kron(kron(I1,a2x),a3),I4)... +kron(kron(kron(I1,a2),a3x),I4)... +kron(kron(kron(I1,I2),a3x),a4)... +kron(kron(kron(I1,I2),a3),a4x)... +kron(kron(kron(a1x,I2),I3),a4)... +kron(kron(kron(a1,I2),I3),a4x))... +(-mu)*kron(kron(kron(n1a,I2),I3),I4)-(dta+mu)*kron(kron(kron(n1b,I2),I3),I4)... +beta*(kron(kron(kron(b1xa,I2),I3),I4)+kron(kron(kron(a1xb,I2),I3),I4))... +(-mu)*kron(kron(kron(I1,n2a),I3),I4)-(dta+mu)*kron(kron(kron(I1,n2b),I3),I4)... +beta*(kron(kron(kron(I1,b2xa),I3),I4)+kron(kron(kron(I1,a2xb),I3),I4))... +(-mu)*kron(kron(kron(I1,I2),n3a),I4)-(dta+mu)*kron(kron(kron(I1,I2),n3b),I4)... +beta*(kron(kron(kron(I1,I2),b3xa),I4)+kron(kron(kron(I1,I2),a3xb),I4))... +(-mu)*kron(kron(kron(I1,I2),I3),n4a)-(dta+mu)*kron(kron(kron(I1,I2),I3),n4b)... +beta*(kron(kron(kron(I1,I2),I3),b4xa)+kron(kron(kron(I1,I2),I3),a4xb))... +z*v*(kron(kron(kron(n1b,n2b),I3),I4)+kron(kron(kron(I1,n2b),n3b),I4)... +kron(kron(kron(I1,I2),n3b),n4b)+kron(kron(kron(n1b,I2),I3),n4b)); h12=-z*k*((kron(kron(kron((a1x+a1),I2),I3),I4))*p2... +(kron(kron(kron(I1,(a2x+a2)),I3),I4))*p1-2*p1*p2*II)... +z*v*(kron(kron(kron(n1b,I2),I3),I4)*k2+kron(kron(kron(I1,n2b),I3),I4)*k1-k1*k2*II); h23=-z*k*((kron(kron(kron(I1,(a2x+a2)),I3),I4))*p3... +(kron(kron(kron(I1,I2),(a3x+a3)),I4))*p2-2*II*p2*p3)... +z*v*(kron(kron(kron(I1,n2b),I3),I4)*k3+kron(kron(kron(I1,I2),n3b),I4)*k2-II*k2*k3); h34=-z*k*((kron(kron(kron(I1,I2),(a3x+a3)),I4))*p4... +(kron(kron(kron(I1,I2),I3),(a4x+a4)))*p3-2*II*p3*p4)... +z*v*(kron(kron(kron(I1,I2),n3b),I4)*k4+kron(kron(kron(I1,I2),I3),n4b)*k3-II*k3*k4); h41=-z*k*((kron(kron(kron(I1,I2),I3),(a4x+a4))*p1... +(kron(kron(kron((a1x+a1),I2),I3),I4)))*p4-2*II*p4*p1)... +z*v*(kron(kron(kron(I1,I2),I3),n4b)*k1+kron(kron(kron(n1b,I2),I3),I4)*k4-II*k4*k1); hm=sparse(h12+h23+h34+h41); 代碼如上,可是計(jì)算量太大,筆記本電腦幾乎要癱了!求高手交給我怎么才能減少計(jì)算量呢?老師說(shuō)可以用一個(gè)小函數(shù)代替那么多直積,可我還是不會(huì)。!跪求大神指點(diǎn)。! |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +7 | zjl的號(hào) 2026-03-16 | 12/600 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +3 | 困于星晨 2026-03-17 | 3/150 |
|
|
[考研] 一志愿985,本科211,0817化學(xué)工程與技術(shù)319求調(diào)劑 +5 | Liwangman 2026-03-15 | 5/250 |
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6 | zlingli 2026-03-13 | 6/300 |
|
|
[考研] 材料工程327求調(diào)劑 +3 | xiaohe12w 2026-03-11 | 3/150 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7 | mysdl 2026-03-11 | 9/450 |
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
[考研] [0860]321分求調(diào)劑,ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
[考研] 一志愿西南交大,材料專碩317求調(diào)劑 +5 | lx8568 2026-03-11 | 5/250 |
|
|
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|
|
[考研] 283求調(diào)劑,材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|