02.08.2013 Views

och reglerteknik - Åbo Akademi

och reglerteknik - Åbo Akademi

och reglerteknik - Åbo Akademi

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

˚ABO AKADEMI<br />

TEKNISKA FACULTY OF<br />

FAKULTETEN TECHNOLOGY<br />

Laboratoriet för Process Control<br />

<strong>reglerteknik</strong> Laboratory<br />

INTRODUKTION TILL<br />

SYSTEM- OCH REGLERTEKNIK<br />

HANNU TOIVONEN<br />

Biskopsgatan 8<br />

FIN–20500 ˚Abo Finland


Kapitel 1<br />

Vad är <strong>reglerteknik</strong>?<br />

M˚alsättningen med denna kurs är att ge en informell introduktion till <strong>reglerteknik</strong>en. För att<br />

svara p˚a fr˚agan vad <strong>reglerteknik</strong> är skall vi börja med att betrakta n˚agra exempel. Det är<br />

mycket vanligt att tekniska system kräver n˚agon form av kontinuerlig manipulering för att<br />

fungera p˚a önskat sätt. N˚agra exempel:<br />

• För att h˚alla en bil (eller annat fordon) p˚a vägen bör den hela tiden styras.<br />

• Fartyg <strong>och</strong> flygplan kräver en pilot för att h˚alla kursen under varierande vindförh˚allanden.<br />

• För att h˚alla konstant hastighet hos en bil bör gaspedalens läge justeras för att kompensera<br />

för varierande förh˚allanden, s˚asom upp- <strong>och</strong> nerförsbackar, vindstyrka, vägunderlag<br />

m.m.<br />

• Uppvärmningen i en byggnad bör kontinuerligt justeras för att h˚alla temperaturen<br />

inne konstant trots förändringar i yttertemperaturen, varierande antal människor i byggnaden<br />

samt andra störningar, s˚asom öppnade dörrar eller fönster.<br />

• Temperaturen i en ugn h˚alls vid önskat värde genom att sl˚a p˚a respektive av uppvärmningen<br />

enligt behov.<br />

• En kamera bör fokuseras för att man skall f˚a en skarp bild. I en videokamera med rörliga<br />

objekt bör denna fokusering utföras kontinuerligt för att kompensera för variationer i<br />

objektets avst˚and.<br />

• I en CD-spelare är det viktigt att skivan roterar med en given konstant hastighet.<br />

Dessutom skall läshuvudet positioneras snabbt <strong>och</strong> exakt vid önskat sp˚ar. Allt detta<br />

skall ske trots olika yttre störningar, s˚asom skakning <strong>och</strong> annan rörelse.<br />

1


• I mobiltelefoni bör den mottagna signalstyrkan vid kommunikation mellan basstationen<br />

<strong>och</strong> de olika mobiltelefonerna h˚allas vid en lämplig niv˚a, oavsett avst˚andet till basstationen.<br />

Detta innebär att signalen bör sändas med större effekt ju längre avst˚andet är.<br />

För att uppn˚a konstant signalstyrka justeras effekten därför kontinuerligt.<br />

• För att h˚alla konstant kvalitet vid framställning av produkter av olika slag krävs vanligen<br />

regler- <strong>och</strong> styr˚atgärder. För att t.ex. producera papper av jämn kvalitet, som har<br />

önskad tjocklek, ytvikt, styrka, färg, m.m., bör en pappersmaskins funktion (massaströmmar,<br />

˚angflöden m.m.) ständigt justeras. Detta bör göras eftersom det alltid finns<br />

variationer i r˚avarans egenskaper, som skulle ge upphov till ojämn papperskvalitet utan<br />

kompenserande mot˚atgärder.<br />

Exemplen ovan är typiska praktiska problem inom <strong>reglerteknik</strong>en. Ett gemensamt drag hos<br />

dessa problem är att man utnyttjar mätningar fr˚an ett system (hastighet hos ett fordon, temperaturen,<br />

bildskärpan, signalstyrkan m.m.) för att p˚averka systemet s˚a att önskad funktion<br />

erh˚alls. Vi kan ge följande definition av begreppet reglering (eng. control; fi. säätö):<br />

Reglering innebär att man p˚a basen av mätningar fr˚an ett system, som ger information om<br />

systemets tillst˚and, aktivt p˚averkar systemet att för att f˚a ett önskat beteende.<br />

Reglerteknik (eng. control engineering; fi. säätötekniikka) är det omr˚ade inom tekniken som<br />

sysslar med reglering. Begreppet styrning (eng. control, fi. ohjaus) är nära besläktat med<br />

reglering, men här antar man inte nödvändigtvis att mätningar av systemets tillst˚and utnyttjas<br />

för styrningen. Automation är ett bredare begrepp som även innefattar de olika implementeringsaspekterna<br />

i samband med reglering, s˚asom sensorer, elektronik, <strong>och</strong> datatekniska<br />

problem. Ett annat närst˚aende omr˚ade är kybernetiken (eng. cybernetics), där man studerar<br />

informationens betydelse för funktionen hos olika typers komplicerade system, s˚asom<br />

maskiner, levande organismer eller samhällen.<br />

En del av de i exemplen beskrivna reglerproblemen kan lösas manuellt, genom enbart<br />

mänsklig inverkan. Detta är fallet för bilstyrningsexemplet. Andra problem är däremot sv˚ara,<br />

eller rentav omöjliga, att lösa manuellt. S˚aledes kräver manövreringen av vissa överljudsflygplan<br />

snabbare reaktionshastighet än vad människan klarar av. Operationen av en CD-spelare<br />

eller justeringen av alla signalstyrkor i ett mobiltelefonisystem kan i praktiken inte heller<br />

göras manuellt. Även en pappersmaskins beteende är s˚a komplicerat (<strong>och</strong> snabbt) att det är<br />

sv˚art att upprätth˚alla en jämn produktkvalitet manuellt. Man har därför sedan länge infört<br />

olika typer av automatiserade metoder för reglering. Metoderna g˚ar tillbaka till James Watts<br />

mekaniska hastighetsregulator för ˚angmaskiner p˚a 1700-talet. Före <strong>och</strong> efter tiden för andra<br />

världskriget började elektriska reglersystem introduceras i större skala. I dag förverkligas (implementeras)<br />

reglering nästan uteslutande med hjälp av datorer, vilket möjliggör användning<br />

av beräkningsintensiva <strong>och</strong> mycket effektiva metoder.<br />

2


Kapitel 2<br />

Grundbegrepp<br />

2.1 Introducerande exempel<br />

För att introducera den problematik <strong>och</strong> de fr˚ageställningar som är aktuella inom <strong>reglerteknik</strong>en<br />

skall vi i det följande betrakta ett par enkla exempel p˚a reglerproblem.<br />

Exempel 2.1 - Farth˚allare.<br />

Betrakta automatisk farth˚allning i en bil, vars avsikt är att h˚alla konstant hastighet. P˚a grund<br />

av ständigt varierande förh˚allanden, s˚asom upp- <strong>och</strong> nerförsbackar, varierande vindstyrka,<br />

vägunderlag m.m. bör gaspedalens läge kontinuerligt justeras för att en konstant hastighet<br />

skall kunna upprätth˚allas. För att undersöka hur detta kan ˚astadkommas bör vi undersöka<br />

hur bilens hastighet y beror av de olika ovan beskrivna faktorerna samt hur vi med hjälp<br />

av gaspedalens läge kan p˚averka hastigheten. För detta behöver vi en matematisk modell<br />

som beskriver sambandet mellan de ing˚aende storheterna. En s˚adan modell kan, ˚atminstone<br />

approximativt, bestämmas med hjälp av enkel mekanik.<br />

Situationen kan illustreras enligt figur 2.1. Enligt Newtons tröghetslag gäller<br />

am = F (2.1)<br />

där a = dy/dt är accelerationen, m är bilens massa <strong>och</strong> F är den totala kraften som p˚averkar<br />

bilen i färdriktningen. Bilen p˚averkas av följande krafter:<br />

• Motorns framdrivande kraft Fd. Vi antar för enkelhets skull att denna kraft är direkt<br />

proportionell mot gaspedalens lägesvinkel u,<br />

Fd(t) = ku(t) (2.2)<br />

Vi antar s˚aledes att motorn reagerar ögonblickligen p˚a gaspedalens läge (vilket givetvis<br />

är en approximation).<br />

3


• Gravitationskraftens komponent Fg i vägen plan (jfr figur 2.1),<br />

där ϕ(t) är vägens lutning (ϕ = 0 motsvarar plan väg).<br />

Fg(t) = −mg sin(ϕ(t)) (2.3)<br />

• Luftmotst˚and (Fluft). Denna kraft ökar med stigande hastighet <strong>och</strong> vi kan som en<br />

relativt god approximation anta att den är direkt proportionell mot skillnaden mellan<br />

hastigheten y <strong>och</strong> vindhastigheten vvind i bilens färdriktning,<br />

där b är en luftmotst˚andskoefficient.<br />

Fluft(t) = −b[y(t) − vvind(t)] (2.4)<br />

• Friktionsmotst˚and fr˚an däck, Ff (t). Vi antar att denna kraft, som är riktad mot bilens<br />

färdrikting <strong>och</strong> därför negativ, beror endast av vägunderlaget.<br />

Eftersom a = dy/dt ger ekvation (2.1) med F = Fluft + Fd + Fg + Ff ,<br />

eller<br />

m dy(t)<br />

dt = −b[y(t) − vvind(t)] + ku(t) + Fg(t) + Ff (t) (2.5)<br />

m dy(t)<br />

+ by(t) = ku(t) + d(t) (2.6)<br />

dt<br />

där d(t) = bvvind(t) + Fg(t) + Ff (t). I modellen (2.6) anger y den variabel som skall regleras<br />

(hastigheten, som skall h˚allas konstant), u anger den variabel som manipuleras för att p˚averka<br />

systemets beteende (gaspedalens läge), <strong>och</strong> d(t) anger en yttre störning som p˚averkar den<br />

reglerade variabeln, <strong>och</strong> som i detta exempel best˚ar av vindens inverkan, gravitationskraften<br />

<strong>och</strong> friktionsmotst˚andet. Vi skall ˚aterkomma till problemet hur automatisk farth˚allning kan<br />

˚astadkommas, men före det skall vi betrakta ytterligare ett exempel.<br />

Modellen (2.6) känns igen som en differentialekvation, närmare bestämt en linjär differentialekvation<br />

av första ordningen. De system som är aktuella inom <strong>reglerteknik</strong>en beskrivs<br />

vanligen just av differentialekvationer. För att illustrera saken betraktar vi följande exempel.<br />

Exempel 2.2 - Temperaturreglering.<br />

Betrakta ett temperaturregleringsproblem enligt figur 2.2. Temperaturen T i ett rum skall<br />

h˚allas konstant trots variationer i yttertemperaturen Ty. Värmeförlusterna genom väggarna<br />

är direkt proportionella mot temperaturskillnaden T − Ty, dvs effektförlusterna ges av<br />

Effekt ut = k(T − Ty) (2.7)<br />

Temperaturen kan regleras med hjälp av effekten P i ett värmeelement. Vi antar för enkelhets<br />

skull att luftens omblandning är god, s˚a att temperaturen kan anses densamma i hela rummet.<br />

4


Figur 2.1: Schematisk illustration av farth˚allningsproblemet.<br />

Om P är mindre än värmeförlusten genom väggarna kommer T att minska, <strong>och</strong> om P är<br />

större än värmeförlusten genom väggarna kommer T att öka. Enligt en enkel energibalans<br />

för rummet är ⎡<br />

⎤<br />

⎣<br />

Ändring av<br />

upplagrad energi<br />

per tidsenhet<br />

⎦ = [Effekt in] − [Effekt ut] (2.8)<br />

Den totala mängden luft i rummet är ρV , där ρ är luftens densitet <strong>och</strong> V är rummets volym.<br />

Ändringen av upplagrad energi per tidsenhet är s˚aledes cρV dT<br />

dt , där c är luftens specifika<br />

värmekapacitet. Vi f˚ar allts˚a<br />

cρV dT<br />

dt = P − k(T − Ty) (2.9)<br />

eller<br />

cρV dT<br />

+ kT = P + kTy<br />

(2.10)<br />

dt<br />

Modellen (2.10) kan jämföras med modellen (2.6) i farth˚allningsproblemet. I modellen (2.10)<br />

anger T den variabel som skall regleras (temperaturen), P anger den variabel som manipuleras<br />

för att p˚averka systemets beteende (effekten till värmeelementet), <strong>och</strong> Ty är en yttre störning<br />

som p˚averkar den reglerade variabeln.<br />

Vi har sett att s˚aväl bilen i exempel 2.1 som temperaturen i exempel 2.2 kan beskrivas<br />

med hjälp av en differentialekvation. Detta är typiskt för s.k. dynamiska system. I de enkla<br />

5


Figur 2.2: Schematisk illustration av temperaturregleringsproblemet.<br />

exemplen ovan fick vi differentialekvationer av första ordningen. I allmänhet brukar systemen<br />

emellertid vara mera komplicerade, <strong>och</strong> man f˚ar differentialekvationer av högre ordning.<br />

Eftersom systemen i exemplen ovan kan beskrivas av samma typer av ekvationer, s˚a kan<br />

reglerproblemen i de b˚ada fallen lösas genom att studera det generella problemet att reglera<br />

system som beskrivs av differentialekvationer. Vi behöver allts˚a inte studera farth˚allningsreglering,<br />

temperaturreglering, osv separat, utan det räcker med att helt generellt studera<br />

regleringen av system som beskrivs av en viss typs differentialekvationer. Däremot är givetvis<br />

den praktiska implementeringen (s˚asom mätapparatur m.m.) problemspecifik.<br />

Ur det ovan sagda följer att <strong>reglerteknik</strong> är en generisk metodvetenskap som inte är bunden<br />

till n˚agon speciell del av tekniken. P˚a engelska talar man om ’enabling technology’, för<br />

att betona att det är fr˚agan om en metodik som gör det möjligt att realisera önskade beteenden<br />

<strong>och</strong> funktioner hos tekniska system. I detta avseende har <strong>reglerteknik</strong>en likheter med<br />

