| 查看: 194 | 回复: 0 | ||
fanbao_keke铁虫 (正式写手)
|
[求助]
求高人帮我修改程序
|
|
求高人帮我修改程序，想作三维图，坐标空间(teta=0:pi/16:pi,F=0:0.1:10), 第三维为E(:,m)=(QuasiEnergy(Index(1:4)))，即每对 teta 与 F 对应四个E，请问 如何作出这种图？帮帮修改下面的程序： clear all omega=30; m=1; n=1; D=1; v=0.2 for teta=0:pi/16:pi for F=0:0.1:10 J=100; T=2*pi/omega; tstep=T/J; U=eye(4,4); for t=tstep:tstep:T H_0=zeros(4,4); H_0(1,2)=-D; H_0(1,4)=-D; H_0(1,3)=-v; H_0(2,1)=-D; H_0(2,3)=-D; H_0(2,4)=-v; H_0(3,1)=-v; H_0(3,2)=-D; H_0(3,4)=-D; H_0(4,1)=-D; H_0(4,2)=-v; H_0(4,3)=-D; H=zeros(4,4); H(1,1)=-(F*cos(teta)*omega*cos(omega*t)); H(2,2)=-(F*sin(teta)*omega*cos(omega*t)); H(3,3)=(F*cos(teta)*omega*cos(omega*t)); H(4,4)=(F*sin(teta)*omega*cos(omega*t)); U=expm(-i*(H_0+H)*tstep)*U; end [C,Ueig]=eig(U); QuasiEnergy=(real(i*log(diag(Ueig))/T)); ExpectH_0=real(C'*(H_0)*C); [ExpSort,Index]=sort(diag(ExpectH_0)); E(:,m,n)=(QuasiEnergy(Index(1:4))); m=m+1; end n=n+1 end [XM,TM]=meshgrid(0:0.1:10,0:pi/16:pi); surf(XM,TM,E(:,XM,TM)) |
回复此楼» 猜你喜欢
关于代码变化问题，想知道的进来
已经有9人回复
-
关于filecode
已经有6人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
-
大家散了吧，后缀研究没有意义，别浪费时间了，过好目前的每一天，不要焦虑
已经有7人回复
-
关于豆爷回答的JTJC与%2F数量
已经有6人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有31人回复
找到一些相关的精华帖子，希望有用哦~
求 摘要 引言 结论翻译 （我已经自己翻译了一遍，希望高人帮忙修改一下）
已经有1人回复
- 编程废又来求助关于哈弗曼树的C++代码……
已经有8人回复
- 求高手帮我分析下 PCR 实验电泳图，怎么改善！！！
已经有6人回复
- 求高手帮我一篇文章的标题修改下
已经有1人回复
- 在线坐等论文摘要，求高手帮助改正！有汉语对应的英语翻译！老师说有些问题！~
已经有2人回复
- 求高手帮忙改一下英文摘要
已经有12人回复
- 【求助】求将c语言改写为C++，求教高人啊！
已经有6人回复
- 求英语高人指点，给老师的英文信，帮忙修改
已经有4人回复
- 修改中译英翻译。急用！请高手帮帮忙
已经有1人回复
- 【求助】请高手帮我看一下：供应商对PVC配方做了什么修改。。
已经有8人回复
- 【求助】求高人帮我修改程序
已经有5人回复
- 【求助】求高手帮忙改写数据实时动态显示程序~
已经有8人回复
- 【求助】matlab程序求助，请高人帮助修改，能运行得不到结果~~~
已经有3人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
20