<?xml version="1.0"  encoding="iso-8859-1" ?>
<?xml-stylesheet type="text/xsl" href="xsl/mathml.xsl"?>
<!DOCTYPE MMLEntities [
<!ENTITY af "&#x02061;">
<!ENTITY alpha "&#x003B1;">
<!ENTITY ap "&#x02048;">
<!ENTITY apos "&#x00027;">
<!ENTITY ApplyFunction "&#x02061;">
<!ENTITY Backslash "&#x02216;">
<!ENTITY beta "&#x003B2;">
<!ENTITY cap "&#x02229;">
<!ENTITY CenterDot "&#x000B7;">
<!ENTITY chi "&#x003C7;">
<!ENTITY compfn "&#x02218;">
<!ENTITY cup "&#x0222A;">
<!ENTITY dd "d">
<!ENTITY Del "&#x02207;">
<!ENTITY Delta "&#x00394;">
<!ENTITY delta "&#x003B4;">
<!ENTITY DifferentialD "d">
<!ENTITY DoubleRightArrow "&#x021D2;">
<!ENTITY ee "e">
<!ENTITY empty "&#x02205;&#x02063;">
<!ENTITY emptyset "&#x02205;&#x02063;">
<!ENTITY epsi "&#x003B5;">
<!ENTITY epsilon "&#x003B5;">
<!ENTITY epsiv "&#x0025B;">
<!ENTITY equiv "&#x02261;">
<!ENTITY eta "&#x003B7;">
<!ENTITY exist "&#x02203;">
<!ENTITY Exists "&#x02203;">
<!ENTITY ExponentialE "e">
<!ENTITY forall "&#x02200;">
<!ENTITY ForAll "&#x02200;">
<!ENTITY Gamma "&#x00393;">
<!ENTITY gamma "&#x003B3;">
<!ENTITY Gammad "&#x003DC;">
<!ENTITY gammad "&#x003DC;">
<!ENTITY ge "&#x02265;">
<!ENTITY gt "&#x0003E;">
<!ENTITY ii "i">
<!ENTITY image "&#x02111;">
<!ENTITY ImaginaryI "i">
<!ENTITY in "&#x02208;">
<!ENTITY infin "&#x0221E;">
<!ENTITY Integral "&#x0222B;">
<!ENTITY InvisibleTimes "&#x02062;">
<!ENTITY iota "&#x003B9;">
<!ENTITY it "&#x02062;">
<!ENTITY kappa "&#x003BA;">
<!ENTITY kappav "&#x003F0;">
<!ENTITY Lambda "&#x0039B;">
<!ENTITY lambda "&#x003BB;">
<!ENTITY lceil "&#x02308;">
<!ENTITY le "&#x02264;">
<!ENTITY LeftBracketingBar "|">
<!ENTITY LeftCeiling "&#x02308;">
<!ENTITY LeftFloor "&#x0230A;">
<!ENTITY lfloor "&#x0230A;">
<!ENTITY LowerRightArrow "&#x02198;">
<!ENTITY mu "&#x003BC;">
<!ENTITY nabla "&#x02207;">
<!ENTITY nbsp "&#x000A0;">
<!ENTITY ne "&#x02260;">
<!ENTITY not "&#x000AC;">
<!ENTITY NotElement "&#x02209;">
<!ENTITY NotEqual "&#x02260;">
<!ENTITY notin "&#x02209;">
<!ENTITY NotSubsetEqual "&#x02288;">
<!ENTITY nsub "&#x02284;">
<!ENTITY nsube "&#x02288;">
<!ENTITY nsubseteq "&#x02288;">
<!ENTITY nu "&#x003BD;">
<!ENTITY Omega "&#x003A9;">
<!ENTITY omega "&#x003C9;">
<!ENTITY OverBar "&#x000AF;">
<!ENTITY part "&#x02202;">
<!ENTITY Phi "&#x003A6;">
<!ENTITY phi "&#x003C6;">
<!ENTITY phiv "&#x003D5;">
<!ENTITY Pi "&#x003A0;">
<!ENTITY pi "&#x003C0;">
<!ENTITY piv "&#x003D6;">
<!ENTITY prime "&#x02032;">
<!ENTITY Product "&#x0220F;">
<!ENTITY Psi "&#x003A8;">
<!ENTITY psi "&#x003C8;">
<!ENTITY rarr "&#x02192;">
<!ENTITY rceil "&#x02309;">
<!ENTITY real "&#x0211C;">
<!ENTITY rfloor "&#x0230B;">
<!ENTITY rho "&#x003C1;">
<!ENTITY rhov "&#x003F1;">
<!ENTITY RightArrow "&#x02192;">
<!ENTITY rightarrow "&#x02192;">
<!ENTITY RightBracketingBar "|">
<!ENTITY RightCeiling "&#x02309;">
<!ENTITY RightFloor "&#x0230B;">
<!ENTITY Sigma "&#x003A3;">
<!ENTITY sigma "&#x003C3;">
<!ENTITY sigmav "&#x003C2;">
<!ENTITY sub "&#x02282;">
<!ENTITY sube "&#x02286;">
<!ENTITY SuchThat "&#x0220B;">
<!ENTITY Sum "&#x02211;">
<!ENTITY tau "&#x003C4;">
<!ENTITY Theta "&#x00398;">
<!ENTITY theta "&#x003B8;">
<!ENTITY thetav "&#x003D1;">
<!ENTITY times "&#x000D7;">
<!ENTITY UpperRightArrow "&#x02197;">
<!ENTITY Upsi "&#x003D2;">
<!ENTITY upsi "&#x003C5;">
<!ENTITY vee "&#x02228;">
<!ENTITY wedge "&#x02227;">
<!ENTITY weierp "&#x02118;">
<!ENTITY xi "&#x003BE;">
<!ENTITY Xi "&#x0039E;">
<!ENTITY zeta "&#x003B6;">
]>
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<title>virtap5.html</title>
<!-- Created by Maple 8.00, IBM INTEL NT -->
<link href="../navi/css/sovell.css" rel="stylesheet" type="text/css"/>
</head>
<body bgcolor="FFFFFF">

