| 查看: 2160 | 回复: 0 | ||
[求助]
MATLAB 中的kron()函数怎么精简？？？新手上路多多关照！！！
|
|
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); 代码如上，可是计算量太大，笔记本电脑几乎要瘫了！求高手交给我怎么才能减少计算量呢？老师说可以用一个小函数代替那么多直积，可我还是不会！！！跪求大神指点！！！ |
回复此楼» 猜你喜欢
散金币祈福
已经有187人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
-
物理学I论文润色/翻译怎么收费?
已经有191人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有114人回复
-
深圳大学2027级光电信息工程/物理学研究生推免
已经有70人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有4人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有257人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有39人回复