数据挖掘第一次实验

import random
from matplotlib import pyplot
#蒙特卡洛算法

class MTKL:
    def __init__(self, n, m):
        self.n = n
        self.m = m

    def MC(self):
        n = self.n
        m = self.m
        for i in range(m):

            # 在(0,1)之间去随机数
            x = random.random()
            y = random.random()

            if ((x * x * x) - y >= 0):
                n += 1

        return (n / m)

    def mc(self, n, m):
        for i in range(m):
            # 在0,1之间去随机数
            x = random.random()
            y = random.random()
            if ((x * x * x) - y >= 0):
                n += 1
        s = n / m

        return (s)

    def draw_MC(self):
        xx = []
        yy = []
        for i in range(self.m):
            x = random.random()
            y = random.random()

            if ((x * x * x) >= y):
                xx.append(x)
                yy.append(y)

        pyplot.plot(xx, yy)
        pyplot.show()

test = MTKL(0, 10000)

test.MC()
print('结果:')
print(test.MC())
print(test.draw_MC())

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值