Beschreiben einer Mathlab-Funktion in C

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

MainzerMath

Forum-Newbie


	Beiträge: 2

	Anmeldedatum: 19.02.14

	Wohnort: ---

	Version: ---

Verfasst am: 19.02.2014, 16:59 Titel: Beschreiben einer Mathlab-Funktion in C

Hallo an alle, ein nicht Matlabkenner braucht hilfe von euch. Ich entwickle inC++. Muss aber ein Mathlab Code von in C++ schreiben.

Es geht hier um ein Kreisfitting. Es geht nur um eine Stelle im Code am Ende.

Code:

function Par = CircleFitByPratt(XY)

%--------------------------------------------------------------------------
%
% Circle fit by Pratt
% V. Pratt, "Direct least-squares fitting of algebraic surfaces",
% Computer Graphics, Vol. 21, pages 145-152 (1987)
%
% Input: XY(n,2) is the array of coordinates of n points x(i)=XY(i,1), y(i)=XY(i,2)
%
% Output: Par = [a b R] is the fitting circle:
% center (a,b) and radius R
%
% Note: this fit does not use built-in matrix functions (except "mean"),
% so it can be easily programmed in any programming language
%
%--------------------------------------------------------------------------

n = size(XY,1); % number of data points

centroid = mean(XY); % the centroid of the data set

% computing moments (note: all moments will be normed, i.e. divided by n)

Mxx=0; Myy=0; Mxy=0; Mxz=0; Myz=0; Mzz=0;

for i=1:n
Xi = XY(i,1) - centroid(1); % centering data
Yi = XY(i,2) - centroid(2); % centering data
Zi = Xi*Xi + Yi*Yi;
Mxy = Mxy + Xi*Yi;
Mxx = Mxx + Xi*Xi;
Myy = Myy + Yi*Yi;
Mxz = Mxz + Xi*Zi;
Myz = Myz + Yi*Zi;
Mzz = Mzz + Zi*Zi;
end

Mxx = Mxx/n;
Myy = Myy/n;
Mxy = Mxy/n;
Mxz = Mxz/n;
Myz = Myz/n;
Mzz = Mzz/n;

% computing the coefficients of the characteristic polynomial

Mz = Mxx + Myy;
Cov_xy = Mxx*Myy - Mxy*Mxy;
Mxz2 = Mxz*Mxz;
Myz2 = Myz*Myz;

A2 = 4*Cov_xy - 3*Mz*Mz - Mzz;
A1 = Mzz*Mz + 4*Cov_xy*Mz - Mxz2 - Myz2 - Mz*Mz*Mz;
A0 = Mxz2*Myy + Myz2*Mxx - Mzz*Cov_xy - 2*Mxz*Myz*Mxy + Mz*Mz*Cov_xy;
A22 = A2 + A2;

epsilon=1e-12;
ynew=1e+20;<=
IterMax=20;
xnew = 0;

% Newton's method starting at x=0

for iter=1:IterMax
yold = ynew;
ynew = A0 + xnew*(A1 + xnew*(A2 + 4.*xnew*xnew));
if (abs(ynew)>abs(yold))
disp('Newton-Pratt goes wrong direction: |ynew| > |yold|');
xnew = 0;
break;
end
Dy = A1 + xnew*(A22 + 16*xnew*xnew);
xold = xnew;
xnew = xold - ynew/Dy;
if (abs((xnew-xold)/xnew) < epsilon), break, end
if (iter >= IterMax)
disp('Newton-Pratt will not converge');
xnew = 0;
end
if (xnew<0.)
fprintf(1,'Newton-Pratt negative root: x=%f\n',xnew);
xnew = 0;
end
end

% computing the circle parameters

DET = xnew*xnew - xnew*Mz + Cov_xy;
Center = [Mxz*(Myy-xnew)-Myz*Mxy , Myz*(Mxx-xnew)-Mxz*Mxy]/DET/2;

Par = [Center+centroid , sqrt(Center*Center'+Mz+2*xnew)];

end % CircleFitByPratt

Funktion ohne Link?

Ich habe eigenlich alles schon in meiner C++ Routine laufen. Der Mittelpunkt wird korrekt berechnet .

Ich verstehe nur die folgende Codezeile nicht:

Code:

sqrt(Center*Center'+Mz+2*xnew)

Funktion ohne Link?

"

Ich benütze in meinen C++ Code für den Vektor "Center" einfach zwei double Variablen dX und dY. Alle andere Variable Mz und xnew sind auch double.

Wie müsste die Zeile sqrt(Center*Center'+Mz+2*xnew) in C aussehen wenn ich anstatt "Center" dX und dY benütze.

Vielen Dank im Voraus

Alex

Sirius3

Forum-Guru


	Beiträge: 441

	Anmeldedatum: 11.11.11

	Wohnort: ---

	Version: ---

Verfasst am: 19.02.2014, 20:51 Titel:

Code:

Center*Center'

Funktion ohne Link?

ist das Skalarprodukt

Code:

(dx*dx + dy*dy)

Funktion ohne Link?

MainzerMath

Themenstarter

Forum-Newbie


	Beiträge: 2

	Anmeldedatum: 19.02.14

	Wohnort: ---

	Version: ---

Verfasst am: 19.02.2014, 23:47 Titel: Vielen Dank

Vielen Dank,

das war es. Jetzt habe ich auch den Radius. ich habe schon mehrere Bekannte gefragt und keiner konnte was mit dem Hochkomma am Ende anfangen. Du hast mir echt weiter geholfen.

VG Alex


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.