ingenjörsmatematiken <strong>och</strong> datatekniken. Reglertekniska problem är viktiga inom alla delar<br />

av tekniken <strong>och</strong> <strong>reglerteknik</strong> är därför ett ämne som studeras inom flera ingenjörsomr˚aden,<br />

s˚asom:<br />

• Elektroteknik.<br />

Reglering av elmotorer, reglering av spänningsaggregat, UPS m.m.<br />

• Robotik.<br />

Reglering av robotar rörelse, automatisk navigation m.m.<br />

6


• Mekanik.<br />

Varvtalsreglering av motorer, aktivfjädring, ABS bromsar m.m.<br />

• Processteknik<br />

Höga kvalitetskrav p˚a framställda produkter, begränsningen av t.ex. r˚amaterialanvändning,<br />

energiförbrukning <strong>och</strong> utsläpp till ett minimum skulle inte kunna uppn˚as utan<br />

l˚angt g˚aende reglering <strong>och</strong> automation av processerna. Reglertekniken utgör ett av de<br />

viktigaste verktygen för att uppn˚a kvalitets- <strong>och</strong> produktivitetkraven inom processindustrin.<br />

• Datateknik.<br />

Reglering förverkligas i praktiken med hjälp av datorer. Regler- <strong>och</strong> styrprogrammen är<br />

realtidssystem <strong>och</strong> dessutom ofta inbyggda system. Reglering <strong>och</strong> automation hör till<br />

de viktigaste tillämpningsomr˚adena av datateknik.<br />

Reglerproblem är, s˚asom vi skall se, ocks˚a av intresse utanför tekniken, t.ex. inom ekonomin<br />

eller medicinen. Mera teoretiska aspekter av reglerproblem studeras dessutom i tillämpad<br />

matematik.<br />

2.2 Signaler <strong>och</strong> system<br />

Vi har i samband med exemplen ovan talat om variabler, s˚asom y(t), u(t) osv, som är funktioner<br />

av tiden. S˚adana variabler kallas signaler, <strong>och</strong> de kan karakteriseras genom att de<br />

inneh˚aller information av olika slag. Signalen y(t) i exempel 2.1 ger t.ex. information om<br />

bilens hastighet som funktion av tiden. Förutom signaler har vi system, som kännetecknas av<br />

den verkan de har p˚a signaler. Bilen i exempel 2.1 är ett system som beskriver hur signalen<br />

y(t) beror av signalerna u(t) <strong>och</strong> d(t).<br />

2.2.1 Blockschema<br />

Man brukar ange sambanden mellan olika signaler <strong>och</strong> system i form av blockscheman. Figur<br />

2.3 visar ett system S med tv˚a insignaler, u(t) <strong>och</strong> d(t), samt en utsignal y(t). Här är u(t) en<br />

styrsignal, som vi kan manipulera för att p˚averka systemet, medan d(t) är en störning, som<br />

vi ej kan manipulera men som p˚averkar systemet. Signalen y(t) är en utsignal fr˚an systemet,<br />

som vi kan mäta.<br />

Exempel 2.3<br />

Bilen i exempel 2.1 är ett system med styrsignalen u(t) (gaspedalens läge) <strong>och</strong> störningen Fg<br />

(gravitationskraften), samt utsignalen v(t) (hastigheten). Själva systemet beskrivs av sambandet<br />

mellan insignalerna <strong>och</strong> utsignalerna, dvs differentialekvationen (2.6).<br />

7


u<br />

✲<br />

S<br />

d<br />

❄<br />

✲ y<br />

Figur 2.3: Ett system S med styrsignalen u, störningen d <strong>och</strong> utsignalen y.<br />

Blockschema är bekväma för att˚ask˚adliggöra strukturen hos sammansatta system. Konstruktionen<br />

av blockschema kan göras med hjälp av elementen i figurerna 2.4–2.6. Figur 2.4 visar<br />

tv˚a seriekopplade system, där utsignalen y1 fr˚an systemet S1 är insignal till systemet S2.<br />

Figur 2.5 visar förgrening av en signal. Observera att signalerna här uppfattas som funktioner<br />

eller informationsflöden, <strong>och</strong> förgreningen skapar s˚aledes tv˚a identiska kopior av signalen<br />

u. Kombination av tv˚a signaler genom addition eller subtraktion symboliseras med en cirkel<br />

enligt figur 2.6. Mera komplexa systemkopplingar kommer att behandlas längre fram.<br />

u<br />

✲ S1<br />

u<br />

y1<br />

✲ S2<br />

✲ y<br />

Figur 2.4: Seriekopplade system.<br />

u<br />

u<br />

✲<br />

✲<br />

Figur 2.5: Förgrening av signal.<br />

8


u1<br />

✲<br />

+<br />

❤<br />

+ ✻<br />

u2<br />

u1 + u2<br />

✲<br />

u1<br />

u1 − u2<br />

✲<br />

+<br />

❤ ✲<br />

− ✻<br />

u2<br />

(a) (b)<br />

Figur 2.6: Summering (a) <strong>och</strong> subtraktion (b) av signaler.<br />

2.2.2 Statiska <strong>och</strong> dynamiska system<br />

Det är viktigt att skilja mellan statiska <strong>och</strong> dynamiska system. Ett statiskt system kännetecknas<br />

av att utsignalen y(t) är beroende av endast insignalens värde u(t) vid samma tidpunkt,<br />

dvs<br />

y(t) = f(u(t)) (2.11)<br />

där f(u) är funktion. Figur 2.7 visar insignalen <strong>och</strong> utsignalen hos ett statiskt system, d˚a det<br />

sker stegvisa förändring i insignalen. Utsignalen följer insignalen ögonblickligen, utan n˚agon<br />

tröghet. I motsats till statiska system har dynamiska system en tröghet som gör att utsignalen<br />

y(t) är beroende av tidigare värden p˚a insignalen u, dvs<br />

y(t) = F (u(τ), τ ≤ t) (2.12)<br />

Dynamiska system system kan vanligen modelleras med hjälp av differentialekvationer, av<br />

vilka vi sett exempel p˚a i exempel 2.1 <strong>och</strong> 2.2.<br />

Problem 2.1<br />

Beskriv en elektrisk krets best˚aende av ett motst˚and med resistansen R som ett system,<br />

där spänningen u(t) över motst˚andet är insignal <strong>och</strong> strömmen i(t) är utsignal. Är systemet<br />

statiskt eller dynamiskt?<br />

9


1.2<br />

1<br />

0.8<br />

0.6<br />

u<br />

0.4<br />

0.2<br />

0<br />

−0.2<br />

0.6<br />

0.5<br />

0.4<br />

0.3<br />

y<br />

0.2<br />

0.1<br />

0<br />

−0.1<br />

−2 0 2 4<br />

tid<br />

6 8 10<br />

Figur 2.7: Responsen hos ett statiskt system.<br />

Problem 2.2<br />

Beskriv en elektrisk krets best˚aende av en spole med induktansen L i serie med ett motst˚and<br />

med resistansen R som ett system, där spänningen u(t) över kretsen är insignal <strong>och</strong> strömmen<br />

i(t) är utsignal. Är systemet statiskt eller dynamiskt?<br />

Exempel 2.4 Enkelt dynamiskt system.<br />

Betrakta ett dynamiskt system<br />

y = Su (2.13)<br />

som beskrivs differentialekvationen<br />

dy(t)<br />

dt<br />

+ ay(t) = bu(t) (2.14)<br />

Figur 2.8 visar insignalen u(t) <strong>och</strong> utsignalen y(t) hos systemet för stegvisa förändringar i<br />

insignalen. Parametervärdena a = 1, b = 1 har använts i figuren. Vi kan i detta skede göra<br />

n˚agra kvalitativa observationer av systemets beteende. Observera att p˚a grund av systemets<br />

tröghet dröjer det en stund innan utsignalen f˚att sitt nya värde efter insignalens förändring.<br />

Om insignalen u är konstant, dvs u(t) = u0 = konstant, kommer utsignalen y att asymp-<br />

totiskt närma sig värdet y = b<br />

a u0 (ty för detta värde är dy/dt = 0 <strong>och</strong> ingen ytterligare<br />

förändring hos y f˚as). Storheten b<br />

a<br />

kallas systemets statiska (eller stationära) förstärkning<br />

. Systemets dynamiska, eller transienta, beteende bestäms ˚a sin sida av parametern a: ju<br />

större positivt värde a har, desto snabbare varierar y(t). Detta kan ses om man definierar<br />

avvikelsen (differensen) ydiff (t) fr˚an det nya statiska värdet efter stegförändringen u0, dvs<br />

10


1.2<br />

1<br />

0.8<br />

0.6<br />

u<br />

0.4<br />

0.2<br />

0<br />

−0.2<br />

1.2<br />

1<br />

0.8<br />

0.6<br />

y<br />

0.4<br />

0.2<br />

0<br />

−0.2<br />

−2 0 2 4<br />

tid<br />

6 8 10<br />

Figur 2.8: Respons hos systemet som beskrivs av ekvation (2.14) för stegformiga förändringar<br />

i insignalen.<br />

ydiff (t) = y(t) − b<br />

a u0. Eftersom dydiff (t)/dt = dy(t)/dt <strong>och</strong> aydiff (t) = ay(t) − bu0 ger<br />

insättning i (2.14) följande differentialekvation för ydiff (t):<br />

dydiff (t)<br />

dt<br />

+ aydiff (t) = 0 (2.15)<br />

Vi ser att för en given avvikelse ydiff (t) gäller att derivatan dydiff (t)/dt = dy(t)/dt är desto<br />

större ju större värde parametern a har. Detta innebär att systemet reagerar desto snabbare<br />

ju större värde parametern har.<br />

Systemets transienta <strong>och</strong> statiska responser kan anges mera explicit i systemekvationen<br />

genom att skriva denna i formen<br />

T dy(t)<br />

dt<br />

+ y(t) = Ku(t) (2.16)<br />

där K = b/a är den statiska förstärkningen <strong>och</strong> T = 1/a kallas systemets tidskonstant, <strong>och</strong><br />

är direkt proportionell mot den tid som det tar för systemet att reagera för en förändring i<br />

insignalen.<br />

Trögheten hos dynamiska system beror vanligen p˚a olika typers energiupplagring eller p˚a<br />

transportfördröjningar. I farth˚allningsexemplet är det bilens upplagrade rörelseenergi som<br />

ger upphov till trögheten, <strong>och</strong> i temperaturregleringexemplet är det den i luften lagrade<br />

värmeenergin. I Problem 2.2 lagras energi i spolens elektromagnetiska fält. Trögheten gör<br />

11


att insignalen u(t) till ett dynamiskt system p˚averkar det framtida förloppet hos systemets<br />

utsignal y(t). För att kunna reglera <strong>och</strong> styra dynamiska system är det därför viktigt att ha<br />

en modell som beskriver det framtida beteendet hos systemet.<br />

2.3 Systemtekniska ämnen<br />

Signaler <strong>och</strong> system är viktiga inte endast inom <strong>reglerteknik</strong>en, utan mera allmänt inom<br />

s.k. systemvetenskaper. Speciellt viktiga är dessa begrepp inom signalbehandling. Medan man<br />

inom <strong>reglerteknik</strong>en använder information i uppmätta signaler för att kunna p˚averka ett<br />

system s˚a att det beter sig p˚a önskat sätt, är man inom signalbehandlingen intresserad av att<br />

manipulera själva signalerna, t.ex. för att filtrera bort brus eller komprimering av data som<br />

en signal inneh˚aller.<br />

Exempel 2.5 - Effektreglering i mobiltelefoni.<br />

Den mottagna signalstyrkorna fr˚an de olika mobiltelefonerna till basstationen h˚alls konstanta<br />

oavsett avst˚andet till basstationen, jfr figur 2.9. Detta sker genom att basstationen sänder<br />

information om den mottagna signalstyrkan till mobiltelefonen, som ändrar effekten p˚a den<br />

utsända signalen enligt behov. Detta är ett reglerproblem: mätningar fr˚an systemet (styrkan<br />

hos mottagen signal) används för att manipulera systemet (den utsända signalens effekt).<br />

Figur 2.9: Mobiltelefonerna justerar den utsända effekten p˚a basen av den mottagna signalstyrkan<br />

vid basstationen.<br />

Exempel 2.6 - Signalfiltrering i mobiltelefoni.<br />

Signalen mellan basstation <strong>och</strong> mobiltelefon p˚averkas p˚a grund av flervägsutbredning, jfr figur<br />

12


2.10. Detta kan beskrivas med hjälp av ett system S som p˚averkar den utsända signalen s:<br />

y = Ss (2.17)<br />

där y är den vid telefonen mottagna signalen. Den ursprungliga signalen s kan rekonstrueras<br />

fr˚an den mottagna signalen y genom att bestämma ett filter F, s˚a att<br />

ˆs = Fy ≈ s (2.18)<br />

För att kunna bestämma filtret F bör systemet S vara känt. Detta är ett signalbehandlingsproblem:<br />

signalen fr˚an systemet filtreras för att ta fram den ursprungliga signalen, men<br />

själva systemet manipuleras inte.<br />

Figur 2.10: Den mottagna signalen y vid mobiltelefonen är förvrängd p˚a grund av flervägsutbredning.<br />

Den ursprungliga signalen s kan rekonstrueras genom filtering av y.<br />

2.4 Begreppet ˚aterkoppling<br />

Vi skall nu titta n˚agot närmare p˚a de principer som används för att lösa reglerproblem<br />

av den typ som diskuterats ovan. Vi betraktar farth˚allningsproblemet i exempel 2.1. Bilens<br />

hastighet kunde (˚atminstone approximativt) beskrivas med differentialekvatioen (2.6), som<br />

vi här ˚aterger i n˚agot modifierad form:<br />

dy(t)<br />

dt + a1y(t) = b1u(t) + c1d(t) (2.19)<br />

13


❄<br />

+<br />

✲ ❡ e ✲ u<br />

y<br />

Gc<br />

✲ Gp<br />

✲<br />

− ✻<br />

d<br />

Figur 2.11: Ett ˚aterkopplat system.<br />

