CAPÍTULO 3 DINÂMICA DOS SISTEMAS BIOLÓGICOS E … · algébricas FCTUC/EngBiomedica/ADC/2007 8...

Post on 08-Nov-2020

0 views 0 download

Transcript of CAPÍTULO 3 DINÂMICA DOS SISTEMAS BIOLÓGICOS E … · algébricas FCTUC/EngBiomedica/ADC/2007 8...

FCTUC/EngBiomedica/ADC/2007 1

CAPÍTULO 3

DINÂMICA DOS SISTEMAS BIOLÓGICOS E FISIOLÓGICOS

FCTUC/EngBiomedica/ADC/2007 2

Objectivos

Modelo: equações diferenciais• Resolução de equações diferencias ?• Transformada de Laplace !

Função de Transferência• Plano complexo• Estabilidade• Resposta transitória e estacionária• Sistemas de primeira e segunda ordem

FCTUC/EngBiomedica/ADC/2007 3

Indice3.1 Modelização: equações diferenciais3.2 Transformada de Laplace

• Transformada de funções elementares• Propriedades da TLP

3.3 Inversão da TLP• Fracções parciais• Resolução equações diferenciais• Polinómio característico• Estabilidade reposta natural

3.4 Função de Transferência• Estabilidade reposta forçada• Pano complexo• Polos: resposta transitória• Zeros

3.5 Resposta temporal de sistemas lineares• Sistemas de primeira ordem• Sistemas de segunda ordem• Sistemas de ordem superior

FCTUC/EngBiomedica/ADC/2007 4

Indice3.1 Modelização: equações diferenciais3.2 Transformada de Laplace

• Transformada de funções elementares• Propriedades da TLP

3.3 Inversão da TLP• Polos e zeros• Polinómio característico• Plano complexo• Estabilidade resposta natural

3.4 Função de Transferência• Estabilidade resposta forçada• Plano complexo• Polos• Zeros

3.5 Resposta temporal de sistemas lineares• Sistemas de primeira ordem• Sistemas de segunda ordem• Sistemas de ordem superior

FCTUC/EngBiomedica/ADC/2007 5

Do capítulo 2 …

….. como ???

FCTUC/EngBiomedica/ADC/2007 6

Indice3.1 Modelização: equações diferenciais

3.2 Transformada de Laplace• Transformada de funções elementares• Propriedades da TLP

3.3 Inversão da TLP• Polos e zeros• Polinómio característico• Plano complexo

FCTUC/EngBiomedica/ADC/2007 7

O problema ?

Resolução de uma equação diferencial

Resolução de uma equação algébrica

Transformada de Laplace

domínio s

domínio t

Solução no domínio s

domínio sManipulações

Transformada Inversa

A Solução !

Solução no domínio temporal

domínio tMétodo normal

algébricas

FCTUC/EngBiomedica/ADC/2007 8

Indice3.1 Modelização: equações diferenciais

3.2 Transformada de Laplace• Transformada de funções elementares• Propriedades da TLP

3.3 Inversão da TLP• Fracções parciais• Resolução equações diferenciais• Polinómio característico• Estabilidade

FCTUC/EngBiomedica/ADC/2007 9

Transforma Laplace

( )f t ( )F sL

∫∞

