c++ lernen. schleifen? wer kann helfen??

[Sebastian.S]

Hallo,
ich bin gerade dabei, c++ zu lernen. Meine erste Schleife soll das Heron-Verfahren anwenden, um eine Wurzel zu ziehen. (Heron ist auf: http://home.t-online.de/home/arndt.brue ... frame2.htm dokumentiert)

Meine Schleife sieht so aus:

Code: Alles auswählen

#include <iostream>
using namespace std;

int main() {

double wurzel_aus=2, startwert=2, quadrat=4;
int durchlauf=10, durchlauf2=0;

cout << "Bitte geben Sie die zu radizierende Zahl ein: _" ;
cin >> wurzel_aus;
cout << "\nStartwert? _";
cin >> startwert;


	 while (durchlauf>0) {
	//while (quadrat = wurzel_aus) {
		startwert=((startwert+wurzel_aus/startwert)/2);
		quadrat=startwert*startwert;
		if (quadrat=wurzel_aus) --durchlauf; else durchlauf=durchlauf;
		++durchlauf2;
		// --durchlauf;
	}
cout << "Ergebnis: " << startwert << endl << (startwert*startwert) << endl << durchlauf2 << endl;

}

Das Ergebnis:

Code: Alles auswählen

Bitte geben Sie die zu radizierende Zahl ein: _159753

Startwert? _1
Ergebnis: 404.478
163602
10

Das Programm macht immer nur zehn Durchläufe, obwohl die Wurzel noch nicht vollständig berechnet ist. Warum wird die if-Schleife ausgeführt??

[feltel]

... moved nach "Softwareentwicklung und -paketierung, Scripting"

[keyem]

Du hast einen typischen C-Fehler in Deinem Programm: in C/C++ wird für den Vergleich der "=="-Operator benutzt, für die Zuweisung das "=".

keyem

[Sebastian.S]

Ups. Wohl ein Anfängerfehler...

Noch eine Frage: Ich habe das Problem jetzt gelöst, aber die Schleifer wird nie verlassen, obwohl das Programm testweise auch bei 65535 Durchläufen die Bedingung Quadrat==wurzel_aus erfüllt hatte. Warum wird die Schleife nicht verlassen?

Schleife wurde angepasst:

Code: Alles auswählen

	while (quadrat!=wurzel_aus) {
		startwert=((startwert+wurzel_aus/startwert)/2);
		quadrat=startwert*startwert;
		//--durchlauf; //if (-1 < quadrat-wurzel_aus < 1) {--durchlauf;}
		 ++durchlauf2;
	}

[jogix]

Äh?! Du hast aber nicht
if (-1 < quadrat-wurzel_aus < 1)
als Abfrage da stehen oder?! Das geht nämlich nicht!

Ansonsten kann ich Dir nur raten, bei numerischen Problemen die numerischen Ausdrücke nicht mit == bzw. != abzufragen, da oft kleine Ungenauigkeiten von z.B. 1*10^-8 in den ausdrücken zu finden sind und somit == nicht erfüllt wird!

[Sebastian.S]

Äh?! Du hast aber nicht
if (-1 < quadrat-wurzel_aus < 1)
als Abfrage da stehen oder?! Das geht nämlich nicht!

Ansonsten kann ich Dir nur raten, bei numerischen Problemen die numerischen Ausdrücke nicht mit == bzw. != abzufragen, da oft kleine Ungenauigkeiten von z.B. 1*10^-8 in den ausdrücken zu finden sind und somit == nicht erfüllt wird!

Das

Code: Alles auswählen

if (-1 < quadrat-wurzel_aus < 1)

war eine Änderung, die ich probiert habe, weil ich das mit den Ungenauigkeiten vermutet hatte. Versehentlich hatte ich also nicht die originale Version ins Posting kopiert. Ist ja aber auch kommentiert.

Wie kann ich denn diese Ungenauigkeiten umgehen? Gibt es einen Ungefähr-Operator?

[Dookie]

Hi,

du könntest einen "ungefähr" Operator mit einer Funktion simulieren.

in Python würd ich den so machen

Code: Alles auswählen

nearly = lambda a,b,precision: abs(a-b) < precision
>>> x = 12345.6789
>>> y = 12345.6788
>>> print nearly(x,y,0.001)
1

Das Prinzip ist ganz einfach, ziehe den wert b vom wert a ab und vergleiche die absolute Differenz mit Deinem Schwellwert.

Gruß

Dookie

[jogix]

Hi,

ich jetzt nicht die Absolutwertberechnung von C++, würde das aber entsprechend so aufrufen:

Code: Alles auswählen

if (abs(quadrat-wurzel_aus) < delta ) ...

Damit bildest Du den Betrag der Differenz, die ja um die Null herum sein sollte. Durch den Betrag ist das ganze immer positiv und somit sinnvoll gegen einen kleinen Wert Delta zu prüfen.
Delta kannst Du als Konstante definieren und dann mit der Genauigkeit spielen: z.B.
delta = 1e-5;

[Sebastian.S]

Vielen Dank für die Antworten...


Leider habe ich da so ein Problem mit der abs-Funktion:

Code: Alles auswählen

test1.c: In function `int main(...)':
test1.c:15: warning: `double' used for argument 1 of `abs(int)'

Also erwartet abs einen integer-Wert, womit die Genauigkeit aber wieder hinfällig wäre.

[keyem]

"fabs" löst Dein Problem!

Viel Erfolg weiterhin,
keyem

[Sebastian.S]

Falls es jemanden interessiert, stelle ich das Programm einfach mal hierher:

Code: Alles auswählen

#include <iostream>
#include <cmath>

 using namespace std;
// g++ -I/usr/include/g++-3/ -o test3 test1.c



 void main(void) {

 double wurzel_aus=2, startwert=2, quadrat=4;
 unsigned int durchlauf=65535, durchlauf2=0;
 double delta;

 cout << "Bitte geben Sie die zu radizierende Zahl ein: _" ;
 cin >> wurzel_aus;
 cout << "\nStartwert? _";
 cin >> startwert;
  cout << "\nGenauigkeit? (in 1*10 ^ -x) _";
 cin >> delta;

 delta=(-1)*fabs(delta);
delta=1*pow(10, delta);
//cout << delta;


	while (fabs(quadrat-wurzel_aus) > delta) {
		startwert=((startwert+wurzel_aus/startwert)/2);
		quadrat=startwert*startwert;
		 ++durchlauf2;
	}

 cout << "\n\nErgebnis: " << startwert << endl << (startwert*startwert) << endl << durchlauf2 << endl<<quadrat<<endl;

 }