där a1 = b/m, b1 = k/m <strong>och</strong> c1 = 1/m. Anta nu att den önskade hastigheten, eller<br />

hastighetens börvärde, är y(t) = r (t.ex. r = 80 km/h). Enligt modellen (2.19) är y(t) = r<br />

om vi väljer insignalen<br />

u(t) = 1<br />

[a1r − c1d(t)] (2.20)<br />

b1<br />

ty d˚a är dy/dt = 0 för det önskade värdet y(t) = r. Insignalvalet (2.20) ger s˚aledes det önskade<br />

statiska värdet hos y. Enligt (2.20) <strong>och</strong> (2.19) ges den transienta responsen, d˚a hastigheten<br />

y(t) är olikt r, av<br />

dy(t)<br />

dt + a1y(t) = a1r (2.21)<br />

Detta är ett system av första ordningen, med en dynamik av den typ som vi hade i exempel<br />

2.4. Om y(t) är olikt r i början, tar det därför en tid innan hastigheten närmat sig r.<br />

Trots att ovan beskrivna förfarande verkar att fungera i teorin har den emellertid n˚agra<br />

uppenbara nackdelar:<br />

• Metoden kräver att störningen d(t) är exakt känd. M.a.o. skall gravitationskraften Fg(t),<br />

luftmotst˚andet bvvind(t) samt friktionsmotst˚andet Ff vara exakt kända för att proceduren<br />

skall kunna tillämpas. Vi kunde först˚as i princip mäta vägplanets lutning ϕ <strong>och</strong><br />

bilens massa m, <strong>och</strong> därmed bestämma Fg(t) = mg sin(ϕ(t)). Även luftmotst˚andet<br />

kan uppskattas genom att noggrant bestämma luftmotst˚andskoefficienten b <strong>och</strong> vindstyrkan<br />

vvind(t). Detta förfarande har den uppenbara nackdelen att den fordrar noggrann<br />

kännedom om alla störningar (Fg(t), vvind(t), Ff osv) <strong>och</strong> andra faktorer (s˚asom<br />

b) som p˚averkar den reglerade utsignalen. I praktiken är det emellertid i allmänhet helt<br />

orealistiskt att ha fullständig kunskap om alla störningar som p˚averkar ett system.<br />

• Även om kännedom om alla faktorer som p˚averkar y(t) funnes, har vi inte p˚averkat<br />

snabbheten hos dynamiken i ekvation (2.21). Svarets snabbhet bestäms av a1 = b/m.<br />

Om m är stort kan förändringen i hastigheten ske mycket l˚angsamt!<br />

14


d<br />

u<br />

y<br />

−1 0 1 2<br />

tid<br />

3 4 5<br />

Figur 2.12: Responsen hos systemet (2.19) för stegformad störning d˚a ingen reglering används.<br />

Begränsningen med det ovan beskrivna förfarandet är att man försöker bestämma styrsignalen<br />

u(t) som en funktion av endast s˚adana variabler som p˚averkar den reglerade signalen y(t),<br />

men utnyttjar ej mätningar av själva y! Genom att ocks˚a utnyttja mätning av y(t) f˚ar vi<br />

information om avvikelsen fr˚an det önskade värdet, r − y(t). Genom att utnyttja denna<br />

mätning kan otillräcklig kunskap om störningarna <strong>och</strong> systemet kompenseras. Dessutom är det<br />

möjligt att göra systemets respons snabbare. Principen att bestämma styrsignalen u(t) som<br />

funktion av den reglerade utsignalen kallas ˚aterkoppling (eng. feedback; fi. takaisinkytkentä).<br />

Situationen illustreras i figur 2.11, som visar hur utsignalen y(t) fr˚an systemet Gp ˚aterkopplas<br />

för att bestämma styrsignalen. Kretsen i figur 2.11 kallas ˚aterkopplad krets eller sluten krets.<br />

Blocket Gc kallas regulator. Signalen r är börvärdet eller ledvärdet för utsignalen y, <strong>och</strong> signalen<br />

e(t) = r − y(t) (2.22)<br />

är regleravvikelsen eller reglerfelet.<br />

Det visar sig att ˚aterkopplingsprincipen i all sin enkelhet är en mycket kraftfull metod för<br />

att p˚averka dynamiska systems beteende. Vi skall illustrera effekten av ˚aterkoppling genom<br />

att undersöka vad som händer om vi använder den enkla reglerlagen<br />

u(t) = Kp [r − y(t)] = Kpe(t) (2.23)<br />

dvs regulatorn är en konstant, Gc = Kp. Reglerlagen (2.23) kallas proportionell regulator eller<br />

helt enkelt P-regulator, eftersom styrsignalen u(t) är direkt proportionell mot regleravvikelsen<br />

15


d<br />

u<br />

y<br />

−1 0 1 2<br />

tid<br />

3 4 5<br />

Figur 2.13: Responsen för stegformad störning d˚a systemet (2.19) regleras med en P-regulator<br />

(2.23).<br />

r − y(t). Substitution av (2.23) i systemekvationen (2.19) ger för den slutna kretsen differentialekvationen<br />

dy(t)<br />

dt + (a1 + b1Kp)y(t) = b1Kpr + c1d(t) (2.24)<br />

Vi ser att för en konstant stegstörning d(t) = d, konvergerar utsignalen mot det statiska<br />

svaret (jfr exempel 2.4)<br />

y(t) → b1Kp<br />

r +<br />

a1 + b1Kp<br />

c1<br />

a1 + b1Kp<br />

d d˚a t → ∞ (2.25)<br />

För stora värdet p˚a regulatorparametern Kp gäller enligt (2.25) y(t) ≈ r d˚a t → ∞. Förutom<br />

det statiska svaret har även den transienta responsen p˚averkats. Eftersom systemets transienta<br />

respons beror av värdet hos storheten a1 + b1Kp, kan den i princip göras godtyckligt<br />

snabb genom att göra Kp stor.<br />

Vi har s˚aledes med den enkla proportionella reglerlagen (2.23) kunnat minska störningarnas<br />

inverkan samt gjort systemets transienta respons snabbare utan att mäta störningarna<br />

d(t). Resultatet beror inte heller av n˚agon kunskap om systemparametern a1. Fullständig<br />

eliminering av en stegstörnings inverkan kräver emellertid ett oändligt stort värde p˚a Kp.<br />

Det visar sig att Kp i praktiken inte kan göras hur stort som helst, eftersom t.ex. sm˚a tidsfördröjningar,<br />

som alltid finns i verkliga system, d˚a gör att den slutna kretsen blir instabil.<br />

(I farth˚allningsexemplet har vi t.ex. försummat motorns dynamik.) De begränsningar som<br />

stabilitetskravet innebär kommer att diskuteras senare.<br />

16


d<br />

u<br />

y<br />

−1 0 1 2<br />

tid<br />

3 4 5<br />

Figur 2.14: Responsen för stegformad störning d˚a systemet (2.19) regleras med en I-regulator<br />

(2.26).<br />

För att undvika statiska regleravvikelser efter en stegstörning behöver vi i stället för Pregulatorn<br />

en reglerlag som kontinuerligt justerar u(t) s˚a länge y(t) är olikt referensvärdet<br />

r. M.a.o. skall u(t) ökas s˚a länge y(t) < r <strong>och</strong> minskas s˚a länge y(t) > r, tills ett s˚adant<br />

värde för u(t) uppn˚atts för vilket y(t) = r. Detta är precis vad en människa skulle göra för<br />

att reglera y(t) till börvärdet. Matematiskt kan en s˚adan reglerlag beskrivas med hjälp av en<br />

integrator i formen<br />

u(t) = Ki<br />

t<br />

τ=0<br />

[r − y(τ)] dτ = Ki<br />

t<br />

τ=0<br />

e(τ)dτ (2.26)<br />

Reglerlagen (2.26) kallas integrerande regulator eller helt enkelt I-regulator. Principen hos en<br />

I-regulator är den, att s˚a länge r − y(t) > 0 ökar integralens värde, varvid u(t) blir större <strong>och</strong><br />

f˚ar värdet hos y(t) att växa. Detta h˚aller p˚a tills y(t) = r, dvs regleravvikelsen har körts till<br />

noll, varefter u(t) har ett konstant värde. En liknande funktion gäller om r − y(t) < 0. P˚a<br />

detta sätt elimineras inverkan av stegstörningar s˚a att ingen statisk avvikelse f˚as. Observera<br />

att detta har uppn˚atts utan att känna till stegstörningens storlek!<br />

En svaghet hos I-regulatorn är att eftersom regleravvikelsen r − y(t) integreras blir Iregulatorns<br />

respons mot reglerfel l˚angsam. I-regulatorns förm˚aga att eliminera statiska regleravvikelser<br />

<strong>och</strong> P-regulatorns snabba respons kan kombineras i PI-regulatorn<br />

u(t) = Kpe(t) + Ki<br />

t<br />

τ=0<br />

e(τ)dτ (2.27)<br />

Observera att PI-regulatorn reagerar först d˚a ett reglerfel orsakad av en störning redan<br />

förekommer, dvs e(t) = r − y(t) = 0. För att förekomma en regleravvikelse redan innan den<br />

17


d<br />

u<br />

y<br />

−1 0 1 2<br />

tid<br />

3 4 5<br />

Figur 2.15: Responsen för stegformad störning d˚a systemet (2.19) regleras med en PI-regulator<br />

(2.27).<br />

hunnit uppst˚a inför man därför vanligen ytterligare en deriverande term i regulatorn. P˚a<br />

detta sätt f˚as en PID-regulator, som ges av<br />

u(t) = Kpe(t) + Ki<br />

t<br />

τ=0<br />

de(t)<br />

e(τ)dτ + Kd<br />

dt<br />

(2.28)<br />

Avsikten med den deriverande termen är att f˚a styrsignalen u(t) att reagera p˚a förändringsriktningen<br />

hos y(t), för att p˚a detta sätt kunna motverka en regleravvikelse redan innan den<br />

hunnit uppst˚a. P˚a detta sätt kan regulatorns respons göras snabbare. I praktiken kan den<br />

deriverande verkan Kd emellertid inte göras hur stor som helst, eftersom regulatorn d˚a blir<br />

alltför känslig mot det brus som man alltid har i praktiken.<br />

De olika regulatortypernas egenskaper illustreras i figurerna 2.12–2.16 som visar styrsignalen<br />

u(t) <strong>och</strong> utsignalen y(t) för systemet (2.19) efter en stegformad störning.<br />

PID-regulatorn är den överlägset vanligaste standardregulatorn för enkla reglerproblem<br />

i industrin <strong>och</strong> andra tekniska system, <strong>och</strong> den finns implementerad i alla processdatorsystem.<br />

För mera komplicerade reglerproblem, där systemet som skall regleras har ett mera<br />

komplicerat dynamiskt beteende, behövs dock mera invecklade reglerstrategier.<br />

En ännu enklare regulatorn är p˚a-av-regulatorn eller ON-OFF-regulatorn. I en p˚a-avregulator<br />

har styrsignalen u endast tv˚a värden: umax <strong>och</strong> umin, <strong>och</strong> värdet bestäms beroende<br />

p˚a om systemets utsignal är större eller mindre än ledvärdet, t.ex.<br />

<br />

umax om y(t) < r<br />

u(t) =<br />

(2.29)<br />

umin om y(t) > r<br />

18


d<br />

u<br />

y<br />

−1 0 1 2<br />

tid<br />

3 4 5<br />

Figur 2.16: Responsen för stegformad störning d˚a systemet (2.19) regleras med en PIDregulator<br />

(2.28).<br />

P˚a-av-regulatorer används i enkla system som ofta har en l˚angsam dynamisk respons, <strong>och</strong> för<br />

vilka det räcker med att y(t) är i närheten av ledvärdet r. Typiska tillämpningar av p˚a-av<br />

reglering är olika typer av temperaturregleringsproblem (jfr exempel 2.2), s˚asom reglering av<br />

temperaturen i hus eller bilar, kylsk˚ap, strykjärn, osv. Här sl˚as uppvärmningen p˚a eller av<br />

beroende p˚a temperaturen. Det diskontinuerliga funktionssättet hos regulatorn f˚ar utsignalen<br />

(temperaturen) att svänga kring ledvärdet, men detta kan mycket väl accepteras i de flesta<br />

temperaturregleringsproblem.<br />

Framkoppling<br />

Trots att ˚aterkoppling av den reglerade signalen är viktig för att kunna kompensera för<br />

inverkan av delvis okända störningar kan information om störningar naturligtvis utnyttjas<br />

för att uppn˚a bättre resultat. I farth˚allningsexemplet kunde t.ex. information om en störning i<br />

form av en kommande uppförsbacke användas för att accelerera redan innan en minskning av<br />

hastigheten f˚as. Denna princip att använda (partiell) kunskap om störningar för att bestämma<br />

den manipulerbara styrsignalen kallas framkoppling (eng. feedforward; fi. myötäkytkentä).<br />

Figur 2.17 visar ett blockschema av ett reglersystem med framkoppling.<br />

2.5 Digital reglering<br />

Regulatorer implementeras idag nästan uteslutande digitalt med hjälp av datorer. Detta<br />

innebär att den ovan beskrivna PID-regulatorn kan förverkligas endast approximativt.<br />

19


Gf ✛<br />

❄<br />

+<br />

✲ ❡ e ✲<br />

❄❡ u<br />

y<br />

Gc<br />

✲ ✲ Gp<br />

✲<br />

− ✻<br />

Figur 2.17: Ett system med ˚aterkoppling <strong>och</strong> framkoppling.<br />

Figur 2.18 visar ett typiskt digitalt reglersystem. Den kontinuerliga utsignalen y(t) diskretiseras<br />

med hjälp av en A/D-omvandlare, som genererar en diskret sekvens,<br />

d<br />

yk = y(kh), k = 0, 1, 2, . . . (2.30)<br />

Här är h samplingstiden (eng. sampling time; fi. näytteenottoväli). Sekvensen yk, k = 0, 1, 2, . . .<br />

processeras sedan digitalt för att generera en diskret styrsignalsekvens uk, k = 0, 1, 2, . . . till<br />

