网页资讯视频图片知道文库贴吧地图采购
进入贴吧全吧搜索

 
 
 
日一二三四五六
       
       
       
       
       
       

签到排名:今日本吧第个签到,

本吧因你更精彩,明天继续来努力!

本吧签到人数:0

一键签到
成为超级会员,使用一键签到
一键签到
本月漏签0次!
0
成为超级会员,赠送8张补签卡
如何使用?
点击日历上漏签日期,即可进行补签。
连续签到:天  累计签到:天
0
超级会员单次开通12个月以上,赠送连续签到卡3张
使用连续签到卡
10月25日漏签0天
mathematica吧 关注:19,908贴子:74,148
  • 看贴

  • 图片

  • 吧主推荐

  • 游戏

  • 6回复贴,共1页
<<返回mathematica吧
>0< 加载中...

请问我为何无法重现这个程序?

  • 只看楼主
  • 收藏

  • 回复
  • temp_no_2
  • 大写小写
    3
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
在英文wiki上看到这张图,想照猫画虎自己写一个却发现一运行就耗尽我E5机子的16G内存然后卡死。
但是看图上好像很快就能跑出结果,请问是我哪里写错了么?附上原图和我写的版本。谢谢了。度盘1bnpNfFx


  • temp_no_2
  • 大写小写
    3
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
贴一下代码吧,不知道会不会被度娘改。
bifurcate[f_, a0_, k0_, k_, while_: (True &)] := NestWhileList[f, a0, while, 2, k0+k-1] // Drop[#, k0] &
logistic = {r, y} \[RightTeeArrow] r y (1-y);
Row@{"Initial value: ", a0 = RandomReal[{0.1, 0.9}]}
Row@{"Points per r: ", density = 10^2}
Row@{"Initial k: ", k0 = 10^4}
Row@{"Time taken: ", Timing[
plotData=ParallelTable[
{ConstantArray[r, density], bifurcate[logistic[r, #] &, a0, k0, density]}\[Transpose],
{r, 3.5, 4, 0.0001}] // Flatten[#, 1] &;][[1]], " s"}
Row@{"Data length: ", plotData // Length}


2025-10-25 01:46:50
广告
不感兴趣
开通SVIP免广告
  • 无影东瓜
  • 小吧主
    12
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
bifurcate[f_, a0_, k0_, k_, while_: (True &)] :=
NestWhileList[f, a0, while, 2, k0 + k - 1] // Drop[#, k0] &
logistic[r_, y_] := r y (1 - y);
Row@{"Initial value: ", a0 = RandomReal[{0.1, 0.9}]}
Row@{"Points per r: ", density = 10^2}
Row@{"Initial k: ", k0 = 10^4}
Row@{"Time taken: ",
Timing[plotData =
ParallelTable[{ConstantArray[r, density],
bifurcate[logistic[r, #] &, a0, k0,
density]}\[Transpose], {r, 3.5, 4, 0.0001}] //
Flatten[#, 1] &;][[1]], " s"}
Row@{"Data length: ", plotData // Length}
ListPlot[plotData, PlotStyle -> {PointSize[0], Opacity[.1]},
ImageSize -> 800, PlotRange -> {All, {0, 1}}]

你的问题出在logistic的定义,截图里的那个 \[RightTeeArrow] 符号在前文有定义,所以在你发的截图里这一行的字母都是绿色的,而你的代码这一行则是蓝色的
把这里直接改成函数的定义,或者在前面加一行
a_ \[RightTeeArrow] b_ := Function[a, b]
都可以,这里为了直观就直接用函数定义了
另外在用了并行的情况下这个点的计算是很快的,只是由于点太多所以ListPlot画图需要花好久


  • 无影东瓜
  • 小吧主
    12
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
多谢来自@xycr的指正
那个符号其实不是\[RightTeeArrow],而是\[Function],是Function的中缀形式,输入方式可以是Esc+fn+Esc或者\[Function]

一长一短,差距还是挺明显的


登录百度账号

扫二维码下载贴吧客户端

下载贴吧APP
看高清直播、视频!
  • 贴吧页面意见反馈
  • 违规贴吧举报反馈通道
  • 贴吧违规信息处理公示
  • 6回复贴,共1页
<<返回mathematica吧
分享到:
©2025 Baidu贴吧协议|隐私政策|吧主制度|意见反馈|网络谣言警示