|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
clear8 t' D8 E7 K* R" O- r
clc
/ M2 e! X) h, X/ Umag=0:0.1:30;
( K' ], D, M" B4 Dfor i=1:length(mag)
3 c9 X$ o& o7 ?+ W/ ]( B, G4 d& P. _tau_v0=4.0d0;
/ y9 u' U3 R3 j S+ n, m% Z+ Jgamma=40e-6;
$ G, W- P& l# d z! J* waRF=1e5; $ P3 w" ^* k/ f! x9 x# {9 u$ m7 g
T=40;
4 c$ P2 j- O0 A- q# v" f! Y, MTnd=150;
; W! f' F+ ]5 U1 W8 JTx=14;2 O/ W8 s/ s( Z5 c, x" P
Txx=33;3 H% H0 z& d: K8 [' q$ T
g=4.57e-4;% u' ~. h. m$ o" h
beta=100;( `( a+ ~ ~6 Y9 `
Tbd=1;: s. [$ |! ]- l- F2 c- [" ?
kb=8.617e-5;
; f0 Q0 s* ~9 M0 I( bg1=4.57e-4;! Z$ C2 M4 ^" i2 [3 D p3 E2 h3 H
mub=0.058e-3;1 g* j) z5 Y# S- |6 @6 |) u. [
delta_s=1.0*mub*mag;2 J; h/ H$ C9 I# P3 Y0 K
delta_v=2.0*mub*mag;4 ~2 m' {( v C3 b! v
delta_a=2.0*mub*mag;
d8 }2 Q: J6 {E_vs1=2*delta_a;%表示两个亮激子之间的能量差7 I: y! T/ t/ Y! o8 c4 J- M
E_vs=4*delta_s;%表示两个双激子之间的能量差
7 g: r& |0 j; @
: W8 v. t) K3 k7 }' C. ys1=(gamma*gamma)./(gamma^2+E_vs1.^2).*(1.0/tau_v0);%激子
# L) x# S: D( cs=(gamma*gamma)./(gamma^2+E_vs.^2).*(1.0/tau_v0);%双激子
0 j9 B# T8 G; V. }0 f& C- A) ~if mag=0
Z3 a6 b) m2 g V8 _# w
! z- f5 g6 U1 |9 `; J( h7 ~我这里想做一个判断,如果mag=0我就可以这样,如果mag不等于0可以那样,那么这段代码可以怎么改呢 |
6 E2 i [/ ?, N9 M* R4 \' Y |
|