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

Fehlerfortpflanzung und lineare interpolationsgleichung

 

Pigmy
Forum-Anfänger

Forum-Anfänger


Beiträge: 10
Anmeldedatum: 22.03.16
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 26.06.2019, 15:12     Titel: Fehlerfortpflanzung und lineare interpolationsgleichung
  Antworten mit Zitat      
Hallo zusammen,

ich habe eine Gleichung:
C = E*Kf/100.
E nimmt dabei konstante fehlerbehaftete Werte an.
Zu dieser Gleichung würde ich jetzt gerne die Fehler abschätzen.
Die Fehlerfortpflanzungsformel sieht ja so aus:
\sigma_c = sqrt{\left(\frac{\partial C}{\partial E}\right)^2*\sigma_E^2 + \left(\frac{\partial C}{\partial Kf}\right)^2*\sigma_{Kf}^2}
Die Berechnung von Kf ist etwas aufwändiger. Hierzu habe ich folgende Stützstellen gegeben:
Code:
xs = [0.7; 0.8; 0.9; 1.0; 1.1; 1.2; 1.3; 1.4; 1.5];
ys = [1.29; 1.16; 1.07; 1.00; 0.94; 0.89; 0.85; 0.81; 0.79];[

Für beliebige xi Werte zwischen den jeweilige Intervallen soll einfach linear interpoliert werden mit
yi = y1 + (y2 - y1) / (x2 - x1) * xi - x1)
also z.B. für xi = 1.45
Code:
x1 = 1.4;
x2 = 1.5;
y1 = 0.81;
y2 = 0.79;


Wie stelle ich das jetzt an wenn ich keine Funktionsgleichung für Kf habe?
Gibt es da eine Möglichkeit?

Ich habe direkt daran gedacht, mir mit polyfit eine angenäherte Polynomfunktion zu berechnen und diese in die Fehlerfortpflanzungsformel einzubauen.
Code:

n = 7;
p = polyfit(xs,ys,n);

Erhalte damit auch ein R² = 0.99999945...
Wenn ich mir jetzt die Gleichung mit den Polynomen aufstelle oder polyval verwende und z.B. an der Stelle xi = 1 auswerte
yi= polyval(p,1); erhalte ich ein yi = 1.000174... und nicht von 1. Dies ist natürlich logisch aber die Genauigkeit reicht mir nicht aus.
Jetzt könnte ich natürlich einfach immer mehr polynome nehmen, allerdings muss ich diese Gleichung ja auch noch ableiten.

Habt Ihr eine Idee?
Danke für eure Unterstützung
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: 26.06.2019, 16:13     Titel: Re: Fehlerfortpflanzung und lineare interpolationsgleichung
  Antworten mit Zitat      
Hallo Pigmy,

Ich zumindest kann dir nicht folgen.

Zitat:
Jetzt könnte ich natürlich einfach immer mehr polynome nehmen, allerdings muss ich diese Gleichung ja auch noch ableiten.

Immer mehr Polynome?

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

Forum-Anfänger

Forum-Anfänger


Beiträge: 10
Anmeldedatum: 22.03.16
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 26.06.2019, 16:27     Titel:
  Antworten mit Zitat      
Vielleicht schlecht ausgedrückt.
Ich könnte ja einfach die polyfit Funktion mit einem höheren Grad verwenden. Dadurch erhalte ich ja auch mehr Polynome (n+1). Dadurch kann ja oftmals die Genauigkeit erhöht werden. Aber ich bekomme eben auch einen zusätzlichen Term.
Aber ich habe die Ableitungen jetzt so weit automatisiert, dass das keine Arbeit mehr macht.
Trotzdem würde ich gerne wissen, wie man in diesem Fall eine genauere Funktionsvorschrift erhält.
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: 26.06.2019, 20:37     Titel:
  Antworten mit Zitat      
Hallo Pigmy,

Zitat:
Ich könnte ja einfach die polyfit Funktion mit einem höheren Grad verwenden. Dadurch erhalte ich ja auch mehr Polynome (n+1).

Nein, polyfit berechnet immer nur ein Polynom. Meinst du: "Dann hat das Polynom mehr Koeffizienten" bzw. "einen höheren Grad"?

Zitat:
Dadurch kann ja oftmals die Genauigkeit erhöht werden.

Das gilt nur in einem sehr begrenzten Bereich. Eine Interpolation mit hohen Polynomengraden neigt zu gigantischem Schwingen in den Randbereichen und ist im Allgemeinen nicht hilfreich.

Zitat:
Trotzdem würde ich gerne wissen, wie man in diesem Fall eine genauere Funktionsvorschrift erhält.

Eine Funktionsvorschrift wofür?

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

Forum-Anfänger

Forum-Anfänger


Beiträge: 10
Anmeldedatum: 22.03.16
Wohnort: ---
Version: ---
     Beitrag Verfasst am: 27.06.2019, 07:48     Titel:
  Antworten mit Zitat      
Ok, ja, mehr Koeffizienten hat das Polynom dann, das meinte ich.

Eine Funktionsvorschrift für Kf brauche ich, weil Kf wiederum von E und einer weiteren Variable, RE, abhängt. Um dann die Fehlerabschätzung machen zu können, brauche ich doch ne Gleichung für Kf, oder nicht?

Vielleicht probiere ich es noch einmal:

Ich möchte den Fehler der Formel C = E*GW*Kf/100 abschätzen. GW hat keinen Fehler. Der Fehler von E ist mir bekannt. Der Fehler von Kf ist mir nicht bekannt. Kf errechnet sich aus den xs und den dazugehörigen ys Werten, bzw. wird dazwischen linear interpoliert. Die xs Werte wiederum ergeben sich aus xs = E/RE. Der Fehler für E ist mir immer noch bekannt und der Fehler für RE auch.

Ist es jetzt vielleicht klarer?
Was ich jetzt gemacht habe ist mir die Polynomsgleichung für Kf mit den Koeffizienten durch polyfit berechnen lassen. Die Koeffizenten kann ich dann ja in die allgemeine Form Kf = p_1x^n+p_2x^{n-1} + ... + p_{n+1}x^0 einsetzen und ich bekomme meine angenäherte Gleichung für Kf.
Diese kann ich ja dann in meine Ausgangsgleichung für E einsetzen und die Fehlerabschätzung durchführen.
Private Nachricht senden Benutzer-Profile anzeigen
 
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 - 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.