<a href="../mpl/virtap5.mws" target="_blank">virtap5.mws</a>

<a name="MapleAutoBookmark1" />
<p align="center">
<b><u><font color="#000000" size="5">Monisilmukkainen vaihtovirtapiiri</font></u></b>
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Otetaan tarkastelun kohteeksi RLC-vaihtovirtapiiri jossa on käämejä, vastuksia ja kondensaattoreita.</font>
</p>
<p align="left">
</p>
<p align="left">
<b><font color="#000000" size="5">Kytkentä 1</font></b>
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Tarkastellaan virtapiiriä, jossa yksinkertaiseen RLC-piiriin on kodensaattorin rinnalle kytketty piiri, jossa on vastus ja toinen kondensaattori. Alla kuva piirin kokoonpanosta. <br />
</font>
</p>
<p align="center">
<img src="images/virtap51.gif" width="244" height="102" alt="[Maple Bitmap]" />
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Piirin komponentit kohdistavat jännitteeseen kolme erilaista pudotusta,</font>
</p>
<p align="left">
</p>
<p align="center">
<img src="images/virtap52.gif" width="122" height="58" alt="E[L] = L*diff(I(t),t)" align="middle" />
<font color="#000000">&nbsp;&nbsp;&nbsp;käämin yli,</font>
</p>
<p align="center">
<font color="#000000">&nbsp;</font>
<img src="images/virtap53.gif" width="89" height="32" alt="E[R] = R*I(t)" align="middle" />
<font color="#000000">&nbsp;&nbsp;&nbsp;&nbsp;vastuksen yli sekä</font>
</p>
<p align="center">
<img src="images/virtap54.gif" width="69" height="58" alt="E[C] = 1/C" align="middle" />
<math xmlns='http://www.w3.org/1998/Math/MathML'>
  <mrow>
    <mo>&Integral;</mo>
    <mrow>
      <mn>&ImaginaryI;</mn>
      <mo>&ApplyFunction;</mo>
      <mfenced>
        <mi>t</mi>
      </mfenced>
    </mrow>
    <mo>&InvisibleTimes;</mo>
    <mrow>
      <mo>&DifferentialD;</mo>
      <mi>t</mi>
    </mrow>
  </mrow>