systemet. För system med en kontinuerlig dynamik bör den diskreta styrsignalen omvandlas<br />

till en kontinuerlig styrsignal u(t) till systemet. Detta utförs med en D/A-omvandlare, som<br />

genererar en styckevis konstant styrsignal enligt<br />

u H(t) = uk, kh ≤ t < kh + h (2.31)<br />

Efter D/A-omvandlaren brukar man ännu ha ett filter H för att utjämna diskontinuiteterna<br />

hos uH(t).<br />

Den digitala regulatorn Gd som beräknar styrsignalen uk ur mätsignalen yk kan bestämmas<br />

p˚a olika sätt. Om samplingstiden är kort är en vanlig metod att helt enkelt diskretisera<br />

den kontinuerliga systemmodellen eller reglerlagen genom att approximera derivator med<br />

finita differenser <strong>och</strong> integraler med summor.<br />

Om vi inför derivata-approximationen<br />

dy(t)<br />

dt<br />

1<br />

≈ [y(t + h) − y(t)] (2.32)<br />

h<br />

f˚as t.ex. för systemet som beskrivs av differentialekvationen (2.19) den diskreta approximationen<br />

1<br />

h [y(kh + h) − y(kh)] + a1y(kh) = b1u(kh) + c1d(kh) (2.33)<br />

eller<br />

y(kh + h) + (ha1 − 1)y(kh) = hb1u(kh) + hc1d(kh) (2.34)<br />

20


k<br />

✲<br />

+<br />

✻<br />

−<br />

ek<br />

❣ ✲<br />

Gd<br />

uk uH(t) ✲ D/A ✲<br />

H<br />

u(t)<br />

✲<br />

Gp<br />

Figur 2.18: Digitalt reglersystem.<br />

y(t)<br />

✲<br />

Införs beteckningarna yk = y(kh), uk = u(kh), dk = d(kh) <strong>och</strong> ad = ha1−1, bd = hb1, cd = hc1<br />

f˚as<br />

(2.35)<br />

yk+1 + adyk = bduk + cddk<br />

vilket är en differensekvation av första ordningen.<br />

De kontinuerliga standardregulatorerna har sina diskreta motsvarigheter, s˚a att en diskret<br />

(digital) P-regulator ges av<br />

uk = Kpek, (2.36)<br />

en digital I-regulator ges av<br />

<strong>och</strong> en digital PID-regualtor ges av<br />

k<br />

uk = Ki en<br />

n=1<br />

k<br />

uk = Kpek + Ki<br />

n=1<br />

A/D<br />

✲<br />

yk<br />

(2.37)<br />

en + Kd[ek − ek−1] (2.38)<br />

Den digitala PID-regulatorn brukar implementeras i en s.k. differensform,<br />

uk = uk−1 + Kp[ek − ek−1] + Kiek + Kd[ek − 2ek−1 + ek−2] (2.39)<br />

där summan i uttrycket (2.38) eliminerats genom att subtrahera uk−1 fr˚an uk.<br />

Exempel 2.7 - Digital effektreglering i mobiltelefoni.<br />

Ett enkelt exempel p˚a hur en digital I-regulator fungerar är effektregleringen i mobiltelefoni,<br />

jfr exempel 2.5 <strong>och</strong> figur 2.9. L˚at uk beteckna effekten hos den utsända signalen under tidsintervalllet<br />

k, <strong>och</strong> l˚at yk vara styrkan hos den mottagna signalen under samma tidsintervall.<br />

Man strävar till att signalstyrkan vid mottagaren har ett givet värde r. Eftersom signalens<br />

dämpning varierar som funktion av avst˚and <strong>och</strong> terräng bör effekten hela tiden ändras s˚a att<br />

önskad mottagen signalstyrka kan upprätth˚allas. Enda sättet att p˚a ett tillfredsställande sätt<br />

21


uppn˚a detta i praktiken är via ˚aterkoppling, s˚a att uk ökas s˚a länge yk < r <strong>och</strong> minskas s˚a<br />

länge yk > r. Detta motsvarar digital I-reglering enligt<br />

k<br />

uk = Ki<br />

n=1<br />

där Ki > 0. Reglerlagen kan även skrivas i formen<br />

[r − yn]<br />

uk = uk−1 + [r − yk]<br />

vilket explicit visar att reglerlagen har de önskade egenskaperna, dvs uk > uk−1 om r−yk > 0<br />

<strong>och</strong> uk < uk−1 om r − yk < 0.<br />

2.6 Exempel p˚a ˚aterkoppling<br />

˚Aterkoppling är en mycket kraftfull metod för att p˚averka systems beteende även i s˚adana<br />

fall d˚a systemets dynamik eller störningarna är endast ofullständigt kända. S˚asom vi sett kan<br />

t.ex. regleravvikelsen pga en okänd stegstörning fullständigt elimineras med en integrerande<br />

regulator. Det är därför naturligt att ˚aterkopplingsprincipen utnyttjas ocks˚a i flera sammanhang<br />

utanför tekniken. I nästan alla processer som kan beskrivas med de mycket generella<br />

begreppen signaler <strong>och</strong> system finner man ˚aterkoppling av signaler. Viktiga exempel finns<br />

bl.a. inom ekonomiska <strong>och</strong> biologiska system.<br />

Exempel 2.8 - ˚Aterkoppling i vardagen.<br />

• När vi g˚ar eller cyklar h˚aller vi oss upprätta genom att använda den information som<br />

v˚ara sinnen ger oss. Utan denna ˚aterkoppling skulle vi mycket snabbt tappa balansen<br />

<strong>och</strong> falla.<br />

• När vi griper ett förem˚al ˚aterkopplar vi visuell information av förem˚alets läge i förh˚allande<br />

till handen. Robotar kan p˚a liknande sätt styras av visuell information fr˚an en<br />

digitalkamera för att gripa förem˚al.<br />

• Vi justerar temperaturen i en dusch genom att känna p˚a temperaturen <strong>och</strong> manipulera<br />

varm- eller kallvattenkranen tills rätt temperatur f˚as.<br />

• En föreläsare f˚ar ˚aterkoppling eller feedback (fi. ’palaute’) fr˚an ˚ahörarna via fr˚agor,<br />

kritik, diskussioner o.dyl. Om han beaktar denna feedback har vi ett sluten krets som<br />

kan konvergera till ett, förhoppningsvis bättre, tillst˚and. Utan s˚adan ˚aterkoppling vet<br />

föreläsaren ej vad ˚ahörarna f˚ar ut av föreläsningen.<br />

• En studerande är missnöjd med sitt tentresultat <strong>och</strong> förbereder sig bättre till följande<br />

tillfälle. Den studerande använder tentresultatet för ˚aterkoppling s˚a att önskat resultat<br />

uppn˚as.<br />

22


• En hund jagar en katt. Detta kan tolkas som ˚aterkoppling där kattens position är<br />

börvärdet för hunden. P˚a liknande sätt söker en m˚alsökande missil sig till ett rörligt<br />

m˚al.<br />

• ’Du skall behandla andra s˚a som du själv vill bli behandlad’.<br />

Inom ekonomin finns flera fenomen som kan beskrivas i form av dynamiska system. Detta<br />

gäller t.ex. nationalekonomin eller aktiemarknaden, vars respons till olika variabler uppvisar<br />

en tröghet. I de ekonomiska systemen finns ocks˚a flera exempel p˚a ˚aterkopplingsmekanismer.<br />

Exempel 2.9 - ˚Aterkoppling i nationalekonomin.<br />

Den ekonomiska aktiviteten kan styras p˚a olika sätt, t.ex. med hjälp av skatteniv˚an. Med en<br />

god politik kan ekonomisk aktivitet stimuleras under en l˚agkonjunktur genom skattesänkningar,<br />

medan överhettning inom ekonomin kan reduceras genom skattehöjningar under en högkonjuktur.<br />

Ekonomin är emellertid ett dynamiskt system, som reagerar p˚a förändringar med en<br />

viss tröghet. Det är därför viktigt att göra förändringar av olika slag vid rätta tidpunkter. Politiker<br />

som fattar beslut om skatter är emellertid ocks˚a dynamiska system i sin beslutsfattning.<br />

Därför är det vanligt att förändringarna görs när det redan är för sent: skattelättnader introduceras<br />

när konjunkturerna redan svängt till högkonjunktur, <strong>och</strong> höjda skatter införs först<br />

under l˚agkonjunkturen. I denna situation gör ˚aterkopplingen fr˚an eknomin till beslutsfattning<br />

konjunktursvängningarna snarare kraftigare än mindre.<br />

Inom biologin finns flera viktiga processer <strong>och</strong> fenomen, som kan först˚as <strong>och</strong> analyseras<br />

inom ramen för teorin om dynamiska system <strong>och</strong> ˚aterkoppling. Detta beror p˚a att livet i<br />

hög grad g˚ar ut p˚a att reagera p˚a omgivningen, <strong>och</strong> är därför p˚a sätt <strong>och</strong> vis en stor reglerprocess.<br />

Viktiga ˚aterkopplingsprocesser inom biologin finns i synnerhet inom fysiologin,<br />

i cellernas funktion samt inom ekosystemen. En viktig skillnad mellan tekniska <strong>och</strong> biologiska<br />

system är, att medan man i tekniska system konstruerar en regulator för att f˚a en<br />

önskad funktion, s˚a har i biologiska system hela den ˚aterkopplade kretsen uppst˚att samtidigt,<br />

genom olika evolutionsmekanismer. Det därför inte alltid självklart vad som är ’regulator’ <strong>och</strong><br />

vad som är ’det reglerade systemet’, utan ˚aterkopplingsmekanismens funktion är snarare att<br />

göra hela systemet mindre känsligt för externa störningar. Teorin för dynamiska system <strong>och</strong><br />

˚aterkoppling kan i dessa fall bidra till att skapa en först˚aelse för systemets funktion. Följande<br />

exempel ger en uppfattning om vilken roll dynamiska system <strong>och</strong> ˚aterkoppling spelar inom<br />

biologin.<br />

Exempel 2.10 - Kroppens temperaturreglering.<br />

Kroppstemperaturen h˚alls mycket konstant oavsett ansenliga variationer i yttre temperatur<br />

<strong>och</strong> den genererade energin p.g.a fysisk aktivitet. Detta ˚astadkoms genom ˚aterkoppling s˚a<br />

att blodets temperatur p˚averkar kroppens uppvärmnings- <strong>och</strong> avkylningsmekanismer, s˚asom<br />

svettning, mekanisk skakning <strong>och</strong> blodcirkulationen i de yttre blodkärlen. Utan ˚aterkoppling<br />

fr˚an temperaturen skulle konstant kropppstemperatur inte kunna upprätth˚allas.<br />

23


Exempel 2.11 - Reglering av blodets glukoshalt.<br />

Det faktum att blodets glukoshalt h˚alls mycket konstant har kunnat förklaras genom en fysiologisk<br />

˚aterkopplingsmekanism. Halten glukos i blodet är s˚a gott som konstant, ca 5 mmol/l,<br />

oavsett stora variationer i tillförd glukos, s˚asom när man äter sötsaker. Detta möjliggörs<br />

genom att en glukoshalt som överstiger 5 mmol/l leder till generering av insulin, som avlägsnar<br />

glukos fr˚an blodet, medan en glukoshalt som understiger 5 mmol/l leder till generering av ett<br />

annat hormon, glukagon, som i sin tur stimulerar avsöndringen av glukos fr˚an muskler <strong>och</strong><br />

övriga organ till blodet. Detta är ett exempel p˚a ˚aterkopplad reglering: den variabel som skall<br />

h˚allas vid ett konstant värde (glukoshalten) p˚averkar variabler (insulin- <strong>och</strong> glukagonhalterna)<br />

som i sin tur reglerar glukoshalten. En detaljerad analys av kroppens glukosreglering<br />

har visat att regleringen fungerar som en I-regulator, <strong>och</strong> kan kompensera mot konstanta<br />

belastningsstörningar som p˚averkar glukoshalten. Defekter i denna ˚aterkopplingsmekanism<br />

har allvarliga följder <strong>och</strong> leder till olika typer av diabetes.<br />

Exempel 2.12 - Reglering av proteinsyntes i cellerna.<br />

Proteinsyntesen i cellerna styrs av cellkärnornas DNA. Syntesen av ett givet protein äger rum<br />

d˚a aktuell gen är aktiv, varvid motsvarande nukleidsyresekvens i kärnans DNA översätts<br />

till proteinets aminosyresekvens. Ett protein framställs i korrekta mängder <strong>och</strong> vid rätta<br />

tidpunkter genom˚aterkoppling: aktiviteten hos generna regleras av proteinkoncentrationerna.<br />

För att först˚a cellernas funktion bör man reda ut ˚aterkopplingsmekanismerna. De dynamiska<br />

processerna i en cell är emellertid synnerligen komplicerade, <strong>och</strong> ett av de viktigaste m˚alen<br />

inom cellbiologin för tillfället är att kartlägga dessa processer. Denna strävan har gett upphov<br />

till ett nytt tvärvetenskapligt omr˚ade, den s.k. systembiologin.<br />

Exempel 2.13 - Reglermekanismer i biosfären.<br />

Biosfären inneh˚aller en del intressanta ˚aterkopplingsmekanismer som är viktiga för att först˚a<br />

globala egenskaper s˚asom jordens klimat eller atmosfärens <strong>och</strong> oceanernas sammansättningar.<br />

Betrakta t.ex. följande fakta:<br />

• Atmosfärens sammansättning har h˚allits praktiskt taget konstant den tid det funnits liv<br />

p˚a land. Syrehalten är t.ex. vid det bekväma värdet 21 vol-%. Om värdet understeg ca 15<br />

vol-% skulle en tändsticka slockna <strong>och</strong> vi skulle kvävas pga syrebrist, <strong>och</strong> om det översteg<br />

25 vol-% skulle vegetationen självantändas. Den konstanta halten är anmärkningsvärd,<br />

eftersom syre reagerar med andra ämnen, <strong>och</strong> utan ett tillskott som exakt kompenserar<br />

den reagerade syremängden skulle halten ej kunna förbli konstant.<br />

• Man vet att oceanernas salthalt har h˚allit sig mycket konstant under den tid det funnits<br />

