dimensions missmatch beim Arbeiten mit function_handle

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Christian St.

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 04.04.2016, 11:34 Titel: dimensions missmatch beim Arbeiten mit function_handle

Hallo miteinander, ich versuchs einfach mal und hoffe dass mir wer helfen kann.

Ich habe ein nichtlineares Gleichungssystem das ich gerne mit einem Newtonverfahren lösen möchte und habe dazu zwei m_files gschrieben (bzw. abgekupfert Confused

) eines das die Schleife mit dem Newtonverfahren enthält:

Code:

function [x,fx,it]=newton(fun,dfun,x0,tol,nmax)
%Parameter:
%fun: Funktion
%dfun: Ableitung der Funktion
%x0: Startwert
%tol: Toleranzwert
%bmax: Maximale Anzahl der Iterationsschritte

%Ausgabe:
%x: angenäherte Nullstelle
%fx: Funktionswert an der Stelle X
%it: Anzahl der benötigten Iterationsschritte

x=x0;
it=0;
a=1+tol;

while a>tol && it < nmax
xtemp=x-(fun(x).*(dfun(x).^-1));
a=norm(xtemp-x);
x=xtemp;
it=it+1;
end
fx=fun(x);
end

Funktion ohne Link?

und das zweite das das Gleichungssystem erstellt:

Code:

fun=@(x) [x(1)-x(2)-x(3); x(3)-x(4)-x(8); x(4)-x(5)-x(6); x(6)+x(5)-x(7); x(7)+x(8)-x(9); x(9)+x(2)-x(10);
mein.^2.*R1 + x(2).^2.*R2 + x(10).^2.*R10 - pges;
x(3).^2.*R3 + x(4).^2.*R4 + x(5).^2.*R5 + x(7).^2.*R7 + x(9).^2.*R9 - x(2).^2.*R2;
x(6).^2.*R6 - x(5).^2.*R5;
x(8).^2.*R8 - x(7).^2.*R7 - x(6).^2.*R6 - x(4).^2.*R4];

dfun= @(x) [1-1-1; 1-1-1; 1-1-1; 1+1-1; 1+1-1; 1+1-1;
2*x(2).*R2 + 2.*x(10).*R10;
2*x(3).*R3 + 2.*x(4).*R4 + 2.*x(5)*R5 + 2.*x(7).*R7 + 2.*x(9)*R9 + 2.*x(2).*R2;
2*x(6).*R6 + 2.*x(5).*R5;
2*x(8).*R8 -2.*x(7).*R7 - 2.*x(6).*R6 - 2.*x(4).*R4];

x0= [1;1;1;1;1;1;1;1;1;1];
tol=10^-9;
nmax=200;
[x,fx,it]=newton(fun,dfun,x0,tol,nmax)

Funktion ohne Link?

ich bekomme immer folgende Fehlermeldung:

error: vertical dimensions mismatch (9x1 vs 1x2)
error: called from
para3> at line -1 column -1
newton at line 20 column 10
para3 at line 42 column 9

kann mir aber nicht erklären wie der Fehler entstehen soll...

Hat wer eine Idee was ich falsch gemacht habe?

SG Christian

Winkow

Moderator


	Beiträge: 3.842

	Anmeldedatum: 04.11.11

	Wohnort: Dresden

	Version: R2014a 2015a

Verfasst am: 04.04.2016, 12:14 Titel:

wenn dimmensionsfehler auftreten ist es meinst einfacher mit dem debugger selber rein zu gucken was passiert an der stelle. setz doch mal die debug marke an der stelle.
das beispiel ist ja an sich so nicht lauffähig. darum kannst du das nur selber machen.
_________________

richtig Fragen

Christian St.

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 04.04.2016, 12:47 Titel:

Ja das ist richtig ich wollte nicht alles mit code zupflastern...

der Debugger hält immer, wie in der Errormeldung, in Zeile 20 bei dem ausdruck:

Code:

xtemp=x-(fun(x).*(dfun(x).^-1));

