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 t≥t1
pochodne(1)=(E-R1*x(1))/L1; //opis funkcji
i0=1;
t0=0; //definicja parametrów początkowych
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);
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.
the_mariov