liv i havet. Salthalten är vid en niv˚a som är lämplig för oceanernas djur- <strong>och</strong> växtliv. Den<br />

konstanta salthalten är anmärkningsvärd eftersom floder ständigt tillför nya mineral till<br />

oceanerna. I kombination med avdunstningen skulle detta leda till en gradvis ökning<br />

av salthalten (jämför Döda Havet). Det finns allts˚a n˚agon mekanism som avskaffar salt<br />

fr˚an oceanerna som kompenserar för salttillförseln.<br />

24


• Jordens medeltemperatur har h˚allit sig mycket konstant de senaste ca 3.5 miljarder<br />

˚aren, med en maximal avvikelse p˚a ca ±5 ◦ C. Under samma period har solens str˚alning<br />

ökat med ca 30%. Man kan räkna ut att denna ökning i solens str˚alning skulle ge upphov<br />

till en betydande höjning av jordens medeltemperatur (flera tiotals grader) utan n˚agon<br />

mekanism som motverkar denna höjning.<br />

Orsakerna till ovan beskrivan förh˚allanden är olika typer av biologisk aktivitet, som inför<br />

˚aterkopplingsmekanismer. Ännu för 30 ˚ar sedan var den allmänt accepterade uppfattningen<br />

den, att biologiska organismer endast passivt adapterar sig till den omgivande miljön, <strong>och</strong><br />

att det är just därför som miljön i varje enskilt fall är lämplig för djur- <strong>och</strong> växtlivet där.<br />

Man har först relativt nyligen, under de senaste 20–30 ˚aren, insett att situationen är mera<br />

komplicerad än s˚a, <strong>och</strong> att den globala miljön även p˚averkas av biosfären genom olika sorters<br />

˚aterkopplingar. Denna insikt har möjliggjorts dels tack vare ett nytt globalt perspektiv p˚a<br />

jorden som helhet, <strong>och</strong> dels tack vare en ökad först˚aelse av dynamiska system. I de ovan<br />

nämnda exemplen har man följande typers ˚aterkoppling fr˚an biosfären:<br />

• Atmosfärens sammansättning h˚alls konstant av biosfären. Syre tillförs atmosfären genom<br />

de gröna växternas fotosyntes. En ökning av syrehalten ökar biologisk aktivitet (s˚asom<br />

förmultning) <strong>och</strong> bränder, vilka förbrukar syre. En minskning av syrehalten ˚ater gör de<br />

syreförbrukande processerna l˚angsammare.<br />

• Den salt som tillförs oceanerna med floder <strong>och</strong> förvittring avlägsnas av mikro-organismer<br />

som binder mineral (bl.a. i sina skal) <strong>och</strong> sedimenteras p˚a havets botten. Ju mera salt,<br />

desto flera mikro-organismer, vilket h˚aller salthalten nere.<br />

• Jordens medeltemperatur regleras via vegetationens inverkan: en höjd temperatur f˚ar<br />

växter som reflekterar infallande str˚alning att öka i antal, eftersom den lokala temperaturen<br />

är n˚agot lägre i närheten av dessa växter; en sjunkande temperatur ˚ater f˚ar växter<br />

som absorberar infallande str˚alning att öka i antal, eftersom den lokala temperaturen<br />

är n˚agot högre vid dessa växter. Slutresultatet blir globalt att temperaturförändringar<br />

pga av variationer i infallande str˚alning motverkas: en ökad str˚alning ger upphov till en<br />

vegetation som reflekterar mera str˚alning, <strong>och</strong> en minskning av str˚alningen ger upphov<br />

till en vegetation som absorberar mera str˚alning.<br />

I vart <strong>och</strong> ett av fallen bidrar s˚aledes biosfären med en ˚aterkopplingsmekanism som motverkar<br />

förändringar helt enkelt genom den inverkan som dessa förändringar har p˚a biosfärens aktivitet.<br />

För att först˚a de globala förh˚allandena p˚a jorden bör dessa ˚aterkopplingar beaktas.<br />

Teorin att biologisk aktivitet p˚averkar atmosfärens <strong>och</strong> oceanernas sammansättningar <strong>och</strong><br />

även klimatet kallas Gaia-teorin (fr˚an Gaia – jordens gudinna i grekisk mytologi). Fenomenet<br />

är analogt med de fysiologiska reglerprocesserna, <strong>och</strong> man har därför ocks˚a infört begreppet<br />

’planetär fysiologi’. Gaia-teorin har intressanta implikationer för diskussionen av miljöproblem,<br />

inte minst för den av människan förorsakade växthuseffekten. Teorin säger˚a ena sidan att<br />

25


iosfären har en förm˚aga att kompensera för inverkan av störningar, inklusive de som människan<br />

˚astadkommer. ˚ A andra sidan säger oss teorin ocks˚a att denna förm˚aga är beroende av<br />

speciella ˚aterkopplingsmekanismer i miljön. Det är därför speciellt viktigt att värna om att<br />

naturens ˚aterkopplingsmekanismer h˚alls intakta. Om dessa mekanismer skadas kan följderna<br />

vara drastiska. Kunskapen om mekanismerna är ännu ganska bristfällig, men klart är att<br />

vegetation, s˚asom tropikernas regnskogar eller oceanernas plankton, spelar en viktig roll.<br />

2.7 Negativ <strong>och</strong> positiv ˚aterkoppling<br />

I samtliga exempel hittills har effekten av ˚aterkopplingen varit att motverka förändringen hos<br />

ett system. Man talar i detta fall om negativ˚aterkoppling. Om vi ändrar tecknet hos regulatorn<br />

Gc i figur 2.11 kommer˚aterkopplingens effekt att vara den omvända: en förändring i utsignalen<br />

p˚averkar insignalen i en riktning som ger upphov till en ännu större utsignalförändring. Detta<br />

kallas positiv ˚aterkoppling.<br />

Medan negativ ˚aterkoppling stabiliserar ett system vid ett önskat tillst˚and, har positiv<br />

˚aterkoppling motsatt verkan: systemet f˚as att snabbt driva ifr˚an starttillst˚andet. Positiv<br />

˚aterkoppling har begränsad användning inom tekniken, men förekommer allmänt i andra<br />

situationer. N˚agra exempel f˚ar belysa funktionen hos positiv ˚aterkoppling.<br />

• Autokatalys i kemiska reaktioner. Slutprodukten fungerar även som en katalysator<br />

för tidigare reaktionssteg, <strong>och</strong> därmed försnabbar reaktionen. Autokatalys förekommer<br />

allmänt i biokemiska reaktioner.<br />

• Inom ekonomin finns en mängd positivt ˚aterkopplade processer. Ett exempel är introduktionen<br />

av nya produkter p˚a marknaden. Försäljningen av t.ex. en mobiltelefon<br />

med nya egenskaper f˚ar fart först d˚a tillräckligt m˚anga människor redan skaffat sig<br />

produkten. Detta är positiv ˚aterkoppling fr˚an produktens utbreddhet till försäljningen.<br />

En liknande ˚aterkoppling inom ekonomin gynnar stora producenter: en hög volym gör<br />

det möjligt att framställa billigare produkter <strong>och</strong> satsa p˚a produktutveckling, med en<br />

ytterligare volymökning till följd.<br />

• ’De rika blir allt rikare, de fattiga allt fattigare’. Ett negativt exempel p˚a positiv˚aterkoppling!<br />

2.8 Litteratur<br />

Det finns ett antal utmärkta läroböcker i <strong>reglerteknik</strong>. Bland svenskspr˚akiga böcker i omr˚adet<br />

ges en praktisk <strong>och</strong> mycket elementär introduktion till reglertekniska idéer av Hägglund<br />

(1997). N˚agot mera teoretiska introduktioner ges av Lennartson (2001), Glad <strong>och</strong> Ljung<br />

(1989), Schmidtbauer (1995) <strong>och</strong> Thomas (1992). Harnefors <strong>och</strong> medarbetare (2004) ger en<br />

26


generellare introduktion till signaler <strong>och</strong> system, som förutom <strong>reglerteknik</strong> ocks˚a behandlar<br />

grunderna i signalbehandling. Bland engelskspr˚akiga böcker kan nämnas Albertos <strong>och</strong> Mareels<br />

(2010), som ger en elementär <strong>och</strong> allmänbildande introduktion till omr˚adet.<br />

De icke-tekniska exemplen p˚a reglerprocesser som i korthet beskrivits i avsnitt 2.6 inbegriper<br />

alltför omfattande problemomr˚aden för att kunna behandlas p˚a grundkursniv˚a. För<br />

den intresserade kan vi nämna t.ex. Khoo (2000), som beskriver reglermekanismer i fysiologiska<br />

system. Gaia-teorin introducerades av den brittiske kemisten James Lovelock (1982).<br />

En kvantitativ behandling av dynamiska system <strong>och</strong> ˚aterkopplingsmekanismer i atmosfären<br />

<strong>och</strong> oceanerna diskuteras bl.a. i Kump <strong>och</strong> medförfattare (2000).<br />

Referenser<br />

Albertos, P. <strong>och</strong> I. Mareels (2010). Feedback and Control for Everyone. Springer.<br />

Harnefors, L., J. Holmberg <strong>och</strong> J. Lundqvist (2004). Signaler <strong>och</strong> system. Liber AB.<br />

Hägglund, T. (1997). Praktisk processreglering. Studentlitteratur.<br />

Khoo, M. (2000). Physiological Control Systems. IEEE Press.<br />

Kump, L. R., J. F. Kasting R. G. Crane (2000). The Earth System. Prentice Hall.<br />

Lennartson, B. (2001). Reglerteknikens grunder. Studentlitteratur.<br />

Ljung, L. <strong>och</strong> T. Glad (1989). Reglerteknik. Grundläggande teori Studentlitteratur.<br />

Lovelock, J. (1982). Gaia: A New Look at Life on Earth. Oxford University Press.<br />

Schmidtbauer, B. (1995). Analog <strong>och</strong> digital <strong>reglerteknik</strong>. Studentlitteratur.<br />

Thomas, B. (1992). Modern <strong>reglerteknik</strong>. Liber AB.<br />

27


Kapitel 3<br />

Dynamiska system<br />

3.1 Enkla systemtyper <strong>och</strong> deras stegsvar<br />

För att kunna konstruera regulatorer för dynamiska system bör systemens egenskaper vara<br />

kända. Innan vi g˚ar vidare till att behandla modeller för dynamiska system skall vi ge en<br />

kortfattad kvalitativ beskrivning av de viktigaste systemtyperna som man bör känna till. För<br />

att ge en kvalitativ inblick i de olika systemtyperna kommer vi att betrakta deras stegsvar,<br />

dvs utsignalen y(t) efter en stegförändring i insignalen u(t).<br />

Första ordningens system.<br />

Stegsvaret hos ett system av första ordningen har getts i exempel 2.4, figur 2.8. Karakteristisk<br />

för ett system av denna typ är att tidsderivatan dy/dt är olikt noll omedelbart efter<br />

stegförändringen. Exemplen i kapitel 2 illustrerar hur system av denna typ kan f˚as genom<br />

fysikalisk modellering.<br />

System med tv˚a tidskonstanter.<br />

Genom att ha tv˚a system av första ordningen i serie f˚as ett system med tv˚a tidskonstanter.<br />

Stegsvaret karakteriseras av en l˚angsammare respons i början, med dy/dt = 0 vid t = 0, jfr<br />

figur 3.1.<br />

System av denna typ f˚as i praktiken p˚a samma sätt som första ordningens system. Ett<br />

system med tv˚a tidskonstanter f˚as t.ex. genom att modellera motorn i exempel 2.1 som ett<br />

första ordningens system.<br />

System med översväng.<br />

Figur 3.2 visar stegsvaret hos ett system med översväng. Detta är typiskt för system mekaniska<br />

system med fjädrande element. N˚agra exempel är en kran med en hängande last, en pendel,<br />

fjädringen i en bil eller varvtalet hos en motor med flexibel koppling. Systemets tendens<br />

till svängningar gör det sv˚arare att reglera. En viktig uppgift för regleringen är att dämpa<br />

svängningarna. För att kunna göra detta bör man ha en tillräckligt god modell av systemet.<br />

28


1.2<br />

1<br />

0.8<br />

0.6<br />

u<br />

0.4<br />

0.2<br />

0<br />

−0.2<br />

1.2<br />

1<br />

0.8<br />

y 0.6<br />

0.4<br />

0.2<br />

0<br />

−2 0 2 4 6 8 10 12 14 16<br />

tid<br />

Figur 3.1: Responsen hos ett system med tv˚a tidskonstanter.<br />

System med dödtid.<br />

Dödtid eller tidsfördröjning (eng. dead time, time delay; fi. kuollut aika, aikaviive) innebär att<br />

det tar en tid L innan en förändring i insignalen syns i utsignalen. Stegsvaret hos ett system<br />

av första ordningen med dödtid illustreras i figur 3.3.<br />

Eftersom styrsignalen inverkan syns först efter en tid är system med dödtid sv˚ara att<br />

reglera. För att kunna vidta korrekt regler˚atgärd p˚a basen av mätningen y(t) vid tiden t<br />

bör regulatorn kunna förutse hur systemet kommer att bete sig mellan tiden t <strong>och</strong> t + L, d˚a<br />

regler˚atgärdens inverkan p˚a utsignalen kan observeras.<br />

Dödtider förorsakas i praktiken vanligen av olika typer av transporttider, s˚asom vid förflyttning<br />

av material eller vid vätske- eller gasflöden. System med dödtid är därför mycket<br />

vanliga i processindustrin.<br />

System med integration.<br />

Om utsignalen bestäms av integralen hos insignalen, enligt<br />

t<br />

y(t) = k<br />

0<br />

(u(τ) − u0) dτ (3.1)<br />

finns det endast ett värde u0 p˚a insignalen för vilket utsignalen y h˚alls konstant. Varje<br />

avvikelse fr˚an u0 f˚ar y antingen att ständigt öka eller ständigt minska. Se figur 3.4. Typiska<br />

exempel p˚a system med integration är mängden av material y som man har i ett lager med<br />

konstant utströmning u0. Mängden y h˚alls konstant endast om inflödet u till lagret är exakt<br />

