II rząd-.docx

(64 KB) Pobierz

Politechnika Lubelska w Lublinie

Laboratorium Metod Numerycznych

Temat: Równania różniczkowe

II rzędu

Nr. 10

Data:  07.01.2013 r.

Głuszczak Paweł

 

Grupa:

EIST 3.1

Ocena:

 

1.Cel ćwiczenia.

Analiza stanów nieustalonych w obwodach elektrycznych.

2.Schemat układu.

3.Dane:

        E=10V                                                          L1=L2=0,1H

        R1=R2=10Ω                                                  t1=2(L1/ R1)=0,02s

5.Dla czasu 0d t=0 do t1.

a) skrypt

 

clc;

clear;

xdel;

 

function [pochodne]=stany(t, x)

pochodne(1)=(E-R1*x(1)-R1*x(2)-R2*x(1))/L2;

pochodne(2)=(E-R1*x(1)-R1*x(2))/L1;                //opis funkcji

endfunction

 

E=10;

R1=10;

R2=10;                                   //definicja parametrów

L1=0.1;

L2=0.1;

 

t0=0;

i20=0;                             //definicja parametrów początkowych

i30=0;

 

t=0:0.0001:0.1;                   //definicja wektora czasu

 

roz=ode([i20;i30;],t0,t,stany);   //polecenie rozwiązania układu równań

 

i2=roz(1,:);

i3=roz(2,:);

uR2=R2*i2;        //wyznaczenie napięć

 

subplot(411);

plot2d(t,i2,2,leg='i2(t)', axesflag=5);      //wykresy

subplot(412);

plot2d(t,i3,2,leg='i3(t)', axesflag=5);

subplot(413);

plot2d(t,[uR2' ],style=[6],leg='uR2(t)',axesflag=5);

 

 

b) wykresy

6.Dla czasu tt1

a) skrypt

clc;

clear;

xdel;

 

function [pochodne]=stany(t, x)

        pochodne(1)=(E-R1*x(1))/L1;          //opis funkcji

endfunction

 

E=10;

R1=10;

R2=10;                                   //definicja parametrów

L1=0.1;

L2=0.1;

 

i0=1;

t0=0;                       //definicja parametrów początkowych

 

t=0:0.0001:0.1;                   //definicja wektora czasu

 

roz=ode([i0],t0,t,stany);   //polecenie rozwiązania układu równań

 

i=roz(1,:);

uR1=R1*i;       //wyznaczenie napięć

uL1=E-uR1;

 

subplot(311);

plot2d(t,i,2,leg='i(t)', axesflag=5);

subplot(312);

plot2d(t,[uR1' uL1'],style=[4,9],leg='uR1(t)@uL1(t)',axesflag=5);

 

b) wykresy

7.Wnioski.

Zastosowanie programu Scilab umożliwia szybkie rozwiązanie złożonych układów RL,RC i RLC.W ćwiczeniu zostały wyznaczone przebiegi wszystkich prądów oraz napięć na poszczególnych elementach.W stanie początkowym przełącznik pierwszy i drugi jest załączony, więc mamy dwa równania, natomiast po upływie czasu t1 przełącznik drugi zostaje rozwarty i analizujemy tylko jedno oczko, zatem mamy trzecie równanie.

Zgłoś jeśli naruszono regulamin