ChrisBatho/Lorenz Plot

From J Wiki
Jump to: navigation, search

Describe Chris-Batho here.

Lorenzcontrol.jpg

This is the control panel [{{#file: "lorenz.ijs"}} Download script: lorenz.ijs ]

NB. The Lorenz Attractor for j601  version 1.02
NB. by [[User:Chris Batho|Chris Batho]], and Oleg Kobchenko
NB. http://www.jsoftware.com/jwiki/Scripts
NB. June 2007  Note Stereo View by Devon McCormick
NB. s r b nominal 10 , 28 , 8%3
Lorenz2.jpg

This is the rotating view

Main function to calc and display lorenz. [{{#file: "lorenz.ijs"}} Download script: lorenz.ijs ]

NB.===================================================
NB. main function to calc and display lorenz
flag=:0
require 'viewmat graph numeric trig plot'

lorenz=: 3 : 0
[A=: 0
[B=: 0
[C=: 0
'X Y Z'=:0.01, 0.01, 0.01   NB. define X Y Z
's r b num_inc' =. y
for. i.num_inc do.
dt=. 0.005
X=: X + (s*(Y - X))*dt
Y=: Y + ((X * r) -(X*Z) - Y)*dt
Z=: Z + ((X*Y) - b*Z)*dt
[A=: A,X
[B=: B,Y
[C=: C,Z
pd 'pensize 1'
pd 'axes 1 1 1'
pd 'type line'
pd 'labels 1'
pd 'viewsize 2'
pd 'title Lorenz'
pd 'color darkblue'
'viewpoint 7.2 1.2 6.4 ;framebackcolor white' plot A; B; C
end.
)

Funcions associated with main. [{{#file: "lorenz.ijs"}} Download script: lorenz.ijs ]

NB.===================================================
NB. associated funcions	
rotz_z_=: ] +/ .* ((cos,sin),:-@sin,cos)@[ , 0 0 1"_

stereo=: 3 : 0
pd 'sub 1 2;new;boxed 0'
pd 'color red'
pd 'viewpoint ',":vp=. _4.6 3.9 1.5
pd 'viewsize 1'
pd A;B;C
pd 'new;boxed 0'
pd 'color red'
NB. pd 'viewpoint ',":1r25p1 rotz_z_ vp
pd 'viewpoint ',":1r50p1 rotz_z_ vp
pd 'viewsize 1'
pd A;B;C
pd 'show'
)

rotator=: 3 : 0
pd 'color green'
pd 'viewpoint ',":vp=. _4.6 3.5 1.4
pd 'viewsize 1'
angle=: 1r24p1
while. angle < 1p1%3 do.
  pd A;B;C
  pd 'viewpoint ',":angle rotz_z_ vp
  pd 'show'
  delay 1
  angle=: angle +(1p1%400)
  end.
'Rotation completed 'printf here=: ''
)

sub=:(<:@] { [)`(([ (<:@] { [) >@{.@]) (<:@] {"1 [) >@{:@])@.(<:@#@$@[ *. #@]>1:)

delay=: 3 : 0
c=: y
c=: 3000 * c  NB. 400 5 sec
while. c > 0 do. c=: c - 1
c
end.
)
Lorenz3.jpg

This is the Stereo View.

GUI input. [{{#file: "lorenz.ijs"}} Download script: lorenz.ijs ]

NB.===================================================
require 'strings printf'

LORENZ=: 0 : 0
pc Lorenz;pn "Lorenz Equation Explorer";
menupop "File";
menu new "&New" "" "" "";
menu open "&Open" "" "" "";
menusep;
menu exit "&Exit" "" "" "";
menupopz;
xywh 165 15 30 11;cc editS edit;
xywh 166 42 30 11;cc editR edit;
xywh 166 67 30 11;cc editB edit;
xywh 161 93 38 11;cc editNumiter edit;
xywh 166 4 32 10;cc static1 static;cn "S nom 10";
xywh 166 31 31 10;cc static2 static;cn "R nom 28";
xywh 167 56 35 10;cc static3 static;cn "B nom 8/3";
xywh 162 82 40 10;cc static4 static;cn "Num. Steps";
xywh 103 95 34 11;cc bnPlot button;cn "Plot";
xywh 9 95 48 11;cc bnRotate button;cn "Rotate 30 deg";
xywh 11 3 34 11;cc bnEndview button;cn "End View";
xywh 11 18 34 11;cc bnSideview button;cn "Side View";
xywh 12 33 34 11;cc bnTopview button;cn "Top View";
xywh 3 54 43 11;cc bnOriginalview button;cn "Original View";
xywh 12 78 42 11;cc bnStereoview button;cn "Stereo View";
xywh 59 1 99 80;cc stDesc static;cn "stDesc";
pas 6 6;pcenter;
rem form end;
)

Lorenz_run=: 3 : 0
wd LORENZ
NB. initialize form here
wd 'set editS 10'
wd 'set editR 28'
wd 'set editB 2.6667'
wd 'set editNumiter 2000'
wd 'set stDesc *',Desc rplc LF;''
wd 'pshow;'
)

Lorenz_close=: 3 : 0
wd'pclose'
)

fire=: 3 : 0
s=: 1>.{.0".editS
r=: 1>.{.0".editR
b=: 1>.{.0".editB
num_inc =. 1>.{.0".editNumiter
lorenz s, r, b, num_inc
)

Lorenz_bnPlot_button=: 3 : 0
fire''
flag=:1
'Plot completed 'printf nhere=: ''
)

Lorenz_bnRotate_button=: 3 : 0
if. flag = 0 do. wdinfo 'Press the Plot button first' return. end.
pd 'pensize 1'
pd 'axes 1 1 1'
pd 'type line'
pd 'labels 1'
pd 'viewsize 2'
pd 'title Lorenz'
pd 'color darkblue'
'viewpoint  _4.6 3.5 1.4 ;framebackcolor white' plot A; B; C
rotator ''
)
Lorenz_bnEndview_button=: 3 : 0
if. flag = 0 do. wdinfo 'Press the Plot button first' return. end.
pd 'pensize 1'
pd 'axes 1 1 1'
pd 'type line'
pd 'labels 1'
pd 'viewsize 2'
pd 'title Lorenz'
pd 'color darkblue'
'viewpoint 1 4 1;framebackcolor white' plot A; B; C
)
Lorenz_bnSideview_button=: 3 : 0
if. flag = 0 do. wdinfo 'Press the Plot button first' return. end.
pd 'pensize 1'
pd 'axes 1 1 1'
pd 'type line'
pd 'labels 1'
pd 'viewsize 2'
pd 'title Lorenz'
pd 'color darkblue'
'viewpoint 4 4 1;framebackcolor white' plot A; B; C
)

Lorenz_bnTopview_button=: 3 : 0
if. flag = 0 do. wdinfo 'Press the Plot button first' return. end.
pd 'pensize 1'
pd 'axes 1 1 1'
pd 'type line'
pd 'labels 1'
pd 'viewsize 2'
pd 'title Lorenz'
pd 'color darkblue'
'viewpoint 1 1 5 ;framebackcolor white' plot A; B; C
)
Lorenz_bnOriginalview_button=: 3 : 0
if. flag = 0 do. wdinfo 'Press the Plot button first' return. end.
pd 'pensize 1'
pd 'axes 1 1 1'
pd 'type line'
pd 'labels 1'
pd 'viewsize 2'
pd 'title Lorenz'
pd 'color darkblue'
'viewpoint 7.2 1.2 6.4 ;framebackcolor white' plot A; B; C
)


Lorenz_bnStereoview_button=: 3 : 0
if. flag = 0 do. wdinfo 'Press the Plot button first' return. end.
pd 'reset'
stereo''
)

Desc=: 0 : 0
A Lorenz equation explorer which demonstrates
 dynamically the unique properties of this strange
attractor.
Different fixed viewpoints, 30 deg rotation,
 and a Stereo view are supported. The three variables
 s r and b are all adjustable for viewing variations in form.
)
Lorenz_run''
NB.lorenz 13 28 2.67 160