In[]:=
γn=17600000;​​Η0=1000;​​Μ0=140;​​μ0=4*π*10^-7;​​δσ=1;​​ω=3.1*10^10;​​k0=150*10^2;​​v0=0;​​σ=2*10^4;​​K=0.5;​​d=10*10^-6;​​γ=0;​​Pmin=0.04;​​Pmax=Pmin;​​P01=Pmin;​​P02=0.0;​​weight=260;​​forma=1;​​tp=20;​​tg=200;​​L=1;​​tmax=545.0;
In[]:=
Y0=*
1+*
μ0*σ
k0^2
(ω-k0*v0)
​​Y1=*
1+*
μ0*σ
k0^2
(ω-k0*v0)
Out[]=
-1.1411+1.51727
Out[]=
-1.1411+1.51727
In[]:=
ωh=γn*Η0​​ωm=4*π*γn*Μ0//N​​ω1=
ωh(ωh+ωm)
//N
Out[]=
17600000000
Out[]=
3.09635×
10
10
Out[]=
2.92356×
10
10
In[]:=
θ0=
1
1-*Y0
(-*ωm^2*Y0+ωm*(1+*Y0)(ωh+ω))​​θ1=
1
1-*Y1
(-*ωm^2*Y1+ωm*(1+*Y1)(ωh-ω))
Out[]=
1.29768×
20
10
-3.06372×
20
10

Out[]=
7.84251×
20
10
+2.67181×
20
10

In[]:=
V0=
θ0*d
2ω
​​V1=
θ1*d
2ω
Out[]=
20930.4-49414.8
Out[]=
126492.+43093.7
In[]:=
χ0=K*V0*k0​​χ1=K*V1*k0
Out[]=
1.56978×
8
10
-3.70611×
8
10

Out[]=
9.48691×
8
10
+3.23203×
8
10

In[]:=
κ0=*V1*k0*
δσ
2
​​κ1=*V0*k0*
δσ
2
Out[]=
-3.23203×
8
10
+9.48691×
8
10

Out[]=
3.70611×
8
10
+1.56978×
8
10

In[]:=
η0=
ω^2-ω1^2-θ0*k0*d
2ω
​​η1=
ω^2-ω1^2-θ1*k0*d
2ω
Out[]=
1.40027×
9
10
+7.41222×
8
10

Out[]=
-1.83161×
8
10
-6.46405×
8
10

In[]:=
ϑ0=
δσ
2
*χ1​​ϑ1=
δσ
2
*χ0
Out[]=
4.74346×
8
10
+1.61601×
8
10

Out[]=
7.84889×
7
10
-1.85306×
8
10

In[]:=
PlotV=ReImPlot
1
2ω
1
1-**
1+*
μ0*σ
k0^2
(ω-k0*v0)
-*ωm^2**
1+*
μ0*σ
k0^2
(ω-k0*v0)
+ωm*1+**
1+*
μ0*σ
k0^2
(ω-k0*v0)
(ωh+ω)*d,
1
2ω
1
1-**
1+*
μ0*σ
k0^2
(ω-k0*v0)
-*ωm^2**
1+*
μ0*σ
k0^2
(ω-k0*v0)
+ωm*1+**
1+*
μ0*σ
k0^2
(ω-k0*v0)
(ωh-ω)*d,{σ,0,20000},PlotLegends{{"
V
0
","
V
-1
"},"ReIm"},ReImLabels"действительная часть","мнимая часть",PlotStyle{Red,Blue},AxesLabel"σ,
1
Ом*м
","V, м/c",PlotRangeAll
Out[]=
V
0
V
-1
действительная часть
мнимая часть
​
In[]:=
ReImPlotK*
1
2ω
1
1-**
1+*
μ0*σ
k0^2
(ω-k0*v0)
-*ωm^2**
1+*
μ0*σ
k0^2
(ω-k0*v0)
+ωm*1+**
1+*
μ0*σ
k0^2
(ω-k0*v0)
(ωh+ω)*d*k0,K*
1
2ω
1
1-**
1+*
μ0*σ
k0^2
(ω-k0*v0)
-*ωm^2**
1+*
μ0*σ
k0^2
(ω-k0*v0)
+ωm*1+**
1+*
μ0*σ
k0^2
(ω-k0*v0)
(ωh-ω)*d*k0,{σ,0,20000},PlotLegends{{"
χ
0
","
χ
-1
"},"ReIm"},ReImLabels"действительная часть","мнимая часть",PlotStyle{Red,Blue},AxesLabel"σ,
1
Ом*м
","χ,1/с",PlotRangeAll
Out[]=
χ
0
χ
-1
действительная часть
мнимая часть
​