| 查看: 2074 | 回复: 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); 代码如上，可是计算量太大，笔记本电脑几乎要瘫了！求高手交给我怎么才能减少计算量呢？老师说可以用一个小函数代替那么多直积，可我还是不会！！！跪求大神指点！！！ |
回复此楼» 猜你喜欢
计算三阶力常数用thirdorder扩胞时参数如何选择？临近原子数怎么取值
已经有7人回复
-
2026年循环经济功能材料国际会议（ICFMCE 2026）
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有291人回复
-
2026年第五届电气、电子与信息工程国际会议（ISEEIE 2026）
已经有0人回复
-
哈尔滨理工大学物理系招收物理学考研调剂
已经有5人回复
-
0702一志愿吉大B区求调剂
已经有5人回复
-
求调剂
已经有0人回复
-
0702一志愿吉大B区求调剂有论文
已经有0人回复
-
请问还有没有用Latex写文章的小伙伴们？
已经有0人回复
-
光学工程学硕调剂信息
已经有8人回复
