HPO:機密日誌

■SYNCモデルにショートカットを加える

ただいま無限ループ中．．．．

無限ループは脱出したものの全然ショートカットの効果が感じられない。なぜ？

' SYNC model SW sim program

' started 2005/05/14

' by Hideki Hirayama

#prompt

'#console

'初期設定

dim table_original[21,21] as single

dim table_temp[21,21] as single

dim retina_model[3,3] as single

dim i,j,k,l,dummy,yes_no as integer

dim x,y,z as single

dim mean_table,var_table as single

dim sw_table[20,20,2]

'table_original：一様乱数初期化サブルーチン

sub randam_table()

dim mean_value, var_value as single '乱数の最大値

input "テーブル初期値の平均値";mean_value

input "テーブル初期値の分散";var_value

for i=0 to 21

for j=0 to 21

if i=0 or j=0 then

table_original[i,j] = 0

else

table_original[i,j] = bell_rand(mean_value,var_value)

endif

next

next

endsub

'table_original：特定の値による初期化サブルーチン（メインで選択）

sub specified_table()

for i=0 to 21

for j=0 to 21

if i>10 and i<13 and j>10 and j<13 then

table_original[i,j] = 3

else

table_original[i,j] = 0

endif

next

next

endsub

' フィルター部分の定義：正方格子

sub fil_initial()

dim periphe_value, center_divident as single 'フィルターの周辺値、中央の分母

input "フィルターの周辺値";periphe_value

input "漸増させる分母";center_divident

for i=1 to 3

for j=1 to 3

retina_model[i,j]=periphe_value

next

next

retina_model[2,2]=center_divident

endsub

'正規分布乱数サブルーチン：引数c1,c2は、一様乱数

'[参照]wikipedia：乱数列:

'http://ja.wikipedia.org/wiki/%E4%B9%B1%E6%95%B0#.E6.AD.A3.E8.A6.8F.E4.B9.B1.E6.95.B0

function bell_rand(mm,vv) as single

'mm=平均値、vv=分散

dim d,e as single

d=Sqr(-2 * Log(1-Rnd())/Log(10))

e=2*3.1415926535897932*Rnd()

bell_rand=d*Cos(e)*vv+mm

't=sqrt(-2.0 * log(1-rnd()));

'u=2*PI*rnd();

'r1=t*cos(u) ; r2=t*sin(u);

end function

'スモールワールドリンク生成サブルーチン?

'つなぎ変えでなく、ショートカットを追加する

'とりあえず最大でもひとつのセルにひとつリンク

'相互リンクとする

'sw_table[x_n,y_n,2]=対象セルとリンクしているセルのｘよｙ（１＝ｘ、２＝ｙ）

sub sw_link(n) as single

'n=生成するリンクの数

dim x_n1,y_n1,x_n2,y_n2 as single

for i=1 to n

do

do

x_n1=Int(Rnd()*20+1)

loop until x_n1>0 and x_n1<21

do

y_n1=Int(Rnd()*20+1)

loop until y_n1>0 and y_n1<21

if sw_table[x_n1,y_n1,1]=0 then

do

do

x_n2=Int(Rnd()*20+1)

loop until x_n2>0 and x_n2<21

do

y_n2=Int(Rnd()*20+1)

loop until y_n2>0 and y_n2<21

loop while x_n1=x_n2 and y_n1=y_n2

print "(";x_n1;",";y_n1;")-(";x_n2;",";y_n2;")"

sw_table[x_n1,y_n1,1]=x_n2

sw_table[x_n1,y_n1,2]=y_n2

sw_table[x_n2,y_n2,1]=x_n1

sw_table[x_n2,y_n2,2]=y_n1

exit do

end if

loop

next

end sub

'table_original：表示サブルーチン

Sub p_table()

cls

for i=1 to 20

for j=1 to 20

x = Int(table_original[i,j])

Locate i*3+5,j+5

color x

print x;

color 7

next

next

endsub

'table_original：時間を進める、フィルターを一度かける

sub dev_table()

dim mean_temp,var_temp as single

mean_table=0

var_table=0

for i=1 to 20

for j=1 to 20

table_temp[i,j]=0

next

next

for i=1 to 20

for j=1 to 20

table_temp[i,j]=table_original[i,j]+(8-table_original[i,j])/retina_model[2,2]+table_temp[i,j]

If table_original[i,j]>=7 then

table_temp[i,j]=-10

table_temp[i-1,j-1]=table_temp[i-1,j-1]+retina_model[1,1]

table_temp[i-1,j]=table_temp[i-1,j]+retina_model[1,2]

table_temp[i-1,j+1]=table_temp[i-1,j+1]+retina_model[1,3]

table_temp[i,j-1]=table_temp[i,j-1]+retina_model[2,1]

table_temp[i,j+1]=table_temp[i,j+1]+retina_model[2,3]

table_temp[i+1,j-1]=table_temp[i+1,j-1]+retina_model[3,1]

table_temp[i+1,j]=table_temp[i+1,j]+retina_model[3,2]

table_temp[i+1,j+1]=table_temp[i+1,j+1]+retina_model[3,3]

if sw_table[i,j,1]>0 then

table_temp[sw_table[i,j,1],sw_table[i,j,2=table_temp[sw_table[i,j,1],sw_table[i,j,2+retina_model[3,3]

' print "(";i;",";j;")-(";sw_table[i,j,1];",";sw_table[i,j,2];")"

end if

End If

next

next

for i=1 to 20

for j=1 to 20

if table_temp[i,j]<0 then

table_original[i,j]=0

else

table_original[i,j]=table_temp[i,j]

endif

mean_temp=mean_temp+table_original[i,j]

next

next

mean_table=mean_temp/400

for i=1 to 20

for j=1 to 20

var_temp=var_temp+(table_original[i,j]-mean_table)^2

next

next

var_table=Sqr(var_temp/400)

endsub

'メイン

'　初期化はどちらか選択

Randomize

do

print "初期値はランダムにとりますか？サブルーチン指定の特定の値にしますか？"

input "（0=ランダム,それ以外特定）";yes_no

if yes_no=0 then

randam_table()

else

specified_table()

end if

do

print "ランダムリンクを加えますか？"

input "（0=生成しない,1以上200未満＝生成するリンクの数）";yes_no

loop until yes_no>=0 and yes_no<200

if yes_no>0 then

sw_link(Int(yes_no))

end if

fil_initial()

p_table()

print

print

input "試行数";k

for l=1 to k

dev_table()

p_table()

print

print

print "試行回数＝";l;"　平均値＝";mean_table;"　分散値＝";var_table

' Sleep(500)

next

input "終了しました。続けますか？（1＝終了）";dummy

loop until dummy=1

end

ツイートする

Permalink | コメント(0) | トラックバック(0)