Polyfit zu Funktion um schreiben

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

martin1234

Forum-Newbie


	Beiträge: 8

	Anmeldedatum: 05.09.11

	Wohnort: ---

	Version: ---

Verfasst am: 05.09.2011, 17:43 Titel: Polyfit zu Funktion um schreiben

Hi...

also ich habe mit dem polyfit- und dem polyval-Befehl, ein Polynom erstellt (6. Grades).
Also besitze ich nun die verschiedenen Koeffizienten meines Polynoms. Nun würde ich gerne eine Funktion erstellen und in Abhängigkeit der 'X-Koordinate' haben, da ich mit meinem Polynom weiter rechnen möchte (ich muss es integrieren). Hoffe mir kann jemand helfen. Denke eigentlich das sollte machbar sein.

Code:

werte_des_warmekapazitat = xlsread('C....);

temp = werte_des_warmekapazitat(:,1);
cp = werte_des_warmekapazitat(:,2);

plot(temp, cp, ':s', 'MarkerSize', 2, 'MarkerFaceColor','b')
title('Wärmekapazität über cp bei 25 MPa')
xlabel('Temeperatur')
ylabel('cp')

p = polyfit(temp, cp, 6);
pp = polyval(p, temp);
hold on
plot(temp, pp, '-r')

legend('cp','Kurve', 'location','NorthWest')

Funktion ohne Link?

_________________

Vielen Dank

Thomas84

Forum-Meister


	Beiträge: 546

	Anmeldedatum: 10.02.10

	Wohnort: ---

	Version: ---

Verfasst am: 06.09.2011, 07:35 Titel:

Code:

% daten generieren
x = 1:10;
y = rand(size(x));

% fit
p = polyfit(x,y,6);

% function handle erstellen
f = @(arg) polyval(p,arg);

% integrieren
int = quad(f,1,4);

%ploten
xvec = 1:0.01:10;
plot(x,y,'r+',xvec,f(xvec));

Funktion ohne Link?

Das Stichwort heißt function handle.

viele Grüße
Thomas

martin1234

Themenstarter

Forum-Newbie


	Beiträge: 8

	Anmeldedatum: 05.09.11

	Wohnort: ---

	Version: ---

Verfasst am: 06.09.2011, 11:29 Titel:

Hey vielen Dank Thomas,

ist es denn noch möglich ein Integral nicht zu einem gewissen wert laufen zu lassen, sondern zu einer variablen.

Code:

int = quad(f,300,x)
% dieses integral würde ich gerne in einen fsolve befehl stecken und dann nach diesem "x" auflösen
fsolve(...=int,x) % kann den genauen fsolve befehl noch nicht

Funktion ohne Link?

_________________

Vielen Dank

Thomas84

Forum-Meister


	Beiträge: 546

	Anmeldedatum: 10.02.10

	Wohnort: ---

	Version: ---

Verfasst am: 06.09.2011, 14:41 Titel:

Code:

% daten generieren
x = 1:10;
y = rand(size(x));

% fit
p = polyfit(x,y,6);

% function handle erstellen
f = @(arg) polyval(p,arg);

% function handle für das Integral
f_int = @(x_end) quad(f,1,x_end);

% Lösung der Gleichung Integral(f) (x_end) = 1
xend = fsolve(@(x_end) f_int(x_end) - 1,1)

% Probe
quad(f,1,xend)

Funktion ohne Link?

Auch das geht mit function handles. Man braucht ein wenig Übung, aber wenn man weiss wie man damit umgehen kann sind sie sehr praktisch.


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 |

RSS

Hosted by:

Copyright © 2007 - 2025 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.