Mein MATLAB Forum - goMatlab.de

Mein MATLAB Forum

 
Gast > Registrieren       Autologin?   
Bücher:

Studierende:
Praktikum - Vergleich von MATLAB und Scilab
Branche: mehrere
IAV GmbH - Berlin, Chemnitz

Abschlussarbeit / Praktikum: Entwicklung Matlab (m/w)
Branche: Informationstechnologie, Elektrotechnik, Elektronik
GIGATRONIK Technologies GmbH - Ulm

Praktikum im Bereich Vorausentwicklung - Aktoren sicherheitsrelevanter Anwendungen
Branche: Fahrzeugtechnik, Fahrzeugteile, Maschinentechnik, Industrietechnik
MAHLE International GmbH - Stuttgart

Masterarbeit - Optimierung von Online erstellten 3D-Modellen
Branche: Luft-, Raumfahrttechnik, Wehrtechnik
MBDA Deutschland - Schrobenhausen

Praktikum / Abschlussarbeit im Bereich Assistenz- und Regelsysteme
Branche: Fahrzeugtechnik, Fahrzeugteile, Luft-, Raumfahrttechnik
Bertrandt Technologie GmbH - Mönsheim bei Weissach

weitere Angebote

Partner:


Vermarktungspartner


Forum
      Option
[Erweitert]
  • Diese Seite per Mail weiterempfehlen
     


Gehe zu:  
Neues Thema eröffnen Neue Antwort erstellen

2 Dimensionale Interpolation

 

Avo
Forum-Anfänger

Forum-Anfänger


Beiträge: 20
Anmeldedatum: 04.12.17
Wohnort: ---
Version: 2016b
     Beitrag Verfasst am: 03.01.2018, 16:39     Titel: 2 Dimensionale Interpolation
  Antworten mit Zitat      
Hallo,

ich habe drei Matrizen gleicher Größe mit unterschiedlich Messwerten, bei denen jede Zelle als ein Messpunkt zu verstehen ist. Konkret: Es gibt die Matrizen H, Q, und P und H(1,1), Q(1,1) und P(1,1) gehören zusammen, da sie zum selben Zeitpunkt aufgenommen wurden. Plottet man jetzt H und P über Q erhält man zwei Kennfelder.

Ich möchte nun einen H und einen Q Wert vorgeben, ohne dass diese gemessen wurden und ein passendes P mithilfe der gemessenen Daten schätzen.

mein bisheriger Code sieht wie folgt aus:
Code:
load('PHQ_Data.mat');
[temp, row] = max(Q.*H);
[~, col] = max(temp);
row = row(col);
q100 = Q(row,col);
h100 = H(row,col);
h0 = h100*0.5;
q75 = q100*0.75;
h75 = (h0)/q100*q75+h0;

rowNo = findRow(Q,3,q75,1);
vectorAtRowNo = interp1([floor(rowNo), floor(rowNo)+1],H(floor(rowNo):floor(rowNo)+1,1:length(Q(1,:))),rowNo);
columnNo = findColumn(vectorAtRowNo,1,h75,1);
temp_oben = interp1([floor(columnNo),floor(columnNo)+1],[P(floor(rowNo),(floor(columnNo))),P(floor(rowNo),(floor(columnNo)+1))],columnNo);
temp_unten = interp1([floor(columnNo),floor(columnNo)+1],[P(floor(rowNo)+1,(floor(columnNo))),P(floor(rowNo)+1,(floor(columnNo)+1))],columnNo);
temp_value = interp1([floor(rowNo),floor(rowNo)+1],[temp_oben temp_unten], rowNo);
P_value = temp_value;


Er basiert dabei stark auf eine vorherige Lösung einer ähnlichen Problemstellung, nur das dort die erste Referenz-Matrix, in jeder Spalte konstant war und somit in findRow nur eine Spalte (z.B. hier die 3) durchsucht werden musste, hier sind die Q-Werte zwar ähnlich, aber nicht konstant. Wie komme ich an eine sinnvolle start-Spalte?

Wäre auch für komplett andere Lösungsvorschläge offen.


Besten Gruß,
Avo

