| 查看: 2167 | 回复: 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); 代码如上，可是计算量太大，笔记本电脑几乎要瘫了！求高手交给我怎么才能减少计算量呢？老师说可以用一个小函数代替那么多直积，可我还是不会！！！跪求大神指点！！！ |
回复此楼» 猜你喜欢
深圳大学应用技术学院招聘凝聚态物理博士后
已经有124人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
物理学I论文润色/翻译怎么收费?
已经有220人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有0人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复