| 查看: 2037 | 回复: 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); 代码如上，可是计算量太大，笔记本电脑几乎要瘫了！求高手交给我怎么才能减少计算量呢？老师说可以用一个小函数代替那么多直积，可我还是不会！！！跪求大神指点！！！ |
回复此楼» 猜你喜欢
反铁磁体中的磁性切换：两种不同的机制已成功可视化
已经有0人回复
-
求标准粉末衍射卡号 ICDD 01-076-1802
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有128人回复
-
新西兰Robinson研究所招收全奖PhD
已经有0人回复
-
石墨烯转移--二氧化硅衬底石墨烯
已经有0人回复
-
笼目材料中量子自旋液体基态的证据
已经有0人回复
-
数学教学论硕士可以读数学物理博士吗？
已经有0人回复
-
德国亥姆霍兹Hereon中心汉堡分部招镁合金腐蚀裂变SCC课题方向2026公派博士生
已经有4人回复
-
澳门大学 应用物理及材料工程研究院 潘晖教授课题组诚招博士后
已经有11人回复