Annäherung an Pi mit Monte carlo, paintbox

[Borgschüler]

hmm in der schule haben wir das immer mit der paintbox gemacht.

randomize wird im button eh nur einmal aufgerufen. wenn ich es bei formcreate reinschreibe, ändert isch leider garnichts. er macht irgendwie immer nur die hälfte von der paintbox. es sei denn du machst die paintbox richtig klein, aber das kanns auch nicht sein!

[Klaus01]

 if (sqrt(x*x+y*y) < 200) ->

Delphi-Quellcode:

if sqrt(sqr(x)+sqr(y)) <= sqrt(sqr(paintbox.width)+(sqr(paintbox.heigth)) then
  //im radius
else
 //asserhalb des Radius

Grüße
Klaus

[einbeliebigername]

Hallo,

Zitat von Klaus01:

if sqrt(sqr(x)+sqr(y)) <= sqrt(sqr(paintbox.width)+(sqr(paintbox.heigth)) then

Das ist aber nicht im Sinne des Erfinders. Denn damit ist der Radius des Kreis so groß wie die Diagonale der PaintBox. Das hat dann die Folge, dass alle Punkte, die man sieht, innerhalb des Kreises liegen. Besser ist das:

Delphi-Quellcode:

  R:= Min(paintbox.width, paintbox.heigth)
  ...
    x:= random(R);
    y:= random(R);
    if (sqrt(x*x+y*y) < R) then
  ...

Einbeliebigername.

[TiGü]

Wie soll es denn richtig aussehen?
Wenn ich deinen Code ausführe, dann kommt das bei rum:

[guinnes]

Zitat von TiGü:

Wie soll es denn richtig aussehen?

Genau so ! : http://de.wikipedia.org/wiki/Monte-Carlo-Simulation

[Borgschüler]

tigü ja genau, so soll es aussehen!
tigü, hast du meinen quelltext imn ersten beitrag eins zu eins übernommen?
hmm, dann probier ich mal, den radius mit paintbox width auszudürcken!

EDIT: Ich werd noch wahnsinnig, jetzt hab ich < paintbox1.width und es geht immer noch nicht!


übrigens:

als zum berechnen von pi dachte ich mir

K/Q = 200²Pi/200²

erklärung: Verhältnis der punkte von Kreis und Quadrat(gesamtes) = verhältnis von Viertelkreisgleichung mit festgelegtem radius und fläche von quadrat

und daraus drück ich mir aus:

K/Q mal 4 = Pi

[Borgschüler]

so so schaut das aus bei mir, siehe anhang

[Valle]

Kann es sein dass die Paintbox auf dem Formular einfach zu klein ist?

Liebe Grüße,
Valentin

[DeddyH]

Vielleicht legt sich kurz ein anderes Fenster drüber. Verschieb das Zeichnen doch mal ins OnPaint der PaintBox, da gehört es ja auch hin. Bei mir wird die Grafik auch korrekt dargestellt.

[edit] Meine Änderungen:

Delphi-Quellcode:

type
  TForm1 = class(TForm)
    PaintBox1: TPaintBox;
    Button1: TButton;
    procedure Button1Click(Sender: TObject);
    procedure PaintBox1Paint(Sender: TObject);
    procedure FormCreate(Sender: TObject);
  private
    { Private-Deklarationen }
    FPainting: Boolean;
  public
    { Public-Deklarationen }
  end;

...

procedure TForm1.Button1Click(Sender: TObject);
begin
  paintbox1.Width := 200;
  paintbox1.Height := 200;
  FPainting := not FPainting;
  PaintBox1.Invalidate;
end;

procedure TForm1.FormCreate(Sender: TObject);
begin
  randomize;
end;

procedure TForm1.PaintBox1Paint(Sender: TObject);
var
  i,x,y : integer;
begin
  if FPainting then
    for i := 1 to 100000 do
      begin
        x:= random(200);
        y:= random(200);
        if (sqrt(x*x+y*y) < 200) then
          paintbox1.canvas.pixels[x,y] := clred
        else
          paintbox1.canvas.pixels[x,y] := clgreen;
      end
  else
    PaintBox1.Canvas.FillRect(Rect(0, 0, PaintBox1.Width, PaintBox1.Height));
end;

[/edit]

[Borgschüler]

Liebe Leute, je nach dem wo die paintbox platziert ist, ist ein anderer teil zu sehn, bzw an einer stelle sieht man wirklich ein perfektes quadrat! aber seht selbst.

Wie kann das bitte sein? geht garnicht in meinen kopf!

ich werde jetzt die pi berechnung durchführen, und dann das ergebnis präsentieren.