[kobiaszu]

Mircy, mierzę sobie dwie wartości wyjściowe z mojego modelu w simulinku. W jaki sposób mogę zmierzyć tłumienie i częstotliwość? Mogę wyoprowadzić sobie przebiegi bloczkiem "to workspace" i tam ogarnąc to jakąś funkcją?

#matlab #automatyka

[KapitanDystrybuanta]

@kobiaszu: Możesz to aproksymować A*exp(-beta*t)*sin(phi + omega*t)?
To wtedy mógłbyś zminimalizować wskaźnik Q(beta,phi, omega) = (A*exp(-beta*t)*sin(phi + omega*t) - pomiar).^2
Służy do tego funkcja lsqnonlin.

[KapitanDystrybuanta]

@kobiaszu: Coś mi się znaki mnożenia usunęły, ale widać o co chodzi :)

[kobiaszu]

@KapitanDystrybuanta: Nie mam zielone pojęcia jak zrobić to o czym mówisz :D
Mógłbyś na to rzucić okiem? Interesuje mnie pomiar tłumienia i częstotliwości dla U,V oraz phi i theta.
http://speedy.sh/xAsJ8/strojenie.m

[KapitanDystrybuanta]

@kobiaszu: Teraz nie mam dostępu do kompa, ale za jakąś godzinę mógłbym zerknąć.

[kobiaszu]

@KapitanDystrybuanta: spoko, dzięki za pomoc bo już nie ogarniam tego.

[KapitanDystrybuanta]

@kobiaszu: Trochę ciężko mi się połapać w skrypcie, a pisałeś coś o pomiarach z Simulinka.
Podsyłam skrypt do identyfikacji parametrów przykładowego pomiaru i funkcję celu.
Skoro pisałeś o tłumieniu i częstotliwości, to zakładam, że przebiegi mają postać jak opisałem wyżej.

http://speedy.sh/kRzjr/main.m
http://speedy.sh/rHqMx/cel.m

[kobiaszu]

Znaczy pomiar mam już zapisane jako timeseries do workspace i nie wiem jak z tego wyciągnąć tłumienie i częstotliwość. Pod parametr podaję nazwę, a pod t i pomiar mój timeseries?

[KapitanDystrybuanta]

@kobiaszu: Jeśli przebiegi mają taką postać jak u mnie to wystarczy tylko, że za t i pomiar z main.m wstawisz dane z Twojego timeseries i powinno działać