lika stort som utströmmen u0. För u > u0 ökar y med konstant hastighet <strong>och</strong> för u < u0<br />

minskar y med konstant hastighet till lagret är tomt. Ett konkret exempel är vätskeinneh˚allet<br />

29


1.2<br />

1<br />

0.8<br />

0.6<br />

u<br />

0.4<br />

0.2<br />

0<br />

−0.2<br />

y<br />

1.5<br />

1<br />

0.5<br />

0<br />

−2 −1 0 1 2 3<br />

tid<br />

4 5 6 7 8<br />

Figur 3.2: Responsen hos ett system med översväng.<br />

y i en beh˚allare med konstant utflöde u0 <strong>och</strong> inflödet u.<br />

Ett annat viktigt exempel p˚a ett system med integration är servomotorer som används<br />

vid positionsreglering. Positionen y styrs med spänningen u till motorn. Rörelsehastigheten<br />

är proportionell mot spänningen, dvs dy/dt = ku. En avvikelse fr˚an u = 0 ger därför upphov<br />

till en ih˚allande förändring i positionen.<br />

System med inverssvar.<br />

Figur 3.5 visar stegsvaret hos en process med inverssvar. Kännetecknande för dessa system<br />

är att svaret startar ˚at motsatta h˚allet innan det närmar sig det statiska värdet. Behovet<br />

att beakta den initiala motsatta verkan av en regler˚atgärd gör att system med inverssvar är<br />

mycket sv˚ara att reglera.<br />

System med inverssvar uppst˚ar t.ex. d˚a tv˚a delprocesser i ett system verkar ˚at motsatt<br />

h˚all. Om en av processerna har ett snabbt svar <strong>och</strong> den andra har ett l˚angsamt svar, kommer<br />

den snabba processen till en början att f˚a utsignalen att g˚a mot ett h˚all innan den l˚angsamma<br />

processen hinner p˚averka utsignalen. Ett exempel p˚a denna typs system är vid förbränning<br />

av fasta bränslen (s˚asom flis). En ökning av bränslemängden f˚ar temperaturen i eldhärden<br />

först att minska eftersom det inkommande bränslet har en lägre temperatur. Efter att förbränningen<br />

av det tillförda bränslet kommit i g˚ang ökar temperaturen. Andra exempel p˚a<br />

system med inverssvar finns t.ex. vid reglering av vissa flygplan. Inom ekonomiska system är<br />

inverssvar vanliga: t.ex. en skattesänkning sänker till en början skatteinkomsterna pga den<br />

lägre skatten, men kan senare resultera i en större skatteintäkt tack vare ökad ekonomisk<br />

aktivitet.<br />

30


1.2<br />

1<br />

0.8<br />

0.6<br />

u<br />

0.4<br />

0.2<br />

0<br />

−0.2<br />

1.2<br />

1<br />

0.8<br />

0.6<br />

y<br />

0.4<br />

0.2<br />

0<br />

−0.2<br />

−2 0 2 4<br />

tid<br />

6 8 10<br />

Figur 3.3: Responsen hos ett system av första ordningen med dödtid (L = 1).<br />

Instabila system.<br />

Instabila system karakteriseras av att de divergerar fr˚an sitt begynnelsetillst˚and om de lämnas<br />

˚at sig själva. Ett enkelt exempel p˚a ett instabilt system är en inverterad pendel, som faller<br />

om den inte kontinuerligt balanseras. Det enda sättet att h˚alla ett instabilt system vid ett<br />

börvärde är genom att använda ˚aterkoppling. Vissa flygplan är instabila <strong>och</strong> behöver därför<br />

ständiga regler˚atgärder för att h˚alla kursen. Exempel p˚a instabila system inom processindustrin<br />

är vissa reaktorer med exoterma reaktioner, som bör kylas tillräckligt för att h˚alla<br />

reaktionen under kontroll. Ett annat exempel p˚a instabilitet är att backa ett fordon med släp.<br />

Systemet är instabilt eftersom den minsta avvikelse i kursen f˚ar släpet att driva ˚at sidan. Det<br />

enda sättet att h˚alla kursen är att ständigt kompensera kursavvikelserna med styrningen. En<br />

regulator som stabiliserar systemet kan backa ett fordon med släp utan sv˚arighet.<br />

3.2 Linjära system<br />

I kapitel 2 har vi sett att dynamiska system beskrivs av differentialekvationer. Vi skall i<br />

detta kapitel närmare introducera modeller av de vanligaste typerna av dynamiska system.<br />

Generellt kan ett linjärt dynamiskt system G med insignalen u <strong>och</strong> utsignalen y beskrivas<br />

med den linjära differentialekvationen<br />

dny(t) +a1<br />

dtn dn−1y(t) dy(t)<br />

+· · ·+an−1<br />

dtn−1 dt +any(t) = b0<br />

där n är systemet ordning. Vanligtvis gäller m < n.<br />

31<br />

dmu(t) du(t)<br />

+· · ·+bm−1<br />

dtm dt +bmu(t) (3.2)


1.2<br />

1<br />

0.8<br />

0.6<br />

u<br />

0.4<br />

0.2<br />

0<br />

−0.2<br />

6<br />

5<br />

4<br />

y 3<br />

2<br />

1<br />

0<br />

−2 0 2 4<br />

tid<br />

6 8 10<br />

Figur 3.4: Responsen hos ett system med integration.<br />

Systemekvationen (3.2) kan skrivas i en bekvämare form genom att introducear beteckningen<br />

p = d<br />

(3.3)<br />

dt<br />

för differentialoperatorn. Eftersom<br />

följer generellt att<br />

d2 <br />

d d<br />

y =<br />

dt2 dt dt y<br />

<br />

= p 2 y (3.4)<br />

d k<br />

dt k y = pk y (3.5)<br />

Ekvation (3.2) kan s˚aledes skrivas i formen<br />

eller<br />

p n y + a1p n−1 y + · · · + an−1py + any = b0p m u + b1p m−1 u + · · · + bm−1pu + bmu (3.6)<br />

där vi introducerat polynomen<br />

A(p)y = B(p)u (3.7)<br />

A(p) = p n + a1p n−1 + · · · + an−1p + an<br />

B(p) = b0p m + b1p m−1 + · · · + bm−1p + bm<br />

32<br />

(3.8)<br />

(3.9)


1.2<br />

1<br />

0.8<br />

0.6<br />

u<br />

0.4<br />

0.2<br />

0<br />

−0.2<br />

1.5<br />

1<br />

y 0.5<br />

Löser vi ut y ur ekvation (3.7) f˚as<br />

eller<br />

där vi definierat<br />

0<br />

−0.5<br />

−2 −1 0 1 2 3<br />

tid<br />

4 5 6 7 8<br />

Figur 3.5: Responsen hos ett system med inverssvar.<br />

y(t) = B(p)<br />

u(t) (3.10)<br />

A(p)<br />

y(t) = G(p)u(t) (3.11)<br />

G(p) = B(p)<br />

A(p)<br />

(3.12)<br />

Eftersom differentialoperatorn p är en operator kan även G(p) uppfattas som en operator som<br />

transformerar signalen u(t) till en annan signal y(t). Operatorn G(p) kallas systemets överföringsoperator.<br />

Alternativt kan man betrakta G(p) som en funktion av operatorn p, varför<br />

man ocks˚a kallar G(p) överföringsfunktion (eng. transfer function; fi. siirtofunktio). Av n˚agon<br />

orsak är den förra termen bruklig i svenskan, medan den senare används i engelskspr˚akig text.<br />

Anmärkning 3.1<br />

Begreppet överföringsoperator kan defineras mera rigoröst via Laplace-transformen som en<br />

operator som opererar p˚a en Laplace-transformerad signal. Denna metod är mera generell<br />

än den som används här, men överföringsoperatorn har samma form i bägge fallen. Laplacetransformen<br />

kommer att behandlas i senare kurser.<br />

System som beskrivs av linjära differentialekvationer av typen (3.2) har ett antal viktiga<br />

egenskaper, vilka gör deras behandling speciellt enkel. Egenskaperna följer direkt ur differentialoperatorns<br />

egenskaper.<br />

33


1.2<br />

1<br />

0.8<br />

0.6<br />

u<br />

0.4<br />

0.2<br />

0<br />

−0.2<br />

100<br />

80<br />

60<br />

y<br />

40<br />

20<br />

0<br />

−2 −1 0 1 2 3 4 5<br />

tid<br />

Figur 3.6: Responsen hos ett instabilt system.<br />

• Superpositionsprincipen: om insignalen u1 till systemet G ger utsignalen y1 = Gu1,<br />

<strong>och</strong> insignalen u2 ger utsignalen y2 = Gu2, gäller att insignalen u1 + u2 ger utsignalen<br />

y1 + y2, dvs<br />

y = G(u1 + u2) = Gu1 + Gu2 = y1 + y2<br />

(3.13)<br />

• Parallellkoppling av tv˚a system med överföringsoperatorerna G1 <strong>och</strong> G2 är ekvivalent<br />

med ett system med överföringsoperatorn G1 + G2, ty utsignalen fr˚an parallellkopplade<br />

system ges av<br />

y = G1u + G2u = (G1 + G2)u (3.14)<br />

Jfr figur 3.7.<br />

• Seriekoppling av tv˚a system med överföringsoperatorerna G1 <strong>och</strong> G2 är ekvivalent med<br />

ett system med överföringsoperatorn G1G2 = G2G1, ty utsignalen fr˚an seriekopplade<br />

system ges av<br />

y = G1 (G2u) = G1G2u = G2G1u (3.15)<br />

Jfr figur 3.8.<br />

Observera speciellt att det ur ovan följer att utsignalen blir densamma oberoende av<br />

systemens ordningsföljd, ty<br />

y(t) = G1(p)G2(p)u(t) = G2(p)G1(p)u(t) (3.16)<br />

34


u<br />

u<br />

u<br />

✲<br />

✲<br />

G1<br />

G2<br />

y1 = G1u<br />

u ✲ y = (G1 + G2)u<br />

✲ G1<br />

y2 = G2u<br />

Figur 3.7: Parallellkopplade system.<br />

y1 = G1u<br />

✲ G2<br />

✲ y = G2y1 = G2G1u<br />

Figur 3.8: Seriekopplade system.<br />

Problem 3.1<br />

Betrakta tv˚a första ordningenens system G1 <strong>och</strong> G2, där G1 definieras av differentialekvationen<br />

dy1(t)<br />

dt + y1(t) = u1(t) (3.17)<br />

<strong>och</strong> G2 definieras av differentialekvationen<br />

dy2(t)<br />

dt + 2y2(t) = 3u2(t) (3.18)<br />

- Bestäm systemens överföringsoperatorer.<br />

- Härled differentialekvationen som beskriver en parallellkoppling av G1 <strong>och</strong> G2. Verifiera att<br />

överföringsoperatorn hos det parallellkopplade systemet ges av (3.14).<br />

- Härled differentialekvationen som beskriver en seriekoppling av G1 <strong>och</strong> G2. Verifiera att<br />

överföringsoperatorn hos det seriekopplade systemet ges av (3.15).<br />

Sambandet (3.11) säger att utsignalen y kan uttryckas genom att multiplicera insignalen<br />

u med överföringsoperatorn G(p). Att detta kan göras följer ur differentialoperatorns linjäritet.<br />

Det faktum att differentialekvationssambandet kan uttryckas i form av en multiplikation<br />

gör det möjligt att härleda differentialekvationer för kopplade system genom rent<br />

algebraiska manipulationer, best˚aende av endast multiplikationer <strong>och</strong> additioner. Betrakta<br />

t.ex. den ˚aterkopplade reglerkretsen i figur 3.9. Genom att uttrycka signalsambanden med<br />

hjälp av överföringsoperatorerna Gp(p) <strong>och</strong> Gc(p) <strong>och</strong> eliminera signalerna e <strong>och</strong> u med hjälp<br />

av algebraiska manipulationer f˚ar vi att sambandet mellan signalen r <strong>och</strong> utsignalen y beskrivs<br />

av överföringsoperatorn<br />

y(t) = G(p)r(t) (3.19)<br />

35


där<br />

r<br />

+<br />

✲ ❡ e ✲ u<br />

Gc<br />

✲ Gp<br />

✲<br />

− ✻<br />

Figur 3.9: ˚ Aterkopplad krets.<br />

G(p) = Gp(p)Gc(p)<br />

1 + Gp(p)Gc(p)<br />

y<br />

(3.20)<br />

Motsvarande differentialekvation kan sedan bestämmas genom att utnyttja sambandet (3.12)<br />

mellan överföringsoperatorn <strong>och</strong> differentialekvationen.<br />

Problem 3.2<br />

Härled sambandet (3.19), (3.20) mellan r <strong>och</strong> y.<br />

Problem 3.3<br />

Betrakta det ˚aterkopplade systemet i figur 3.9, <strong>och</strong> antag att systemet Gp beskrivs av differentialekvationen<br />

dy(t)<br />

+ y(t) = u(t) (3.21)<br />

dt<br />

<strong>och</strong> regulatorn Gc beskrivs av differentialekvationen<br />

du(t)<br />

dt<br />

+ 2u(t) = 2e(t) (3.22)<br />

Bestäm differentialekvationen som beskriver sambandet mellan signalerna r <strong>och</strong> y.<br />

36


3.3 Modellering av enkla standardsystem<br />

I detta avsnitt skall vi mera detaljerat studera de enkla standardsystemtyperna som beskrevs<br />

i avsnitt 3.1. Speciellt skall vi visa hur enkla differentialekvationsmodeller ger upphov till<br />

de olika stegsvaren. Vi skall allts˚a bestämma den utsignal y(t) som f˚as, d˚a insignalen är<br />

stegformad, enligt<br />

samt y(t) = 0, t < 0.<br />

u(t) =<br />

3.3.1 Första ordningens system<br />

0 , för t < 0<br />

usteg , för t ≥ 0<br />

Ett stabilt system av första ordningen beskrivs av differentialekvationen<br />

dy(t)<br />

dt<br />

där a > 0. Man brukar ofta skriva systemekvationen i formen<br />

T dy(t)<br />

