[pstricks] Re : Function of the current in a coil

Thu Feb 24 20:25:13 CET 2011

Dear Thomas,

Concerning your problem about the circuits RL, this study which I had made some years ago must correspond, perhaps, to what you search?

(in french)

http://manuel.luque.perso.neuf.fr/RL/

 best regards,

M.L

-----E-mail d'origine-----
De : Thomas Söll 1und1 <Thomas.Soell at onlinehome.de>
A : pstricks at tug.org
Envoyé le : Mercredi, 23 Février 2011 22:42
Sujet : [pstricks] Function of the current in a coil

Hallo,

I want to draw the current through a coil. The example shows my result. It seems to be ok.
The problem comes, if I want to draw the current for parameters where the maximal current is not reached, as shown in the last example.
Then the parts of the function are not continual.

The endpoints should be calculated and this gives the startpoint from the next curve and so on.

Does anybody knows a solution?

Thanks for help

Thomas

\documentclass[a4paper,12pt]{article}

\usepackage[T1]{fontenc}

\usepackage[latin1]{inputenc}

\usepackage{ngerman}

\usepackage{fancyhdr}

\usepackage{amsmath,amsthm}

\usepackage{enumerate}

\usepackage{alphalph}

\usepackage{calc}         % From LaTeX distribution

\usepackage{multido}      % >From PSTricks

\usepackage{ifthen}       % From LaTeX distribution

\usepackage{array}

\usepackage[dvips]{graphics}

%\usepackage[dvips,usenames,dvipsnames]{xcolor} %% Farben sind im Dokument xcolor.pdf definiert

\usepackage[dvips]{changebar}

\usepackage{pstricks}

\usepackage{pst-plot}

\usepackage{pst-node}

\usepackage{pst-tree}

\usepackage{pst-eps}

\usepackage{pst-ghsb}

\usepackage[tiling]{pst-fill}

\usepackage{pst-text}

\usepackage{pst-blur}

%\usepackage{pst-all}

\usepackage{pst-osci,pst-circ}

\usepackage{pst-poly}     % From pstricks/contrib/pst-poly

\usepackage{pst-optic}

\usepackage{pst-eucl}

\usepackage{pst-math}

\usepackage{pst-func}

%\usepackage{pst-3dplot}

\usepackage{pst-electricfield}

\usepackage{pstricks-add}

%\input{dokdef0}

\begin{document}

\section*{Selbstinduktion}

%
\begin{pspicture}(-6,-1)(4,5)

\newpsstyle{Red}{plotstyle=line,linecolor=red,linewidth=0.15,plotpoints=2500}

\newpsstyle{Blau}{plotstyle=line,linecolor=blue,linewidth=0.15,plotpoints=2500}

\newpsstyle{Green}{plotstyle=line,linecolor=green,linewidth=0.15,plotpoints=2500}

%\psgrid[subgriddiv=5,gridcolor=lightgray,griddots=10,gridlabels=10pt](-4,-1)(4,5)

\pnode(0,0){A}

\pnode(0,4){B}

\pnode(4,4){C}

\pnode(4,0){D}

\pnode(-3,0){E}

\pnode(-3,4){F}

% Dipole node connection

\transformer[primarylabel=$n_1$,secondarylabel=$n_1$,arrows=-o,arrowscale=1.7](B)(A)(C)(D){}

% Wire to complete circuit

\wire[arrows=o-,arrowscale=1.7](F)(B)%

\resistor[labeloffset=0,arrows=o-,arrowscale=1.7](E)(1,4|E){$R$}%

\rput(4.5,1.9){%

 \psscalebox{0.2}{%

  \Oscillo[amplitude1=-3.5,phase1=90,period1= 2E30,damping1=0.02,plotstyle1=Green]}}%

\rput(-0.6,1.4){%

 \psscalebox{0.15}{%

  \Oscillo[amplitude1=0,offset2=3,amplitude2=-3,period2= 2E31,damping2=0.02,phase2=90,plotstyle2=Blau]}}

