SetOptions[Plot,PlotStyle->{Black,AbsoluteThickness[2]},BaseStyle->{16}];

This is the admittance of a series LRC circuit.  Parameters chosen to make it fairly symmetric.

R=500;
L=.140;
C0=10^(-7);
ω0=1/Sqrt[L*C0]
f0=ω0/2/π
β=R/2/L
Y[ω_]:=(ω/L)/Sqrt[(4*β^2*ω^2+(ω0^2-ω^2)^2)]
ϕ[ω_]:=π/2+ArcTan[(ω0^2-ω^2),-R*ω/L]
FreeI[t_]:=Exp[-t*R/2/L]*Sin[ω0*t]
Plot[Y[ω], {ω,0, 2 ω0},PlotRange->All,AxesLabel->{"angular freq, ω","|Y|"},Ticks->{{{ω0,(Subscript[ω, res])},{2ω0, },{ω0/6, },{ω0/3, },{ω0/2, },{2*ω0/3, },{5*ω0/6, },{7*ω0/6, },{4*ω0/3, },{3*ω0/2, },{5*ω0/3, },{11*ω0/6, }},None}]
Plot[ϕ[ω], {ω,0, 2 ω0},PlotRange->{-π/2,π/2},AxesLabel->{"angular freq, ω","(Subscript[ϕ, I])"},Ticks->{{{ω0,(Subscript[ω, res])},{2ω0, },{ω0/6, },{ω0/3, },{ω0/2, },{2*ω0/3, },{5*ω0/6, },{7*ω0/6, },{4*ω0/3, },{3*ω0/2, },{5*ω0/3, },{11*ω0/6, }},{-π/2,0,π/2}}]
Plot[FreeI[t], {t,0, 10*2*π/ ω0},PlotRange->All,AxesLabel->{"time, t","Current"},Ticks->{{{2*π/ ω0,""}},None}]

Now apply a signal that is cos (2wnt) at low, mid and high freqs rel to resonance.  Then apply all together 2,1,10

ω1=2*π*700;
ω2=2*π*1400;
ω3=2*π*2100;
V1=2;
V2=1;
V3=10;
max=.02;
maxV=V1+V2+V3;
Voltage[t_]:=V1*Cos[ω1*t]+V2*Cos[ω2*t]+V3*Cos[ω3*t]
Plot[Voltage[t], {t,0, 4*π/ ω1},PlotRange->{-maxV,maxV},AxesLabel->{"time, t","Voltage, applied"},Ticks->{{{2*π/ ω1,"2π/ω1"}},None}]
Current1[t_]:=V1*Cos[ω1*t+ϕ[ω1]]*Y[ω1]
Plot[Current1[t], {t,0, 4*π/ ω1},PlotRange->{-max,max},AxesLabel->{"time, t","Current, omega"},Ticks->{{{2*π/ ω1,"2π/ω1"}},None}]
Current2[t_]:=+V2*Cos[ω2*t+ϕ[ω2]]*Y[ω2]
Plot[Current2[t], {t,0, 4*π/ ω1},PlotRange->{-max,max},AxesLabel->{"time, t","Current, 2omega"},Ticks->{{{2*π/ ω1,"2π/ω1"}},None}]
Current3[t_]:=V3*Cos[ω3*t+ϕ[ω3]]*Y[ω3]
Plot[Current3[t],{t,0, 4*π/ ω1},PlotRange->{-max,max},AxesLabel->{"time, t","Current, 3omega"},Ticks->{{{2*π/ ω1,"2π/ω1"}},None}]
Current[t_]:=V1*Cos[ω1*t+ϕ[ω1]]*Y[ω1]+V2*Cos[ω2*t+ϕ[ω2]]*Y[ω2]+V3*Cos[ω3*t+ϕ[ω3]]*Y[ω3]
Plot[Current[t], {t,0, 4*π/ ω1},PlotRange->{-max,max},AxesLabel->{"time, t","Current, total"},Ticks->{{{2*π/ ω1,"2π/ω1"}},None}]
Plot[{Voltage[t],R*Current[t]}, {t,0, 4*π/ ω1},PlotRange->{-maxV,maxV},AxesLabel->{"time, t","V&I, total"},Ticks->{{{2*π/ ω1,"2π/ω1"}},None}]