求教下列程序及图什么意思?初学者看不太懂
mathematica吧
全部回复
仅看楼主
level 3
forest353 楼主
(*Parameters, J. Bonse, JAP 2005*)
ClearAll["Global`*"];
s = 10;(*Shape factor, which could be any value above 0 when \
\[Theta]=0, referred to J. E. Sipe - PRB 1983*)
f = 0.1;(*Filling factor, which could be any value above 0 when \
\[Theta]=0, referred to J. E. Sipe - PRB 1983*)
(*\[Epsilon]=(3.69+I*0.0065)^2;(*Dielectric permittivity, unexcited \
state of Si*)*)
\[Epsilon] = -110 +
I*8;(*Dielectric permittivity, obtained from Drude model*)
\[Theta] = 0;(*Incident angle*)
(*Equations referred to J. Bonse - JAP 2005, which is concluded by \
Bonse from J. E. Sipe - PRB 1983*)
rc = (\[Epsilon] - 1)/(\[Epsilon] + 1);
fc = Sqrt[(s)^2 + 1] - s;
gc = (Sqrt[(s)^2 + 4] + s)/2 - Sqrt[(s)^2 + 1];
\[Eta] = 2*Pi*Abs[\[Nu]1 + Conjugate[\[Nu]2]];
(*s-polarization*)
\[Nu]1 = (hss1*(k1y)^2 + hkk1*(k1x)^2)*\[Gamma]t*(Abs[ts])^2;
\[Nu]2 = (hss2*(k2y)^2 + hkk2*(k2x)^2)*\[Gamma]t*(Abs[ts])^2;
(k1y) = (Sin[\[Theta]] + \[Kappa]y)/\[Kappa]1;
(k2y) = (Sin[\[Theta]] - \[Kappa]y)/\[Kappa]2;
(k1x) = \[Kappa]x/\[Kappa]1;
(k2x) = \[Kappa]x/\[Kappa]2;
\[Kappa]1 = Sqrt[(\[Kappa]x)^2 + (Sin[\[Theta]] + \[Kappa]y)^2];
\[Kappa]2 = Sqrt[(\[Kappa]x)^2 + (Sin[\[Theta]] - \[Kappa]y)^2];
hss1 = 2*I/(Sqrt[1 - (\[Kappa]1)^2] +
Sqrt[\[Epsilon] - (\[Kappa]1)^2]);
hss2 = 2*I/(Sqrt[1 - (\[Kappa]2)^2] +
Sqrt[\[Epsilon] - (\[Kappa]2)^2]);
hkk1 = 2*I*
Sqrt[(\[Epsilon] - (\[Kappa]1)^2)*(1 - (\[Kappa]1)^2)]/(\[Epsilon]*
Sqrt[1 - (\[Kappa]1)^2] + Sqrt[\[Epsilon] - (\[Kappa]1)^2]);
hkk2 = 2*I*
Sqrt[(\[Epsilon] - (\[Kappa]2)^2)*(1 - (\[Kappa]2)^2)]/(\[Epsilon]*
Sqrt[1 - (\[Kappa]2)^2] + Sqrt[\[Epsilon] - (\[Kappa]2)^2]);
ts = 2*Abs[
Cos[\[Theta]]]/(Abs[Cos[\[Theta]]] +
Sqrt[\[Epsilon] - (Sin[\[Theta]])^2]);
\[Gamma]t = (\[Epsilon] - 1)/(4*
Pi*(1 + (1 - f)*(\[Epsilon] - 1)*(fc - rc*gc)/2));
(*Equations finished*)
(*Plotting*)
DensityPlot[\[Eta], {\[Kappa]x, -4, 4}, {\[Kappa]y, -4, 4},
PlotPoints -> 100]
Export["03.xls",
Table[\[Eta], {\[Kappa]x, -4, 4, 0.1}, {\[Kappa]y, -4, 4, 0.1}]]
(*
Plot3D[\[Eta],{\[Kappa]x,-4,4},{\[Kappa]y,-4,4},PlotPoints\[Rule]100]
*)
2017年09月11日 14点09分 1
level 3
forest353 楼主
这个是截图,求教各位大佬这是什么意思?
2017年09月11日 14点09分 2
level 3
forest353 楼主
刚接触不久,老师给了一串代码,望各位大佬指教
2017年09月11日 14点09分 3
level 10
代码本身很普通;难点应该要结合你的专业理论才能看懂,不懂你的专业我也看不懂代表什么意思。
2017年09月11日 22点09分 4
那我想算他的傅里叶变换要加什么代码啊?谢谢
2017年09月14日 23点09分
@forest353 参考FourierTransform
2017年09月14日 23点09分
@奥斯马登 我能再问一下,我作出他的傅里叶变换,怎么画出他的图呢?
2017年09月18日 12点09分
@forest353 Plot[Evaluate@ FourierTransform[BesselJ[2, t], t, \[Omega]], {\[Omega], -Pi, Pi}]
2017年09月18日 14点09分
1