代码如下
e={1.131,1.33,1.594,1.33,1.131,1.33,1.594,1.33,1.131,1.33,1.594,1.33,1.131,1.33,1.594,1.33,0.605,0.9468,2.668,1.222};(*各膜层厚度倍率*)
n={2.3,1.63,1.38,1.63,2.3,1.63,1.38,1.63,2.3,1.63,1.38,1.63,2.3,1.63,1.38,1.63,2.3,1.63,1.38,1.63};(*各膜层折射率*)
\[Theta]=Table[ArcCos[(1-1.52*1.52*(Sin[Pi/4]/n[[j]])^2)^0.5]/Degree,{j,Length[n]}];(*各膜层入射角*)
d=Table[e[[j]]*550*n[[j]]/4,{j,Length[n]}];(*各膜层厚度*)
\[Eta]p=Table[n[[j]]/Cos[\[Theta][[j]] Degree],{j,Length[n]}];(*p偏振光的\[Eta]*)
\[Eta]s=Table[n[[j]]*Cos[\[Theta][[j]] Degree],{j,Length[n]}];(*s偏振光的\[Eta]*)
\[CapitalDelta]=Table[2*Pi*n[[j]]*d[[j]]*Cos[\[Theta][[j]] Degree]/\[Lambda]/Degree,{j,Length[n]}];(*各膜层相位厚度*)
A=Table[{{Cos[\[CapitalDelta][[j]]],I*Sin[\[CapitalDelta][[j]]Degree]/\[Eta]p[[j]]},{I*\[Eta]p[[j]]*Sin[\[CapitalDelta][[j]]Degree],Cos[\[CapitalDelta][[j]]Degree]}},{j,Length[n]}];(*\[Eta]p传输矩阵*)
as=Table[{{Cos[\[CapitalDelta][[j]]],I*Sin[\[CapitalDelta][[j]]Degree]/\[Eta]s[[j]]},{I*\[Eta]s[[j]]*Sin[\[CapitalDelta][[j]]Degree],Cos[\[CapitalDelta][[j]]Degree]}},{j,Length[n]}];(*\[Eta]s传输矩阵*)
K=Product[A[[j]],{j,1,20}]//FullSimplify;(*\[Eta]p传输矩阵allX*)
F={1,1.52/Cos[\[Theta][[20]] Degree]}//FullSimplify;
K.F;
{Bp,Cp}=K.F;
\[Eta]p0=1.52/Cos[45 Degree];
Tp[\[Lambda]_]:=(4 \[Eta]p0*1.52/Cos[\[Theta][[20]] Degree])/((\[Eta]p0*Bp+Cp)*Conjugate[(\[Eta]p0*Bp+Cp)]);(*p偏振的透过率*)
Plot[Tp[\[Lambda]]//N,{\[Lambda],500,600},PlotRange->{{500,600},{0,1}}](*p偏振的透过率作图*)
K=Product[as[[j]],{j,1,20}]//FullSimplify;(*\[Eta]s传输矩阵allX*)
Fs={1,1.52Cos[\[Theta][[20]] Degree]}//FullSimplify;(*偏振s末矩阵*)
K.Fs;
{Bs,Cs}=K.Fs;
\[Eta]s0=1.52Cos[45 Degree];(*偏振s初态*)
Ts[\[Lambda]_]:=(4 \[Eta]s0*1.52*Cos[\[Theta][[20]] Degree])/((\[Eta]s0*Bs+Cs)*Conjugate[(\[Eta]s0*Bs+Cs)]);(*s偏振的透过率*)
Plot[Ts[\[Lambda]]//N,{\[Lambda],500,600},PlotRange->{{500,600},{0,1}}](*s偏振的透过率作图*)
Ts[550]//N

这是在Ts[550] // N的输出

各位大佬帮忙看看，小弟感激不尽

你lambda没有值