WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   

Partner:




Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

Fertige Funktion für PT2-Glied?

 

fenderbender
Forum-Anfänger

Forum-Anfänger


Beiträge: 49
Anmeldedatum: 17.06.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 29.10.2010, 08:52     Titel: Fertige Funktion für PT2-Glied?
  Antworten mit Zitat      
Hallo Leute,

ich habe vor, einen statischen Datenwert durch eine Näherung zu ersetzen, die von Sprunghöhe und Zeit abhängig ist, also ein klassisches PT2-Glied.
Zum Verständnis folgender Code:

Code:

t=0:10^-15:250*10^-15; K=1; d=0.2; T=10*10^-15;
a = K - (K.exp(-d.*/T)/sqrt(1-d^2)).*sin((1-d^2).t/T+arctan(sqrt(1-d^2)/d))
b = sum(a)/length(a);
 

Falls Tippfehler drin sind bitte einfach ignorieren. Plottet man beides, sieht man eine Kurve mit Anstiegszeit und Einschwingen und b ergibt den Mittelwert, also eine konstante Linie etwas unterhalb von K.
Der Wert b soll also einen gemessenen Punkt (K ist der statisch gemessene Wert) ersetzen. Das Problem ist, es sind mehrere tausend solcher Punkte und das ganze in einer Doppelschleife (1. Schleife für alle Punkte, 2.Schleife für das PT2-Glied mit length(t) Werten) gibt elend lange Zeiten. Daher die Frage, gibt es eine fertige PT2-Funktion, die mir die Schleife über t=0:10^-15:250*10^-15; schneller macht?
Würde also z.B. nach so etwas suchen (von mir ausgedacht... gibts nicht wirklich)
Code:

b = sum(pt2(d,T,t,K))/length()
 


Jemand nen Tipp?

Danke
Private Nachricht senden Benutzer-Profile anzeigen


Andreas Goser
Forum-Meister

Forum-Meister


Beiträge: 3.654
Anmeldedatum: 04.12.08
Wohnort: Ismaning
Version: 1.0
     Beitrag Verfasst am: 29.10.2010, 09:00     Titel:
  Antworten mit Zitat      
Mir ist nicht klar ob die Control System Toolbox vorhanden ist.

Andreas
Private Nachricht senden Benutzer-Profile anzeigen E-Mail senden
 
Jan S
Moderator

Moderator


Beiträge: 11.057
Anmeldedatum: 08.07.10
Wohnort: Heidelberg
Version: 2009a, 2016b
     Beitrag Verfasst am: 29.10.2010, 10:50     Titel: Re: Fertige Funktion für PT2-Glied?
  Antworten mit Zitat      
Hallo fenderbender,

Dein Code-Beispiel ist wirklich etwas merkwürdig: "K.exp(-d.*/T)"??
Ich sehe auch keine Schleife über t.
Macht "T=10*10^-15;" wirklich genau das, was Du erwartest?!

Möchtest Du die Laufzeit verkürzen? Dann poste doch bitte eine lauffähige Version des Programms. Ich vermute, das würde sich deutlich beschleunigen lassen, wahrscheinlich auch noch vektorisieren. Ich kann aber auf keinen Fall Code "beschleunigen", der nicht der Matlab Syntax genügt. Es bringt ja wenig, nicht-laufen Code in schneller nicht-laufenden Code zu verwandeln...

Gruß, Jan
Private Nachricht senden Benutzer-Profile anzeigen
 
fenderbender
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 49
Anmeldedatum: 17.06.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 29.10.2010, 16:22     Titel:
  Antworten mit Zitat      
Wie gesagt, Tippfehler bitte ignorieren Smile

Gemeint ist die Formel für den aperiodischen Fall (d<1), siehe z.B. Wikipedia:
http://de.wikipedia.org/wiki/PT2-Glied

Die Formel hab ich richtig im Matlabprogramm, nur hier wohl vertippt. Jedefalls ist die Sprungantwort das was ich erwarte.

Ich glaube kaum, dass du mit dem ganzen code etwas anfangen kannst.
Private Nachricht senden Benutzer-Profile anzeigen
 
Jan S
Moderator

Moderator


Beiträge: 11.057
Anmeldedatum: 08.07.10
Wohnort: Heidelberg
Version: 2009a, 2016b
     Beitrag Verfasst am: 29.10.2010, 21:38     Titel:
  Antworten mit Zitat      
Hallo fenderbender,
Zitat:
Ich glaube kaum, dass du mit dem ganzen code etwas anfangen kannst.

Das wird wohl stimmen. Mit dem partiellen Code kann ich aber auch nichts anfangen.
Ich glaube, die Funktion liesse sich deutlich beschleunigen. Ausprobieren kann ich es aber natürlich nicht, da ich ja keine lauffähige Version der Funktion kenne.

Viel Erfolg, Jan
Private Nachricht senden Benutzer-Profile anzeigen
 
Thomas84
Forum-Meister

Forum-Meister


Beiträge: 546
Anmeldedatum: 10.02.10
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 01.11.2010, 07:46     Titel:
  Antworten mit Zitat      
Das ist ja eine lineare inhomogene Differentialgleichung. Die kannst du mit lsim lösen (vorher natürlich in DGL-Sytem erster Ordnung umwandeln)

viele Grüße
Thomas
Private Nachricht senden Benutzer-Profile anzeigen
 
Andreas Goser
Forum-Meister

Forum-Meister


Beiträge: 3.654
Anmeldedatum: 04.12.08
Wohnort: Ismaning
Version: 1.0
     Beitrag Verfasst am: 02.11.2010, 10:11     Titel:
  Antworten mit Zitat      
Nur für den Fall, dass sich jemand interessiert warum ich nach der Control System Toolbox gefragt habe Very Happy Dort wäre alles schon implementiert.

Andreas
Private Nachricht senden Benutzer-Profile anzeigen E-Mail senden
 
Neues Thema eröffnen Neue Antwort erstellen



Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:

Du kannst Beiträge in dieses Forum schreiben.
Du kannst auf Beiträge in diesem Forum antworten.
Du kannst deine Beiträge in diesem Forum nicht bearbeiten.
Du kannst deine Beiträge in diesem Forum nicht löschen.
Du kannst an Umfragen in diesem Forum nicht mitmachen.
Du kannst Dateien in diesem Forum posten
Du kannst Dateien in diesem Forum herunterladen
.





 Impressum  | Nutzungsbedingungen  | Datenschutz | FAQ | goMatlab RSS Button RSS

Hosted by:


Copyright © 2007 - 2025 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.