\rput(-2.8,2.0){%

 \psscalebox{0.15}{%

  \Oscillo[ Wave1= \RectangleA,period1=25,amplitude1=1,offset1=1,plotstyle1=Red]}}

%\wire(A)(B)

\qdisk(A){2pt}

\qdisk(-2,0){2pt}

\pnode(-2,0){M}

\pnode(4.18,1.0){G} \ncangles[angleA=0, armA=1.7cm,angleB=-90,linearc=0.2]{-o}{C}{G}

\pnode(4,1.0){H} \nccurve[angleA=90,angleB=-90]{-o}{D}{H}

\pnode(-0.85,0.73){I} \nccurve[angleA=90,angleB=-90,arrowscale=.85]{-o}{A}{I}

\pnode(-1.05,0.73){J} \nccurve[angleA=90,angleB=-90,arrowscale=.85]{-o}{M}{J}

\pnode(-3.3,1.33){K} \ncangles[angleA=180, armA=.8cm,angleB=-90,linearc=0.2,arrowscale=.85]{-o}{F}{K}

\pnode(-3.1,1.33){L} \nccurve[angleA=90,angleB=-90,arrowscale=.85]{-o}{E}{L}

\end{pspicture}

\begin{center}

\begin{pspicture}(-1,-1)(12,5)

 \psframe*[linecolor=green!20](-0.8,-0.8)(12.5,4.6)

 \psframe*[linecolor=white](0,0)(12,4.0)

 \psgrid[subgriddiv=0,gridcolor=lightgray,griddots=10,gridlabels=0,xunit=1,yunit=1.333](0,0)(11.5,3.0)

\psset{xAxisLabel=$t$, yAxisLabel=$U$,arrowscale=2}

\begin{psgraph}{->}(0,-0.5)(12,1.0){12cm}{4cm}

\makeatletter

\psplot[plotpoints=20,linecolor=red,linewidth=2pt]{0}{11.5}

  [ /yMax 0.5 def /T0 2 def /T1 T0 def ]% <=== this is new !!!!

  { x T1 gt {

      T1 yMax \tx at ScreenCoor

      2 copy 4 2 roll L neg L

      /yMax yMax neg def /T1 T1 T0 add def } if yMax }

  \end{psgraph}

\end{pspicture}

\end{center}
\begin{center}

\begin{pspicture}(-1,-0.5)(12,4)

 \psframe*[linecolor=green!20](-0.8,-0.5)(12.5,4.6)%

 \psframe*[linecolor=white](0,0)(12,4)

 \psgrid[subgriddiv=0,gridcolor=lightgray,griddots=10,gridlabels=0](0,0)(12,4)

\psset{xAxisLabel=$t$, yAxisLabel=$I$,arrowscale=2}

\def\R{120}%    Widerstand

\def\L{0.16}%   Induktivität

\def\t{0.004}%    Zeitkonstante

\def\A{1}%     Amplitude

\begin{psgraph}{->}(0,0)(0,-.5)(12,1.5){12cm}{4cm}

\psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=blue,yMaxValue=5,yMinValue=-2]{0}{2}

  {\A*(1-2.7182818^(-x*\t*\R/\L))}

\psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=blue,yMaxValue=5,yMinValue=-2]{4}{6}

  {\A*(1-2.7182818^(-(x-4)*\t*\R/\L))}

\psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=blue,yMaxValue=5,yMinValue=-2]{8}{10}

  {\A*(1-2.7182818^(-(x-8)*\t*\R/\L))}

\psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=blue,yMaxValue=5,yMinValue=-2]{2}{4}

  {\A*2.7182818^(-(x-2)*\t*\R/\L)}

\psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=blue,yMaxValue=5,yMinValue=-2]{6}{8}

  {\A*2.7182818^(-(x-6)*\t*\R/\L)}

  \psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=blue,yMaxValue=5,yMinValue=-2]{10}{11.5}

  {\A*2.7182818^(-(x-10)*\t*\R/\L)}

  \end{psgraph}

\end{pspicture}

\end{center}

\begin{center}

