regressionsgerade < Maple < Mathe-Software < Mathe < Vorhilfe
|
Status: |
(Frage) überfällig | Datum: | 21:42 Do 01.03.2007 | Autor: | keko |
Aufgabe | Erstelle eine Regressionfunktion zu den vorgegeben Messwerten eines Abkühlungsprozesses!
|
Ich will für meine Messwerte in einem Abkühlungsversuch eine entsprechende Regressionsfunktion errechnen mit deren Hilfe ich nicht nur meine Messpunkte auf den Graph zeichnen, sondern auch eine Funktion auf den Graph zeichnen kann! Jedoch haben meine Versuche bisher immer fehl geschlagen!(ich arbeite mit maple6)
Dies ist mein Versuch:
> restart:
> with(stats):with(fit):with(statplots):
> datax:=[0,2,4,6,8,10,12,14,16,18,20,25,30,35,40,45,50,60,70]:
> datay:=[37,34.5,32.5,30,28.3,26.8,25.1,24.1,23,21.8,20.8,18.8,17,15,14,12.9,12,10,8.9]:
> dataxy:=zip((x,y)->[x,y],datax,datay):
> fy:=x->a0+a1*exp(x):
> f:=LeastSquares(datax,datay,x,curve=fy):
> p1:=scatterplot(datax,datay,color=blue):
> p2:=plot(f(x),x=-1..71,y=7..38):
> with(plots):
> display([p1,p2]);
</task>
ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt!
wäre nett wenn mir jemand helfen könnte!
|
|
|
|
Status: |
(Mitteilung) Reaktion unnötig | Datum: | 02:29 Fr 02.03.2007 | Autor: | leduart |
Hallo
diese Daten kannst du doch sicher nicht mit f=a+b*exp(x) irgendwie modellieren!
1. exp(x) waechst, deine fkt faellt
2. wenn es was exponentielles sein soll dann [mm] f=a+bx+e^{-cx}
[/mm]
von deinem programm versteh ich sonst zu wenig.
Gruss leduart
|
|
|
|
|
Status: |
(Mitteilung) Reaktion unnötig | Datum: | 11:57 Fr 02.03.2007 | Autor: | keko |
Vielen dank für die Hilfe.
Doch leider Funktioniert es auch jetzt nicht. der Fehler muss wohl an einer anderen Stelle liegen.
Ich bekomme bei dem Befehl with(plots) eine folgende Fehlermeldung:
"Warning, the name changecoords has been redefined"
Liegt mein Problem bei diesem Befehl?
|
|
|
|
|
Status: |
(Mitteilung) Reaktion unnötig | Datum: | 12:19 Fr 02.03.2007 | Autor: | Martin243 |
Hallo,
das ist nur eine Warnung und unwichtig. Aber du bringst hier Packages durcheinander:
LeastSquares ist im Package CurveFitting enthalten und nicht in fit (das es gar nicht gibt???)
Außerdem erwartet LeastSquares für fy keine Funktion, ondern nur den Ausdruck (also ohne "x->").
Ersetze einfach mal Testweise die Doppelpunkte durch Semikolons, dann siehst du, dass einige Ausdrücke gar nicht berechnet werden...
Ansonsten käme evtl. noch Statistics[ExponentialFit] in Frage.
Gruß
Martin
|
|
|
|
|
Status: |
(Mitteilung) Reaktion unnötig | Datum: | 15:22 Fr 02.03.2007 | Autor: | keko |
vielen dank, jedoch besteht mein problem weiterhin.
Maple zeichnet mir meine Messpunkte ein, zeichnet jedoch keien Kurve.
kann ich das noch irgendwie auf eine andere Art mit Maple6 bewältigen?
> restart:
> with(stats);
[anova, describe, fit, importdata, random, statevalf, statplots,
transform]
> with(statplots);
[boxplot, histogram, scatterplot, xscale, xshift, xyexchange,
xzexchange, yscale, yshift, yzexchange, zscale, zshift]
> datax:=[0,2,4,6,8,10,12,14,16,18,20,25,30,35,40,45,50,60,70]:
> datay:=[37,34.5,32.5,30,28.3,26.8,25.1,24.1,23,21.8,20.8,18.8,17,15,14,12.9,12,10,8.9]:
> dataxy:=zip((x,y)->[x,y],datax,datay):
> fy:=a0+a1*exp(-(a2*x)):
> f:=LeastSquares(datax,datay,x,curve=fy):
> p1:=scatterplot(datax,datay,color=blue):
> p2:=plot(f(x),x=-1..71,y=7..38):
> with(plots):
Warning, the name changecoords has been redefined
> display([p1,p2]);
|
|
|
|
|
Status: |
(Mitteilung) Reaktion unnötig | Datum: | 18:45 Fr 02.03.2007 | Autor: | Martin243 |
Hallo,
nochmal:
Die Funktion LeastSquares ist hier gar nicht verfügbar (einfach mal dahinter ein Semikolon setzen)!!! Entweder nimmst du dazu noch das Package CurveFitting (falls vorhanden) oder du benutzst die Funktion ExponentialFit im Package Statistics (falls vorhanden) anstelle von LeastSquares.
Gruß
Martin
|
|
|
|
|
Status: |
(Mitteilung) Reaktion unnötig | Datum: | 21:22 Fr 02.03.2007 | Autor: | keko |
nochmals vielen dank, doch leider funktioniert dies auch wieder nicht:
> restart:
> with(Statistics);
Error, (in pacman:-pexports) Statistics is not a package
> with(stats);
[anova, describe, fit, importdata, random, statevalf, statplots,
transform]
> with(statplots);
[boxplot, histogram, scatterplot, xscale, xshift, xyexchange,
xzexchange, yscale, yshift, yzexchange, zscale, zshift]
> datax:=[0,2,4,6,8,10,12,14,16,18,20,25,30,35,40,45,50,60,70]:
> datay:=[37,34.5,32.5,30,28.3,26.8,25.1,24.1,23,21.8,20.8,18.8,17,15,14,12.9,12,10,8.9]:
> dataxy:=zip((x,y)->[x,y],datax,datay):
> fy:=a0+a1*exp(-a0*x);
fy := a0 + a1 exp(-a0 x)
> f:=ExponentialFit(datax,datay,x,curve=fy);
f := ExponentialFit([0, 2, 4, 6, 8, 10, 12, 14, 16, 18, 20, 25, 30,
35, 40, 45, 50, 60, 70], [37, 34.5, 32.5, 30, 28.3, 26.8,
25.1, 24.1, 23, 21.8, 20.8, 18.8, 17, 15, 14, 12.9, 12, 10,
8.9], x, curve = a0 + a1 exp(-a0 x))
> p1:=scatterplot(datax,datay,color=blue):
> p2:=plot(f(x),x=-1..71,y=7..38);
p2 := PLOT(AXESLABELS("x", "y"), VIEW(-1. .. 71., 7. .. 38.))
> with(plots):
Warning, the name changecoords has been redefined
> display([p1,p2]);
Maple zeigt mir immer nur die Messpunkte.
mfg, keko
|
|
|
|
|
Status: |
(Mitteilung) Reaktion unnötig | Datum: | 22:48 Fr 02.03.2007 | Autor: | Martin243 |
Hmmm,
das bedeutet, du hast dieses Package nicht. Hast du es denn mit dem Package CurveFitting ausprobiert? Falls du das auch nicht hast, weiß ich im Moment auch nicht weiter. Da kannst du es höchstens selber programmieren...
Gruß
Martin
|
|
|
|
|
Status: |
(Mitteilung) Reaktion unnötig | Datum: | 13:23 Sa 03.03.2007 | Autor: | keko |
Ja ich habe es auch mit dem Package CurveFitting probiert.
Dann werd ich es mal versuchen selbst zu programieren.
vielen dank für die hilfe!
mfg keko
|
|
|
|
|
Status: |
(Mitteilung) Reaktion unnötig | Datum: | 12:36 Fr 02.03.2007 | Autor: | luis52 |
Moin keko,
moechtest du die Graphik erstellt haben (mit den zugehoeriegen Werten)
oder einen Rat zu maple6? Im letzteren Fall muss ich passen.
|
|
|
|
|
Status: |
(Mitteilung) Reaktion unnötig | Datum: | 15:25 Fr 02.03.2007 | Autor: | keko |
hey luis52,
ich möchte eine Kurve erstellen, die die Funktion der Messpunkte wiederspiegelt.
|
|
|
|
|
Status: |
(Mitteilung) Reaktion unnötig | Datum: | 20:47 Fr 02.03.2007 | Autor: | luis52 |
Hab's doch noch entdeckt, wie man Dateien anhaengen kann.
Dateianhänge: Anhang Nr. 1 (Typ: pdf) [nicht öffentlich] Anhang Nr. 2 (Typ: pdf) [nicht öffentlich]
|
|
|
|
|
Status: |
(Mitteilung) Reaktion unnötig | Datum: | 21:18 Fr 02.03.2007 | Autor: | keko |
vielen dank, aber wahrscheinlich hab ich mich falsch ausgedrückt. ich brauche zwar die aufzeichnung der funktion, aber muss ich das mit hilfe von maple bewerkstelligen.
gruß keko
|
|
|
|
|
Status: |
(Mitteilung) Reaktion unnötig | Datum: | 22:20 Sa 03.03.2007 | Autor: | matux |
$MATUXTEXT(ueberfaellige_frage)
|
|
|
|