24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 262  |  回复: 1
当前主题已经存档。

lovein

金虫 (小有名气)

[交流] 【求助】麻烦帮忙看下一个“随机游走”的程序,为什么运行时老出错!!!

整个程序如下:
(*----------------------------------------------------------------------*)
BeginPackage["RW2D`"]
  RW2D::usage = "RW2D[maxstep, ktry, range, startarea, datadirectory]
    reads the two - dimensional digital image file (text file) and
    exports the mean - square displacement of random walkers as a text
    file, 'msdlist2d.out'.
     'maxstep' is the final value of the time step of the random walker.
     'ktry' is the number of average (i.e., number of walkers).
     'range' is the lower and upper limit of the voxel intensity of
    the pore space.
     'startarea' is the dimension of a square as the percentage to
    the system size.
     The random walker starts to diffuse from inside the square which is
    located at the center of the image.
     'datadirectory' is the directory or folder name where a single
    two-dimensional digital image file (text file) is located. "
  ShowTrajectory::usage = "ShowTrajectory[kk] is for the graphic display
    of the random walk trajectory of the 'kk' - th trial. "
Begin["`Private`"]
(* Reading the two - dimensional data in 'datadirectory' *)
readdata2d[range_List, datadirectory_] :=
    Module[{dataname, datasize, x, m2d, rmin = range[[1]], rmax = \
range[[2]]},
      dataname = First[FileNames["*", {datadirectory}]];
      datasize = Round[Sqrt[Length[ReadList[ToString[dataname], Real]]]];
      m2d =
        Partition[
          Map[If[rmin
(* random walk *)
(* 'xylist' is the trajectory of 'maxstep' - step from 'o2d' *)
tra2d[maxstep_, o2d_List, m2d_List] :=
    Module[{n2d = {{1, 0}, {-1, 0}, {0, 1}, {0, -1}}, now = o2d,
        xylist = {o2d}, outoffield = 0, xmax = ymax = Length[m2d],
        step, new, x, y},
      ran2d := Random[Integer, {1, 4}];
      For[i = 1, i
(* random walk of 'ktry' - walkers *)
(* 'tra2dktry' is list of the trajectory of 'ktry' - walkers *)
rw2d[maxstep_, ktry_, startarea_, m2d_List] :=
    Module[{xdatasize = ydatasize = Length[m2d], xylist, outoffield,
        tra2dktry = {}, failure = 0},
      rano2d[lx_, ly_, sa_] :=
        {Random[Integer, {(1 - sa/100)/2 lx, (1 + sa/100)/2lx}],
         Random[Integer, {(1 - sa/100)/2 ly, (1 + sa/100)/2ly}]};
      o2d := (For[i = 1; {x0, y0} = rano2d[xdatasize, ydatasize, startarea],
            m2d[[x0, y0]] != 0, i++,
            {x0, y0} = rano2d[xdatasize, ydatasize, startarea]]; {x0, y0});
      Do[{xylist, outoffield} = tra2d[maxstep, o2d, m2d];
        AppendTo[tra2dktry, xylist]; AddTo[failure, outoffield], {i, 1, \
ktry}];
      Return[{tra2dktry, failure}];
    ];
(* calculation of the mean - square displacement *)
msd2d[data_List] :=
    Module[{ktry = Dimensions[data][[1]], maxstep = Dimensions[data][[2]] - \
1,
         d1, d2, d3, dd, ii},
      d1 = Transpose[data, {2, 1}];
      d2 = Map[(# - First[d1]) &, d1];
      d3 = Transpose[Map[(#^2) &, d2]];
      dd = Apply[Plus, Transpose[Apply[Plus, d3]/ktry]];
      ii = Range[0, maxstep];
      msdlist2d = Transpose[{ii, dd}];
      Return[msdlist2d];
    ];
RW2D[maxstep_, ktry_, range_List, startarea_, datadirectory_] :=
    Module[{failure, msdlist2d},
      (* reading the two - dimensional data in 'datadirectory' *)
      m2d = readdata2d[range, datadirectory];
      (* execution of the random walk *)
      {tra2dktry, failure} = rw2d[maxstep, ktry, startarea, m2d];
      Print["Failure = ", failure, " / Success = ", ktry];
      (* calculation of the mean - square displacement of 'ktry' walkers *)
      msdlist2d = msd2d[tra2dktry];
      (* file output of the mean - square displacement *)
      (* The file name is 'msdlist2d.out' *)
      Export["msdlist2d.out", N[msdlist2d], "Table"];
      (* graphic output of the mean - square displacement *)
      ListPlot[msdlist2d, PlotJoined -> True,
          AxesLabel -> {"time \n step", "mean-square \n displacement"}];
    ];
(* graphic output of the trajectory *)
(* Show the trajectory of a single walker as an example. *)
(* RGBColor[] of 0, 1, 2 are colors of the pore, obstacle, and trajectory,
  respectively. You can change the color as you like. *)
gtra2d[tra2d_List, m2d_List] :=
    Module[{color0 = RGBColor[0.9, 0.9, 0.9],
            color1 = RGBColor[0.5, 0.5, 0.5],
            color2 = RGBColor[0.8, 0.8, 0.8],
            plmin = 1 - 2, plmax = Length[m2d] + 1},
      tra2dshift = Map[(# - {0.5, 0.5}) &, tra2d];
      tra2dm2d = (aa2 = m2d; Map[(aa2[[#[[1]], #[[2]]]] = 2) &, tra2d]; aa2);
      Show[
        Graphics[
          RasterArray[tra2dm2d /. {0 -> color0, 1 -> color1, 2 -> color2}]],
        Graphics[{Line[Map[Reverse, tra2dshift]], AbsolutePointSize[5.0],
            Point[Reverse[First[tra2dshift]]],
            Point[Reverse[Last[tra2dshift]]]}],
        PlotRange -> {{plmin, plmax}, {plmin, plmax}}, AspectRatio -> 1];
    ];
(* Only for the 'kk' - th trial, the trajectory of the random walk is shown
  graphically *)
ShowTrajectory[kk_] := gtra2d[tra2dktry[[kk]], m2d];
End[];
EndPackage[];
(*----------------------------------------------------------------------*)
(* example of usage *)
(*----------------------------------------------------------------------*)
(* As noted below, first of all, you put the command '<< RW2D.m'. *)
(* Then you can put the command set,
  'RW2D[500, 200, {0, 34000}, 10, "DATAsingle"];' and 'ShowTrajectory[2];'. \
*)
(* 'RW2D[maxstep, ktry, range, startarea, datadirectory]' is the command for
  the performance of the random walk and the calculation of mean - square
  displacement. *)
(* The command 'ShowTrajectory[kk]' is for the graphic display of the
  trajectory of the 'kk' - th trial. *)
(* In the example below, a trajectory of the second trial will be shown
  graphically. *)
(* In[1] := << RW2D.m *)
(* In[2] := RW2D[500, 200, {0, 34000}, 10, "DATAsingle"]; ShowTrajectory[2]; \
*)
不知道为啥,运行的时候老是出错,我又看不大懂这程序,请求帮忙

[ Last edited by lovein on 2009-10-26 at 13:35 ]
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lovein

金虫 (小有名气)

附件是程序中需要调用的图像,即the two - dimensional digital image file (text file)

太长了,只能用附件发。还望帮忙。
灰常灰常感谢!!
2楼2009-10-26 13:39:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lovein 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 11408软件工程求调剂 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 21:50 by zhq0425
[考研] 083000学硕274求调剂 +8 Li李鱼 2026-03-26 8/400 2026-03-28 20:33 by 加油向未来啊
[考研] 食品工程专硕一志愿中海洋309求调剂 +4 小张zxy张 2026-03-26 8/400 2026-03-28 19:25 by lbsjt
[考研] 286求调剂 +12 PolarBear11 2026-03-26 12/600 2026-03-28 12:14 by zllcz
[材料工程] 一志愿C9材料与化工专业总分300求调剂 +8 曼111 2026-03-24 9/450 2026-03-28 07:58 by YYYYX1234
[考研] 331环境科学与工程求调剂 +3 熠然好运气 2026-03-27 3/150 2026-03-28 04:11 by fmesaito
[考研] 266分求材料化工冶金矿业等专业的调剂 +4 哇呼哼呼哼 2026-03-26 4/200 2026-03-27 17:02 by zhyzzh
[考研] 一志愿南师大0703化学 275求调剂 +4 Ripcord上岸 2026-03-27 4/200 2026-03-27 17:00 by zhyzzh
[考研] 283求调剂(080500) +4 A child 2026-03-27 4/200 2026-03-27 15:34 by XPU李庆
[考研] 求调剂 +3 刘柯@ 2026-03-24 4/200 2026-03-27 11:28 by shangxh
[考研] 276求调剂。有半年电池和半年高分子实习经历 +10 材料学257求调剂 2026-03-23 11/550 2026-03-27 10:13 by YCIT- LHL
[考研] 材料调剂 +8 匹克i 2026-03-23 8/400 2026-03-27 08:11 by hypershenger
[考研] 化学工程085602 305分求调剂 +17 RichLi_ 2026-03-25 17/850 2026-03-26 19:44 by plmuchong
[考研] 085602化学工程求调剂。 +4 平乐乐乐 2026-03-26 4/200 2026-03-26 17:57 by fmesaito
[考研] 材料考研求调剂 +3 Dendel 2026-03-23 6/300 2026-03-26 17:51 by fmesaito
[考研] 085602 289分求调剂 +8 WWW西西弗斯 2026-03-24 8/400 2026-03-26 16:33 by 不吃魚的貓
[考研] 求b区院校调剂 +4 周56 2026-03-24 5/250 2026-03-25 17:12 by yishunmin
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[考研] 344求调剂 +3 desto 2026-03-24 3/150 2026-03-24 10:09 by 搏击518
[考研] 石河子大学(211、双一流)硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
信息提示
请填处理意见