\begin{pspicture}(-1,-1)(12,6)

 \psframe*[linecolor=green!20](-0.8,-0.4)(12.5,6.8)

 \psframe*[linecolor=white](0,0)(12,6)

 \psgrid[subgriddiv=0,gridcolor=lightgray,griddots=10,gridlabels=0](0,0)(12,6)

\psset{xAxisLabel=$t$, yAxisLabel=$U_{\text{ind}}$,arrowscale=2}

\def\R{120}%    Widerstand

\def\L{0.16}%   Induktivität

\def\t{0.004}%  Zeitkonstante

\def\A{1}%     Amplitude

\def\B{1}%     Amplitude

\begin{psgraph}{->}(0,0)(0,-1.5)(12,1.5){12cm}{6cm}

\psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=green,yMaxValue=5,yMinValue=-2]{0}{2}

  {\A*(-2.7182818^(-x*\t*\R/\L))}

\psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=green,yMaxValue=5,yMinValue=-2]{4}{6}

  {\A*(-2.7182818^(-(x-4)*\t*\R/\L))}

\psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=green,yMaxValue=5,yMinValue=-2]{8}{10}

  {\A*(-2.7182818^(-(x-8)*\t*\R/\L))}

\psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=green,yMaxValue=1.5,yMinValue=-2]{2}{4}

  {\B*2.7182818^(-(x-2)*\t*\R/\L)}

\psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=green,yMaxValue=1.5,yMinValue=-2]{6}{8}

  {\B*2.7182818^(-(x-6)*\t*\R/\L)}

  \psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=green,yMaxValue=1.5,yMinValue=-2]{10}{11.5}

  {\B*2.7182818^(-(x-10)*\t*\R/\L)}

%  \psline[linewidth=1pt,linecolor=blue](2,0)(2,1.5)%

%  \psline[linewidth=1pt,linecolor=blue](6,0)(6,1.5)%

%  \psline[linewidth=1pt,linecolor=blue](10,0)(10,1.5)%

  \end{psgraph}

\end{pspicture}

\end{center}

\begin{center}

\begin{pspicture}(-1,-0.5)(12,4)

 \psframe*[linecolor=green!20](-0.8,-0.5)(12.5,4.6)%

 \psframe*[linecolor=white](0,0)(12,4)

 \psgrid[subgriddiv=0,gridcolor=lightgray,griddots=10,gridlabels=0](0,0)(12,4)

\psset{xAxisLabel=$t$, yAxisLabel=$I$,arrowscale=2}

\def\R{120}%    Widerstand

\def\L{0.5}%   Induktivität

\def\t{0.004}%    Zeitkonstante

\def\A{1}%     Amplitude

\begin{psgraph}{->}(0,0)(0,-.5)(12,1.5){12cm}{4cm}

\psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=blue,yMaxValue=5,yMinValue=-2]{0}{2}

  {\A*(1-2.7182818^(-x*\t*\R/\L))}

\psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=blue,yMaxValue=5,yMinValue=-2]{4}{6}

  {\A*(1-2.7182818^(-(x-4)*\t*\R/\L))}

\psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=blue,yMaxValue=5,yMinValue=-2]{8}{10}

  {\A*(1-2.7182818^(-(x-8)*\t*\R/\L))}

\psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=blue,yMaxValue=5,yMinValue=-2]{2}{4}

  {\A*2.7182818^(-(x-2)*\t*\R/\L)}

\psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=blue,yMaxValue=5,yMinValue=-2]{6}{8}

  {\A*2.7182818^(-(x-6)*\t*\R/\L)}

  \psplot[linewidth=1pt,algebraic=true,plotpoints=500,linecolor=blue,yMaxValue=5,yMinValue=-2]{10}{11.5}

  {\A*2.7182818^(-(x-10)*\t*\R/\L)}

  \end{psgraph}

\end{pspicture}

\end{center}

\end{document} 

_______________________________________________

PSTricks mailing list

PSTricks at tug.org

http://tug.org/mailman/listinfo/pstricks

archive: http://www.tug.org/pipermail/pstricks/

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://tug.org/pipermail/pstricks/attachments/20110224/b6707eda/attachment.html>