dt<br />

där T = 1/a <strong>och</strong> K = b/a. Systemets överföringsoperator är<br />

<strong>och</strong> det har stegsvaret<br />

(3.23)<br />

+ ay(t) = bu(t) (3.24)<br />

+ y(t) = Ku(t) (3.25)<br />

G(p) = b K<br />

=<br />

p + a T p + 1<br />

ysteg(t) = K(1 − e −t/T )usteg<br />

Jfr figur 2.8. Systemets statiska förstärkning är lika med K, dvs<br />

(3.26)<br />

(3.27)<br />

ysteg(t) → Kusteg, d˚a t → ∞ (3.28)<br />

Parametern T kallas systemets tidskonstant (eng. time constant; fi. aikavakio), <strong>och</strong> är ett<br />

m˚att p˚a hur snabbt systemet reagerar; ju mindre T (> 0), desto snabbare respons. Speciellt<br />

gäller, att vid tiden t = T har (1 − e −1 ) × 100% = 63.2% av den totala förändringen n˚atts.<br />

Problem 3.4<br />

Verifiera att systemet (3.25) har stegsvaret (3.27).<br />

37


3.3.2 System med tv˚a tidskonstanter<br />

Ett system med tv˚a tidskonstanter best˚ar av tv˚a seriekopplade första ordningens system. L˚at<br />

det första systemet, G1, beskrivas av<br />

dy1(t)<br />

T1<br />

dt + y1(t) = K1u(t) (3.29)<br />

med utsignalen y1(t), som är insignal till det andra systemet, G2, som beskrivs av<br />

dy(t)<br />

T2<br />

dt + y(t) = K2y1(t) (3.30)<br />

Enligt tidigare har vi att det seriekopplade systemet har överföringsoperatorn (jfr (3.15))<br />

G(p) = G2(p)G1(p) =<br />

K<br />

(T2p + 1)(T1p + 1)<br />

(3.31)<br />

där K = K1K2. Detta är ett andra ordningens system.<br />

Stegsvaret hos det seriekopplade systemet G2G1 kan bestämmas genom att observera att<br />

om T1 = T2 kan G(p) skrivas i formen<br />

G(p) = KT1/(T1 − T2)<br />

T1p + 1<br />

+ KT2/(T2 − T1)<br />

T2p + 1<br />

(3.32)<br />

Detta karakteriserar G(p) i form av tv˚a parallellkopplade system av första ordningen. Enligt<br />

(3.14) ges stegsvaret av summan av de enskilda systemens stegsvar i (3.32), vilket enligt (3.27)<br />

är<br />

ysteg(t) =<br />

= K<br />

KT1<br />

T1 − T2<br />

<br />

1 − e −t/T1<br />

<br />

+ KT2<br />

<br />

T2 − T1<br />

<br />

1 − T1<br />

e<br />

T1 − T2<br />

−t/T1<br />

T2<br />

− e<br />

T2 − T1<br />

−t/T2<br />

Man kan visa att om T1 = T2 = T ges stegsvaret av<br />

ysteg(t) = K<br />

<br />

1 − (1 + t<br />

T )e−t/T<br />

<br />

usteg<br />

1 − e −t/T2<br />

<br />

usteg<br />

usteg<br />

(3.33)<br />

(3.34)<br />

Figur 3.1 visar stegsvaret för ett system med tv˚a tidskonstanter. Jämfört med ett första<br />

ordningens system har stegsvaret hos ett system med tv˚a tidskonstanter en kontinuerligt<br />

varierande derivata, s˚a att dy(t)/dt = 0 vid t = 0. I analogi med ett första ordningens system<br />

sker förändringen monotont fr˚an begynnelsevärdet till det statiska värdet, <strong>och</strong> stegsvaret<br />

saknar allts˚a översväng.<br />

P˚a motsvarande sätt kan man bilda system med flera tidskonstanter genom seriekoppling<br />

av flera första ordningens system.<br />

Problem 3.5<br />

Verifiera att systemet (3.31) har stegsvaret (3.34) i fallet T1 = T2 = T .<br />

38


3.3.3 System med översväng<br />

Det seriekopplade systemet i ekvation (3.31) är ett system av andra ordningen med tv˚a reella<br />

tidskonstanter. Överföringsoperatorns nämnarpolynom A(p) har i detta fall de tv˚a reella<br />

nollställena −1/T1 <strong>och</strong> −1/T2. System med översväng beskrivs av differentialekvationer med<br />

den egenskapen att polynomet A(p) har komplexkonjugerade nollställen. Ett generellt system<br />

av andra ordningen beskrivs av<br />

d2y(t) dy(t)<br />

+ a1<br />

dt2 dt + a2y(t)<br />

du(t)<br />

= b1<br />

dt + b2u(t) (3.35)<br />

<strong>och</strong> motsvarande överföringsoperator<br />

b1p + b2<br />

G(p) =<br />

p2 (3.36)<br />

+ a1p + a2<br />

För att belysa de karakteristiska egenskaperna hos system av andra ordningen räcker det<br />

med att betrakta fallet d˚a b1 = 0. Systemets statiska förstärkning, den transienta responsens<br />

snabbhet samt översvängens storlek bestäms d˚a entydigt av systemparametrarna. Det visar<br />

sig att dessa tre systemegenskaper kan karakteriseras explicit genom att skriva systemekvationen<br />

i formen<br />

d2y(t) dy(t)<br />

+ 2ζωn<br />

dt2 dt + ω2 ny(t) = Kω 2 nu(t) (3.37)<br />

<strong>och</strong> motsvarande överföringsoperator<br />

G(p) =<br />

p2 + 2ζωnp + ω2 n<br />

Här är K systemets statiska förstärkning, ζ kallas relativ dämpning, <strong>och</strong> ωn är systemets<br />

odämpade egenfrekvens. Alla stabila system av andra ordningen kan skrivas i formen (3.37)<br />

med ζ > 0 <strong>och</strong> ωn > 0. Den relativa dämpningen ζ bestämmer formen hos systemets stegsvar.<br />

I fallet ζ > 1 sägs systemet vara överdämpat. Polynomet A(p) = p2 + 2ζωnp + ω2 n har härvid<br />

tv˚a negativa reella nollställen <strong>och</strong> systemet är ekvivalent med tv˚a seriekopplade system av<br />

första ordningen (jfr avsnitt 3.3.2). I fallet ζ < 1 sägs systemet vara underdämpat. Polynomet<br />

A(p) har i detta fall tv˚a komplexkonjugerade nollställen med negativa realdelar, <strong>och</strong><br />

systemets stegsvar har en översväng, jfr figur 3.2. I gränsfallet ζ = 1 kallas systemet kritiskt<br />

dämpat. Polynomet A(p) har d˚a ett dubbelt nollställe som är reell <strong>och</strong> negativ. Den odämpade<br />

egenfrekvensen bestämmer tidskalan hos systemets respons; ju större ωn, desto snabbare<br />

respons.<br />

Stegsvaren för överdämpade <strong>och</strong> kritiskt dämpade system ges av (3.33) respektive (3.34).<br />

Kω 2 n<br />

I det underdämpade fallet (ζ < 1) ges stegsvaret av<br />

<br />

ysteg(t) = K 1 − 1<br />

β e−ζωnt <br />

[ζ sin(βωnt) + β cos(βωnt)]<br />

<br />

= K 1 − 1<br />

β e−ζωnt <br />

sin(βωnt + ϕ) usteg<br />

39<br />

usteg<br />

(3.38)<br />

(3.39)


där<br />

<strong>och</strong><br />

β =<br />

<br />

1 − ζ 2 (3.40)<br />

ϕ = arccos(ζ) (3.41)<br />

Man kan visa att översvängens maximala storlek M ges av (angiven i % av statiska värdet)<br />

<strong>och</strong> tiden tM vid vilken den maximala översvängen f˚as är<br />

<br />

−ζπ<br />

M = exp × 100% (3.42)<br />

β<br />

tM = π<br />

ωnβ<br />

(3.43)<br />

Figur 3.10 illustrerar stegsvar hos överdämpade, kritiskt dämpade <strong>och</strong> underdämpade<br />

system av andra ordningen.<br />

Problem 3.6<br />

Verifiera att systemet (3.37) har stegsvaret (3.39) om ζ < 1.<br />

y<br />

1.8<br />

1.6<br />

1.4<br />

1.2<br />

1<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

0<br />

−0.2<br />

0 1 2 3 4 5<br />

tid<br />

6 7 8 9 10<br />

Figur 3.10: Responsen hos system av andra ordningen med relativa dämpningarna ζ = 2, 1, 0.5<br />

<strong>och</strong> 0.1 (nerifr˚an räknat) samt odämpade egenfrekvensen ωn = 1 <strong>och</strong> statiska förstärkningen<br />

K = 1.<br />

I fysikaliska system uppst˚ar ett dynamiskt beteende med översväng typiskt av att systemet<br />

oskillerar mellan tv˚a tillst˚and, s˚asom följande exempel illustrerar.<br />

40


k<br />

m<br />

u(t)<br />

Figur 3.11: Fjäderupphängd massa med dämpning.<br />

Exempel 3.1 - Fjädringssystem med dämpning.<br />

Figur 3.11 visar en massa m upphängd i en fjäder. Positionsavvikelsen y(t) fr˚an ett jämviktsläge<br />

p˚averkas av kraften u(t). Vid jämviktsläge gäller y(t) = 0 <strong>och</strong> u(t) = 0. Newtons andra<br />

lag ger<br />

m d2 y(t)<br />

dt 2 = u(t) + Ff (t) + Fd(t)<br />

där Ff (t) är fjäderkraften, som är proportionell mot avvikelsen y(t) <strong>och</strong> riktad i motsatt<br />

riktning,<br />

Ff (t) = −ky(t)<br />

<strong>och</strong> Fd(t) är den dämpande kraften i dämparen, som är proporionell mot hastigheten dy(t)/dt<br />

<strong>och</strong> riktad i motsatt riktning,<br />

Fd(t) = −b dy<br />

dt<br />

Kombination av ekvationerna ger<br />

b<br />

y(t)<br />

m d2y(t) = u(t) − ky(t) − bdy<br />

dt2 dt<br />

eller<br />

d2y(t) b dy k 1<br />

+ + y(t) =<br />

dt2 m dt m m u(t)<br />

Detta är ett system av andra ordningen med överföringsoperatorrepresentationen<br />

y(t) =<br />

1/m<br />

p2 + b k<br />

mp + m<br />

41<br />

u(t)


Jämförelse med (3.38) visar att systemets odämpade egenfrekvens är<br />

<br />

k<br />

ωn =<br />

m<br />

den statiska förstärkningen ges av<br />

K = 1/m<br />

ω 2 n<br />

= 1<br />

k ,<br />

<strong>och</strong> relativa dämpningen är<br />

ζ = b/m<br />

=<br />

2ωn<br />

1 b<br />

√<br />

2 mk<br />

Sambandet visar att en stor dämpning (b) <strong>och</strong> liten fjäderkonstant k gör systemet överdämpat,<br />

medan ett system med liten dämpning <strong>och</strong> stor fjäderkonstant är underdämpad.<br />

Vi skall ännu visa hur stegsvaret hos ett system (3.35) av andra ordningen med b1 = 0<br />

kan härledas. Överföringsoperatorn (3.36) kan skrivas<br />

där<br />

<strong>och</strong><br />

G(p) = G1(p) + G2(p) (3.44)<br />

b1p<br />

G1(p) =<br />

p2 + a1p + a2<br />

(3.45)<br />

G2(p) =<br />

p2 (3.46)<br />

+ a1p + a2<br />

Svaret hos G(p) f˚as s˚aledes som en summa av svaren fr˚an tv˚a parallellkopplade system med<br />

överföringsoperatorerna G1(p) <strong>och</strong> G2(p). Stegsvaret för G2 har givits ovan. Stegsvaret hos<br />

G1 kan bestämmas genom att skriva<br />

G1(p) = p<br />

p2 (3.47)<br />

+ a1p + a2<br />

Detta representerar G1 i form av en seriekoppling av ett system vars stegsvar ges av (3.33),<br />

(3.34) eller (3.39) ˚atföljd av differentialoperatorn p = d<br />

dt . Stegsvaret y1(t) hos G1 kan därför<br />

bestämmas genom att först bestämma stegsvaret hos systemet<br />

˜y1(t) =<br />

p2 u(t)<br />

+ a1p + a2<br />

vilket ges av uttrycket (3.33) om systemet är överdämpat, uttrycket (3.34) om systemet är<br />

kritiskt dämpat, eller uttrycket (3.39) om systemet är underdämpat, varefter y1(t) ges av<br />

b1<br />

b2<br />

b1<br />

y1(t) = d˜y1(t)<br />

dt<br />

42


3.3.4 System med dödtid<br />

Ett system med dödtid har en tidsfördröjning som fördröjer insignalens effekt p˚a utsignalen.<br />

Ett system av första ordningen med dödtiden L beskrivs s˚aledes av differentialekvationen<br />

dy(t)<br />

dt<br />

+ ay(t) = bu(t − L) (3.48)<br />

Vi kan härleda ett uttryck för överföringsoperatorn hos en tidsfördröjning genom att uttrycka<br />

= p. En Taylor-serieutveckling av u(t) ger<br />

denna med hjälp av differentialoperatorn d<br />

dt<br />

d2u(t) hk<br />

+ · · · +<br />

dt2 k!<br />

d2 hk d<br />

+ · · · +<br />

dt2 k!<br />

k<br />

+ · · ·<br />

dtk u(t + h) = u(t) + h du(t) h2<br />

+<br />

<br />

dt 2!<br />

= 1 + h d h2<br />

+<br />

dt 2!<br />

<br />

= 1 + hp + (hp)2 (hp)k<br />

+ · · · +<br />

2! k!<br />

+ · · ·<br />

dku(t) + · · ·<br />

dtk <br />

<br />

u(t)<br />

u(t)<br />

= e hp u(t) (3.49)<br />

där vi i det sista steget utnyttjat Taylor-serieutvecklingen av exponentialfunktionen e x . För<br />

h = −L ger sambandet (3.49) ett uttryck för tidsfördröjningen med hjälp av differentialoperatorn,<br />

