Symbolische Funktionen zu numerischen Funktionen konvertiere

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

skyforce

Forum-Fortgeschrittener


	Beiträge: 82

	Anmeldedatum: 04.09.12

	Wohnort: ---

	Version: ---

Verfasst am: 19.11.2013, 14:05 Titel: Symbolische Funktionen zu numerischen Funktionen konvertiere

Hall zusammen,

ich stehe seit einigen Stunden vor folgendem Problem:

Ich möchte die part. Ableitungen einer Funktion symbolisch berechnen und anschließend numerisch auswerten:

Code:

% Quadratische Funktion der Form: x1^2 + a*x2^2
a = 1;
x = sym('x',[2,1]);
f = x(1)^2 + a*x(2)^2;

% konvertiere die Funktion f und seine part. Ableitungen von symbolisch zu numerisch.
g1 = matlabFunction((diff(f,x(1))));
g2 = matlabFunction((diff(f,x(2))));
f = matlabFunction(f);

Funktion ohne Link?

Leider liefert die matlabFunction folgende Vorschrift:

f = @(x1,x2) x1^2 + x2^2;
g1 = @(x1)x1.*2.0;
g2 = @(x2)x2.*2.0;

Ich möchte jedoch, dass die Funktion f sowie seine part. Ableitungen folgende Form haben:

f = @(x) x(1)^2 + x(2)^2;
g1 = @(x) x(1).*2.0;
g2 = @(x) x(2).*2.0;

Damit kann ich dann beim Auswerten einfach einen Vektor einsetzen. Ich bitte um Eure Hilfe.

Vielen Dank im Voraus!

skyforce

Themenstarter

Forum-Fortgeschrittener


	Beiträge: 82

	Anmeldedatum: 04.09.12

	Wohnort: ---

	Version: ---

Verfasst am: 19.11.2013, 15:14 Titel:

Ich habe einfach die ganzen funktionen in Strings umgewandelt und anschließend eine reform_fkt geschrieben, welche meine gewünschte Form liefert:

Code:

function [ergebnis] = reform_fkt(f)
% Diese funktion form die Funktionsvorschrift wie folgt um:
% @(x1,x2,...)x1^2 + x2*x1 + .... => @(x) x(1)^2 + x(2)*x(1)

fktstr = func2str(f);
ergebnis = fktstr;
dim = nargin(f);
help1 = '';

% @(x1,x2,...) => @(x)
for i=1:dim
help1 = [help1,'x',int2str(i),','];
end
help1 = help1(1:length(help1)-1);

ergebnis = strrep(ergebnis,help1,'x');

% xi => x(i)
for i=1:dim
help = ['x',int2str(i)];
if (strfind(fktstr,help) ~= 0)
ergebnis = strrep(ergebnis,help,['x(',int2str(i),')']);
end
end
end

Funktion ohne Link?

Es wäre aber natürlich besser, wenn die Fkt "matlabFunction" das könnte


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.