Springen naar inhoud

[numeriek]interpolatie


  • Log in om te kunnen reageren

#1

jhnbk

    jhnbk


  • >5k berichten
  • 6905 berichten
  • VIP

Geplaatst op 02 december 2008 - 15:22

Ik heb twee vectoren, w en d, met eental 6 punten in. Bestaat er een mogelijkheid om met octave de coŽfficiŽnten te weten te komen van de interpolatie (cubic spline) voor de twee laatste deeltjes?
Ik zou namelijk het maximum van de punten moeten interpoleren en dat zou zich ergens in die delen bevinden.
Het vel van de beer kunnen verkopen vraagt moeite tenzij deze dood voor je neervalt. Die kans is echter klein dus moeten we zelf moeite doen.

Dit forum kan gratis blijven vanwege banners als deze. Door te registeren zal de onderstaande banner overigens verdwijnen.

#2

EvilBro

    EvilBro


  • >5k berichten
  • 6703 berichten
  • VIP

Geplaatst op 03 december 2008 - 08:00

Ik snap de vraag niet. Kun je duidelijker maken wat je wilt? (plaatje of zo?)

#3

jhnbk

    jhnbk


  • >5k berichten
  • 6905 berichten
  • VIP

Geplaatst op 03 december 2008 - 08:49

:D Ik heb dus een serie punten (metingen) waarvan ik het maximum moet bepalen. Ik weet zeker dat dit maximum zich tussen de laatste 3 punten moet zitten. Nu is de vraag ik rechtstreeks uit Octave de coŽfficiŽnten van de veelterm kan halen door cubic spline interpolatie te gebruiken. (Indien niet zal ik dus zelf aan het werk moeten om deze te berekenen)
Het vel van de beer kunnen verkopen vraagt moeite tenzij deze dood voor je neervalt. Die kans is echter klein dus moeten we zelf moeite doen.

#4

jhnbk

    jhnbk


  • >5k berichten
  • 6905 berichten
  • VIP

Geplaatst op 03 december 2008 - 09:09

Ik vond net de oplossing door als methode 'pp' te zetten krijg ik: interp1(w,d,'pp'):
{
  x =

	  0.73000
	  4.16000
	  6.77000
	  9.07000
	 10.82000
	 13.88000

  P =

	-0.012828   1.698000
	 0.023755   1.654000
	 0.010000   1.716000
	 0.012000   1.739000
	-0.037582   1.760000

  n =  5
  k =  2
  d =  1
}

Nu moet ik nog uitvinden wat die waarden zijn. Iemand enig idee?
Het vel van de beer kunnen verkopen vraagt moeite tenzij deze dood voor je neervalt. Die kans is echter klein dus moeten we zelf moeite doen.

#5

EvilBro

    EvilBro


  • >5k berichten
  • 6703 berichten
  • VIP

Geplaatst op 03 december 2008 - 09:32

Ik heb dus een serie punten (metingen) waarvan ik het maximum moet bepalen.

Om dit te doen heb je geen interpolatie nodig. Ik vermoed daarom dat je aan een continu verband hebt gemeten en dat je op basis van deze metingen een functie opstelt (die je op de een of andere manier kunt verantwoorden). Deze functie is in jouw geval 'spline'. Let op: zonder verantwoording voor de functie is je zoektocht naar een maximum betekenisloos.

Nu is de vraag of ik rechtstreeks uit Octave de coŽfficiŽnten van de veelterm kan halen door cubic spline interpolatie te gebruiken.

Dat kan volgens mij wel. Bekijk de help van spline, mkpp en ppval. Is het echter niet veel makkelijker om gewoon te interpoleren met een kleine stapgrootte en dan het maximum te zoeken?

#6

jhnbk

    jhnbk


  • >5k berichten
  • 6905 berichten
  • VIP

Geplaatst op 03 december 2008 - 14:10

Het gaat 'm inderdaad om het maximum van het continu verband dat er zou moeten zijn.

Interpolatie met kleine stapgrootte levert mij inderdaad een resultaat op maar het probleem is dat ik met meetfouten moet werken. Ik ga het er op houden dat mijn voorlaatste punt eigenlijk bijna het maximum is.
Het vel van de beer kunnen verkopen vraagt moeite tenzij deze dood voor je neervalt. Die kans is echter klein dus moeten we zelf moeite doen.

#7

EvilBro

    EvilBro


  • >5k berichten
  • 6703 berichten
  • VIP

Geplaatst op 03 december 2008 - 14:18

Ik ga het er op houden dat mijn voorlaatste punt eigenlijk bijna het maximum is.

Heb je de meting zelf uitgevoerd? Als dat zo is, dan had je je maximum moeten opzoeken (dat is namelijk een interessant punt). Als dat niet zo is, dan is de persoon die gemeten heeft een prutser. :D

Je kunt natuurlijk op basis van je punten en een gegeven verband de foutgrenzen wel doorrekenen, maar daar word je vast niet gelukkig van...

#8

jhnbk

    jhnbk


  • >5k berichten
  • 6905 berichten
  • VIP

Geplaatst op 03 december 2008 - 14:27

Heb je de meting zelf uitgevoerd? Als dat zo is, dan had je je maximum moeten opzoeken (dat is namelijk een interessant punt). Als dat niet zo is, dan is de persoon die gemeten heeft een prutser. :D

De cursus schreef de metingen voor dus meer mochten we niet doen aangezien het een nogal krap tijdsbestek was. Het uitzoeken van het maximum was inderdaad de juiste werkwijze geweest.

Je kunt natuurlijk op basis van je punten en een gegeven verband de foutgrenzen wel doorrekenen, maar daar word je vast niet gelukkig van...

Neen niet bepaald dus ga ik zoals boven geschreven oplossen.
Het vel van de beer kunnen verkopen vraagt moeite tenzij deze dood voor je neervalt. Die kans is echter klein dus moeten we zelf moeite doen.

#9

jhnbk

    jhnbk


  • >5k berichten
  • 6905 berichten
  • VIP

Geplaatst op 10 december 2008 - 15:53

Weet iemand bij toeval hoe ik in Octave de kleuren kan aanpassen van mijn plot?

plot(w,d,'x',wi,di,w,ds);
Ik zou de 2 eerste reeksen dezelfde kleur moeten geven maar ik geraak er via de handleiding niet uit.
Het vel van de beer kunnen verkopen vraagt moeite tenzij deze dood voor je neervalt. Die kans is echter klein dus moeten we zelf moeite doen.

#10

Bart

    Bart


  • >5k berichten
  • 7224 berichten
  • VIP

Geplaatst op 10 december 2008 - 19:13

Als het goed is, is Octave 99.9% Matlab compatible en zou dit moeten werken:

plot(x1, y1, 'b', x2, y2, 'b', x3, y3, 'r');
waar b = blauw en r = rood

Anders kun je dit eens proberen:
plot(x1, y1, 'b');
hold on;
plot(x2, y2, 'b');
plot(x3, y3, 'r');
hold off;
If I have seen further it is by standing on the shoulders of giants.-- Isaac Newton

#11

jhnbk

    jhnbk


  • >5k berichten
  • 6905 berichten
  • VIP

Geplaatst op 11 december 2008 - 17:36

het werkt! Thx
Het vel van de beer kunnen verkopen vraagt moeite tenzij deze dood voor je neervalt. Die kans is echter klein dus moeten we zelf moeite doen.





0 gebruiker(s) lezen dit onderwerp

0 leden, 0 bezoekers, 0 anonieme gebruikers

Ook adverteren op onze website? Lees hier meer!

Gesponsorde vacatures

Vacatures