u(t − L) = e −Lp u(t) (3.50)<br />

Det följer att ett tidsfördröjningen har överföringsoperatorn<br />

GL(p) = e −Lp<br />

(3.51)<br />

Systemet (3.48) best˚ar av en seriekoppling av en tidsfördröjning <strong>och</strong> ett system av första<br />

ordning, <strong>och</strong> kan s˚aledes uttryckas med hjälp av överföringsoperatorer (dvs med hjälp av<br />

differentialoperatorn) i formen<br />

y(t) = b<br />

p + a e−Lp u(t) (3.52)<br />

3.3.5 System med integration<br />

Ett system med integration inneh˚aller ett delsystem av formen<br />

dy(t)<br />

dt = k [u(t) − u0] (3.53)<br />

jfr ekvation (3.1). Utsignalen y förändras s˚a länge u(t) = u0. I praktiken ing˚ar systemet (3.53)<br />

vanligen som ett led i en seriekoppling med flera delsystem.<br />

43


3.3.6 System med inverssvar<br />

System med inverssvar kan beskrivas med hjälp av en parallellkoppling av delsystem vars<br />

statiska förstärkningar har olika tecken. För att se hur ett system kan ge upphov till inverssvar<br />

av den typ som visas i figur 3.5, betrakta tv˚a parallelkopplade system av första ordningen,<br />

s˚a att<br />

y(t) = y1(t) + y2(t)<br />

där y1(t) ges av<br />

<strong>och</strong> y2(t) ges av<br />

y1(t) = b1<br />

u(t)<br />

p + a1<br />

y2(t) = b2<br />

u(t)<br />

p + a2<br />

För att det parallellkopplade systemet skall ha ett inverssvar av den typ som visas i figur 3.5<br />

bör utsignalen vara positiv för stora t,<br />

y(t) = y1(t) + y2(t) > 0, d˚a t → ∞<br />

medan utsignalen först skall bli negativ, dvs<br />

<br />

dy(t) <br />

<br />

dt =<br />

t=0<br />

dy1(t)<br />

<br />

<br />

<br />

dt <br />

t=0<br />

+ dy2(t)<br />

dt<br />

<br />

<br />

<br />

< 0<br />

t=0<br />

Det statiska värdet efter en stegformad insignal hos det parallellkopplade systemet ges av<br />

y(t) = y1(t) + y2(t) →<br />

b1<br />

a1<br />

+ b2<br />

a2<br />

medan utsignalens tidsderivata vid t = 0 är<br />

<br />

dy(t) <br />

= (b1 + b2) usteg<br />

dt<br />

Detta implicerar olikheterna<br />

t=0<br />

a1b2 + a2b1 > 0<br />

b1 + b2 < 0<br />

<br />

usteg, d˚a t → ∞<br />

Dessa olikheter kan är ekvivalenta med ett enkelt villkor hos det parallellkopplade systemets<br />

överföringsoperator. Vi noterar nämligen att det parallellkopplade systemets överföringsoperator<br />

Gpar(p) =<br />

b1<br />

p + a1<br />

+ b2<br />

p + a2<br />

= (b1 + b2)p + a1b2 + a2b1<br />

(p + a1)(p + a2)<br />

44<br />

(3.54)


<strong>och</strong> de ovan givna olikheterna implicerar d˚a att överföringsoperatorns täljarpolynom<br />

har ett positivt nollställe,<br />

(b1 + b2)p + a1b2 + a2b1<br />

p = − a1b2 + a2b1<br />

> 0<br />

b1 + b2<br />

Detta resultat kan generaliseras, <strong>och</strong> man kan visa att system med inverssvar generellt har<br />

en överföringsoperator vars täljarpolynom har minst ett nollställe med en positiv realdel.<br />

Som ett exempel visar figur 3.5 stegsvaret hos ett system med överföringsoperatorn<br />

Ginv(p) = 2 (−5)<br />

+<br />

p + 1 p + 5<br />

Genom att skriva uttrycket i (3.55) p˚a gemensam nämnare f˚as<br />

Ginv(p) =<br />

−3p + 5<br />

p 2 + 6p + 5<br />

(3.55)<br />

(3.56)<br />

Detta är ett system av andra ordningen. Observera att täljarpolynomet B(p) = −3p + 5 hos<br />

överföringsoperatorn (3.56) har ett positivt reellt nollställe.<br />

3.3.7 System av högre ordning<br />

De ovan härledda stegsvaren kan användas för att bestämma stegsvar hos system av högre<br />

ordning. Det generella fallet kan helt enkelt behandlas genom att överföringsoperatorn (3.12)<br />

partialbr˚akuppdelas enligt<br />

G(p) = B(p)<br />

A(p) =<br />

r<br />

Gk(p) (3.57)<br />

Detta karakteriserar G(p) som en parallellkoppling av r stycken system med överföringsoperatorerna<br />

Gk(p). Utsignalen hos G kan härvid bestämmas som summan av de enskilda systemens<br />

utsignaler. Om polynomet A(p) saknar multipla nollställen kan partialbr˚aksuppdelningen<br />

(3.57) göras s˚a att överföringsoperatorerna Gk(p) är av första eller andra ordningen,<br />

<strong>och</strong> deras stegsvar kan s˚aledes beräknas med de formler som givits i detta kapitel.<br />

3.3.8 Instabila system<br />

Vi har ovan antagit att de studerade systemen är stabila, utan att kvantitativt ange när ett<br />

system är stabilt. Ett system är instabilt om det finns variabler som divergerar (dvs vars<br />

värden växer obegränsat) trots att insignalen är begränsad. För de system som diskuterats<br />

i detta kapitel innebär detta att stegsvaret divergerar. Ett allmänt villkor för att ett system<br />

skall vara stabilt kan härledas genom följande observationer:<br />

45<br />

k


• För ett system av första ordningen enligt i ekvation (3.24) divergerar stegsvaret (3.27)<br />

om T < 0, eller ekvivalent om a < 0. Stabiliteten beror s˚aledes av nollstället hos<br />

överföringsoperatorns nämnare i ekvation (3.26), dvs lösningen till<br />

p + a = 0 (3.58)<br />

Systemet (3.24) är allts˚a instabilt om lösningen till (3.58) är positiv.<br />

• För ett system med tv˚a tidskonstanter T1 <strong>och</strong> T2 divergerar stegsvaret om T1 < 0 eller<br />

T2 < 0, jfr ekvation (3.33), (3.34). Detta är ekvivalent med att nämnaren hos systemets<br />

överföringsoperator i (3.31) har ett positivt nollställe, dvs ekvationen<br />

har en positiv lösning.<br />

(T2p + 1)(T1p + 1) = 0<br />

• För ett underdämpat system som beskrivs av (3.37) divergerar stegsvaret (3.39) om<br />

ζωn < 0. Detta är ekvivalent att nämnaren hos systemets överföringsoperator (3.38)<br />

har ett nollställe med positiv realdel, ty ekvationen<br />

har lösningarna<br />

p 2 + 2ζωnp + ω 2 n = 0<br />

<br />

p1,2 = −ζωn ± (ζωn) 2 − ω2 n<br />

Dessa resultat kan generaliseras till en det allmänna fallet genom att observera att en<br />

överföringsoperator av typen (3.12) alltid kan beskrivas som en seriekoppling enligt<br />

G(p) = B(p)<br />

A(p) = G1(p) · G2(p) · . . . · GM(p)<br />

där faktorerna Gi(p) är öveföringsoperatorer av första eller andra ordningen, vars nollställen<br />

är rötterna till ekvationen<br />

A(p) = 0 (3.59)<br />

Systemet är stabilt om <strong>och</strong> endast om alla delsystem Gi(p) i seriekopplingen är stabila. Fr˚an<br />

resonemangen ovan följer s˚aledes att ett system med överföringsoperatorn G(p) = B(p)/A(p)<br />

är instabilt om nämnarpolynomet A(p) har en eller flera nollställen med positiv realdel.<br />

46


3.4 N˚agot om regulatordesign<br />

Som vi sett kan ett systems beteende p˚averkas med ˚aterkopplad reglering. En vanlig m˚alsättning<br />

är att kompensera för störningars inverkan samt att f˚a de reglerade signalerna i<br />

systemet att följa givna referensvärden. Förutom de redan behandlade problemen belyser<br />

följande exempel n˚agra ytterligare typiska problem där ˚aterkopplad reglering är viktig.<br />

Exempel 3.2 Stabilisering av instabila system.<br />

Betrakta ett instabilt system, t.ex.<br />

dy(t)<br />

dt<br />

+ ay(t) = bu(t)<br />

där a < 0. Det enda sättet att h˚alla utsignalen y(t) fr˚an ett instabilt system vid en önskad<br />

niv˚a är med hjälp av ˚aterkopplad reglering. Om systemet ovan regleras med regulatorn<br />

ges den slutna kretsen av<br />

u(t) = Kp (r(t) − y(t))<br />

dy(t)<br />

dt + (a + bKp) y(t) = bKpr(t)<br />

som är stabilt om Kp väljs s˚a att a + bKp > 0. Det finns tekniska system som kan göras<br />

effektivare om de konstrueras s˚a att de är instabila utan reglering. Flygplan är t.ex. instabila<br />

(i motsats till luftballonger eller luftskepp) <strong>och</strong> kräver ˚aterkopplad reglering.<br />

Exempel 3.3 Modifiering ett systems beteende för enklare manuell styrning.<br />

Ett system som annars skulle vara sv˚art eller omöjligt att styra manuellt kan genom ˚aterkoppling<br />

modifieras s˚a att manuell styrning är möjlig. Ett flygplan vars dynamik är instabil<br />

<strong>och</strong> för snabb för manuell styrning kan t.ex. först stabiliseras med ˚aterkopplad reglering.<br />

Den manuella styrningen ger d˚a börvärdet r(t) för regulatorn, medan reglersignalen u(t) till<br />

systemet bestäms av regulatorn, jfr figur 3.9.<br />

3.4.1 Begränsningar vid regulatordesign<br />

Metoder för design av regulatorer behandlas inte i denna kurs, utan här kommer vi bara att<br />

med ett exempel belysa en typisk problemställning vid regulatordesign. En viktig egenskap hos<br />

regulatordesign är att den alltid best˚ar av en kompromiss mellan ett antal olika m˚alsättningar.<br />

Speciellt har man ofta en kompromiss mellan reglerprestanda, s˚asom snabb respons, versus<br />

stabilitet hos det reglerade systemet. Detta belyses av följande exempel.<br />

47


Figur 3.12: Temperaturreglering med dödtid.<br />

Exempel 3.4 Stabilitet vid ˚aterkoppling.<br />

En vätska som strömmar i en rörledning skall uppvärmas genom tillförsel av ˚anga enligt figur<br />

3.12. Temperaturen hos den inkommande strömmen är Td(t), <strong>och</strong> den önskade temperaturen<br />

y(t) hos den utg˚aende strömmen är r.<br />

˚Angflödet regleras med hjälp av en reglerventil, vars läge u(t) bestäms av en regulator<br />

Gc. Det antas att temperaturökningen i vätskan är direkt proportionell mot u(t),<br />

y − Td = bu(t)<br />

där b är en konstant. Temperaturen y(t) mäts emellertid vid ett ställe efter ˚angtillförseln dit<br />

det tar L = 10 sekunder för vätskan att strömma (jfr figur 3.12). Sambandet mellan u(t) <strong>och</strong><br />

y(t) är därför<br />

y(t) = bu(t − L) + Td(t) (3.60)<br />

L˚at systemet för enkelhets skull regleras med en P-regulator, s˚a att<br />

u(t) = Kp (r − y(t)) (3.61)<br />

Vi kan ocks˚a anta att alla signaler är definierade som avvikelser fr˚an ett fortfarighetstillst˚and,<br />

s˚a att r = 0. Det reglerade systemet (3.60), (3.61) beskrivs d˚a av<br />

y(t) + bKpy(t − L) = Td(t) (3.62)<br />

48


Betrakta sedan vad som händer d˚a den inkommande strömmens temperatur ändras stegformat<br />

fr˚an 0 till 1 vid tiden t = 0. D˚a gäller<br />

y(t) = 0, t ≤ L<br />

y(t) = 1, L < t ≤ 2L<br />

y(t) = 1 − bKp, 2L < t ≤ 3L<br />

y(t) = 1 − bKp(1 − bKp), 3L < t ≤ 4L<br />

.<br />

y(t) = −bKpy(t − L) + 1<br />

Den sista ekvationen innebär att om systemet konvergerar, s˚a att y(t) = y(t − L), s˚a gäller<br />

y(t) →<br />

1<br />

1 + bKp<br />

d˚a t → ∞<br />

Regleravvikelsen kunde d˚a i princip göras godtyckligt liten genom att välja en stor regulatorförstärkning<br />

Kp. Denna analys beaktar emelllertid inte systemets stabilitet. Ur ovan följer<br />

nämligen ocks˚a<br />

y(t) = −bKpy(t − L) + 1<br />

.<br />

= −bKp (−bKpy(t − 2L) + 1) + 1 = (−bKp) 2 y(t − 2L) − bKp + 1<br />

= (−bKp) k y(t − kL) + (−bKp) k−1 + · · · 1<br />

vilket visar, att y(t) → ±∞ om |bKp| > 1. Regulatorförstärkningen Kp kan därför inte väljas<br />

hur stort som helst eftersom systemet d˚a blir instabilt.<br />

Exemplet visar att regulatordesign är en kompromiss mellan reglerprestanda (liten statisk<br />

regleravvikelse) <strong>och</strong> stabilitet hos den slutna kretsen. Den statiska regleravvikelsen kan givetvis<br />

elimineras med hjälp av en PI-regulator med integrerande verkan, men även d˚a kommer kravet<br />

p˚a stabilitet att begränsa den reglerprestanda som kan uppn˚as.<br />

Det finns systematiska metoder för att konstruera regulatorer som ger det ˚aterkopplade<br />

systemet önskade egenskaper med beaktande av stabilitetsbegränsningar. Dessa metoder behandlas<br />

i senare kurser i <strong>reglerteknik</strong>.<br />

49

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!