猫史档案馆

在 Π之城招人啦（）（）中回复

π π

谁说打不出来？

2022-01-08T14:47:38 点赞:0

在【Kitten 4.9.8 新物理介绍】@技术喵怎么还不出新物理到源码图鉴？中回复

听君一席话，如吃一席饭

2022-01-08T14:48:35 点赞:0

在怎么做屏幕滑动中回复

你问问题以前把话说清楚，不然就是有人会也被你整不会了

2022-01-11T17:35:29 点赞:0

在【社区风纪委员】又来啦～中回复

好！

2022-01-11T19:04:34 点赞:0

在 BCM的真正大佬【选出你心目中的大佬】中回复

小小神射手

2022-01-13T19:17:32 点赞:1

在【征稿】征集大佬们上首页的经验与方法中回复

多练

2022-01-13T19:17:58 点赞:0

在编程一渣渣！寻师启示！中回复

你找师傅的目的是什么？如果是为了提升编程能力，方便问不懂的问题的话，我可以做你师父

2022-01-14T12:32:26 点赞:0

在社区观察者#4 中回复

你这就是玩文字游戏，强词夺理了。。。。

2022-01-17T14:12:57 点赞:0

在关修改编程猫不良风气的建议中回复

要多想

2022-01-20T10:39:46 点赞:0

在【小项目】编程猫消息通知器中回复

冷鱼与热猫：那我走？

2022-01-20T20:06:00 点赞:2

在形容你楼下是萌新、大佬还是巨佬（灵感来源：用一句话形容一下你楼下的头像）中回复

萌新

2022-01-21T14:17:40 点赞:0

在【python】混沌系统中回复

常见混沌系统—Lorenz模型

by loa喵ing... （CSDN）

import matplotlib.pyplot as plt
#绘制三维图像
import mpl_toolkits.mplot3d as p3d

'''
Lorenz吸引子生成函数
参数为三个初始坐标，三个初始参数,迭代次数
返回三个一维list
'''
def Lorenz(x0,y0,z0,p,q,r,T):
  #微分迭代步长
  h=0.01
  x=[]
  y=[]
  z=[]
  for t in range(T):
    xt=x0+h*p*(y0-x0)
    yt=y0+h*(q*x0-y0-x0*z0)
    zt=z0+h*(x0*y0-r*z0)
    #x0、y0、z0统一更新
    x0,y0,z0=xt,yt,zt
    x.append(x0)
    y.append(y0)
    z.append(z0)
  return x,y,z

def main():
  #设定参数
  p=10
  q=28
  r=8/3
  #迭代次数
  T=10000
  #设初值
  x0=-16
  y0=-21
  z0=33
  # fig=plt.figure()
  # ax=p3d.Axes3D(fig)
  x,y,z=Lorenz(x0,y0,z0,p,q,r,T)
  ax=plt.subplot(121,projection="3d")
  ax.scatter(x,y,z,s=5)
  ax.set_xlabel('x(t)')
  ax.set_ylabel('y(t)')
  ax.set_zlabel('z(t)')
  ax.set_title('x0=-16 y0=-21 z0=33')
  # plt.axis('off')
  #消除网格

  ax.grid(False)
  #初值微小的变化
  x0=-16
  y0=-21
  z0=33.00001
  xx,yy,zz=Lorenz(x0,y0,z0,p,q,r,T)
  ax=plt.subplot(122,projection="3d")
  ax.scatter(xx,yy,zz,s=5)
  ax.set_xlabel('x(t)')
  ax.set_ylabel('y(t)')
  ax.set_zlabel('z(t)')
  ax.set_title('x0=-16 y0=-21 z0=33.00001')
  ax.grid(False)
  plt.show()

  t=np.arange(0,T)
  plt.scatter(t,x,s=1)
  plt.scatter(t,xx,s=1)
  plt.show()


if __name__=='__main__':
  main()

2022-01-21T18:27:43 点赞:0

在形容你楼下是萌新、大佬还是巨佬（灵感来源：用一句话形容一下你楼下的头像）中回复

神

2022-02-07T14:01:18 点赞:0

在 [编程猫创作者协会]Kitten中锁定视角下的角色运动效果的一种实现方法（360°）中回复

别问我为什么是py代码，问就是冷鱼还没修好上传图片的bug（

2022-02-13T20:46:44 点赞:1

在 [编程猫创作者协会]Kitten中锁定视角下的角色运动效果的一种实现方法（360°）中回复

dddd

2022-02-13T21:21:57 点赞:0

在 [编程猫创作者协会]线性动态规划基本模型1 中回复

ddd