−==0

)()()]([ dtetfsFtfL st

FCTUC/EngBiomedica/ADC/2007 10

FCTUC/EngBiomedica/ADC/2007 11

FCTUC/EngBiomedica/ADC/2007 12

FCTUC/EngBiomedica/ADC/2007 13

FCTUC/EngBiomedica/ADC/2007 14

Algumas transformadas

FCTUC/EngBiomedica/ADC/2007 15

Indice3.1 Modelização: equações diferenciais

3.2 Transformada de Laplace• Transformada de funções elementares

• Propriedades da TLP• Linearidade• Derivação • Integração• Deslocamento• Teorema Valor inicial e final

3.3 Inversão da TLP• Fracções parciais• Resolução equações diferenciais• Polinómio característico• Estabilidade

FCTUC/EngBiomedica/ADC/2007 16

Linearidade

FCTUC/EngBiomedica/ADC/2007 17

Derivação

FCTUC/EngBiomedica/ADC/2007 18

Derivação

FCTUC/EngBiomedica/ADC/2007 19

Integração

FCTUC/EngBiomedica/ADC/2007 20

Integração

FCTUC/EngBiomedica/ADC/2007 21

Deslocamento

FCTUC/EngBiomedica/ADC/2007 22

Deslocamento temporal

FCTUC/EngBiomedica/ADC/2007 23

Teoremas valor inicial e final

FCTUC/EngBiomedica/ADC/2007 24

Utilidade da Transformada?

Se a entrada (força) for igual a f=5 qual a saída (posição) ?

Resolução de uma equação diferencial !

FCTUC/EngBiomedica/ADC/2007 25

Utilidade da Transformada?

FCTUC/EngBiomedica/ADC/2007 26

Utilidade da Transformada?

FCTUC/EngBiomedica/ADC/2007 27

Utilidade da Transformada?

Está calculada a saída !

FCTUC/EngBiomedica/ADC/2007 28

Algumas transformadas

FCTUC/EngBiomedica/ADC/2007 29

ConcluindoConcluindo, a determinação da saída torna-se num problema de fácil resolução • a resolução de uma eq. diferencial transforma-se na

resolução de uma eq. algébrica

No entanto obteve-se a solução em Y(s)Em termos práticos apenas nos interessa y(t)

Será possível usando Y(s), calcular y(t) ?

Para tal recorre-se à inversa da transformada de Laplace.

FCTUC/EngBiomedica/ADC/2007 30

O problema ?

Resolução de uma equação diferencial

Resolução de uma equação algébrica

Transformada de Laplace

domínio s

domínio t

Solução no domínio s

domínio sManipulações

Transformada Inversa

A Solução !

Solução no domínio temporal

domínio tMétodo normal

algébricas

FCTUC/EngBiomedica/ADC/2007 31

Indice3.2 Transformada de Laplace• Transformada de funções elementares• Propriedades da TLP

3.3 Inversão da TLP• Resolução equações diferenciais• Fracções parciais• Polinómio característico• Estabilidade

FCTUC/EngBiomedica/ADC/2007 32

Indice3.2 Transformada de Laplace• Transformada de funções elementares• Propriedades da TLP

3.3 Inversão da TLP• Resolução equações diferenciais• Fracções parciais• Polinómio característico• Estabilidade

FCTUC/EngBiomedica/ADC/2007 33

Resolução equações diferenciais

FCTUC/EngBiomedica/ADC/2007 34

Resolução equações diferenciais

FCTUC/EngBiomedica/ADC/2007 35

Resolução equações diferenciais

FCTUC/EngBiomedica/ADC/2007 36

Y =Yzi + Yzs

Efeito das condições iniciais

Efeito da entrada externa

Resolução equações diferenciais

Função de Transferência

FCTUC/EngBiomedica/ADC/2007 37

Indice3.2 Transformada de Laplace• Transformada de funções elementares• Propriedades da TLP

3.3 Inversão da TLP• Resolução equações diferenciais

• Fracções parciais• Polinómio característico• Estabilidade

FCTUC/EngBiomedica/ADC/2007 38

FCTUC/EngBiomedica/ADC/2007 39

Problema da inversa resolvido !Uso de transformada inversa !

FCTUC/EngBiomedica/ADC/2007 40

Problema resume-se a decomposição em fracções parciais

Cálculo dos residuos (raizes denominador)

• Raizes simples

• Raizes Multiplas

• Raizes complexas conjugadas

Resolução equações diferenciais

FCTUC/EngBiomedica/ADC/2007 41

Resolução equações diferenciais

Inversa

FCTUC/EngBiomedica/ADC/2007 42

1

t

Resolução equações diferenciais

FCTUC/EngBiomedica/ADC/2007 43

1

t

Raizes Simples

FCTUC/EngBiomedica/ADC/2007 44

1

t

Raizes Simples

FCTUC/EngBiomedica/ADC/2007 45

1

t

Raizes de multiplicidade r

FCTUC/EngBiomedica/ADC/2007 46

1

t

Raizes de multiplicidade r

FCTUC/EngBiomedica/ADC/2007 47

Raizes de multiplicidade r

FCTUC/EngBiomedica/ADC/2007 48

Raizes complexas

FCTUC/EngBiomedica/ADC/2007 49

Raizes complexas

FCTUC/EngBiomedica/ADC/2007 50

Raizes complexas

FCTUC/EngBiomedica/ADC/2007 51

Matlab: residuesExemplo de Scripts em MatLab

%---------- factores simples - estritamente propria-----num=1;den=[1 3 2]; %---- den=(s+1)(s+2)[A,P,Ao]=residue(num,den);

%---------------------------------- output ------A=[ -1 1 ]P=[ -2 -1 ]Ao=[ ]

FCTUC/EngBiomedica/ADC/2007 52

Matlab: residues

%--------- factores simples - bipropria ----------------num=[1 0 -9];den=[1 0 -1]; %---- den=(s+1)(s-1)[A,P,Ao]=residue(num,den);

%---------------------------------- output ------A=[ 4 -4]P=[ -1 1]Ao=1

FCTUC/EngBiomedica/ADC/2007 53

Matlab: residues

%---------- factor de multiplicidade 3 --------------num=1;den=[1 3 3 1 ]; %---- den=(s+1)3

[A,P,Ao]=residue(num,den);

%---------------------------------- output ------A=[ 1 -2 -8]P=[ -1 -1 -1]Ao=[ ]

FCTUC/EngBiomedica/ADC/2007 54

Indice3.2 Transformada de Laplace• Transformada de funções elementares• Propriedades da TLP

3.3 Inversão da TLP• Fracções parciais• Resolução equações diferenciais

• Polinómio característico • Estabilidade da resposta natural

FCTUC/EngBiomedica/ADC/2007 55

Polinómio característico e modos do sistema

Se entrada nula, Y = Yzi, resposta não-forçada ou resposta natural

Polinómio característico

Modos (característicos)

FCTUC/EngBiomedica/ADC/2007 56

A resposta natural é uma combinação linear dos modos

Os coeficientes da combinação linear dependem das condições iniciais

FCTUC/EngBiomedica/ADC/2007 57

Indice3.2 Transformada de Laplace• Transformada de funções elementares• Propriedades da TLP

3.3 Inversão da TLP• Fracções parciais• Resolução equações diferenciais• Polinómio característico

• Estabilidade resposta natural

FCTUC/EngBiomedica/ADC/2007 58

Estabilidade (resposta natural)

Modos (característicos)

Se um dos modos tender para infinito, a resposta natural tende para infinito, o sistema diz-se instável

Condição necessária e suficiente de estabilidade:

• Se as raízes do polinómio forem todas menores que zero

• Se as raízes do polinómio forem menores que zero excepto uma (e só uma) que pode ser zero