|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
clear
5 g$ V# S4 U" z3 O$ zclc- K" H3 K- r" B8 h# M" N
mag=0:0.1:30;/ S, V5 `0 h/ S7 ]
for i=1:length(mag)! A# Q1 P& g1 k! P0 n% D b3 X) V
tau_v0=4.0d0;
, w) |; q- x$ q, Q9 X1 ~+ _gamma=40e-6;
, X: m* V7 q( E- R5 saRF=1e5; 9 b( `9 I0 \6 W* w4 J
T=40;& p( n* {6 S, x4 A
Tnd=150;4 [% ^: H7 c4 n1 O: U" ?
Tx=14;7 a5 o2 K$ W0 } B: {6 z) \& ]3 q
Txx=33;
3 E$ P- W* t/ K( K: zg=4.57e-4;
W; B0 c k$ x" _- r. _; o3 T' `beta=100;
^9 U+ r, G& W; L$ HTbd=1;0 C) N+ Q1 m5 j& q7 w
kb=8.617e-5;% ~/ U; m6 D9 q' m
g1=4.57e-4;
% t3 Y( g G( L5 e! j8 _mub=0.058e-3;( S8 t1 ~- w6 B* ~1 f
delta_s=1.0*mub*mag;
* T$ t, V- A+ U3 }9 H4 ydelta_v=2.0*mub*mag;
: |2 ?5 E8 C5 `: R5 Tdelta_a=2.0*mub*mag;$ p9 f, Y4 o: @# R2 x, R% s# \
E_vs1=2*delta_a;%表示两个亮激子之间的能量差% }7 o, y- | C) P2 ]
E_vs=4*delta_s;%表示两个双激子之间的能量差" E2 c7 w' l- q3 ?% }
z! }6 x8 R% @# v
s1=(gamma*gamma)./(gamma^2+E_vs1.^2).*(1.0/tau_v0);%激子1 V- G- E' |5 j O
s=(gamma*gamma)./(gamma^2+E_vs.^2).*(1.0/tau_v0);%双激子1 g) N( `8 {& X& V
if mag=0
( G3 u3 O- T0 k3 ]3 E: @* j! x* F+ G! V+ u8 W
我这里想做一个判断，如果mag=0我就可以这样，如果mag不等于0可以那样，那么这段代码可以怎么改呢 |
* \$ J1 }. I# M% k2 k
|
|
/1 
发表于 2022-7-15 10:14
收藏
淘帖
支持！
反对！
