Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 Klasyfikacja obwodów elektrycznych i metod ich badania Badanie obwodu elektrycznego polega na wyznaczenie pewnych wielkoÅ›ci niewiadomych charakteryzujÄ…cych dany obwód. Rozróżniamy nastÄ™pujÄ…ce typy zagadnieÅ„: 1). Analiza obwodu polega na ocenie jego odpowiedzi na różne sygnaÅ‚y wejÅ›ciowe i na okreÅ›leniu jego wÅ‚aÅ›ciwoÅ›ci 2). Synteza obwodu sprowadza siÄ™ do okreÅ›lenia jego struktury w zależnoÅ›ci od przyjÄ™tych, realizowanych charakterystyk 3). Zagadnienie aproksymacji polega na przybliżeniu wymaganych charakterystyk projektowych przez charakterystyki fizycznie realizowalne Metody badania obwodu można podzielić na: 1). Sieciowe, gdy dana jest peÅ‚na struktura obwodu i parametry jego elementów a poszukujemy rozpÅ‚ywu prÄ…du i rozkÅ‚adów napięć na poszczególnych elementach sieci elektrycznej. 2). Zaciskowe, gdy obwód jest traktowany jako dwójnik, czwórnik lub ogólnie n-wrotnik; wówczas okreÅ›la siÄ™ zależnoÅ›ci pomiÄ™dzy wielkoÅ›ciami zwiÄ…zanymi z zaciskami obwodu bez wnikania w jego strukturÄ™ wewnÄ™trznÄ…. Obwód albo sieć jest kombinacjÄ… elementów poÅ‚Ä…czonych z zewnÄ™trznymi zródÅ‚ami. yródÅ‚a wytwarzajÄ… w sieci sygnaÅ‚y wejÅ›ciowe lub wymuszenia. WynikajÄ…ce stÄ…d napiÄ™cia i prÄ…dy w różnych miejscach sieci sÄ… jej sygnaÅ‚ami wyjÅ›ciowymi albo odpowiedziami. Można powiedzieć, że analiza sieci polega na wyznaczeniu odpowiedzi danej sieci na zadane sygnaÅ‚y wejÅ›ciowe, zaÅ› synteza sieci polega na takim jej zaprojektowaniu, aby uzyskać żądane odpowiedzi na zadane sygnaÅ‚y. Graficznym opisem struktury sieci oraz rodzaju i funkcji jej elementów jest schemat ideowy. Natomiast graficznym opisem wÅ‚aÅ›ciwoÅ›ci zaciskowych sieci jest schemat blokowy. i y(t)=Lx (t) x(t)=i(t) y(t)=L i (t) u L u L wyjÅ›cie wejÅ›cie Schemat ideowy cewki idealnej (bezrezystancyjnej) Schemat blokowy i wymuszenie (prÄ…d) SygnaÅ‚ wyjÅ›ciowy jest proporcjonalny uL odpowiedz (napiÄ™cie samoindukcji) do pochodnej sygnaÅ‚u wejÅ›ciowego di uL = L dt SygnaÅ‚ wyjÅ›ciowy jest proporcjonalny do pochodnej sygnaÅ‚u wejÅ›ciowego Ten sam schemat blokowy może przedstawiać kondensator jeÅ›li: x(t)=u(t), zaÅ› y(t)=C(du/dt); schemat ideowy bÄ™dzie inny. 1 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 W schemacie ideowym sygnaÅ‚y: wejÅ›ciowy i wyjÅ›ciowy sÄ… wielkoÅ›ciami fizycznymi tj. napiÄ™ciami i prÄ…dami, a schemat pokazuje nie tylko ich wzajemnÄ… zależność, ale także typ i rodzaj użytych elementów. Schematy blokowe nie dajÄ… żadnej informacji na temat struktury i rodzaju elementów sieci, a ich wielkoÅ›ci wejÅ›ciowe i wyjÅ›ciowe nie muszÄ… być wielkoÅ›ciami fizycznymi. Obwody elektryczne i ich elementy dzielimy na: 1). Liniowe speÅ‚niajÄ… zasadÄ™ superpozycji: odpowiedz obwodu liniowego na jednoczesne dziaÅ‚anie kilku wymuszeÅ„ jest równa sumie algebraicznej odpowiedzi na każde wymuszenie z osobna. 2). Nieliniowe nie speÅ‚niajÄ… zasady superpozycji, opisujÄ… je równania nieliniowe, w których wielkoÅ›ci elektryczne i ich pochodne wystÄ™pujÄ… w potÄ™dze różnej od jeden. 3). Stacjonarne zÅ‚ożone z elementów o wartoÅ›ciach niezmiennych w czasie. 4). Niestacjonarne nazywane parametrycznymi. 5). Pasywne zÅ‚ożone z elementów pasywnych R, L, C (energia pobrana przez te elementy jest wiÄ™ksza lub równa 0). 6). Aktywne jeÅ›li w skÅ‚ad obwodu wchodzi chociaż jeden element aktywny np. zródÅ‚o, dioda tunelowa, tranzystor, wzmacniacz operacyjny. 7). O parametrach skupionych gdy mogÄ… być pominiÄ™te zjawiska falowe przy przepÅ‚ywie sygnałów. 8). O parametrach rozÅ‚ożonych gdy potrzebny jest skoÅ„czony czas na przeniesienie sygnaÅ‚u z jednego koÅ„ca obwodu do drugiego (np. linia dÅ‚uga zwana transmisyjnÄ…). 9). Odwracalne majÄ… takie same wÅ‚aÅ›ciwoÅ›ci niezależnie od sposobu poÅ‚Ä…czenia i od biegunowoÅ›ci przyÅ‚ożonego napiÄ™cia. 10). Nieodwracalne np. dioda, tranzystor. 11). Dwójniki majÄ… dwa zaciski. 12). Wielobiegunniki majÄ… n zacisków (np. tranzystor ma 3 zaciski). 13). Elementy idealne charakteryzuje jeden rodzaj procesów: a) wytwarzane energii (zródÅ‚a); b) rozpraszanie (rezystory); c) akumulacja (w cewkach i kondensatorach) 14). Elementy rzeczywiste charakteryzujÄ… dwa lub trzy rodzaje procesów z a), b), c). Klasyfikacja sygnałów i obwodów elektrycznych SygnaÅ‚y. Funkcje opisujÄ…ce zmienność w czasie wielkoÅ›ci fizycznych bÄ™dziemy nazywać przebiegami czasowymi tych wielkoÅ›ci lub sygnaÅ‚ami. Możemy zatem powiedzieć, że napiÄ™cie okreÅ›lone wzorem: u(t) = Um sin(É t +Èu ) ma przebieg sinusoidalny lub że ta powyższa zależność okreÅ›la sygnaÅ‚ sinusoidalny. (nazwy sinusoidalny używamy w odniesieniu do powyższego sygnaÅ‚u, niezależnie od tego czy w zapisie posÅ‚ugujemy siÄ™ funkcjÄ… sinus, czy cosinus, tzn. niezależnie od wartoÅ›ci fazy poczÄ…tkowej Èu) 2 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 Nazwy sygnaÅ‚ używamy zwÅ‚aszcza wtedy, kiedy chodzi o podkreÅ›lenie rodzaju zmiennoÅ›ci w czasie, a nie rodzaju wielkoÅ›ci fizycznej. Mówimy np.: sygnaÅ‚y staÅ‚e, sygnaÅ‚y okresowo zmienne, sygnaÅ‚y o skoÅ„czonej energii itp. Nie precyzujÄ…c najczęściej czy chodzi tu o sygnaÅ‚y napięć, prÄ…dów, siÅ‚ elektromotorycznych itp. Istnieje wiele rodzajów sygnałów, np.: " sygnaÅ‚ radiowy, " sygnaÅ‚ optyczny, " sygnaÅ‚ ultradzwiÄ™kowy, " sygnaÅ‚ elektryczny. SygnaÅ‚ami elektrycznymi sÄ… różne napiÄ™cia i prÄ…dy w sieci elektrycznej nazywanej obwodem lub ukÅ‚adem elektrycznym. PojÄ™cie ukÅ‚ad jest ogólne i można je odnieść do wielu dziedzin nauki i techniki, np.: " ukÅ‚ad elektroniczny, " ukÅ‚ad transportowy, " ukÅ‚ad biologiczny, " ukÅ‚ad planetarny. W elektrotechnice ukÅ‚ad jest prostym lub zÅ‚ożonym obwodem elektrycznym skÅ‚adajÄ…cym siÄ™ z oporników (rezystorów), cewek, kondensatorów i zródeÅ‚ energii. SygnaÅ‚y elektryczne sÄ… funkcjami czasu zwiÄ…zanymi zbiorem równaÅ„ wynikajÄ…cych z praw fizycznych (praw Kirchhoffa). Zjawiska fizyczne z dziedziny elektrycznoÅ›ci i magnetyzmu opisujÄ… zależnoÅ›ci matematyczne z dość dużą dokÅ‚adnoÅ›ciÄ…, dlatego rozważania teoretyczne majÄ… w elektrotechnice dużo wiÄ™ksze znaczenie niż w innych dziedzinach techniki. Ponieważ oporniki (rezystory) sÄ… elementami mnożącymi, cewki i kondensatory elementami różniczkujÄ…cymi i caÅ‚kujÄ…cymi, to obwód elektryczny można traktować jako ukÅ‚ad realizujÄ…cy wymienione operacje matematyczne. SygnaÅ‚y zaÅ› sÄ… dowolnymi funkcjami powiÄ…zanymi równaniami uwzglÄ™dniajÄ…cymi wzajemne poÅ‚Ä…czenia elementów. Można wiÄ™c powiedzieć, że obwody elektryczne sÄ… wykorzystywane do przetwarzania różnych sygnałów. Obwody elektryczne sÄ… ukÅ‚adami analogowymi i czÄ™sto mogÄ… być zastÄ…pione przez ukÅ‚ady cyfrowe (komputery), w których sygnaÅ‚y wejÅ›ciowy i wyjÅ›ciowy sÄ… ciÄ…gami liczb czyli sygnaÅ‚ami dyskretnymi. SygnaÅ‚ noÅ›nik wiadomoÅ›ci umożliwiajÄ…cy jej przesyÅ‚anie na odlegÅ‚ość lub w czasie (rejestracja); może mieć postać umownego znaku (np. rysunku, liter) lub przebiegu wielkoÅ›ci fizycznej, którego co najmniej jeden parametr (np. ksztaÅ‚t, czÄ™stotliwość, amplituda) zależy od przesyÅ‚anej nim wiadomoÅ›ci; rozróżnia siÄ™ sygnaÅ‚y elektryczne (np. zmieniajÄ…ce siÄ™ napiÄ™cie lub natężenie prÄ…du), akustyczne (zmieniajÄ…ca siÄ™ czÄ™stotliwość dzwiÄ™ku), optyczne (zmieniajÄ…ce siÄ™ natężenie lub barwa Å›wiatÅ‚a); stosuje siÄ™ też podziaÅ‚ sygnałów na analogowe (ciÄ…gÅ‚e), w przypadku których wielkość reprezentujÄ…ca wiadomość może przyjmować dowolne wartoÅ›ci ze zbioru nieskoÅ„czenie wielu wartoÅ›ci, i dyskretne (nieciÄ…gÅ‚e) 3 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 wielkość może przyjmować tylko okreÅ›lone wartoÅ›ci ze skoÅ„czonego ich zbioru, np. zakodowane w postaci cyfr (sygnaÅ‚ cyfrowy), najczęściej zer i jedynek (sygnaÅ‚ binarny). yródÅ‚o: Encyklopedia Multimedialna PWN 98; opublikowano w: J. Izydorczyk, G. PÅ‚onka, G. Tyma: Teoria Sygnałów, Wydawnictwo HELION, 1999 lub krótko SygnaÅ‚ proces zmian w czasie stanu fizycznego dowolnego obiektu, sÅ‚użący do wizualizacji, rejestracji i przesyÅ‚ania wiadomoÅ›ci (informacji). WÅ‚asnoÅ›ci sygnałów rozważa teoria sygnałów. ZwiÄ…zek pomiÄ™dzy naturÄ… fizycznÄ… sygnałów i zawartÄ… w nich informacjÄ… rozważa teoria informacji. 4 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 5 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 6 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 7 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 UkÅ‚ady UkÅ‚ad to inaczej: a) system o Å›ciÅ›le powiÄ…zanych ze sobÄ… elementach (zbiór elementów wraz z ich relacjami); b) relacja wzajemna, zależność różnych elementów (wielkoÅ›ci, przedmiotów, zjawisk); c) porzÄ…dek rozmieszczenia poszczególnych elementów wzglÄ™dem siebie PodziaÅ‚ ukÅ‚adów UkÅ‚ad otwarty to ukÅ‚ad, na który mogÄ… wpÅ‚ywać zdarzenia spoza ukÅ‚adu. UkÅ‚ad zamkniÄ™ty to ukÅ‚ad, na który zewnÄ™trzne zdarzenia nie majÄ… wpÅ‚ywu. W praktyce spotyka siÄ™ ukÅ‚ady bÄ™dÄ…ce poÅ‚Ä…czeniem ukÅ‚adu otwartego i ukÅ‚adu zamkniÄ™tego. UkÅ‚ad dynamiczny zawiera elementy i/lub przepÅ‚ywy zmienne w czasie. Wyróżniamy tu ukÅ‚ady: a) stabilny odsuniÄ™ty od stanu równowagi dynamicznej wraca do niego samorzutnie, np. wahadÅ‚o (zob. ujemne sprzężenie zwrotne) b) labilny (niestabilny) odchylony od stanu poczÄ…tkowego już do niego nie wraca, ale oddala siÄ™ od niego coraz dalej (np. Å›nieg na stoku i lawina) UkÅ‚ad statyczny nie zmienia siÄ™ w czasie. UkÅ‚ad analogowy UkÅ‚ad analogowy to ukÅ‚ad, w którym zwiÄ…zki pomiÄ™dzy sygnaÅ‚em wejÅ›ciowym a sygnaÅ‚em wyjÅ›ciowym możemy zapisać poprzez równanie różniczkowe lub tzw. transmitancjÄ™ Laplace a uzyskanÄ… na podstawie równaÅ„ różniczkowych. S X(S) Y(S) H(S) 8 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 UkÅ‚ad cyfrowy UkÅ‚ad cyfrowy to ukÅ‚ad, w którym zwiÄ…zki pomiÄ™dzy sygnaÅ‚em wejÅ›ciowym a sygnaÅ‚em wyjÅ›ciowym możemy zapisać za pomocÄ… równaÅ„ różnicowych lub za pomocÄ… tzw. transmitancji Z, którÄ… można wyprowadzić na podstawie tych równaÅ„. .... 0101 1110 1001 1110 1100 1011.... Z X(Z) Y(Z) H(Z) Elementy WystÄ™pujÄ…ce w ukÅ‚adach elementy możemy scharakteryzować za pomocÄ…: równaÅ„ różniczkowych, transmitancji operatorowej, transmitancji widmowej, charakterystyk widmowych, a także poprzez odpowiedzi tych elementów na wymuszenia. Podstawowe elementy to element proporcjonalny, element różniczkujÄ…cy, i element caÅ‚kujÄ…cy. Klasyfikacja ukÅ‚adów UkÅ‚ad liniowy to ukÅ‚ad, w którym wystÄ™pujÄ…ce elementy sÄ… liniowe (idealny rezystor, cewka, kondensator). SygnaÅ‚y przechodzÄ…ce przez te elementy poddawane sÄ… liniowym operacjom matematycznym takim jak: mnożenie sygnaÅ‚u przez staÅ‚y czynnik, różniczkowanie oraz caÅ‚kowanie. W ukÅ‚adach liniowych obowiÄ…zuje zasada superpozycji, zgodnie z którÄ… sygnaÅ‚ na wyjÅ›ciu ukÅ‚adu można wyznaczyć jako sumÄ™ sygnałów wyjÅ›ciowych pochodzÄ…cych od wszystkich sygnałów wejÅ›ciowych. x1 y1 x2 y2 UkÅ‚ad liniowy x3 y3 x1 x2 y=y1+y2+y3 UkÅ‚ad liniowy x3 UkÅ‚ad nieliniowy to ukÅ‚ad, w którym co najmniej jeden element jest nieliniowy. W ukÅ‚adach nieliniowych nie jest speÅ‚niona zasada superpozycji - nie można sygnaÅ‚u wyjÅ›ciowego rozdzielić na skÅ‚adniki pochodzÄ…ce od różnych sygnałów wejÅ›ciowych. 9 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 SygnaÅ‚y sinusoidalnie zmienne WÅ›ród sygnałów elektrycznych zmiennych w czasie duże znaczenie praktyczne majÄ… sygnaÅ‚y przemienne okresowe. Warunek okresowoÅ›ci funkcji (sygnaÅ‚u) można wyrazić zależnoÅ›ciÄ…: f (t + T ) = f (t) Jeżeli warunek okreÅ›lony równaniem (2.1) nie jest speÅ‚niony, to sygnaÅ‚ jest nieokresowy lub aperiodyczny. SygnaÅ‚ okresowy nazywamy przemiennym, jeżeli pole powierzchni ograniczonej przebiegiem sygnaÅ‚u w ciÄ…gu okresu T jest równe zeru, tzn. jeÅ›li jest speÅ‚niony warunek T f (t) d t = 0 +" 0 Szczególne miejsce w elektrotechnice zajmujÄ… sinusoidalne (harmoniczne) sygnaÅ‚y prÄ…du i napiÄ™cia. Głównym tego powodem jest ich naturalna powszechność w przyrodzie oraz Å‚atwość wytwarzania wynikajÄ…ca ze Å›cisÅ‚ego zwiÄ…zku ruchu obrotowego z funkcjami trygonometrycznymi sinus i cosinus. u u(t)=Umsin(Ét+Èu) Um t Èu 0 T Rys. SygnaÅ‚ sinusoidalny napiÄ™ciowy Energia elektryczna docierajÄ…ca do naszych domów otrzymywana jest w generatorach synchronicznych, które sÄ… maszynami elektrycznymi wirujÄ…cymi. Generatory te wytwarzajÄ… napiÄ™cia sinusoidalnie zmienne w czasie o unormowanych parametrach. Ogólna postać napiÄ™cia harmonicznego u(t) = Um sin(É t +È ) [V] u przy czym: t czas w sekundach [s], u(t) wartość chwilowa napiÄ™cia w woltach [V], Um wartość szczytowa napiÄ™cia (amplituda) w woltach [V], Èu faza poczÄ…tkowa napiÄ™cia w radianach [rad], É t +È faza napiÄ™cia w chwili t w radianach [rad], u É = 2 Ä„ / T = 2 Ä„ f pulsacja w radianach na sekundÄ™ [rad/s], f =1/ T czÄ™stotliwość w hercach [Hz], T okres w sekundach [s]. 10 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 SygnaÅ‚y okresowe niesinusoidalnie Rozważmy ukÅ‚ad liniowy, na wejÅ›cie którego podano sygnaÅ‚ okresowy niesinusoidalny. Elementy liniowe tworzÄ…ce ten ukÅ‚ad bÄ™dÄ… oddziaÅ‚ywać na sygnaÅ‚ wejÅ›ciowy, dokonujÄ…c na nim operacji matematycznych charakteryzujÄ…cych siÄ™ liniowoÅ›ciÄ…. Do operacji tych należy zaliczyć mnożenie sygnaÅ‚u przez staÅ‚y czynnik, różniczkowanie oraz caÅ‚kowanie. UkÅ‚ad liniowy można wiÄ™c traktować jako zbiór liniowych przetworników sygnaÅ‚u, bÄ™dÄ…cych kombinacjÄ… podukÅ‚adów proporcjonalnych, różniczkujÄ…cych i caÅ‚kujÄ…cych. DziÄ™ki tym wÅ‚asnoÅ›ciom odpowiedz ukÅ‚adu liniowego na sygnaÅ‚ okresowy niesinusoidalny jest również okresowa. wymuszenie odpowiedz UKAAD okresowe okresowa LINIOWY niesinusoidalne niesinusoidalna Rys. UkÅ‚ad liniowy jako przetwornik sygnaÅ‚u okresowego niesinusoidalnego Ä… f( ),f(t) A t[s] T/2 T nr 1 Ä„ Ä„ 2 Ä… [rad] -A Rys. PrzykÅ‚ad sygnaÅ‚u antysymetrycznego ZwiÄ…zek pomiÄ™dzy sygnaÅ‚ami wejÅ›ciowymi i wyjÅ›ciowymi w ukÅ‚adach liniowych UkÅ‚ady analogowe przetwarzajÄ… wejÅ›ciowe sygnaÅ‚y analogowe dajÄ…c na wyjÅ›ciu również sygnaÅ‚ analogowy zależny od: " sygnaÅ‚u wejÅ›ciowego " parametrów ukÅ‚adu liniowego Duże znaczenie praktyczne majÄ… ukÅ‚ady liniowe, dla których sygnaÅ‚ wejÅ›ciowy x(t) i wyjÅ›ciowy y(t) sÄ… zwiÄ…zane ze sobÄ… równaniem różniczkowym: dn y(t) dn-1y(t) dy(t) dk x(t) dk -1x(t) dx(t) an n + an-1 n-1 + ...+ a1 + a0 y(t) = ak k + ak-1 -1 + ...+ a1 + a0x(t) dt dt dt dt dt dtk Równanie różniczkowe opisuje zwiÄ…zek pomiÄ™dzy sygnaÅ‚em wejÅ›ciowym i jego pochodnymi oraz sygnaÅ‚em wyjÅ›ciowym i jego pochodnymi. Jeżeli znamy postać sygnaÅ‚u wejÅ›ciowego x(t) to znamy również jego pochodne i równanie różniczkowe możemy przeksztaÅ‚cić do postaci: dn y(t) dn-1y(t) dy(t) an n + an-1 n-1 + ..........+ a1 + a0 y(t) = f (t) 11 dt dt dt Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 Współczynniki an,an-1,...,a0 sÄ… bezpoÅ›rednio zwiÄ…zane z parametrami opisujÄ…cymi ukÅ‚ad. Wyznaczanie odpowiedzi y(t) W celu wyznaczenia odpowiedzi y(t) należy rozwiÄ…zać równanie różniczkowe. RozwiÄ…zanie ogólne równania różniczkowego niejednorodnego skÅ‚ada siÄ™ z sumy dwóch rozwiÄ…zaÅ„: " RozwiÄ…zania szczególnego (caÅ‚ki szczególnej) równania różniczkowego niejednorodnego (czyli skÅ‚adowej ustalonej yu(t) lub inaczej nazywanej skÅ‚adowej wymuszonej yw(t) ) " RozwiÄ…zania ogólnego (caÅ‚ki szczególnej) równania różniczkowego jednorodnego (czyli skÅ‚adowej przejÅ›ciowej yp(t) lub inaczej nazywanej skÅ‚adowej swobodnej ys(t) ) dn y(t) dn-1y(t) dy(t) an + an-1 +... + a1 + a0 y(t) = 0 dtn dtn-1 dt 12 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 Reakcja elementów R, L, C na skok jednostkowy 1 dla t e" 0 NapiÄ™cie na rezystorze nie powoduje przesuniÄ™cia wzglÄ™dem Å„Å‚ µ (t) = òÅ‚0 dla t < 0 prÄ…du, wiÄ™c ksztaÅ‚t obydwu sygnałów jest taki sam z ół 1 dla t > 0 Å„Å‚ dokÅ‚adnoÅ›ciÄ… do czynnika skalujÄ…cego R ôÅ‚1 1 µ (t) = dla t = 0 òÅ‚ 2 ôÅ‚ Cewka bÄ™dzie przeciwdziaÅ‚ać zmianom prÄ…du w obwodzie, ół0 dla t < 1 wytwarzajÄ…c siÅ‚Ä™ elektromotorycznÄ… samoindukcji przeciwnie 1/2 skierowanÄ… do wzrastajÄ…cego napiÄ™cia t NienaÅ‚adowany idealny kondensator po podÅ‚Ä…czeniu zasilania można traktować w zasadzie jak zwarcie, teoretycznie prÄ…d zmienia siÄ™ skokowo od zera do nieskoÅ„czonoÅ›ci (w praktyce rezystor ogranicza wartość tego prÄ…du, a jeÅ›li jest on bardzo maÅ‚y to może siÄ™ zdarzyć, że wyparujÄ… przewody) JeÅ›li elementy R, L i C sÄ… poÅ‚Ä…czone szeregowo to ksztaÅ‚t prÄ…du w takim obwodzie, bÄ™dzie zależaÅ‚ od wartoÅ›ci poszczególnych elementów. Stany nieustalone SygnaÅ‚y elektryczne, przy okreÅ›lonej strukturze obwodu i zródÅ‚ach, nie uzyskujÄ… natychmiast swoich wartoÅ›ci ustalonych. MiÄ™dzy jednym stanem stabilnym a drugim wystÄ™pujÄ… wahania napięć i prÄ…dów. PrzyczynÄ… wystÄ™powania zjawisk przejÅ›ciowych w obwodzie (trwajÄ…cych od chwili wystÄ…pienia zakłóceÅ„, aż do chwili ustalenia siÄ™ zjawisk) jest każda zmiana struktury poÅ‚Ä…czeÅ„ lub parametrów elementów, wchodzÄ…cych w skÅ‚ad obwodu. W wielu przypadkach stany nieustalone sÄ… zjawiskami niepożądanymi (np. niepożądane sÄ… zjawiska przejÅ›ciowe wystÄ™pujÄ…ce przy zwarciach i przy wÅ‚Ä…czaniu napięć w obwodach elektrycznych). W innych zaÅ› przypadkach stany nieustalone sÄ… normalnym stanem pracy urzÄ…dzeÅ„ (np. ukÅ‚ady automatycznej regulacji). CechÄ… charakterystycznÄ… zjawisk w obwodach elektrycznych prÄ…du staÅ‚ego lub sinusoidalnie zmiennego jest to, że generatory zasilajÄ…ce te obwody narzucajÄ… sposób zmiennoÅ›ci czasowej prÄ…dów i napięć. W przypadku generatorów prÄ…du staÅ‚ego, napiÄ™cia i prÄ…dy w obwodzie sÄ… wielkoÅ›ciami staÅ‚ymi, a w przypadku generatorów prÄ…du sinusoidalnego, napiÄ™cia i prÄ…dy zmieniajÄ… siÄ™ sinusoidalnie. Tego rodzaju stan obwodu nazywamy ustalonym lub stacjonarnym. W obwodach elektrycznych spotyka siÄ™ także zjawiska spowodowane zmianÄ… dokonanÄ… w obwodzie, jak np. wÅ‚Ä…czenie zródÅ‚a energii do obwodu, zwarcie części obwodu itp. Tego 13 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 rodzaju zmiany w obwodzie wywoÅ‚ane sÄ… dziaÅ‚aniem czynników zewnÄ™trznych, obwód zostaje wyprowadzony ze stanu równowagi, a w obwodzie wytwarza siÄ™ stan zwany przejÅ›ciowym lub nieustalonym. W stanie nieustalonym prÄ…dy i napiÄ™cia w obwodzie zmieniajÄ… siÄ™ inaczej niż siÅ‚y elektromotoryczne generatorów zasilajÄ…cych ukÅ‚ad. Teoretycznie stan przejÅ›ciowy trwa nieskoÅ„czenie dÅ‚ugo, jednak praktycznie po upÅ‚ywie dostatecznie dÅ‚ugiego czasu zjawiska przejÅ›ciowe zanikajÄ… i obwód osiÄ…ga stan ustalony.Przy analizie stanów nieustalonych w obwodach elektrycznych liniowych zarówno napiÄ™cie u, jak i prÄ…d i przedstawiane sÄ… w postaci sumy dwóch skÅ‚adników, a mianowicie skÅ‚adowej ustalonej i skÅ‚adowej przejÅ›ciowej. WyrażajÄ… to zależnoÅ›ci u = uu + u oraz i = iu + ip p SkÅ‚adowe ustalone (uu, iu) zwiÄ…zane sÄ… ze stanem ustalonym, zaÅ› skÅ‚adowe przejÅ›ciowe (up, ip) ze stanem przejÅ›ciowym obwodu. OkreÅ›lanie warunków poczÄ…tkowych Do analizy zjawisk w stanie nieustalonym konieczna jest znajomość stanu poczÄ…tkowego obwodu. WartoÅ›ci wybranych napięć i prÄ…dów w stanie poczÄ…tkowym nazywamy warunkami poczÄ…tkowymi. Warunki poczÄ…tkowe mogÄ… być zerowe lub niezerowe. Jeżeli sÄ… zerowe, to obwód byÅ‚ na poczÄ…tku w stanie bezenergetycznym. Ponieważ gromadzić energiÄ™ mogÄ… tylko cewki i kondensatory, to oprócz zródeÅ‚, napiÄ™cia na wszystkich pojemnoÅ›ciach i prÄ…dy pÅ‚ynÄ…ce przez wszystkie indukcyjnoÅ›ci w chwili t=0s decydujÄ… o zachowaniu siÄ™ obwodu w czasie pózniejszym. Prawa komutacji Rozważamy tu chwilÄ™ t=0s, gdyż zakÅ‚adamy, że wÅ‚aÅ›nie (1) (2) R wtedy nastÄ…piÅ‚a zmiana topologii ukÅ‚adu (np. doÅ‚Ä…czono t = 0 s zródÅ‚o, zwarto element itp.). L Zmiana topologii możliwa jest dziÄ™ki zamykaniu lub e2(t) E1 otwieraniu Å‚Ä…czników. ZakÅ‚adamy, że sÄ… one idealne, to i(t) znaczy, że zamykajÄ… siÄ™ lub otwierajÄ… natychmiast, w (0) stanie zamkniÄ™cia majÄ… zerowÄ… rezystancjÄ™, a w stanie otwarcia nieskoÅ„czonÄ… oraz, że nie wystÄ™puje w nich zjawisko Å‚uku elektrycznego. Zjawiskiem komutacji, a wiÄ™c procesem zmiany struktury ukÅ‚adu rzÄ…dzÄ… prawa komutacji. WynikajÄ… one z ważnego prawa fizyki, prawa zachowania energii. Prawa komutacji można wypowiedzieć nastÄ™pujÄ…co: Energia w polu magnetycznym cewki nie może zmienić siÄ™ skokowo (I prawo komutacji) oraz Energia w polu elektrycznym kondensatora nie może zmienić siÄ™ skokowo (II prawo komutacji) Z matematycznego punktu widzenia funkcje energii cewki i kondensatora 14 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 2 2 iL (t) uC (t) wL (t) = L oraz wC (t) = C muszÄ… być ciÄ…gÅ‚e wzglÄ™dem czasu. 2 2 Ze wzorów na wartoÅ›ci chwilowe energii w cewce i w kondensatorze wynika również ciÄ…gÅ‚ość funkcji prÄ…du iL(t) oraz napiÄ™cia uC(t). Ponadto, ze wzorów na wartoÅ›ci chwilowe strumienia magnetycznego w cewce ¨(t) = L Å" iL (t) oraz Å‚adunku elektrycznego w kondensatorze q(t) = C Å" uC (t) wnioskujemy o ciÄ…gÅ‚oÅ›ci strumienia ¨(t) oraz Å‚adunku q(t) PodsumowujÄ…c powyższe rozważania, prawa komutacji możemy dla chwili opisać nastÄ™pujÄ…cymi równoÅ›ciami: WL (0- ) = WL (0+ ), ¨(0- ) = ¨(0+ ), iL (0- ) = iL (0+ ) I prawo komutacji WC (0- ) = WC (0+ ), q(0- ) = q(0+ ), uC (0- ) = uC (0+ ) II prawo komutacji Dodajmy jeszcze, że chociaż prÄ…d pÅ‚ynÄ…cy przez cewkÄ™ nie może zmieniać siÄ™ skokowo, to napiÄ™cie na cewce może siÄ™ zmieniać skokowo. Podobnie jest z prÄ…dem pÅ‚ynÄ…cym przez kondensator - on również może zmieniać siÄ™ skokowo, ale napiÄ™cie na kondensatorze musi być ciÄ…gÅ‚Ä… funkcjÄ… czasu. Dla rezystora możliwa jest skokowa zmiana zarówno prÄ…du jak i napiÄ™cia, chyba że rezystor poÅ‚Ä…czony jest szeregowo z cewkÄ… lub równolegle z pojemnoÅ›ciÄ…. Metody analizy stanów nieustalonych 1. Metoda klasyczna zwiÄ…zana z klasycznymi metodami rozwiÄ…zywania równaÅ„ różniczkowych i różniczkowo-caÅ‚kowych 2. Metoda operatorowa oparta o przeksztaÅ‚cenie Laplace a 3. Metoda zmiennych stanu zwiÄ…zana z zastosowaniem funkcji macierzy Równanie charakterystyczne Na podstawie równania dn x(t) dn-1x(t) d x(t) + an-1 n-1 + L + a1 + a0 x(t) = f (t) n d t d t d t piszemy tzw. równanie charakterystyczne, które jest równaniem algebraicznym wzglÄ™dem zmiennej pomocniczej s. Równanie to ma postać sn + an-1 sn-1 + L + a1 s + a0 = 0 Równanie to ma n jedno- lub wielokrotnych pierwiastków rzeczywistych lub zespolonych sprzężonych. Z każdym pierwiastkiem si lub z grupÄ… pierwiastków zwiÄ…zany jest kolejny skÅ‚adnik rozwiÄ…zania dla xp(t). Zależność skÅ‚adnika od wÅ‚asnoÅ›ci danego zestawiono w tabeli. 15 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 Tabela 1. SkÅ‚adniki rozwiÄ…zania xs(t) w zależnoÅ›ci od rodzaju pierwiastków si funkcje wchodzÄ…ce w skÅ‚ad xs(t) wÅ‚asnoÅ›ci pierwiastków si - pojedynczy, rzeczywisty xi (t) = C1 esit si - podwójny, rzeczywisty xi (t) = (C1 + C2 t) esit si - potrójny, rzeczywisty xi (t) = (C1 + C2 t + C3 t2 ) esi t ... itd. si,j = Ä… Ä… j ² xij (t) = eÄ… t (C1 cos(² t) + C2 sin(² t)) dwa pierwiastki zespolone sprzężone, pojedyncze si,j = Ä… Ä… j ² xij (t) = eÄ… t [ (C1 + C2 t) cos(² t) dwa pierwiastki zespolone + (C3 + C4 t) sin(² t) ] sprzężone, podwójne ... itd. UkÅ‚ady II rzÄ™du Załóżmy równanie charakterystyczne drugiego stopnia (ukÅ‚ad RLC) R 1 s2 + s + = 0 L LC a pierwiastki równania charakterystycznego dane sÄ… wzorami 2 R R 1 ëÅ‚ öÅ‚ s1,2 = - Ä… ìÅ‚ ÷Å‚ - = -Ä… Ä… ² 2L 2L LC íÅ‚ Å‚Å‚ 2 R R 1 1 ëÅ‚ öÅ‚ 2 przy czym Ä… = , ² = ìÅ‚ ÷Å‚ - = Ä… - 2L 2L LC LC íÅ‚ Å‚Å‚ Przy użyciu tych oznaczeÅ„ wyrażamy poszczególne pierwiastki równania w sposób s1 = nastÄ™pujÄ…cy -Ä… + ² , s2 = -Ä… - ² Jeżeli zaÅ‚ożymy, że indukcyjność L i pojemność C sÄ… staÅ‚e, to rezystancjÄ™ R można dobrać tak, że wyróżnik równania charakterystycznego może być dodatni, ujemny, a w przypadku granicznym - staje siÄ™ zerem. Zależnie od wartoÅ›ci rezystancji rozróżniamy wiÄ™c trzy przypadki. L R > 2 1. Przy C wielkość b przedstawia liczbÄ™ rzeczywistÄ…, przy czym wobec Ä…>² - obydwa pierwiastki sÄ… rzeczywiste i ujemne, s1= -Ä…+²<0, s2= - Ä…-²<0. Fizycznie odpowiada temu Å‚adowanie kondensatora ze zródÅ‚a napiÄ™cia staÅ‚ego poprzez rezystancjÄ™ R i indukcyjność L, majÄ…ce charakter aperiodyczny (nieokresowy). 16 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 i , u c , u L E u c i t 0 t 1 u L Rys. Przebiegi prÄ…du i napięć w obwodzie szeregowym RLC L R = 2 2. W przypadku granicznym przy C wielkość ² staje siÄ™ zerem, pierwiastki równania charakterystycznego sÄ… sobie równe i tworzÄ… jeden pierwiastek podwójny, s1=s2= -Ä…, rzeczywisty i ujemny. Fizycznie odpowiada temu Å‚adowanie kondensatora ze zródÅ‚a napiÄ™cia staÅ‚ego poprzez rezystancjÄ™ R i indukcyjność L, majÄ…ce charakter aperiodyczny krytyczny (nieokresowy krytyczny) L R < 2 .3. Przy C wielkość ² przedstawia liczbÄ™ urojonÄ…. Wprowadzamy oznaczenie ² = jÉ gdzie É jest już liczbÄ… rzeczywistÄ…, która speÅ‚nia równanie 1 Ä…2 + É2 = LC Obydwa pierwiastki równania charakterystycznego sÄ… zespolone sprzężone, równe odpowiednio s1 = -Ä… + jÉ, s2 = -Ä… - jÉ Fizycznie odpowiada temu Å‚adowanie kondensatora ze zródÅ‚a napiÄ™cia staÅ‚ego poprzez rezystancjÄ™ R i indukcyjność L takie, że przebiegi napiÄ™cia na kondensatorze i prÄ…du w funkcji czasu sÄ… oscylacyjne tÅ‚umione, w szczególnoÅ›ci sinusoidalne tÅ‚umione. Rys. Przebiegi oscylacyjne tÅ‚umione prÄ…du i napiÄ™cia na kondensatorze w ukÅ‚adzie szeregowym RLC 17 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 We wszystkich trzech przypadkach pierwiastki równania charakterystycznego leżą w lewej półpÅ‚aszczyznie, Re s1<0 oraz Re s2<0; w zwiÄ…zku z tym skÅ‚adowa przejÅ›ciowa odpowiedzi uCp(t) maleje do zera dla czasu dążącego do nieskoÅ„czonoÅ›ci. W przypadku 1 obydwa pierwiastki leżą na ujemnej części osi rzeczywistej, symetrycznie wzglÄ™dem punktu -Ä…, a w przypadku 2 obydwa pierwiastki tworzÄ… jeden podwójny równy -Ä…; w każdym z tych przypadków skÅ‚adowa przejÅ›ciowa uCp(t) maleje asymptotycznie do zera. Wreszcie w przypadku 3 pierwiastki sl i s2 sÄ… zespolone sprzężone, w konsekwencji skÅ‚adowa przejÅ›ciowa uCp(t) maleje oscylacyjnie do zera. Stany nieustalone Metoda klasyczna przykÅ‚ady obliczeniowe (rozwiÄ…zane podczas wykÅ‚adu) Zadanie 1. Dla obwodu przedstawionego na rys. 1 oblicz i narysuj przebiegi czasowe prÄ…du pÅ‚ynÄ…cego przez cewkÄ™ iL(t) oraz napięć uR (t) i uL(t) . Zbadaj wpÅ‚yw wartoÅ›ci i znaku warunku poczÄ…tkowego na cewce na odpowiedzi czasowe obwodu. Jak wartoÅ›ci elementów R oraz L wpÅ‚ywajÄ… na czas ustalania siÄ™ przebiegów? Dane: U = 20 V , R = 4 &! , L = 20 mH . uR (t) uL (t) a) b) (1) (2) R L t = 0s R L VU U iL(t) (0) Rys. 1. a) obwód do zadania 1, b) przystosowany do analizy w PSpice Rys. 1.1. Wykresy czasowe sygnałów przy zerowym warunku poczÄ…tkowym Rys. 1.2. Wykresy czasowe przy dodatnim warunku poczÄ…tkowym I0 > U / R 18 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 Rys. 1.3. Wykresy czasowe przy ujemnym warunku poczÄ…tkowym Rys. 1.4. WpÅ‚yw indukcyjnoÅ›ci na szybkość ustalania siÄ™ przebiegów ( I0 = 0A ) Rys. 1.5. WpÅ‚yw rezystancji na szybkość ustalania siÄ™ przebiegów ( I0 = 0A ) Zadanie 2. W ukÅ‚adzie z rys. 2a w chwili t = 0 s otwarto klucz K, przez co odÅ‚Ä…czono zasilanie z rzeczywistego zródÅ‚a napiÄ™cia staÅ‚ego. Znajdz wzory, jakimi opisane sÄ… po komutacji: napiÄ™cie na kondensatorze i prÄ…d w cewce. SporzÄ…dz wykresy czasowe tych sygnałów przejÅ›ciowych. Dane: R1 = 60 &!, R2 = 40 &! , L = 10 H, C = 25 mF, E = 50 V . 19 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 uR2(t) a) b) K (1) (2) R2 t = 0 s (1) R2 (2) E uC (t) uC (t) L C C L uL(t) R1 i(t) i(t) (0) (0) Rys. 2. Obwód do zadania 2, a) przed komutacjÄ…, b) po komutacji Rys. 2.1. Wykresy dla obwodu z zadania 2 Zadanie 3. W ukÅ‚adzie przedstawionym na rys. 3a w chwili t = 0 s otwarto Å‚Ä…cznik, przez co odÅ‚Ä…czono zasilanie sinusoidalne. Wyznacz przebieg prÄ…du w cewce oraz napiÄ™cia na kondensatorze w stanie nieustalonym. Dane: e(t) = 70sin(10t + 30o ) V, R1 =14 &!, R2 = 2 &!, L = 0,5 H, C = 20 mF iL(t) iL(t) R1 (1) a) b) (1) t = 0s i2(t) L i1(t) L R2 R2 uC (t) uC (t) (2) (2) e(t) C C (0) (0) Rys. 3. Obwód do zadania 3, a) przed komutacjÄ…, b) po komutacji Rys. 3.1. Wykres prÄ…du pÅ‚ynÄ…cego przez cewkÄ™ z rys. 3 20 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 Rys. 3.2. Wykres napiÄ™cia na kondensatorze z rys. 3 PrzeksztaÅ‚cenie Laplace a Wprowadzenie WÅ›ród metod czÄ™stotliwoÅ›ciowych badania ukÅ‚adów analogowych najczęściej znajduje zastosowanie metoda przeksztaÅ‚cenia Laplace'a. PodstawowÄ… cechÄ… tej metody jest algebraizacja obliczeÅ„ stanów dynamicznych. Algebraizacja polega na zastÄ…pieniu dziaÅ‚ania różniczkowania funkcji czasu przez pomnożenie funkcji zmiennej zespolonej zwanej transformatÄ… przez parametr zespolony s, oraz na zastÄ…pieniu caÅ‚kowania funkcji czasu w granicach od 0 do t przez podzielenie transformaty przez tenże parametr. MetodÄ™ przeksztaÅ‚cenia Laplace'a zaliczamy do metod operatorowych, a zespół twierdzeÅ„ i reguÅ‚ zwiÄ…zanych z zastosowaniem przeksztaÅ‚cenia Laplace'a nazywamy rachunkiem operatorowym. Zalety rachunku operatorowego 1. Prostota dokonywania operacji na równaniach algebraicznych 2. Wprowadzenie warunków poczÄ…tkowych wprost i na samym poczÄ…tku 3. Możliwość rozwiÄ…zywania przypadków nie posiadajÄ…cych rozwiÄ…zania przy metodzie klasycznej 4. Możliwość korzystania z tablic 5. PrzeksztaÅ‚cenie Laplace a dla funkcji spotykanych w technice jest wzajemnie jednoznaczne PrzeksztaÅ‚cenie Laplace a W przeksztaÅ‚ceniu Laplace'a, zwanym też transformacjÄ… Laplace'a, rozpatruje siÄ™ dwie funkcje: 1. funkcjÄ™ f(t) argumentu rzeczywistego (zmiennej rzeczywistej) t; funkcjÄ™ f(t) nazywamy funkcjÄ… oryginalnÄ…, oryginaÅ‚em lub też funkcjÄ… czasu; w elektrotechnice argument t oznacza zazwyczaj czas, 2. funkcjÄ™ F(s) argumentu zespolonego (zmiennej zespolonej) zwanego też parametrem zespolonym, okreÅ›lonÄ… wzorem 21 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 " F(s) = (t) e-s t d t +"f 0 zwanÄ… transformatÄ… funkcji czasu, jej funkcjÄ… przeksztaÅ‚conÄ… lub obrazem. W odniesieniu do funkcji czasu, o której mówimy, że jest transformowalna wedÅ‚ug Laplace'a, czynimy nastÄ™pujÄ…ce zaÅ‚ożenia: 1. znika dla argumentów ujemnych, tzn. f(t) =0 dla t<0, 2. jest jednoznacznie okreÅ›lona w caÅ‚ym przedziale od 0 do +" oraz jest w tym przedziale ciÄ…gÅ‚a, z wyjÄ…tkiem co najwyżej skoÅ„czonej liczby punktów nieciÄ…gÅ‚oÅ›ci pierwszego rodzaju, tzn. takich, w których nastÄ™puje skok funkcji o skoÅ„czonÄ… wartość, 3. wzrasta co do wartoÅ›ci bezwzglÄ™dnej nie szybciej niż funkcja wykÅ‚adnicza, tzn. dla danej funkcji f(t) można dobrać takÄ… liczbÄ™ dodatniÄ… M oraz takÄ… staÅ‚Ä… Ä… nieujemnÄ…, że dla wszelkich wartoÅ›ci argumentu t zachodzi f (t) < M eÄ… t Rozpatrzmy przykÅ‚ad wyznaczania transformaty Laplace a funkcji wykÅ‚adniczej Ä… t f (t) = e przy staÅ‚ej a>0 " " -1 1 F(s) = e-s t d t = e(Ä… -s) t = (*) +"eÄ… t s - Ä… 0 s - Ä… 0 Podany wynik otrzymuje siÄ™ przy zaÅ‚ożeniu Re(s) > Ä… , a wiÄ™c gdy à jest wiÄ™ksza od Ä…. Obliczana caÅ‚ka jest wówczas zbieżna. Gdy natomiast à < Ä… , caÅ‚ka ta jest rozbieżna. Wielkość Ä… nazywamy odciÄ™tÄ… zbieżnoÅ›ci transformaty. Gdy zatem parametr zespolony à jest poÅ‚ożony na pÅ‚aszczyznie zmiennej zespolonej na prawo od prostej à = Ä… równolegÅ‚ej do osi urojonej, wówczas transformata istnieje i jest okreÅ›lona wzorem (*) PrzeksztaÅ‚cenie okreÅ›lone wzorem " F(s) = (t) e-s t d t +"f 0 oznaczamy symbolem ! F(s) = !{f(t)} i nazywamy przeksztaÅ‚ceniem prostym, bo sÅ‚uży ono do wyznaczania transformaty danej funkcji czasu. Jeżeli natomiast dana jest transformata F(s), a szukamy funkcji czasu f(t), wówczas piszemy zależność odwrotnÄ… -1 f(t) = ! {F(s)} stanowiÄ…cÄ… zapis przeksztaÅ‚cenia odwrotnego Laplace'a. 22 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 Można wykazać, że przeksztaÅ‚cenie odwrotne, napisane w postaci caÅ‚ki wzglÄ™dem zmiennej zespolonej s, ma postać c+ j" 1 f (t) = F(s) es t d s +" 2 Ä„ j c- j" gdzie c jest liczbÄ… rzeczywistÄ… dodatniÄ…, nie mniejszÄ… od odciÄ™tej zbieżnoÅ›ci transformaty, a caÅ‚kowanie przebiega wzdÅ‚uż prostej równolegÅ‚ej do osi urojonej ukÅ‚adu współrzÄ™dnych zgodnie z kierunkiem wzrostu argumentu urojonego. Podstawowe wÅ‚asnoÅ›ci PrzeksztaÅ‚cenia Laplace a 1. liniowość PodstawowÄ… wÅ‚asnoÅ›ciÄ… przeksztaÅ‚cenia Laplace'a jest jego liniowość; innymi sÅ‚owy przeksztaÅ‚cenie Laplace'a speÅ‚nia zasadÄ™ superpozycji. W odniesieniu do dwóch funkcji czasu, przy zaÅ‚ożeniu, że 1 i 2 sÄ… skalarami, wÅ‚asność liniowoÅ›ci przeksztaÅ‚cenia Laplace a możemy wyrazić nastÄ™pujÄ…co !{1 f1(t) + 2 f2 (t)} = 1 F1(s) + 2 F2 (s) w której 1 i 2 sÄ… skalarami. 2. splot Jednym z podstawowych pojęć rachunku operatorowego jest splot dwóch funkcji czasu. DefiniujÄ… go równoważnie poniższe caÅ‚ki oznaczone niewÅ‚aÅ›ciwe " " f (t) = f1(t) *f (t) = f1(Ä )f2 (t -Ä )dÄ = f1(t -Ä )f2 (Ä )dÄ +" +" 2 -" -" Ponieważ w ramach przeksztaÅ‚cenia Laplace'a rozważamy funkcje, które znikajÄ… dla chwil t ujemnych, to można dla takich funkcji zawÄ™zić przedziaÅ‚ caÅ‚kowania zmiennej we wzorze i okreÅ›lić splot inaczej t t f (t) = f1(t) *f (t) = f1(Ä )f2 (t -Ä )dÄ = f1(t -Ä )f2 (Ä )dÄ +" +" 2 0 0 Do transformaty splotu odnosi siÄ™ twierdzenie Borela wyrażone zależnoÅ›ciÄ… !{f1(t) * f2 (t)} = F1(s)F2 (s) Twierdzenie to, którego dowód pomijamy, możemy wysÅ‚owić w nastÄ™pujÄ…cy sposób: Transformata splotu dwóch funkcji czasu równa siÄ™ iloczynowi transformat tych funkcji. Transformaty Laplace'a typowych sygnałów majÄ… postać ilorazu bÄ™dÄ…cego funkcjÄ… wymiernÄ… wzglÄ™dem parametru s L(s) F(s) = N(s) W najprostszym przypadku zakÅ‚adamy, że mianownik N(s) nie ma pierwiastków wielokrotnych, a stopieÅ„ licznika L(s) jest mniejszy od stopnia mianownika N(s). Zgodnie z twierdzeniem o rozkÅ‚adzie, rozpatrywanej transformacie odpowiada funkcja czasu 23 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 n L(sk ) f (t) = exp(sk t) " N'(sk ) k=1 gdzie n oznacza stopieÅ„ wielomianu N(s), a sk sÄ… pierwiastkami równania N(s) = 0 3. Pochodna Pochodna splotu 4. CaÅ‚ka F(s) d îÅ‚ d îÅ‚ L [+" f (t)]= L f (t)Å‚Å‚ = sF(s) - f (0) L {f1(t) " f2(t)}Å‚Å‚ = sF1(s) Å"F2(s) ïÅ‚dt śł ïÅ‚dt śł s ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ Wybrane transformaty Laplace a 24 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 25 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 26 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 Zasilanie impulsowe obwodu W badaniach sygnałów duże znaczenie majÄ… transformaty impulsów jednorazowych. a) Pojedynczy impuls prostokÄ…tny W najprostszym przypadku mamy do czynienia z jednorazowym impulsem prostokÄ…tnym, który powstaje w chwili t=a, znika w chwili t=b, a przez czas trwania impulsu wynoszÄ…cy (b-a) ma staÅ‚Ä… amplitudÄ™ A. Impuls taki możemy traktować jako sumÄ™ algebraicznÄ… dwóch funkcji jednostkowych pomnożonych przez amplitudÄ™ A i opóznionych odpowiednio o a oraz b jednostek czasu. Powyższe rozważania nasuwajÄ… zapis jednorazowego impulsu prostokÄ…tnego nastÄ™pujÄ…cÄ… funkcjÄ… czasu f (t) = A[µ(t - a) - µ (t - b)] a) b) c) f(t) f(t) f(t) A·µ(t-a) A·µ(t-a) A A A t t 0 t 0 a b 0 a b a b -A·µ(t-b) -A -A f(t) f(t) A·µ(t-a) A A 0 t 0 t a b a b -A·µ(t-b) -A f (t) = A[µ (t - a) - µ (t - b)] Wspomniano, że jednorazowy impuls prostokÄ…tny zapisany jest funkcjÄ… czasu f (t) = A[µ (t - a) - µ(t - b)] W oparciu o przytoczone uprzednio wzory stwierdzamy, że transformata tego impulsu wynosi A F(s) = [exp(-as) - exp(-bs)] s 27 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 b) Pojedynczy impuls sinusoidalny Analogicznie, w stosunku do pojedynczego impulsu prostokÄ…tnego, można wyznaczyć transformatÄ™ jednorazowego impulsu sinusoidalnego półfalowego, który daje siÄ™ wyrazić sumÄ… dwóch funkcji sinusoidalnych: podstawowej i opóznionej o pół okresu. f(t) f(t) A A A·sin(Ét)·µ(t) 0 t 0 a a 2a 3a t -A f(t) f(t) A A A·sin(Ét)·µ(t) 0 t 0 a a 2a 3a t -A A·sin(É(t-a))·µ(t-a) f(t) f(t) A A A·sin(Ét)·µ(t) 0 t 0 a a 2a 3a t -A A·sin(É(t-a))·µ(t-a) f (t) = sin(Ét)µ (t) + sin(É(t - a))µ (t - a) Jednorazowy impuls sinusoidalny można zapisać w postaci f (t) = sin(Ét)µ (t) + sin(É(t - a))µ (t - a) W oparciu o przytoczone uprzednio wzory stwierdzamy, że transformata tego impulsu wynosi É F(s) = [1+ exp(-as)] 2 s2 + É Metoda operatorowa Patrz - zadania przedstawione na wykÅ‚adzie 28 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 Metoda zmiennych stanu Jeżeli badany jest ukÅ‚ad liniowy z wymuszeniami uporzÄ…dkowanymi w wektor wymuszeÅ„ u1(t) îÅ‚ Å‚Å‚ ïÅ‚u (t)śł 2 ïÅ‚ śł u(t) = ïÅ‚ śł : ïÅ‚ śł (t)śł ïÅ‚u p ðÅ‚ ûÅ‚ oraz odpowiedziami uporzÄ…dkowanymi w wektor odpowiedzi y1(t) îÅ‚ Å‚Å‚ ïÅ‚ y2 (t)śł ïÅ‚ śł y(t) = ïÅ‚ śł : ïÅ‚ śł yq (t)śł ïÅ‚ ðÅ‚ ûÅ‚ to wektor odpowiedzi można wyznaczyć z zależnoÅ›ci y(t) = Cx(t) + Du(t) (MZS 1) w której wektor x1(t) îÅ‚ Å‚Å‚ ïÅ‚x (t)śł 2 ïÅ‚ śł x(t) = ïÅ‚ śł : ïÅ‚ śł (t)ûÅ‚ ðÅ‚xn jest nazywany wektorem stanu ukÅ‚adu, a jego elementy - zmiennymi stanu. Zmiennymi stanu nazywamy wielkoÅ›ci x1(t), x2(t), ..., xn(t), które należy zadać w chwili poczÄ…tkowej t = t0, aby przy zadanych wymuszeniach u1(t), u2(t), ..., up(t) okreÅ›lić jednoznacznie zachowanie siÄ™ ukÅ‚adu dla t > t0. W obwodach elektrycznych zmiennymi stanu sÄ… prÄ…dy pÅ‚ynÄ…ce przez cewki i napiÄ™cia na kondensatorach - wielkoÅ›ci te muszÄ… speÅ‚niać warunki poczÄ…tkowe wynikajÄ…ce z praw komutacji. Wektor stanu jest rozwiÄ…zaniem równania stanu d x(t) = Ax(t) + Bu(t) (MZS 2) dt które jako równanie różniczkowe pierwszego rzÄ™du musi mieć zadany wektor stanu poczÄ…tkowego x(t0 ) = x0 (MZS 3) Macierze w równaniach (MZS 1) i (MZS 2) okreÅ›la siÄ™ jako: A - macierz ukÅ‚adu, wymiar n × n, B - macierz wymuszeÅ„, wymiar n × p, C - macierz odpowiedzi, wymiar q × n, D - macierz transmisyjna, wymiar q × p. Różniczkowe równanie macierzowe (MZS 2) z warunkiem poczÄ…tkowym (MZS 3) dla chwili czasu t = 0 ma rozwiÄ…zanie w postaci 29 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 t x(t) = exp(At)x0 + exp(A(t -Ä ))Bu(Ä )dÄ (MZS 4) +" 0 przy czym funkcja wykÅ‚adnicza macierzy ukÅ‚adu zdefiniowana jest zależnoÅ›ciÄ… " (At)k exp(At) = (MZS 5) " k! k=0 RozwiÄ…zanie równania stanu można wyznaczyć również w postaci sumy skÅ‚adowej ustalonej i przejÅ›ciowej x(t) = xu (t) + xp (t) W prostych obwodach elektrycznych wyznaczenie rozwiÄ…zania dla skÅ‚adowej ustalonej nie nastrÄ™cza zwykle wiÄ™kszych trudnoÅ›ci. Z rozwiÄ…zania dla skÅ‚adowej ustalonej wynika również warunek poczÄ…tkowy dla tej skÅ‚adowej oraz dla skÅ‚adowej przejÅ›ciowej xp (0) = x(0) - xu (0) Dla skÅ‚adowej przejÅ›ciowej równanie stanu w postaci macierzowej można zapisać jako d xp (t) - Axp (t) = 0 dt i jego rozwiÄ…zanie ma postać xp (t) = exp(At)xp (0) Wyznaczanie funkcji exp(At) z zastosowaniem twierdzenia Sylvestera Rozpatrujemy macierz ukÅ‚adu A stopnia n, której elementy sÄ… znane, macierz jednostkowÄ… oznaczonÄ… przez 1 oraz skalar bÄ™dÄ…cy liczbÄ… rzeczywistÄ… lub zespolonÄ…. RóżnicÄ™ macierzy 1 - A nazywamy macierzÄ… charakterystycznÄ… macierzy kwadratowej A, a jej wyznacznik det(1 - A), który jest wielomianem stopnia n wzglÄ™dem , wielomianem charakterystycznym macierzy A. PrzyrównujÄ…c do zera ten wielomian, otrzymujemy równanie stopnia n wzglÄ™dem det(1 - A) = 0 zwane równaniem charakterystycznym macierzy kwadratowej A; równanie to można zapisać krótko Õ() = 0 Pierwiastki równania Õ() = 0 nazywamy pierwiastkami charakterystycznymi macierzy kwadratowej A lub częściej jej wartoÅ›ciami wÅ‚asnymi. Twierdzenie Sylvestera, w swojej podstawowej postaci, pozwala wyrazić dowolny wielomian macierzy kwadratowej A stopnia n w postaci wielomianu stopnia n-1 wzglÄ™dem macierzy A. Jest ono w zakresie algebry macierzy odpowiednikiem znanego z algebry wzoru interpolacyjnego Lagrange'a, który pozwala wyrazić wartość funkcji badanej w pewnym przedziale w zależnoÅ›ci od n znanych wartoÅ›ci tej funkcji w n punktach tego przedziaÅ‚u. Twierdzenie Sylvestera możemy stosować nie tylko do wielomianu wzglÄ™dem macierzy A, ale również do funkcji przestÄ™pnych, rozwijalnych w szereg nieskoÅ„czony. Do takich funkcji należy funkcja wykÅ‚adnicza macierzy A, okreÅ›lona zależnoÅ›ciÄ… (MZS 5). Twierdzenie Sylvestera pozwala na wyrażenie takich funkcji w postaci zamkniÄ™tej, tzn. prostszej niż za pomocÄ… szeregu nieskoÅ„czonego. Dla funkcji wykÅ‚adniczej macierzy kwadratowej A stopnia 30 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 n, majÄ…cej n pierwiastków charakterystycznych różnych od siebie, z twierdzenia Sylvestera wynika wzór "( 1 - A) s n exp(At) = (MZS 6) "exp( t) s`"r - r ) r "( s r=0 s`"r W przypadku szczególnym, w którym macierz kwadratowa A jest stopnia n = 2, wzór (MZS 6) przybiera prostÄ… postać (21 - A) e1t (11 - A) e2t exp(At) = + (MZS 7) 2 - 1 1 - 2 PrzykÅ‚ad rozgaÅ‚Ä™zionego obwodu RLC z niezerowym wymuszeniem [8] StosujÄ…c metodÄ™ zmiennych stanu, obliczyć przebiegi prÄ…du w cewce iL(t) oraz napiÄ™cia na kondensatorze uC(t) po zamkniÄ™ciu wyÅ‚Ä…cznika w ukÅ‚adzie jak na rysunku 1. Dane liczbowe: a) E = 14 V, R1 = 5 &!, R2 = 9 &!, C = 0.1 F, L = 1 H b) E = 13 V, R1 = 5 &!, R2 = 0.2 &!, C = 0.1 F, L = 0.1 H c) E = 200 V, R1 = 200 &!, R2 = 0.05 &!, C = 0.01 F, L = 0.1 H t=0 i1 L R1 iC i2 E C R2 Rys. 1. Schemat obwodu RLC z niezerowym wymuszeniem RozwiÄ…zanie Przypadek a) uC (t) = x1 Å„Å‚ Oznaczamy zmienne stanu òÅ‚i (t) = x2 ół L 0 îÅ‚ Å‚Å‚ Jak Å‚atwo zauważyć, warunki poczÄ…tkowe (dla t = 0): X0 = ïÅ‚0śł ðÅ‚ ûÅ‚ E îÅ‚ R2 Å‚Å‚ ïÅ‚ R1 + R2 śł 9 îÅ‚ Å‚Å‚ ïÅ‚ śł Stan ustalony obliczony metodÄ… klasycznÄ… Xu = = ïÅ‚1śł ïÅ‚ E śł ðÅ‚ ûÅ‚ ïÅ‚ R1 + R2 śł ðÅ‚ ûÅ‚ uC Å„Å‚ E R1i1 + uC = E Ò! i1 = - + ôÅ‚ R1 R1 ôÅ‚ ôÅ‚ di2 UkÅ‚ad równaÅ„ otrzymany z praw Kirchhoffa L + R2i2 = uC òÅ‚ dt ôÅ‚ ôÅ‚ duC i1 = i2 + i3 = C + i2 ôÅ‚ dt ół 31 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 duC i1 uC Å„Å‚ E = - - + ôÅ‚ dt C R1C R1C ôÅ‚ przeksztaÅ‚camy w równania stanu òÅ‚ R2 uC ôÅ‚di2 = - i2 + ôÅ‚ ół dt L L UwzglÄ™dniajÄ…c przyjÄ™te oznaczenia zmiennych stanu, możemy równania zapisać ostatecznie jako dx1 Å„Å‚ 1 1 E = - x1 - x2 + ôÅ‚ dt R1C C R1C ôÅ‚ òÅ‚ 1 R2 ôÅ‚dx2 = x1 - x2 ôÅ‚ ół dt L L îÅ‚- 1 1 Å‚Å‚ - ïÅ‚ śłîÅ‚x1(t) Å‚Å‚ îÅ‚ 1 0Å‚Å‚îÅ‚ x1(t) îÅ‚ Å‚Å‚ d R1C C ïÅ‚ śłïÅ‚EÅ‚Å‚ lub w postaci macierzowej + śłïÅ‚ R1C ïÅ‚x = ïÅ‚ (t)śł 0 0ûÅ‚ðÅ‚0 śł dt (t)śł ïÅ‚ 1 R2 śłðÅ‚x2 ûÅ‚ ïÅ‚ śł ðÅ‚ 2 ûÅ‚ ûÅ‚ - ðÅ‚ ïÅ‚ śł ðÅ‚ L L ûÅ‚ Odpowiedz czasowa ma dwie skÅ‚adowe: ustalonÄ… i przejÅ›ciowÄ… x(t) = xu (t) + xp (t) Dla skÅ‚adowej przejÅ›ciowej równanie stanu w postaci macierzowej można zapisać jako îÅ‚- 1 1 Å‚Å‚ - x1p (t) ïÅ‚ śłîÅ‚ Å‚Å‚ x1p (t) îÅ‚ Å‚Å‚ d R1C C ïÅ‚x (t)śł = ïÅ‚ 1 R2 śłïÅ‚ (t)śł dt 2p ïÅ‚ śł ðÅ‚ ûÅ‚ ðÅ‚x2p ûÅ‚ - ïÅ‚ śł ðÅ‚ L L ûÅ‚ x1p (t) x1p (t) îÅ‚ Å‚Å‚ Å‚Å‚ d îÅ‚- 2 -10 îÅ‚ Å‚Å‚ Dla danych zadania powyższe równanie przybiera postać ïÅ‚x (t)śł = ïÅ‚x (t)śł ïÅ‚ śł dt 1 9 2p 2p ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ SkÅ‚adowÄ… przejÅ›ciowÄ… policzymy nastÄ™pujÄ…co: xp (t) = exp(At)xp (0) 0 9 îÅ‚ Å‚Å‚ îÅ‚ Å‚Å‚ îÅ‚- 9 Å‚Å‚ gdzie skÅ‚adowa przejÅ›ciowa w zerze ma postać xp (0) = x(0) - xu (0) = - = ïÅ‚0śł ïÅ‚1śł ïÅ‚-1śł ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ Obliczenie macierzy exp(At) metodÄ… Sylvestera WartoÅ›ci wÅ‚asne macierzy 1, 2 obliczymy z równania charakterystycznego Å„Å‚ 0 üÅ‚ + 2 10 îÅ‚ Å‚Å‚ îÅ‚- 2 -10 Å‚Å‚ îÅ‚ Å‚Å‚ det[1 - A]= detòÅ‚ïÅ‚ - ïÅ‚ śłżł = detïÅ‚ -1 + 9śł = 1 - 9 ðÅ‚0 śł ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ ûÅ‚ ół þÅ‚ + 2 10 = = ( + 2)( + 9) + 10 = 2 + 11 + 28 -1 + 9 Jako rozwiÄ…zanie równania charakterystycznego otrzymujemy dwa pierwiastki rzeczywiste ujemne 2 + 11 + 28 = 0 " = 9, " = 3 1 = -4 , 2 = -7 Wskazuje to na aperiodyczny charakter przebiegów. UwzglÄ™dniajÄ…c wzór Sylvestera (MZS 7), otrzymujemy 32 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 îÅ‚- 7 0 2 -10 Å‚Å‚ e-4t Å„Å‚ Å‚Å‚ îÅ‚- üÅ‚ exp(At) = - Å" - òÅ‚ïÅ‚ śłżł + 3 0 - 7śł ïÅ‚ 1 - 9 ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ ół þÅ‚ îÅ‚- 4 0 2 -10 Å‚Å‚ e-4t îÅ‚- 5 10 Å‚Å‚ îÅ‚- 2 10 e-7t Å„Å‚ Å‚Å‚ îÅ‚- üÅ‚ e-7t Å‚Å‚ + Å" - = - + = òÅ‚ żł ïÅ‚ śł ïÅ‚-1 2 śł ïÅ‚-1 5 śł 3 0 - 4śł ïÅ‚ 1 - 9 3 3 ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ ół þÅ‚ 5 2 10 îÅ‚ Å‚Å‚ ïÅ‚3 e-4t - 3 e-7t 3 (e-7t - e-4t )śł = ïÅ‚ śł 1 5 2 ïÅ‚ - e-7t ) e-7t - e-4t śł (e-4t ðÅ‚ 3 3 3 ûÅ‚ Na podstawie wzoru xp (t) = exp(At)xp (0) otrzymujemy postać skÅ‚adowej przejÅ›ciowej 5 2 10 8 35 îÅ‚ îÅ‚ Å‚Å‚ ïÅ‚3 e-4t - 3 e-7t 3 (e-7t - e-4t )Å‚Å‚îÅ‚- 9 ïÅ‚3 e-7t - 3 e-4t śł śł Å‚Å‚ xp (t) = = ïÅ‚ śł ïÅ‚ śł ïÅ‚ śł 1 5 2 4 7 ïÅ‚ - e-7t ) e-7t - e-4t e-7t śłðÅ‚-1ûÅ‚ ïÅ‚ - e-4t śł (e-4t ðÅ‚ 3 3 3 ûÅ‚ ðÅ‚ 3 3 ûÅ‚ UwzglÄ™dniajÄ…c obie skÅ‚adowe, możemy zapisać (korzystajÄ…c ze wzoru x(t) = xu (t) + xp (t) ) 8 35 1 îÅ‚ Å‚Å‚ îÅ‚9 + (8e-7t - 35e-4t )Å‚Å‚ ïÅ‚3 e-7t - 3 e-4t śł 9 ïÅ‚ śł îÅ‚ Å‚Å‚ 3 x(t) = + = ïÅ‚ śł ïÅ‚ śł ïÅ‚1śł 4 7 1 ïÅ‚ - e-4t ðÅ‚ ûÅ‚ 1 + (4e-7t - 7e-4t ) śł ïÅ‚ śł e-7t ðÅ‚ 3 3 ûÅ‚ ðÅ‚ 3 ûÅ‚ RozwiÄ…zanie zadania stanowiÄ… przebiegi: 1 uC (t) = [9 + (8e-7t - 35e-4t )] V 3 1 iL (t) = [1 + (4e-7t - 7e-4t )] A 3 Przebiegi czasowe napiÄ™cia uC(t) i prÄ…du iL(t) pokazano na rys. 1.1. i rys. 1.2. OtrzymaliÅ›my przebiegi o charakterze aperiodycznym. 10 y ( t ) 0 9 0 0 t 2 Rys. 1.1. Przebieg napiÄ™cia na kondensatorze uC(t) 1.2 y ( t ) 1 1 0 0 t 2 Rys. 1.2. Przebieg prÄ…du w cewce iL(t) Przypadek b) PostÄ™pujemy analogicznie jak w przypadku a). Przyjmujemy takie same oznaczenia zmiennych stanu. Wektor stanu poczÄ…tkowego jest także taki sam. Po uwzglÄ™dnieniu danych zadania (rozważanego przypadku b) możemy zapisać wektor stanu ustalonego 33 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 E îÅ‚ R2 Å‚Å‚ ïÅ‚ îÅ‚ Å‚Å‚ R1 + R2 śł = 0.5 X = ïÅ‚ śł u ïÅ‚2.5śł E ïÅ‚ śł ðÅ‚ ûÅ‚ ïÅ‚ R1 + R2 śł ðÅ‚ ûÅ‚ SkÅ‚adowÄ… przejÅ›ciowÄ… policzymy analogicznie jak w przypadku a). Możemy zapisać równanie stanu dla tejże skÅ‚adowej îÅ‚- 1 1 Å‚Å‚ - x1p (t) ïÅ‚ śłîÅ‚ Å‚Å‚ x1p (t) îÅ‚ Å‚Å‚ d R1C C ïÅ‚x (t)śł = ïÅ‚ 1 R2 śłïÅ‚ (t)śł dt 2p ïÅ‚ śł ðÅ‚ ûÅ‚ ðÅ‚x2p ûÅ‚ - ïÅ‚ śł ðÅ‚ L L ûÅ‚ i dla danych przypadku b): x1p (t) x1p (t) îÅ‚ Å‚Å‚ Å‚Å‚ d îÅ‚- 2 -10 îÅ‚ Å‚Å‚ ïÅ‚x (t)śł = ïÅ‚x (t)śł ïÅ‚10 - 2 śł dt 2p 2p ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ Obliczamy wartoÅ›ci wÅ‚asne macierzy A, korzystajÄ…c z równania charakterystycznego Å„Å‚ 0 üÅ‚ + 2 10 îÅ‚ Å‚Å‚ îÅ‚- 2 -10 Å‚Å‚ îÅ‚ Å‚Å‚ det[1 - A]= detòÅ‚ïÅ‚ - ïÅ‚10 - 2 śłżł = detïÅ‚ -10 + 2śł = ðÅ‚0 śł ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ ûÅ‚ ół þÅ‚ + 2 10 = = ( + 2)( + 2) + 100 = 2 + 4 +104 -10 + 2 W wyniku rozwiÄ…zania równania charakterystycznego otrzymujemy dwa pierwiastki zespolone 2 + 4 +104 = 0 " = -400; " = Ä… j20 1 = -2 + j10 ; 2 = -2 - j10 Wskazuje to na oscylacyjny charakter przebiegów. KorzystajÄ…c ze wzoru Sylvestera, możemy zapisać Å„Å‚ 2 îÅ‚- Å‚Å‚ îÅ‚- 2 -10 üÅ‚ Å‚Å‚ e(-2+ j10)t - j10 0 exp(At) = - òÅ‚ïÅ‚ 0 - 2 - j10śł - ïÅ‚10 - 2 śłżł + j20 ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ ół þÅ‚ Å„Å‚ Å‚Å‚ îÅ‚- 2 -10 Å‚Å‚ e(-2- j10)t îÅ‚- 2 + j10 0 üÅ‚ + - = òÅ‚ żł ïÅ‚ śł ïÅ‚10 - 2 śł j20 0 - 2 + j10 ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ ół þÅ‚ Å‚Å‚ e(-2- j10)t îÅ‚ Å‚Å‚ e(-2+ j10)t îÅ‚- j10 10 j10 10 = - + = ïÅ‚ śł ïÅ‚-10 j10śł j20 -10 - j10 j20 ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ j10t Å‚Å‚ e-2t îÅ‚j10(e j10t j10t + e- ) 10(e- j10t - e ) = = ïÅ‚ śł j10t j10t j20 10(e - e-j10t ) j10(e + e- j10t )ûÅ‚ ïÅ‚ śł ðÅ‚ îÅ‚e j10t + e-j10t e j10t - e-j10t Å‚Å‚ - ïÅ‚ śł = e-2t ïÅ‚ j10t 2 j10t j10t j2 j10t śł ïÅ‚e - e- e + e- śł ïÅ‚ śł j2 2 ïÅ‚ śł ðÅ‚ ûÅ‚ KorzystajÄ…c z zależnoÅ›ci 34 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 jÉt jÉt jÉt e - e-jÉt e - e sinÉt = ; cosÉt = 2j 2 możemy ostatecznie zapisać îÅ‚e-2t cos10t - e-2tsin10tÅ‚Å‚ exp(At) = ïÅ‚ śł -2t ïÅ‚ ðÅ‚e sin10t e-2tcos10t śł ûÅ‚ UwzglÄ™dniajÄ…c wzór xp (t) = exp(At)xp (0) oraz fakt, że 0 0.5 îÅ‚ Å‚Å‚ îÅ‚ Å‚Å‚ îÅ‚- 0.5 Å‚Å‚ xp (0) = x0 (0) - xu (0) = - = ïÅ‚0śł ïÅ‚2.5śł ïÅ‚- 2.5śł ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ możemy zapisać skÅ‚adowÄ… przejÅ›ciowÄ… îÅ‚e-2t cos10t - e-2t sin10tÅ‚Å‚îÅ‚- 0.5 îÅ‚ Å‚Å‚ Å‚Å‚ e-2t (2.5sin10t - 0.5cos10t) xp (t) = ïÅ‚ śłïÅ‚ = ïÅ‚ śł -2t ïÅ‚ ûÅ‚ śł ðÅ‚e sin10t - e-2tcos10t śłðÅ‚- 2.5śł ïÅ‚- e-2t (0.5sin10t + 2.5cos10t)ûÅ‚ ûÅ‚ ðÅ‚ KorzystajÄ…c ze wzoru x(t) = xu (t) + xp (t) , możemy zapisać x(t) = xp (t) + xu (t) = îÅ‚ Å‚Å‚ îÅ‚0.5 + e-2t (2.5sin10t - 0.5cos10t)Å‚Å‚ 0.5 îÅ‚ Å‚Å‚ e-2t (2.5sin10t - 0.5cos10t) = + = ïÅ‚ śł ïÅ‚ śł ïÅ‚2.5śł ïÅ‚- e-2t (0.5sin10t + 2.5cos10t)ûÅ‚ śł ðÅ‚ ûÅ‚ ïÅ‚ - e-2t (0.5sin10t + 2.5cos10t)ûÅ‚ śł ðÅ‚ ðÅ‚2.5 Ostatecznie można zapisać poszukiwane przebiegi jako uC (t) = [0.5 + e-2t (2.5sin10t - 0.5cos10t)] V iL (t) = [2.5 - e-2t (0.5sin10t + 2.5cos10t)] A Przebiegi te zilustrowano na rys. 1.3 i rys. 1.4. OtrzymaliÅ›my przebiegi o charakterze oscylacyjnym. TÅ‚umienie jest stosunkowo duże z pulsacjÄ… É wynoszÄ…cÄ… 10 rad/s . u [V] c 2,5 2 1,5 1 0,5 0 0 0,2 0,4 0,6 0,8 1 1,2 1,4 1,6 1,8 2 t [s] -0,5 -1 Rys. 1.3. Przebieg napiÄ™cia na kondensatorze uC(t) 35 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 iL [A] 4,5 4 3,5 3 2,5 2 1,5 1 0,5 0 0 0,2 0,4 0,6 0,8 1 1,2 1,4 1,6 1,8 2 t [s] Rys. 1.4. Przebieg prÄ…du w cewce iL(t) Przypadek c) RozwiÄ…zujÄ…c, postÄ™pujemy analogicznie jak w przypadku b). Po uwzglÄ™dnieniu danych zadania (rozważanego przypadku) możemy zapisać wektor stanu ustalonego E îÅ‚ R2 Å‚Å‚ ïÅ‚ îÅ‚ Å‚Å‚ R1 + R2 śł E" 0.05 X = ïÅ‚ śł u ïÅ‚ śł E 1 ïÅ‚ śł ðÅ‚ ûÅ‚ ïÅ‚ R1 + R2 śł ðÅ‚ ûÅ‚ SkÅ‚adowÄ… przejÅ›ciowÄ… obliczymy analogicznie jak w przypadkach a) oraz b) i dla danych przypadku c) x1p (t) x1p (t) îÅ‚ Å‚Å‚ Å‚Å‚ d îÅ‚- 0.5 -100 îÅ‚ Å‚Å‚ ïÅ‚x (t)śł = ïÅ‚x (t)śł ïÅ‚ dt 10 - 0.5śłðÅ‚ 2p ûÅ‚ 2p ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ Obliczamy wartoÅ›ci wÅ‚asne macierzy A, korzystajÄ…c z równania charakterystycznego Å„Å‚ 0 üÅ‚ + 0.5 100 îÅ‚ Å‚Å‚ îÅ‚- 0.5 -100 Å‚Å‚ îÅ‚ Å‚Å‚ det[1 - A]= detòÅ‚ïÅ‚ - = detïÅ‚ = ïÅ‚ 10 - 0.5śłżł ðÅ‚ -10 + 0.5śł ðÅ‚0 śł ðÅ‚ ûÅ‚ ûÅ‚ ûÅ‚ ół þÅ‚ + 0.5 100 = = ( + 0.5)( + 0.5) + 1000 = 2 + + 1000.25 -10 + 0.5 Z rozwiÄ…zania równania charakterystycznego otrzymujemy parÄ™ zespolonych, sprzężonych wartoÅ›ci wÅ‚asnych 2 + + 1000.25 = 0 " = -4000, " = Ä… j20 10 1 = -0.5 + j10 10 , 2 = -0.5 - j10 10 KorzystajÄ…c ze wzoru Sylvestera, możemy zapisać: 36 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 Å‚Å‚ e(-0.5+ j10 10)t îÅ‚- j10 10 100 exp(At) = - ïÅ‚ śł j20 10 ïÅ‚ -10 - j10 10ûÅ‚ śł ðÅ‚ Å‚Å‚ e(-0.5-j10 10)t îÅ‚j10 10 100 + = ïÅ‚ śł j20 10 -10 j10 10ûÅ‚ ïÅ‚ śł ðÅ‚ îÅ‚ j10 10t j10 10t j10 10t e + e- e - e- j10 10t Å‚Å‚ ïÅ‚ śł - 10 2 j2 śł = e-0.5t ïÅ‚ ïÅ‚ śł j10 10t j10 10t e - e- j10 10t e + e- j10 10t śł ïÅ‚ 10 ïÅ‚ śł j2 2 ðÅ‚ ûÅ‚ Na podstawie znanych zależnoÅ›ci możemy ostatecznie zapisać: îÅ‚ e-0.5t cos10 10t - 10e-0.5tsin10 10tÅ‚Å‚ exp(At) = ïÅ‚ śł 10e-0.5tsin10 10t e-0.5tcos10 10t ïÅ‚ śł ðÅ‚ ûÅ‚ UwzglÄ™dniajÄ…c wzór xp (t) = exp(At)xp (0) oraz to, że 0 0.05 îÅ‚ Å‚Å‚ îÅ‚ Å‚Å‚ îÅ‚- 0.05 Å‚Å‚ xp (0) = x0 (0) - xu (0) = - = ïÅ‚0śł ïÅ‚ śł ïÅ‚ śł 1 -1 ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ możemy zapisać skÅ‚adowÄ… przejÅ›ciowÄ… îÅ‚ îÅ‚- Å‚Å‚ e-0.5 tcos10 10t - 10e-0.5 t sin10 10tÅ‚Å‚ 0.05 xp (t) = ïÅ‚ śłïÅ‚ śł = -1 10e-0.5 tsin10 10t e-0.5 t cos10 10t ïÅ‚ śłðÅ‚ ûÅ‚ ðÅ‚ ûÅ‚ îÅ‚ Å‚Å‚ e-0.5 t ( 10sin10 10t - 0.05cos10 10t) = ïÅ‚ śł ïÅ‚- e-0.5 t (0.05 10sin10 10t + cos10 10t)ûÅ‚ śł ðÅ‚ KorzystajÄ…c ze wzoru x(t) = xu (t) + xp (t) możemy zapisać: îÅ‚ Å‚Å‚ 0.05 îÅ‚ Å‚Å‚ e-0.5 t ( 10sin10 10t - 0.05cos10 10t) x(t) = xp (t) + xu (t) = + = ïÅ‚ śł ïÅ‚ śł ïÅ‚- e-0.5 t (0.05 10sin10 10t + cos10 10t)ûÅ‚ 1 śł ðÅ‚ ûÅ‚ ðÅ‚ îÅ‚0.05 + e-0.5 t ( 10sin10 10t - 0.05cos10 10t)Å‚Å‚ = ïÅ‚ śł 1 ïÅ‚ - e-0.5 t (0.05 10sin10 10t + cos10 10t) śł ðÅ‚ ûÅ‚ Ostatecznie można zapisać poszukiwane przebiegi jako uC (t) = [0.05 + e-0.5t ( 10sin10 10t - 0.05cos10 10t)] V iL (t) = [1 - e-0.5t (0.05sin10 10t + cos10 10t)] A i zilustrować je na rys. 1.5 dla napiÄ™cia na kondensatorze oraz na rys. 1.6 dla prÄ…du pÅ‚ynÄ…cego przez cewkÄ™. Otrzymane przebiegi majÄ… charakter oscylacyjny. RóżniÄ… siÄ™ od otrzymanych w przypadku b) czterokrotnie mniejszym współczynnikiem tÅ‚umienia oraz pulsacjÄ… blisko trzy razy wiÄ™kszÄ…, wynoszÄ…cÄ… okoÅ‚o 32 rad/s . 37 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 uc [V] 4 3 2 1 0 0 0,2 0,4 0,6 0,8 1 1,2 1,4 1,6 1,8 2 t [s] -1 -2 -3 Rys. 1.5. Przebieg napiÄ™cia na kondensatorze uC(t) iL [A] 2 1,5 1 0,5 0 0 0,2 0,4 0,6 0,8 1 1,2 1,4 1,6 1,8 2 t [s] Rys. 1.6. Przebieg prÄ…du w cewce iL(t) Analiza czÄ™stotliwoÅ›ciowa sygnałów Na wykÅ‚adzie przedstawiono wszystkie treÅ›ci zawarte w: A. SzczepaÅ„ski, M. Trojnar: Obwody elektryczne. Symulacja komputerowa wybranych zagadnieÅ„. Oficyna Wydawnicza Politechniki Rzeszowskiej, Rzeszów 2006, na str. 119-129 oraz na str. 55-71 PozostaÅ‚e treÅ›ci prezentowane podczas wykÅ‚adu: Podstawowe wÅ‚asnoÅ›ci PrzeksztaÅ‚cenia Fouriera 7. Przeskalowanie 1 jÉ x(at) "! X ( ) dla a > 0 a a Przeskalowanie sygnaÅ‚u w osi czasu prowadzi do odwrotnego przeskalowania jego widma, tzn. Å›ciÅ›niÄ™cie (krótszy czas trwania) sygnaÅ‚u (a>1) prowadzi do rozszerzenia jego widma, rozciÄ…gniÄ™cie (dÅ‚uższy czas trwania) zaÅ› sygnaÅ‚u (a<1) do zwężenia widma. 38 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 Rysunek [3]: PrzykÅ‚adowy sygnaÅ‚ x(t) oraz jego wersja przeskalowana x(2t) oraz widma obu sygnałów Transmitancja czÄ™stotliwoÅ›ciowa ukÅ‚adów liniowych TransmitancjÄ™ czÄ™stotliwoÅ›ciowÄ… ukÅ‚adu (czasami mówi siÄ™ o charakterystyce amplitudowo-fazowej lub też o funkcji przenoszenia ukÅ‚adu) można zdefiniować wykorzystujÄ…c odpowiedz impulsowÄ… ukÅ‚adu: +" H ( jÉ) = +"k(t)e jÉ tdt -" Czyli transformata odwrotna charakterystyki amplitudowo-fazowej ukÅ‚adu jest równa odpowiedzi impulsowej ukÅ‚adu (odpowiedzi na impuls Diraca) +" 1 k (t ) = +"H ( jÉ)e jÉ t dÉ 2Ä„ -" Transmitancja czÄ™stotliwoÅ›ciowa determinuje sposób przetwarzania elementarnych sygnałów harmonicznych, czyli jeÅ›li ukÅ‚ad jest pobudzony sygnaÅ‚em harmonicznym o danej pulsacji w i amplitudzie zespolonej Uwe(jw), to amplituda zespolona odpowiedzi na to pobudzenie jest równa U (jÉ) = H (jÉ)U (jÉ) wy we Bardzo czÄ™sto, szczególnie w obliczeniach inżynierskich funkcjÄ™ H(jÉ) wygodniej jest przedstawić w nastÄ™pujÄ…cej postaci: jÕ (É) H (jÉ) = H (É)e gdzie H(É) = H (jÉ) 39 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 H (É) - charakterystyka amplitudowa ukÅ‚adu Õ(É) - charakterystyka fazowa ukÅ‚adu Charakterystyki czÄ™stotliwoÅ›ciowe sygnałów TransformatÄ™ Fouriera F( jÉ) można zapisać w postaci wykÅ‚adniczej F(jÉ) = F(É)exp[jÕ(É)] FunkcjÄ™ F(É) nazywa siÄ™ charakterystykÄ… czÄ™stotliwoÅ›ciowÄ… moduÅ‚u lub widmem amplitudowym sygnaÅ‚u, a funkcjÄ™ Õ(É) charakterystykÄ… czÄ™stotliwoÅ›ciowÄ… fazy lub widmem fazowym sygnaÅ‚u opisanego funkcjÄ… czasu f (t) . SamÄ… funkcjÄ™ F( jÉ) nazywa siÄ™ po prostu charakterystykÄ… czÄ™stotliwoÅ›ciowÄ… lub widmem sygnaÅ‚u. W celu lepszego zbadania widma sygnaÅ‚u sporzÄ…dza siÄ™ wykresy funkcji F(É) i Õ(É) . Wprowadzonych powyżej pojęć widma amplitudowego oraz fazowego nie należy mylić z charakterystykami czÄ™stotliwoÅ›ciowymi transmitancji ukÅ‚adu, które nazywa siÄ™ także wykresami Bodego, gdyż znaczenie ich jest inne. Z charakterystyk czÄ™stotliwoÅ›ciowych transmitancji ukÅ‚adu dowiadujemy siÄ™ jakim zmianom, po przejÅ›ciu przez ukÅ‚ad, podlegajÄ… sygnaÅ‚y harmoniczne o różnych czÄ™stotliwoÅ›ciach podawane na wejÅ›cie ukÅ‚adu; tak wiÄ™c z charakterystyk transmitancji uzyskujemy informacjÄ™ o zmianie amplitudy i różnicy faz miÄ™dzy sygnaÅ‚em na wejÅ›ciu i na wyjÅ›ciu. Z charakterystyk czÄ™stotliwoÅ›ciowych danego sygnaÅ‚u możemy natomiast uzyskać informacje o jego skÅ‚adzie czÄ™stotliwoÅ›ciowym, tzn. o amplitudach i fazach poczÄ…tkowych sygnałów harmonicznych w nim zawartych. PrzykÅ‚ad 1 str. 127-128 w książce [2] 40 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 41 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 PrzykÅ‚ad 2 str. 128-129 w książce [2] 42 Dr inż. Mariusz Trojnar WykÅ‚ady nr 1,2,3,4,5 Obwody i SygnaÅ‚y 2 PrzykÅ‚ad 3 str. 78-90 w książce [2] Wykorzystano nastÄ™pujÄ…ce materiaÅ‚y: 1. G. MasÅ‚owski, WykÅ‚ady z przedmiotu SygnaÅ‚y i Systemy dla Studentów kierunku Informatyka na Wydziale Elektrotechniki i Informatyki Politechniki Rzeszowskiej. 2. A. SzczepaÅ„ski, M. Trojnar, Obwody elektryczne. Komputerowa symulacja wybranych zagadnieÅ„, Oficyna Wydawnicza Politechniki Rzeszowskiej, Rzeszów, 2006. 3. T. ZieliÅ„ski, Cyfrowe przetwarzanie sygnałów, Wydawnictwa Komunikacji i AÄ…cznoÅ›ci, Warszawa, 2005. 4. A. SzczepaÅ„ski, M. Trojnar, Obwody i SygnaÅ‚y, Oficyna Wydawnicza Politechniki Rzeszowskiej, Rzeszów, 2006. 5. J. Osiowski, J. Szabatin, Podstawy teorii obwodów, tom. I, Wydawnictwa Naukowo- Techniczne, Warszawa, 1992. 6. S. Bolkowski, Teoria obwodów elektrycznych, WNT, Warszawa, 1995. 7. K. Rzepka, WykÅ‚ady z przedmiotu Obwody i sygnaÅ‚y dla Studentów WydziaÅ‚u Elektrotechniki i Informatyki Politechniki Rzeszowskiej. 8. J. Bajorek, G. DraÅ‚us, Podstawy elektrotechniki III, Oficyna Wydawnicza Politechniki Rzeszowskiej, Rzeszów, 2005. 9. K. Snopek, WykÅ‚ady z przedmiotu SygnaÅ‚y i Systemy dla Studentów Politechniki Warszawskiej. 10. R. Kurdziel, Podstawy elektrotechniki, Wydawnictwo Naukowo-Techniczne, Warszawa, 1973. 11. Wikipedia; http://www.wikipedia.pl/ 43