findColumn.m
 Beschreibung:
Code für findRow function

Download
 Dateiname:  findColumn.m
 Dateigröße:  1.01 KB
 Heruntergeladen:  34 mal
findRow.m
 Beschreibung:
Code für findColumn functionr

Download
 Dateiname:  findRow.m
 Dateigröße:  1.03 KB
 Heruntergeladen:  32 mal
PHQ_Data.mat
 Beschreibung:
Test Datensatz

Download
 Dateiname:  PHQ_Data.mat
 Dateigröße:  4.7 KB
 Heruntergeladen:  33 mal
PQ_Kennfeld.fig
 Beschreibung:

Download
 Dateiname:  PQ_Kennfeld.fig
 Dateigröße:  26.7 KB
 Heruntergeladen:  33 mal
HQ_Kennfeld.fig
 Beschreibung:
der mittlere Punkt der geplotteten Gerade soll geschätzt werden

Download
 Dateiname:  HQ_Kennfeld.fig
 Dateigröße:  25.05 KB
 Heruntergeladen:  30 mal
Private Nachricht senden Benutzer-Profile anzeigen


Harald
Forum-Meister

Forum-Meister


Beiträge: 18.215
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2014a
     Beitrag Verfasst am: 03.01.2018, 17:33     Titel:
  Antworten mit Zitat      
Hallo,

das klingt eigentlich nach
Code:
Pneu = interp2(H, Q, P, Hneu, Qneu)


Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
Avo
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 20
Anmeldedatum: 04.12.17
Wohnort: ---
Version: 2016b
     Beitrag Verfasst am: 03.01.2018, 17:49     Titel:
  Antworten mit Zitat      
Da sagt er mir: Input grid is not a valid MESHGRID.

Was macht denn ein gültiges Meshgrid aus? Von dem was ich glaube verstanden zu haben, müssen die Vektoren mit denen es erzeugt wird (teilweise?) monoton und mit gleichbleibenden Abständen sein.
Ich erzeuge hier zwar das Grid nicht aus Vektoren, aber gleichbleibende Abstände habe ich schonmal nicht und wenn man sich die Kennfelder anguckt, sieht man das teilweise mehrere y-Werte einem x-Wert zugeordnet sind.
Da ich davon ausgehen muss, dass die Daten in einer solchen Form zu Verarbeitung ankommen: Irgendeine Idee wie ich sie in gültige Meshgrid Form bringen kann, ohne sie signifikant zu ändern?
Private Nachricht senden Benutzer-Profile anzeigen
 
Harald
Forum-Meister

Forum-Meister


Beiträge: 18.215
Anmeldedatum: 26.03.09
Wohnort: Nähe München
Version: ab 2014a
     Beitrag Verfasst am: 03.01.2018, 18:20     Titel:
  Antworten mit Zitat      
Hallo,

ich hab das zwar jetzt ebenso wenig ausprobiert, aber versuch mal griddata statt interp2.

Grüße,
Harald
Private Nachricht senden Benutzer-Profile anzeigen
 
Avo
Themenstarter

Forum-Anfänger

Forum-Anfänger


Beiträge: 20
Anmeldedatum: 04.12.17
Wohnort: ---
Version: 2016b
     Beitrag Verfasst am: 03.01.2018, 18:35     Titel:
  Antworten mit Zitat      
Hallo,

das scheint zu funktionieren, der herauskommende Wert wirkt zunächst realistisch, muss sich aber noch bei weiteren Tests beweisen.
Danke!

Besten Gruß,
Avo

PS: hatte in der Zwischenzeit die Daten nach Q sortiert, da ich auf Q abgebildet hatte und somit das Problem mit den doppelten y-Werten 'behoben'. Scheint aber für griddata keinen Unterschied zu machen.
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
.


goMatlab ist ein Teil des goForen-Labels
goForen.de goMATLAB.de goLaTeX.de goPCB.de


 Impressum  | Nutzungsbedingungen  | Datenschutz  | Werbung/Mediadaten | Studentenversion | FAQ | goMatlab RSS Button RSS


Copyright © 2007 - 2018 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks
Partner: LabVIEWforum.de

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.