|
Megatron |
Forum-Newbie
|
|
Beiträge: 5
|
|
|
|
Anmeldedatum: 06.04.08
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 06.04.2008, 11:19
Titel: Funktion aus Messwerten
|
|
|
|
|
Hallo, ich bin ein Matlab-Neuling und stehe vor folgendem Problem:
Ich habe in Diagram mit empirisch ermittelten Werten eines Wirkungsgrades in Abhängigkeit der Geschwindigkeit.
w0=[0 2 4.4 6.5 8.6 10.8 12.8 15.2 17.2 19.2 21.6];
Eta=[0 0.15 0.3 0.38 0.415 0.485 0.56 0.62 0.66 0.69 0.7];
Ich habe die Wertepaare dann in zwei Spaltenvektoren eingetragen und sie mit der spline Funktion interpoliert.
w=0:0.1:22;
EtaP=spline(w0,Eta,w);
Wenn ich den Graphen plotte sieht er auch genau so aus wie das Diagramm. Das ist schon mal super.
Was ich aber jetzt eigentlich suche ist eine Funktion Eta(w0) bzw. EtaP(w) die ich dann in andere Gleichung für Eta einsetzten kann.
Habe ich überhaupt die richtige Herangehensweise? Bzw. ist das so überhaupt möglich?
Und wo ich schon mal dabei bin, hab ich noch eine wahrscheinlich ziemlich simple Frage. Ich habe eine Funktion bzw. Gleichung: zur Vereinfachung z.B.
P=2*H
Gibt es einen Befehl zur Umstellung nach H? Und wie lasse ich mir einen Wert für H(P=2) ausgeben?
Ich hoffe es kann mir jemand helfen, dass Problem brennt mir schon das ganze Wochenende unter den Fingern.
Vielen Danke!!
|
|
|
|
|
nschlange |
Ehrenmitglied
|
|
Beiträge: 1.320
|
|
|
|
Anmeldedatum: 06.09.07
|
|
|
|
Wohnort: NRW
|
|
|
|
Version: R2007b
|
|
|
|
|
|
Verfasst am: 06.04.2008, 11:40
Titel:
|
|
Hi,
wenn Du die Daten plottest (plot(w0,Eta)) kannst Du Tools->Basic fitting versuchen, einfache Funktionen an die Daten anzufitten. Welche Klasse von Funktionen allerdings dafür geeignet ist bestimmt der physikalische Zusammenhang. Wenn Du die Curve Fitting Toolbox hast kannst Du auch eigene Funktionen angeben. Oder Du durchsuchst meine Postings, da habe ich mal die Adresse zu einer freien Fitting Toolbox angegeben.
Da ein Wirkungsgrad wohl nicht beliebig steigen kann würde ich vielleicht eine Logarithmus- oder Sättigungsfunktion nehmen.
_________________
Viele Grüße
nschlange
"Chuck Norris ejakuliert fluessigen Stahl!"
|
|
|
Megatron |
Themenstarter
Forum-Newbie
|
|
Beiträge: 5
|
|
|
|
Anmeldedatum: 06.04.08
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 06.04.2008, 15:38
Titel:
|
|
Wow, das funktioniert echt gut, ich habe ein Polynom 7.Ordnung gewählt und das passt ziemlich genau.
Vielen Dank!!!
Könntest du mir denn auch noch sagen wie ich mit Hilfe von Matlab aus diesen beiden Gleichungen
H=[0:0.1:11]
roh=1.225*(1-0.02256*H).^4.2559
Pvp=FG*sqrt((2*cw^2*sigma)./(ca^3*roh));
H von Pvp=z.B.4000 herausbekommen kann?
Vielen Dank!!!!
|
|
|
nschlange |
Ehrenmitglied
|
|
Beiträge: 1.320
|
|
|
|
Anmeldedatum: 06.09.07
|
|
|
|
Wohnort: NRW
|
|
|
|
Version: R2007b
|
|
|
|
|
|
Verfasst am: 06.04.2008, 20:27
Titel:
|
|
Das mit Polynom 7. Ordnung würde ich mir noch mal überlegen...das könnte Dir als 'unphysikalisch' um die Ohren gehauen werden
Vielleicht wäre z.B. folgendes Modell passender:
f(x) = a*(1-exp(-b*x)) mit a = 0.8223 und b = 0.0905
Dann hätte man auch ein Grenze des Wirkungsgrades bei ca. 0.8 und nicht wie bei einem Polynom +- inf.
Zur zweiten Frage:
Du willst ein Pvp finden, so dass H(Pvp) = 4000 ist?
Könnte man so machen. Man sucht den Index des Vektors, deren Eintrag die Bedingung erfüllt. Eine Toleranz muss wegen der Diskretheit der Werte eingeführt werden.
_________________
Viele Grüße
nschlange
"Chuck Norris ejakuliert fluessigen Stahl!"
|
|
|
Megatron |
Themenstarter
Forum-Newbie
|
|
Beiträge: 5
|
|
|
|
Anmeldedatum: 06.04.08
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 10.04.2008, 21:04
Titel:
|
|
Ist ja abgefahren das funktioniert ja wirklich!!!!
Ich hab etwas gebraucht um zu verstehen was ich da genau suche und welchen Wert Matlab mir dann ausgibt aber jetzt hab ichs verstanden.
Danke für die Hilfe!!!
Kann ich diese Suchfunktion auch nutzen um den Schnittpunkt von 2 Graphen heraus zu bekommen?
Ich hab es mal mit nem solve Befehl versucht, bin aber gescheitert.
|
|
|
Megatron |
Themenstarter
Forum-Newbie
|
|
Beiträge: 5
|
|
|
|
Anmeldedatum: 06.04.08
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 10.04.2008, 21:56
Titel:
|
|
Wow ist das abgefahren,
ich hab es gerade mal ausprobiert mit dem find-Befehl, den Schnittpunkt der beiden Graphen zu finden, und es hat echt funktioniert.
Vielen Dank noch mal!!!!
Hab echt was dazugelernt
Edit von Bijick: Code-Formatierung eingefügt. Bitte beim nächsten Mal selbst machen!
|
|
|
|
|
Einstellungen und Berechtigungen
|
|
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
| RSS
Hosted by:
Copyright © 2007 - 2024
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.
|
|