Zachodniopomorski Uniwersytet Technologiczny
Modelowanie Procesów Przemysłowych
Temat: Model głośnika elektrodynamicznego
Paweł Jan Nowak	Prowadzący: Dr Przemysław Orłowski

Kod z programu:

s=tf('s')

r=8

l=0.005

m=0.008

b=1/2

k=500

c=5

A=[-r/l 0 -b/l; 0 0 1; b/m -k/m -c/m]

B=[1/l; 0; 0]

i=eye(3)

zwe=([1 0 0]*(s*i-A)^(-1)*B)^(-1)

Transfer function:

s^3 + 2225 s^2 + 1.069e006 s + 1e008

------------------------------------

200 s^2 + 1.25e005 s + 1.25e007

bode(zwe) (1)

tw=[0 0 1]*(s*i-A)^(-1)*B

Transfer function:

-3.469e-016 s^3 - 7.648e-013 s^2 + 1.25e004 s - 3.425e-008

----------------------------------------------------------

s^3 + 2225 s^2 + 1.069e006 s + 1e008

bode(tw) (2)

C=i

D=[ 0; 0; 0]

t=C*(i*s-A)^(-1)*B+D

Transfer function from input to output...

200 s^2 + 1.25e005 s + 1.25e007

#1: ------------------------------------

s^3 + 2225 s^2 + 1.069e006 s + 1e008

3.903e-023 s^3 + 8.684e-020 s^2 - 4.74e-015 s + 1.25e004

#2: --------------------------------------------------------

s^3 + 2225 s^2 + 1.069e006 s + 1e008

-3.469e-016 s^3 - 7.648e-013 s^2 + 1.25e004 s - 3.425e-008

#3: ----------------------------------------------------------

s^3 + 2225 s^2 + 1.069e006 s + 1e008

impulse(t) (3)

step(t) (4)

1. Charakterystyka częstotliwościowa Bodego funkcji przejścia Zwe(s)

1. Charakterystyka częstotliwościowa Bodego funkcji przejścia Tw(s)

1. Charakterystyki impulsowe

1. Charakterystyki skokowe