level 2
伤CAT心
楼主
求助,各位大佬,程序中间有没有错误,现在的情况是没有办法出图
代码如下
中间要对特征值进行处理,最后画图,原本的问题是特征值排序的问题,我用了一位大佬解答的程序片段,现在的情况是没有办法出图,不清楚是哪里出了问题,结果是后面的图片
p = K0^2 \[Omega] - Ke^2 \[Omega] + Kf^2 \[Omega] -
Kv^2 \[Omega] + \[Epsilon] \[Omega]^2 - \[Omega]^3 -
2 Ke Kf Subscript[l, B] + 2 K0 Kv Subscript[l, B] + \[Epsilon]
\!\(\*SubsuperscriptBox[\(l\), \(B\), \(2\)]\) - \[Omega]
\!\(\*SubsuperscriptBox[\(l\), \(B\), \(2\)]\) + \[Omega]^2 Subscript[
w, A] + \!\(
\*SubsuperscriptBox[\(l\), \(B\), \(2\)]\
\*SubscriptBox[\(w\), \(A\)]\) - K0^2 Subscript[w, B] +
Ke^2 Subscript[w, B] - Kf^2 Subscript[w, B] +
Kv^2 Subscript[w, B] - 2 \[Epsilon] \[Omega] Subscript[w, B] +
2 \[Omega]^2 Subscript[w, B] -
2 \[Omega] Subscript[w, A] Subscript[w, B] + \[Epsilon]
\!\(\*SubsuperscriptBox[\(w\), \(B\), \(2\)]\) - \[Omega]
\!\(\*SubsuperscriptBox[\(w\), \(B\), \(2\)]\) + Subscript[w, A]
\!\(\*SubsuperscriptBox[\(w\), \(B\), \(2\)]\);
\[Omega]tsol = Solve[p == 0, \[Omega]];
\[Omega]1 = Evaluate[ \[Omega] /. \[Omega]tsol[[1]]];
Subscript[g, s] = (
2 Ke Kf -
2 K0 Kv - \[Epsilon] Subscript[l, B] + \[Omega]1 Subscript[l, B] -
Subscript[l, B] Subscript[w, A])/(\[Omega]1 - Subscript[w, B]);
K0 = 1; Kv = 1; Subscript[l, B] = 1.14; Subscript[w, B] = 1.97; Ke = \
0; Kf = 0; Subscript[w, 0] = 2 Subscript[w, B]/3;
Subscript[w, A] = 0;
matrix[\[Epsilon]_] = ( {
{Subscript[w, A] + \[Epsilon] + I*(Subscript[g, s]) , (K0 + Ke) +
I*(Kv + Kf)},
{ (K0 - Ke) + I*(Kv - Kf), Subscript[w, B] - I* Subscript[l, B]}
});
bb = Solve[Det[matrix[#] - \[Omega]*IdentityMatrix[2]] == 0, \[Omega]];
cc = \[Omega] /. bb;
AWlist = Range[-1, 1, 0.01];
dataRe = Map[Transpose@{AWlist, #} &,
Re@Transpose@Map[SortBy[cc, Re] &, AWlist]];
ListLinePlot[dataRe, PlotRange -> Full]
AWlist2 = Range[-1, 1, 0.01]; dataIm =
Map[Transpose@{AWlist2, #} &,
Im@Transpose@Map[SortBy[cc, Im] &, AWlist2]];
ListLinePlot[dataIm, PlotRange -> Full]

2024年11月25日 12点11分
1
代码如下
中间要对特征值进行处理,最后画图,原本的问题是特征值排序的问题,我用了一位大佬解答的程序片段,现在的情况是没有办法出图,不清楚是哪里出了问题,结果是后面的图片
p = K0^2 \[Omega] - Ke^2 \[Omega] + Kf^2 \[Omega] -
Kv^2 \[Omega] + \[Epsilon] \[Omega]^2 - \[Omega]^3 -
2 Ke Kf Subscript[l, B] + 2 K0 Kv Subscript[l, B] + \[Epsilon]
\!\(\*SubsuperscriptBox[\(l\), \(B\), \(2\)]\) - \[Omega]
\!\(\*SubsuperscriptBox[\(l\), \(B\), \(2\)]\) + \[Omega]^2 Subscript[
w, A] + \!\(
\*SubsuperscriptBox[\(l\), \(B\), \(2\)]\
\*SubscriptBox[\(w\), \(A\)]\) - K0^2 Subscript[w, B] +
Ke^2 Subscript[w, B] - Kf^2 Subscript[w, B] +
Kv^2 Subscript[w, B] - 2 \[Epsilon] \[Omega] Subscript[w, B] +
2 \[Omega]^2 Subscript[w, B] -
2 \[Omega] Subscript[w, A] Subscript[w, B] + \[Epsilon]
\!\(\*SubsuperscriptBox[\(w\), \(B\), \(2\)]\) - \[Omega]
\!\(\*SubsuperscriptBox[\(w\), \(B\), \(2\)]\) + Subscript[w, A]
\!\(\*SubsuperscriptBox[\(w\), \(B\), \(2\)]\);
\[Omega]tsol = Solve[p == 0, \[Omega]];
\[Omega]1 = Evaluate[ \[Omega] /. \[Omega]tsol[[1]]];
Subscript[g, s] = (
2 Ke Kf -
2 K0 Kv - \[Epsilon] Subscript[l, B] + \[Omega]1 Subscript[l, B] -
Subscript[l, B] Subscript[w, A])/(\[Omega]1 - Subscript[w, B]);
K0 = 1; Kv = 1; Subscript[l, B] = 1.14; Subscript[w, B] = 1.97; Ke = \
0; Kf = 0; Subscript[w, 0] = 2 Subscript[w, B]/3;
Subscript[w, A] = 0;
matrix[\[Epsilon]_] = ( {
{Subscript[w, A] + \[Epsilon] + I*(Subscript[g, s]) , (K0 + Ke) +
I*(Kv + Kf)},
{ (K0 - Ke) + I*(Kv - Kf), Subscript[w, B] - I* Subscript[l, B]}
});
bb = Solve[Det[matrix[#] - \[Omega]*IdentityMatrix[2]] == 0, \[Omega]];
cc = \[Omega] /. bb;
AWlist = Range[-1, 1, 0.01];
dataRe = Map[Transpose@{AWlist, #} &,
Re@Transpose@Map[SortBy[cc, Re] &, AWlist]];
ListLinePlot[dataRe, PlotRange -> Full]
AWlist2 = Range[-1, 1, 0.01]; dataIm =
Map[Transpose@{AWlist2, #} &,
Im@Transpose@Map[SortBy[cc, Im] &, AWlist2]];
ListLinePlot[dataIm, PlotRange -> Full]
