博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
随机漫步
阅读量:6835 次
发布时间:2019-06-26

本文共 1871 字,大约阅读时间需要 6 分钟。

hot3.png

随机漫步

from random_walk import RandomWalk

#!/usr/bin/env python# -*- coding:utf-8 -*-from random import choicefrom matplotlib import pyplot as pltclass RandoWalk():    '''一个生成随机漫步数据的类'''    def __init__(self, num_points=5000):        '''        初始化随机漫步的属性        :param num_points:        '''        self.num_points = num_points        # 所有随机漫步都始于(0,0)        self.x_values = [0]        self.y_values = [0]    def fill_walk(self):        '''计算随机漫步包含的所有点'''        # 不断漫步,直到列表达式到指定长度        while len(self.x_values) < self.num_points:            # 决定前进方向以及沿这个方向前进的距离            x_direction = choice([1, -1])            x_distance = choice([0, 1, 2, 3, 4])            x_step = x_direction * x_distance            y_direction = choice([1, -1])            y_distance = choice([0, 1, 2, 3, 4])            y_step = y_direction * y_distance            # 拒绝原地踏步            if x_step == 0 and y_step == 0:                continue            # 计算下一个点的x和y值            next_x = self.x_values[-1] + x_step            next_y = self.y_values[-1] + y_step            self.x_values.append(next_x)            self.y_values.append(next_y)while True:    # 创建一个RandomWalk实例,并将其包含的点都绘制出来    rw = RandoWalk(20000)    rw.fill_walk()    # 设置绘图窗口尺寸    # plt.figure(figsize=(10, 6))  # 单位英寸    plt.figure(dpi=128, figsize=(10, 6))  # 分辨率    point_numbers = list(range(rw.num_points))    plt.scatter(rw.x_values, rw.y_values, c=point_numbers, cmap=plt.cm.Blues, edgecolors='none', s=5)    # 突出起点和终点    plt.scatter(0, 0, c='green', edgecolors='none', s=100)    plt.scatter(rw.x_values[-1], rw.y_values[-1], c='red', edgecolors='none', s=100)    # 隐藏坐标轴    plt.axes().get_xaxis().set_visible(False)    plt.axes().get_yaxis().set_visible(False)    plt.show()    keep_running = raw_input("Make another walk?(y/n):")    if keep_running == 'n':        break

转载于:https://my.oschina.net/u/3646781/blog/1586436

你可能感兴趣的文章
zabbix专题:第十二章 zabbix proxy分布式监控配置
查看>>
为什么总觉得自己不适合搞IT?
查看>>
vmware克隆server2008R2造成SID冲突
查看>>
python调用zabbix api接口实时展示数据
查看>>
VMware下Windows2003R2虚拟机磁盘扩容方法
查看>>
运维经验分享(六)-- 深究crontab不能正确执行Shell脚本的问题(二)
查看>>
利用Linux的文件命名规范在Windows中建立“高权限”文件
查看>>
失败者的共同特点
查看>>
Tokyo Tyrant基本规范(4)--协议
查看>>
【Go语言】【14】GO语言的接口类型
查看>>
配置CAS应用客户端
查看>>
摘抄--apache工作模式详解
查看>>
更改sybase下设备名
查看>>
不少朋友在安装IDES 4.71的过程中都遇到了下面的出错提示:
查看>>
企业的人性和狼性
查看>>
mySQL教程 第10章 事务和锁
查看>>
Hello, Kafka World
查看>>
Exchange 2010和Exchange 2016共存部署-10:配置多域名证书
查看>>
SFB 项目经验-03-共存迁移-Lync 2013-TO-SFB 2015-完成
查看>>
F5 配置手册 -F5 BIG-IP 10.1-2-配置-基本参数
查看>>