!Page-2投稿者：SECOND 投稿日：2011年 9月 1日(木)21時00分38秒 |
|
!Page-2 の始め !------------------------------------------ u0 == IDN( 0 階空間の(実,虚,原)各・行vector) ! m0 u0 m0 == f0() ! |re(ta0),im(ta0), 0, 0| |1, 0, 0, 0| ta0 = f0()１階写像空間の実軸単位・行vector ! |re(tb0),im(tb0), 0, 0| |0, 1, 0, 0| tb0 = f0() 〃 〃 〃 虚軸単位・行vector ! | 0 , 0 , X, 0| |0, 0, X, 0| ! |re(to0),im(to0), 0, 1| |0, 0, 0, 1| to0 = f0() 〃 〃 〃 原点位置・行vector ! !----------------------------------------------- mat u=t0*u ! t0 u u ! |re(ta0),im(ta0), 0, 0| |re(a) ,im(a) , 0, 0| |re(a) ,im(a) , 0, 0| ! |re(tb0),im(tb0), 0, 0|*|re(b) ,im(b) , 0, 0|=|re(b) ,im(b) , 0, 0| ! | 0 , 0 , X, 0| | 0 , 0 , X, 0| | 0 , 0 , X, 0| ! |re(to0),im(to0), 0, 1| |re(zo),im(zo), 0, 1| |re(zo),im(zo), 0, 1| ! !----------------------------------------------- draw picture(zo0,zr0,zi0) with u ! picture in 0 u picture in u ! |re(zr0),im(zr0), 0, 1| |re(a) ,im(a) , 0, 0| |re(zr),im(zr), 0, 1| ! |re(zi0),im(zi0), 0, 1|*|re(b) ,im(b) , 0, 0|=|re(zi),im(zi), 0, 1| ! |re(zo0),im(zo0), 0, 1| | 0 , 0 , 0, 0| |re(zo),im(zo), 0, 1| ! |re(zo),im(zo), 0, 1| ! !------------------------------------------------------------------ !４）0 階座標系を、Ｎ階 写像して、その中に、0 階の座標値で絵を描く。 ! ! 各階で、次の階の z0,a,b を、０階～１階の相対座標 to_,ta_,tb_ で送る。 !------------------------------------------------------------------ SUB F23x(k, zo,a,b) ! 元空間の、zo=原点 a=実軸単位ベクトル b=虚軸単位ベクトル IF 0< k THEN !次階空間の、to_=原点 ta_=実軸単位ベクトル tb_=虚軸単位ベクトル LET r$=r$& "1" CALL F23x(k-1, re(to0)*a+im(to0)*b+zo, re(ta0)*a+im(ta0)*b, re(tb0)*a+im(tb0)*b) !次階のzo,a,b LET r$=r$& "2" CALL F23x(k-1, re(to1)*a+im(to1)*b+zo, re(ta1)*a+im(ta1)*b, re(tb1)*a+im(tb1)*b) IF o1=2 THEN LET r$=r$& "3" CALL F23x(k-1, re(to2)*a+im(to2)*b+zo, re(ta2)*a+im(ta2)*b, re(tb2)*a+im(tb2)*b) END IF ELSE LET col=INT(cx)+1 !N 階の、原点・座標軸で、0 階の絵を描く DRAW f012_n( re(zo0)*a+im(zo0)*b+zo, re(zr0)*a+im(zr0)*b+zo, re(zi0)*a+im(zi0)*b+zo ) IF o1=2 THEN LET cx=cx+3^(3-N) ELSE LET cx=cx+2^(5-N) LET no=no+1 END IF LET r$=r$(1:LEN(r$)-1) END SUB !-------------------------------------------------------------- !３）0 階絵(zo0,zr0,zi0)を、N 階写像した座標値(zo,zr,zi)で描く。 !-------------------------------------------------------------- SUB F23(k, zo,zr,zi) !初期値( zo,zr,zi ) は 0 階集合 IF 0< k THEN LET r$=r$& "1" CALL F23(k-1, f0(zo),f0(zr),f0(zi)) !1～N 階目集合 の呼出し LET r$=r$& "2" CALL F23(k-1, f1(zo),f1(zr),f1(zi)) !２分岐型は ココまで IF o1=2 THEN LET r$=r$& "3" CALL F23(k-1, f2(zo),f2(zr),f2(zi)) !３分岐型は ココまで END IF ELSE IF o1=2 THEN LET col=MOD(cx,3^3)+1 ELSE LET col=MOD(cx,2^5)+1 DRAW f012_n( zo,zr,zi ) !N 階集合 LET cx=cx+1 LET no=no+1 END IF LET r$=r$(1:LEN(r$)-1) END SUB !----------- ３角形 △(zo,zr,zi) を描く。zo に * の印と、no#route(path) を付ける。 ! プロッター !----------- PICTURE f012_n(zo,zr,zi) IF Nmax<=N THEN SET POINT COLOR col PLOT POINTS: zo ELSEIF 4<=N THEN SET LINE COLOR col PLOT LINES: zo; zr; zi; zo ELSE SET LINE COLOR N+1 IF N=Ne THEN SET LINE width 2 PLOT LINES: zo; zr; zi SET LINE width 1 !--- SET LINE STYLE 3 PLOT LINES: zo; zi SET LINE STYLE 1 !--- SET POINT STYLE 3 SET POINT COLOR N+1 PLOT POINTS: zo SET POINT STYLE 1 !--- IF N=Ne THEN SET TEXT font "", 9*bms PLOT label,AT zo:STR$(no)& r$ SET TEXT font "",18*bms END IF END IF END PICTURE END !--------------- ! key pannel !--------------- MODULE m MODULE OPTION ARITHMETIC COMPLEX PUBLIC STRING fn$(60) PUBLIC NUMERIC mlb,mrb, bms, o(5) SHARE NUMERIC bx0(5),bxw(5),by0(5),byw(5),Hn(5),Tn(5), bcol(60), chE ! ASK bitmap SIZE i,j !(ask bitmap size)=(ask pixel size)-1 LET bms=MIN(i,j)/500 ! ! ch.1 DATA END , "N=N-1" , "N=N+1" , コッホの曲線 DATA レヴィのＣ曲線, ドラゴン集合 , ペアノの曲線 , 葉脈曲線 DATA カニの行列 , クリヌキ正五角形, ブロッコリー , 高木関数 DATA Extra-1 雲 , Extra-2 雷 , Tragon トラゴン, シルピンスキーの… !ガスケット DATA "" ! ch.2 DATA "N=0","N=1", 2, 3 DATA 4, 5 , 6, 7 DATA 8, 9 , 10, 11 DATA 12, 13 , 14, 15 DATA 16, 17 DATA "" ! ch.3 DATA call F23, call F23x, draw D23, call dp23, call F23mx DATA "" ! LET i=1 DO READ IF MISSING THEN EXIT DO: fn$(i) IF fn$(i)="" THEN EXIT DO LET chE=chE+1 LET o(chE)=i !top pointer( channel) DO LET i=i+1 READ fn$(i) IF fn$(i)="" THEN EXIT DO LOOP LOOP ! CALL m.init( 1, .01,.05, 4,16, .98/ 4,.14/4) !(ch, x0,y0, H,T, xw,yw) CALL m.init( 2, .01,.01, 18,18, .98/18,.14/4) !組ch, top_cell位置(x0,y0), 横数H,総数T, cell幅(xw,yw) CALL m.init( 3, .01,.96, 5, 5, .98/ 5,.14/4) !stop EXTERNAL SUB keyin( ky,col) ASK VIEWPORT a,b,c,d ASK WINDOW e,f,g,h SET VIEWPORT 0,1,0,1 SET WINDOW 0,1,0,1 LET ky=0 DO WHILE mlb=1 CALL sensor(i,col) IF i<>ky THEN IF 0< ky THEN CALL colkey(IP(ky/100), MOD(ky,100), bcol(o(IP(ky/100))+MOD(ky,100))) IF 0< i THEN CALL colkey(IP(i/100) , MOD(i,100) , col) LET ky=i END IF WAIT DELAY 0 LOOP WAIT DELAY .1 IF 0< ky THEN CALL colkey(IP(ky/100), MOD(ky,100), bcol(o(IP(ky/100))+MOD(ky,100))) SET VIEWPORT a,b,c,d SET WINDOW e,f,g,h END SUB EXTERNAL SUB sensor(i,col) MOUSE POLL mx,my,mlb,mrb FOR ch=1 TO chE LET j=INT((mx-bx0(ch))/bxw(ch)) LET i=INT((my-by0(ch))/byw(ch)) *Hn(ch) +j IF 0<=j AND j< Hn(ch) AND 0<=i AND i< Tn(ch) THEN LET i=i+ch*100 EXIT SUB END IF NEXT ch LET i=0 END SUB EXTERNAL SUB colkey(ch, k, col) ASK TEXT HEIGHT w0 LET w0=pixely(w0)-pixely(w0/17) !11pt.w0/9~23 18pt.w0/13~38 ASK TEXT JUSTIFY w1$,w2$ !--- SET TEXT font "MS UI Gothic",11*bms SET TEXT JUSTIFY "center","half" SET TEXT background "transparent" SET AREA COLOR col DRAW box(ch, k) WITH SHIFT( MOD(k,Hn(ch))*bxw(ch)+bx0(ch), INT(k/Hn(ch))*byw(ch)+by0(ch)) SET TEXT background "opaque" SET TEXT JUSTIFY w1$,w2$ SET TEXT font "",w0 !w0 は save/restore の繰返しで、小さめ→消滅、大きめ→発散。 END SUB ! ├ bxw ┤ ! ─┏━━━┓ ! byw┃ ┃ ! by0→┗━━━┛ ! ↑ ! bx0 EXTERNAL PICTURE box(ch, k ) PLOT AREA :0,0; bxw(ch),0; bxw(ch),byw(ch); 0,byw(ch) PLOT LINES :0,0; bxw(ch),0; bxw(ch),byw(ch); 0,byw(ch); 0,0 PLOT TEXT,AT .5*bxw(ch), .5*byw(ch) :fn$(o(ch)+k) END PICTURE EXTERNAL SUB BG_col(ch, k, col) ASK VIEWPORT a,b,c,d ASK WINDOW e,f,g,h SET VIEWPORT 0,1,0,1 SET WINDOW 0,1,0,1 LET bcol(o(ch)+k)=col CALL colkey(ch, k, col) SET VIEWPORT a,b,c,d SET WINDOW e,f,g,h END SUB EXTERNAL SUB init(ch, x0,y0, H,T, xw,yw) LET Hn(ch)=H LET Tn(ch)=T LET bx0(ch)=worldy(pixelx(x0)) LET by0(ch)=worldy(pixely(y0)) LET bxw(ch)=worldy(pixelx(xw))-worldy(pixelx(0)) LET byw(ch)=worldy(pixely(yw))-worldy(pixely(0)) FOR k=0 TO Tn(ch)-1 LET bcol(o(ch)+k)=0 CALL colkey(ch, k, 0) NEXT k END SUB EXTERNAL SUB clear2 ASK VIEWPORT a,b,c,d ASK WINDOW e,f,g,h SET VIEWPORT 0,1,0,1 SET WINDOW 0,1,0,1 CLEAR FOR ch=1 TO chE FOR k=0 TO Tn(ch)-1 LET bcol(o(ch)+k)=0 CALL colkey(ch, k, 0) NEXT k NEXT ch SET VIEWPORT a,b,c,d SET WINDOW e,f,g,h END SUB END MODULE 
|