</math>
<font color="#000000">&nbsp;&nbsp;&nbsp;&nbsp;kondensaattorin yli.</font>
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Lisäksi piiriin syötetään sinimuotoista vaihtovirtaa </font>
<img src="images/virtap55.gif" width="133" height="32" alt="E(t) = E[0]*sin(omega*t)" align="middle" />
<font color="#000000">. Muodostetaan yllä olevalle piirille differentiaaliyhtälöt kiertämällä piiriä virtasilmukoiden </font>
<img src="images/virtap56.gif" width="29" height="32" alt="I[1]" align="middle" />
<font color="#000000">&nbsp;ja </font>
<img src="images/virtap57.gif" width="29" height="32" alt="I[2]" align="middle" />
<font color="#000000">&nbsp;muodostamia silmukoita myöden. Kirchhoffin lain mukaan molempien virtojen kulkiessa komponentin yli tapahtuu jännitteen muutos. Saadaan yhtälöt </font>
</p>
<p align="left">
</p>
<p align="center">
<img src="images/virtap58.gif" width="109" height="58" alt="R[1]*I[1](t)+1/C[1]" align="middle" />
<img src="images/virtap59.gif" width="317" height="58" alt="int(I[1](t)-I[2](t),t)+L[1]*diff(I[1](t),t) = E[0]*sin(omega*t)" align="middle" />
<font color="#000000">,</font>
</p>
<p align="left">
</p>
<p align="center">
<img src="images/virtap510.gif" width="39" height="58" alt="1/C[1]" align="middle" />
<img src="images/virtap511.gif" width="235" height="58" alt="int(I[2](t)-I[1](t),t)+R[2]*I[2](t)+1/C[2]" align="middle" />
<img src="images/virtap512.gif" width="101" height="58" alt="int(I[2](t),t) = 0" align="middle" />
<font color="#000000">.</font>
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Jotta integraalitermeistä päästään eroon, yhtälöt derivoidaan: </font>
</p>
<p align="center">
</p>
<p align="center">
<img src="images/virtap513.gif" width="425" height="76" alt="L[1]*diff(I[1](t),`$`(t,2))+R[1]*diff(I[1](t),t)+(I[1](t)-I[2](t))/C[1] = E[0]*omega*cos(omega*t)" align="middle" />
<font color="#000000">,</font>
</p>
<p align="center">
<img src="images/virtap514.gif" width="282" height="72" alt="(I[2](t)-I[1](t))/C[1]+R[2]*diff(I[2](t),t)+I[2](t)/C[2] = 0" align="middle" />
<font color="#000000">.</font>
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Tämä on differentiaaliyhtälöryhmä pakotetulle värähtelylle. Ratkaistaan yhtälöt ja siten virrat virtasilmukoissa 1 ja 2. <br />
<br />
Laskujen aluksi on syytä hävittää mahdollisista aiemmista laskuista jääneet muuttujat: <br />
</font>
</p>
<p align="left">
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tr>
<td valign="top"><tt>&gt; &nbsp;&nbsp;</tt></td>
<td width="100%" valign="top"><b><font color="#FF0000">restart;</font></b>
</td></tr>
</table>
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Määritellään yhtälöryhmä ja sen tuntemattomat.</font>
</p>
<p align="left">
</p>
<p align="left">
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tr>
<td valign="top"><tt>&gt; &nbsp;&nbsp;</tt></td>
<td width="100%" valign="top"><b><font color="#FF0000">ryhma:= L[1]*diff(i[1](t), t$2)+R[1]*diff(i[1](t), t)+(i[1](t)-i[2](t))/C[1]=E[0]*omega*cos(omega*t), i[2](t)/C[2]+(-i[1](t)+i[2](t))/C[1]+R[2]*diff(i[2](t), t)=0;</font></b>
</td></tr>
</table>
</p>
<p align="center">
<img src="images/virtap515.gif" width="761" height="63" alt="ryhma := L[1]*diff(i[1](t),`$`(t,2))+R[1]*diff(i[1](t),t)+(i[1](t)-i[2](t))/C[1] = E[0]*omega*cos(omega*t), i[2](t)/C[2]+(-i[1](t)+i[2](t))/C[1]+R[2]*diff(i[2](t),t) = 0" />
</p>
<p align="left">
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tr>
<td valign="top"><tt>&gt; &nbsp;&nbsp;</tt></td>
<td width="100%" valign="top"><b><font color="#FF0000">tuntemattomat:= {i[1](t), i[2](t)};</font></b>
</td></tr>
</table>
</p>
<p align="center">
<img src="images/virtap516.gif" width="221" height="27" alt="tuntemattomat := {i[2](t), i[1](t)}" />
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Sijoitetaan vakioille arvot ja asetetaan alkuehto. Kytkettävä lähde on verkkovirtalähde ja sitä ennen piirissä ei ollut virtaa. </font>
</p>
<p align="left">
</p>
<p align="left">
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tr>
<td valign="top"><tt>&gt; &nbsp;&nbsp;</tt></td>
<td width="100%" valign="top"><b><font color="#FF0000">ryhma1:= subs({L[1]=0.1, R[1]=5, C[1]=0.001, R[2]=3, C[2]=0.0005, E[0]=230, omega=2*Pi*50}, {ryhma});</font></b>
</td></tr>
</table>
</p>
<p align="left">
<img src="images/virtap517.gif" width="85" height="20" alt="ryhma1 := {.1*diff(i[1](t),`$`(t,2))+5*diff(i[1](t),t)+.1e4*i[1](t)-.1e4*i[2](t) = 23000*Pi*cos(100*Pi*t), 3000.000000*i[2](t)-.1e4*i[1](t)+3*diff(i[2](t),t) = 0}" />
<br />
<img src="images/virtap518.gif" width="892" height="63" alt="ryhma1 := {.1*diff(i[1](t),`$`(t,2))+5*diff(i[1](t),t)+.1e4*i[1](t)-.1e4*i[2](t) = 23000*Pi*cos(100*Pi*t), 3000.000000*i[2](t)-.1e4*i[1](t)+3*diff(i[2](t),t) = 0}" />
</p>
<p align="left">
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tr>
<td valign="top"><tt>&gt; &nbsp;&nbsp;</tt></td>
<td width="100%" valign="top"><b><font color="#FF0000">alkuehto:= i[1](0)=0, D(i[1])(0)=0, i[2](0)=0;</font></b>
</td></tr>
</table>
</p>
<p align="center">
<img src="images/virtap519.gif" width="310" height="27" alt="alkuehto := i[1](0) = 0, D(i[1])(0) = 0, i[2](0) = 0" />
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Ratkaistaan yhtälöryhmä numeerisin laskentamenetelmin.</font>
</p>
<p align="left">
</p>
<p align="left">
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tr>
<td valign="top"><tt>&gt; &nbsp;&nbsp;</tt></td>
<td width="100%" valign="top"><b><font color="#FF0000">rtk:= dsolve({ryhma1[], alkuehto}, tuntemattomat, type=numeric, output=listprocedure);</font></b>
</td></tr>
</table>
</p>
<p align="center">
<img src="images/virtap520.gif" width="850" height="45" alt="rtk := [t = proc (t) option `Copyright (c) 1993 by the University of Waterloo. All rights reserved.`; t end proc, i[1](t) = proc (t) local res, solnproc, outpoint, ndsol; option `Copyright (c) 2000 by ..." />
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Poimitaan ratkaisusta virtaa kuvaavat termit.</font>
</p>
<p align="left">
</p>
<p align="left">
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tr>
<td valign="top"><tt>&gt; &nbsp;&nbsp;</tt></td>
<td width="100%" valign="top"><b><font color="#FF0000">virrat:= subs(rtk, [i[1](t), i[2](t)]):</font></b>
</td></tr>
</table>
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Piirretään kuvaaja virralle piirin eri osissa. Sininen käyrä on jaetun kondensaattorin </font>
<img src="images/virtap521.gif" width="35" height="32" alt="C[1]" align="middle" />
<font color="#000000">&nbsp;läpi kulkeva virta, punainen käyrä on virtakierroksen 1 virta ja vihreä virtakierroksen 2 virta.</font>
</p>
<p align="left">
</p>
<p align="left">
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tr>
<td valign="top"><tt>&gt; &nbsp;&nbsp;</tt></td>
<td width="100%" valign="top"><b><font color="#FF0000">with(plots):</font></b>
</td></tr>
</table>
</p>
<p align="left">
<tt><pre><font color="#0000FF" size="2">Warning, the name changecoords has been redefined<br />
</font></pre></tt>
</p>
<p align="left">
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tr>
<td valign="top"><tt>&gt; &nbsp;&nbsp;</tt></td>
<td width="100%" valign="top"><b><font color="#FF0000">display(plot(virrat[1]-virrat[2], 0..0.2, color=blue),<br />
plot(virrat[1], 0..0.2, color=red), <br />
plot(virrat[2], 0..0.2, color=green)); <br />
</font></b>
</td></tr>
</table>
</p>
<p align="center">
<img src="images/virtap522.gif" width="448" height="448" alt="[Maple Plot]" />
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Vaihtovirta näyttää jakautuvan hieman eri vaiheisena rinnakkaisten kondensaattorien kesken. Alkuehdon aikaansaama epävakaisuus näyttää vaimenevan varsin nopeasti. &nbsp;</font>
</p>
<p align="left">
</p>
<p align="left">
<b><font color="#000000" size="5">Kytkentä 2</font></b>
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Tarkastellaan seuraavaksi RLC-virtapiiriä, jossa käämin rinnalle on kytketty sama vastus ja kondensaattori kuin edellä. Alla kuva piirin kokoonpanosta. </font>
</p>
<p align="left">
</p>
<p align="center">
<img src="images/virtap523.gif" width="266" height="117" alt="[Maple Bitmap]" />
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Huomaa virtakierrosten valinta! Pyrimme pitämään laskennan helpottamiseksi yhtälöryhmän kertaluvun mahdollisimman alhaalla virtakierrosvalinnoilla. Kirchhoffin lain perusteella saatavat yhtälöt ovat derivoidussa muodossa.</font>
</p>
<p align="left">
</p>
<p align="center">
<img src="images/virtap524.gif" width="480" height="76" alt="L[1]*diff(I(t),`$`(t,2))+R[1]*diff(I[1](t)+I[2](t),t)+(I[1](t)+I[2](t))/C[1] = E[0]*omega*cos(omega*t)" align="middle" />
<font color="#000000">,</font>
</p>
<p align="center">
<img src="images/virtap525.gif" width="532" height="72" alt="R[1]*diff(I[1](t)+I[2](t),t)+(I[1](t)+I[2](t))/C[1]+R[2]*diff(I[2](t),t)+I[2](t)/C[2] = E[0]*omega*cos(omega*t)" align="middle" />
<font color="#000000">.</font>
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Muodostetaan yhtälöryhmä ja tuntemattomat: </font>
</p>
<p align="left">
</p>
<p align="left">
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tr>
<td valign="top"><tt>&gt; &nbsp;&nbsp;</tt></td>
<td width="100%" valign="top"><b><font color="#FF0000">ryhma:= L[1]*diff(i[1](t), t$2)+R[1]*(diff(i[1](t), t)+diff(i[2](t), t))+(i[1](t)+i[2](t))/C[1]=E[0]*omega*cos(omega*t), R[1]*(diff(i[2](t), t)+diff(i[2](t), t))+(i[1](t)+i[2](t))/C[1]+R[2]*diff(i[2](t), t)+i[2](t)/C[2]=E[0]*omega*cos(omega*t);</font></b>
</td></tr>
</table>
</p>
<p align="left">
<img src="images/virtap526.gif" width="583" height="63" alt="ryhma := L[1]*diff(i[1](t),`$`(t,2))+R[1]*(diff(i[1](t),t)+diff(i[2](t),t))+(i[1](t)+i[2](t))/C[1] = E[0]*omega*cos(omega*t), 2*R[1]*diff(i[2](t),t)+(i[1](t)+i[2](t))/C[1]+R[2]*diff(i[2](t),t)+i[2](t)/..." />
<br />
<img src="images/virtap527.gif" width="476" height="59" alt="ryhma := L[1]*diff(i[1](t),`$`(t,2))+R[1]*(diff(i[1](t),t)+diff(i[2](t),t))+(i[1](t)+i[2](t))/C[1] = E[0]*omega*cos(omega*t), 2*R[1]*diff(i[2](t),t)+(i[1](t)+i[2](t))/C[1]+R[2]*diff(i[2](t),t)+i[2](t)/..." />
</p>
<p align="left">
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tr>
<td valign="top"><tt>&gt; &nbsp;&nbsp;</tt></td>
<td width="100%" valign="top"><b><font color="#FF0000">tuntemattomat2:= {i[1](t), j[2](t)};</font></b>
</td></tr>
</table>
</p>
<p align="center">
<img src="images/virtap528.gif" width="230" height="27" alt="tuntemattomat2 := {j[2](t), i[1](t)}" />
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Muodostetaan alkuehto ja annetaan vakioille arvot. Kytkettävä lähde on verkkovirtalähde ja sitä ennen piirissä ei ollut virtaa. Vakiot ovat samat kuin edellä. </font>
</p>
<p align="left">
</p>
<p align="left">
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tr>
<td valign="top"><tt>&gt; &nbsp;&nbsp;</tt></td>
<td width="100%" valign="top"><b><font color="#FF0000">ryhma2:= subs({L[1]=0.1, R[1]=5, C[1]=0.001, R[2]=3, C[2]=0.0005, E[0]=230, omega=2*Pi*50}, {ryhma});</font></b>
</td></tr>
</table>
</p>
<p align="left">
<img src="images/virtap529.gif" width="564" height="63" alt="ryhma2 := {13*diff(i[2](t),t)+.1e4*i[1](t)+3000.000000*i[2](t) = 23000*Pi*cos(100*Pi*t), .1*diff(i[1](t),`$`(t,2))+5*diff(i[1](t),t)+5*diff(i[2](t),t)+.1e4*i[1](t)+.1e4*i[2](t) = 23000*Pi*cos(100*Pi*t)..." />
<br />
<img src="images/virtap530.gif" width="642" height="63" alt="ryhma2 := {13*diff(i[2](t),t)+.1e4*i[1](t)+3000.000000*i[2](t) = 23000*Pi*cos(100*Pi*t), .1*diff(i[1](t),`$`(t,2))+5*diff(i[1](t),t)+5*diff(i[2](t),t)+.1e4*i[1](t)+.1e4*i[2](t) = 23000*Pi*cos(100*Pi*t)..." />
</p>
<p align="left">
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tr>
<td valign="top"><tt>&gt; &nbsp;&nbsp;</tt></td>
<td width="100%" valign="top"><b><font color="#FF0000">alkuehto:= i[1](0)=0, D(i[1])(0)=0, i[2](0)=0;</font></b>
</td></tr>
</table>
</p>
<p align="center">
<img src="images/virtap531.gif" width="310" height="27" alt="alkuehto := i[1](0) = 0, D(i[1])(0) = 0, i[2](0) = 0" />
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Ratkaistaan yhtälöryhmä numeerisesti.</font>
</p>
<p align="left">
</p>
<p align="left">
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tr>
<td valign="top"><tt>&gt; &nbsp;&nbsp;</tt></td>
<td width="100%" valign="top"><b><font color="#FF0000">rtk2:= dsolve({ryhma2[], alkuehto}, tuntemattomat, type=numeric, output=listprocedure);</font></b>
</td></tr>
</table>
</p>
<p align="center">
<img src="images/virtap532.gif" width="858" height="45" alt="rtk2 := [t = proc (t) option `Copyright (c) 1993 by the University of Waterloo. All rights reserved.`; t end proc, i[1](t) = proc (t) local res, solnproc, outpoint, ndsol; option `Copyright (c) 2000 by..." />
</p>
<p align="left">
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tr>
<td valign="top"><tt>&gt; &nbsp;&nbsp;</tt></td>
<td width="100%" valign="top"><b><font color="#FF0000">virrat:= subs(rtk2, [i[1](t), i[2](t)]):</font></b>
</td></tr>
</table>
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Yhtälöryhmän ratkaiseminen </font>
<b><font color="#000000">dsolve</font></b>
<font color="#000000">-komennolla symbolisesti on myös mahdollista, mutta tämän on hyvin hidasta ja sieventämisessä joudutaan helposti vaikeuksiin.</font>
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Piirretään kuvaaja virralle piirin eri osissa. Punainen käyrä on jaetun käämin </font>
<img src="images/virtap533.gif" width="33" height="32" alt="L[1]" align="middle" />
<font color="#000000">&nbsp;läpi kulkeva virta, vihreä virtakierroksen 2 virta ja sininen on näiden summavirta.</font>
</p>
<p align="left">
</p>
<p align="left">
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tr>
<td valign="top"><tt>&gt; &nbsp;&nbsp;</tt></td>
<td width="100%" valign="top"><b><font color="#FF0000">display(plot(virrat[1], 0..0.2, color=red), <br />
plot(virrat[2], 0..0.2, color=green), <br />
plot(virrat[1]+virrat[2], 0..0.2, color=blue));</font></b>
</td></tr>
</table>
</p>
<p align="center">
<img src="images/virtap534.gif" width="448" height="448" alt="[Maple Plot]" />
</p>
<p align="left">
<table width="100%" border="0" cellpadding="0" cellspacing="0">
<tr>
<td valign="top"><tt>&gt; &nbsp;&nbsp;</tt></td>
<td width="100%" valign="top"></td></tr>
</table>
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Tuloksesta näkyy, miten käämi pyrkii vastustamaan virran muutosta ja sen läpi kulkeva virta on hyvin vähäinen. Vastaavasti se indusoi jopa primääripiiriä voimakkaamman värähtelyn lisäyspiiriin.</font>
</p>
<p align="left">
</p>
<p align="left">
<i><b><font color="#000000">Tehtävä</font></b></i>
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Tarkastellaan vielä alla kuvattua piiriä, jossa on samanaikaisesti molemmat edellä tarkastellut lisäkytkennät. <br />
</font>
</p>
<p align="center">
<img src="images/virtap535.gif" width="227" height="210" alt="[Maple Bitmap]" />
</p>
<p align="left">
</p>
<p align="left">
<font color="#000000">Muodosta kolmen differentiaaliyhtälön ryhmä yllä esitetylle piirille. Käytä apuna kuvaan piirrettyjä virtakierroksia. Piirrä kuvaajat virran kululle eri komponenteissa. Vihje: Olet oikeilla jäljillä, mikäli virtakierron 3 värähtely on vastakkaisvaiheista kiertojen 1 ja 2 kanssa.</font>
</p>
<p align="left">
</p>

<p><b>Linkkejä</b></p>


 <a href="../xml_mpl/virtap3.xml">vaihtovirtapiirin pakotettu värähtely</a><br/>
<a href="../xml/algvke.xml">epähomogeeninen vakiokertoiminen lineaariyhtälö</a><br/>
 <a href="../xml/ryhyle.xml">differentiaaliyhtälöryhmä</a><br/>

<p align="left">
<i><font color="#000000">JP &amp; SKK &amp; MS 12.07.2001</font></i>
</p>



</body>
</html>
