|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
clear5 P0 H0 [4 I/ i/ O0 d# @; Q
clc
0 k, g, n" L( \7 i4 pmag=0:0.1:30;& r) A) D, M, ]
for i=1:length(mag)0 T6 l( b' [/ V7 i' h
tau_v0=4.0d0;3 I0 [7 ?" J% V- @
gamma=40e-6;
* `9 N5 z8 J) x/ B" D5 a8 e2 i' iaRF=1e5;
7 Z. s; z: l1 P$ cT=40;
; _: y8 Y d/ ZTnd=150;4 u* J: D0 h7 `
Tx=14;# x/ \4 X+ p5 q$ l' N
Txx=33;
6 C/ r0 n) D5 Y) x: U( gg=4.57e-4;6 z; E( N" _3 n( z. \" ~
beta=100;
: g( X2 @( c/ Z- O; UTbd=1;
5 ]- {, c7 L, D. [0 Ekb=8.617e-5;
0 c- q( L1 p" X3 | V" Qg1=4.57e-4;
5 D: \; W. k. {( A( dmub=0.058e-3;
' G# l' {! o4 H$ [/ A$ }delta_s=1.0*mub*mag;
! x6 t% k3 ?7 L( c* F% Z" Z, Qdelta_v=2.0*mub*mag;6 M, _6 _' i: t: o! ~2 z
delta_a=2.0*mub*mag;, g h; z8 v8 {3 T
E_vs1=2*delta_a;%表示两个亮激子之间的能量差
' H3 ]7 d' M( X* S* B4 aE_vs=4*delta_s;%表示两个双激子之间的能量差
/ ?2 J0 V! L3 a& f G) r# M( t2 x! ~ g
s1=(gamma*gamma)./(gamma^2+E_vs1.^2).*(1.0/tau_v0);%激子5 g" w( ^- q1 x
s=(gamma*gamma)./(gamma^2+E_vs.^2).*(1.0/tau_v0);%双激子
5 }# w6 b# S1 _8 u$ W; G- [' J0 fif mag=0
, U% x5 `8 ?! m2 }. G1 q
3 ]& T" ^# ]' |- q+ s我这里想做一个判断,如果mag=0我就可以这样,如果mag不等于0可以那样,那么这段代码可以怎么改呢 |
" X2 o* y( ^3 B |
|