|
|
INF - Infinity lösen, Fakultät programmieren |
|
sun-J |
Forum-Newbie
|
|
Beiträge: 7
|
|
|
|
Anmeldedatum: 08.10.22
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 08.10.2022, 09:46
Titel: INF - Infinity lösen, Fakultät programmieren
|
|
Hallo,
das ist mein erster Beitrag in einem Forum, aber ich stehe momentan mit meiner Suche im www an. Für eine Studi-Hausaufgabe muss ich die Fakultät mit so wenigigen MATlab-eigenen Funktionen programmieren, wie möglich. Ausgedacht habe ich mir folgende Funktion:
function f=fakultaet(n)
f=1;
for k=1:1:n
f=f*k;
end
Mit dem Funktionsaufruf >> fakultaet(170) bekomme ich ab Werten >170 als Ergebnis "INF". Soweit ich weiß, handelt es sich um einen Overflow. Matlab rechnet in double-Werten und die bits und bytes reichen nicht mehr aus, um so große Zahlen darzustellen. Soweit die Theorie.
Ich frage mich nun: Rechnet Matlab im Hintergrund mit den realen Werten und kann es nur nicht mehr am Bildschirm richtig ausgeben, oder muss ich meine Variablen umformatieren? Habt ihr eine Idee, wie ich das umsetzen kann, dass ich nicht mehr INF bekomme?
Danke schonmal für's Lesen bis hierher.
LG Sonja
|
|
|
|
|
Harald |
Forum-Meister
|
|
Beiträge: 24.492
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
Verfasst am: 08.10.2022, 22:31
Titel:
|
|
Hallo,
MATLAB rechnet mit Double, d.h. es gibt eine höchste darstellbare Zahl.
Da ist es egal, ob du deine Funktion aufrufst oder die fertige Funktion factorial.
Mir ist keine Möglichkeit bekannt, wie man mit größeren Zahlen arbeiten kann.
Die Frage ist: gibt es Anwendungen, für die solche Zahlen benötigt werden? Ich kenne jedenfalls keine.
Es kann höchstens sein, dass als Zwischenergebnis solche Zahlen auftreten können, z.B. bei Binomialkoeffizienten. Dann ist es wichtig, die Operationen in geschickter Reihenfolge auszuführen. Das macht MATLAB:
Grüße,
Harald
_________________
1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
|
|
|
sun-J |
Themenstarter
Forum-Newbie
|
|
Beiträge: 7
|
|
|
|
Anmeldedatum: 08.10.22
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 09.10.2022, 17:20
Titel:
|
|
Hi,
Danke für deine schnelle Antwort. Genau für einen Binomialkoeffizienten brauche ich die Fakultät. Ich schreib' mal meinen Lösungsansatz:
%Berechne Fakultaet
function f=fakultaet(n)
f=1;
double n;
for k=1:1:n
f=f*k;
end
%Berechne Binomialkoeffizient
function b=binom(n,k)
h=n-k;
b=fakultaet(n)/(fakultaet(h)*fakultaet(k));
end
Und diese beide Funktionen brauche ich widerum für meine Hauptaufgabe, die Wahrscheinlichkeit P(X>=5) auszurechnen. Nur leider habe ich Werte, z.B. n=1000 und k=5, die mir schnell mal den Wert INF ausgeben. Vielleicht kann man die Zwischenergebnisse auch logarithmieren, damit MATlab das in komprimierter Form darstellen kann? Ansonsten muss ich wohl auf die Funktion nchoosek ausweichen.
Grüße, Sonja
|
|
|
Harald |
Forum-Meister
|
|
Beiträge: 24.492
|
|
|
|
Anmeldedatum: 26.03.09
|
|
|
|
Wohnort: Nähe München
|
|
|
|
Version: ab 2017b
|
|
|
|
|
|
Verfasst am: 09.10.2022, 20:58
Titel:
|
|
Hallo,
das Problem liegt darin, dass du Zähler und Nenner des Binomialkoeffizienten ausmultiplizierst. Du kannst dir im Code von nchoosek eine Möglichkeit ansehen, das zu vermeiden.
Gut möglich, dass es auch andere Strategien gibt.
Grüße,
Harald
_________________
1.) Ask MATLAB Documentation
2.) Search gomatlab.de, google.de or MATLAB Answers
3.) Ask Technical Support of MathWorks
4.) Go mad, your problem is unsolvable ;)
|
|
|
sun-J |
Themenstarter
Forum-Newbie
|
|
Beiträge: 7
|
|
|
|
Anmeldedatum: 08.10.22
|
|
|
|
Wohnort: ---
|
|
|
|
Version: ---
|
|
|
|
|
|
Verfasst am: 10.10.2022, 09:55
Titel:
|
|
Hallo Harald,
danke für deine Hilfe. Damit sollte ich weiter kommen.
Schönen Tag noch!
LG
|
|
|
|
|
Einstellungen und Berechtigungen
|
|
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.
|
|