代码:
\[CapitalOmega] = 2 Pi/365;
\[Omega] = 2 Pi;
\[CapitalPsi] = 0.40666 Pi;
\[Phi][t_] :=
ArcCos[Sqrt[
Sin[\[CapitalOmega] t]^2/((Tan[\[CapitalPsi]] \
Cos[\[CapitalOmega] t])^2 + 1)] - \[Omega] t];
\[Theta][t_] := \[CapitalPsi] - Pi/2 +
ArcTan[Sqrt[Cos[\[CapitalOmega] t]^2 + Cot[\[CapitalPsi]]^2]/
Sin[\[CapitalOmega] t]];
Plot[\[Phi][t_], {t, 1, 365}]
Plot[\[Theta][t_], {t, 1, 365}]
然后错误提示是这样:
Plot::exclul: {Im[Sin[2/365 \[Pi] Pattern[<<2>>]]^2/(1+10.9687 Power[<<2>>])]-0,Im[-2 \[Pi] t_+Sqrt[Sin[<<1>>]^2/Plus[<<2>>]]]-0} must be a list of equalities or real-valued functions. >>
截个图:
谢谢!
\[CapitalOmega] = 2 Pi/365;
\[Omega] = 2 Pi;
\[CapitalPsi] = 0.40666 Pi;
\[Phi][t_] :=
ArcCos[Sqrt[
Sin[\[CapitalOmega] t]^2/((Tan[\[CapitalPsi]] \
Cos[\[CapitalOmega] t])^2 + 1)] - \[Omega] t];
\[Theta][t_] := \[CapitalPsi] - Pi/2 +
ArcTan[Sqrt[Cos[\[CapitalOmega] t]^2 + Cot[\[CapitalPsi]]^2]/
Sin[\[CapitalOmega] t]];
Plot[\[Phi][t_], {t, 1, 365}]
Plot[\[Theta][t_], {t, 1, 365}]
然后错误提示是这样:
Plot::exclul: {Im[Sin[2/365 \[Pi] Pattern[<<2>>]]^2/(1+10.9687 Power[<<2>>])]-0,Im[-2 \[Pi] t_+Sqrt[Sin[<<1>>]^2/Plus[<<2>>]]]-0} must be a list of equalities or real-valued functions. >>
截个图:
谢谢!