Funktion ohne Link?

überraschenderweise haben sich nun die Dimensionen verändert, (op1 is 10x1, op2 is 0x0)

ich verstehe aber nicht wieso, eigentlich sollten beide 10x1 sein

hier ist der komplette Code:

Code:

clc;
clear all;
debug_on_error (1)
# Erstellen der Parameter:
mein=10; pges=7*10^5;
# Mediumeigenschaften:
P=965; v=0.326*10^-6;

# Durchmesser:
d1 = 0.2073 ; d3 = 0.2073 ; d9 = 0.2073 ; d10 = 0.2073 ;
d8 = 0.1071 ; d2 = 0.1071 ;
d4 = 0.0825 ; d7 = 0.0825 ;
d5 = 0.0703 ; d6 = 0.0703 ;

# Längen:
l1=120 ; l10=120 ; l2=60 ; l3=50 ; l9=50 ; l4=25 ; l7=25 ; l5=15 ; l6=15 ; l8=75 ;

# Lambda:
A200=0.016 ; A100=0.018 ; A80=0.020 ; A65=0.022 ;

# Teilstreckenwiderstände
R1=A200*l1*8/(d1*pi.^2*d1.^4*P); R2=A100*l2*8/(d2*pi.^2*d2.^4*P);
R3=A200*l3*8/(d3*pi.^2*d3.^4*P); R4=A80*l4*8/(d4*pi.^2*d4.^4*P);
R5=A65*l5*8/(d5*pi.^2*d5.^4*P); R6=A65*l6*8/(d6*pi.^2*d6.^4*P);
R7=A80*l7*8/(d7*pi.^2*d7.^4*P); R8=A100*l8*8/(d8*pi.^2*d8.^4*P);
R9=A200*l9*8/(d9*pi.^2*d9.^4*P); R10=A200*l10*8/(d10*pi.^2*d10.^4*P);

fun=@(x) [x(1) - x(2) - x(3); x(3) - x(4) - x(8);
x(4)-x(5)-x(6); x(6)+x(5)-x(7);
x(7) + x(8) - x(9); x(9) + x(2) - x(10);
mein.^2.*R1 + x(2).^2.*R2 + x(10).^2.*R10 - pges;
x(3).^2.*R3 + x(4).^2.*R4 + x(5).^2.*R5 + x(7).^2.*R7 + x(9).^2.*R9 - x(2).^2.*R2;
x(6).^2.*R6 - x(5).^2.*R5; x(8).^2.*R8 - x(7).^2.*R7 - x(6).^2.*R6 - x(4).^2.*R4];

dfun= @(x) [1-1-1; 1-1-1;
1-1-1; 1+1-1;
1+1-1; 1+1-1;
2*x(2).*R2 + 2.*x(10).*R10;
2*x(3).*R3 + 2.*x(4).*R4 + 2.*x(5)*R5 + 2.*x(7).*R7 + 2.*x(9)*R9 + 2.*x(2).*R2;
2*x(6).*R6 + 2.*x(5).*R5;
2*x(8).*R8 -2.*x(7).*R7 - 2.*x(6).*R6 - 2.*x(4).*R4];

x0= [1;1;1;1;1;1;1;1;1;1];
tol=10^-9;
nmax=200;
[x,fx,it]=newton(fun,dfun,x0,tol,nmax)

Funktion ohne Link?

Winkow

Moderator


	Beiträge: 3.842

	Anmeldedatum: 04.11.11

	Wohnort: Dresden

	Version: R2014a 2015a

Verfasst am: 04.04.2016, 15:32 Titel:

in der letzten zeile von dfun hast du ein leerzeichen vergessen.
_________________

richtig Fragen

Christian St.

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 05.04.2016, 07:25 Titel:

Aha,

Ich wusste nicht dass die Leerzeichen auch ausgewertet werden,
Das Gleichungssystem scheint zwar nicht richtig zu sein, aber er rechnet jetzt zumindest.

Danke für die Hilfe!
SG Christian


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 - 2024 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.