memo

2016年7月18日月曜日

t f B! P L
# coding: utf-8
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt

# データ数
N = 800*2

# サンプリング間隔
dt = 1/800000

# 信号の周波数
f = 1000

t = np.linspace(1,N,N)*dt - dt
y = np.sin(2*np.pi*f*(t))*800

for i in range(y.size):
    if y[i] >100:
        y[i] = 100
    elif y[i] <-100:
        y[i] = -100

noise  =  np.random.randn(y.size)

y=y+noise*10

plt.plot(t,y,"-o", markersize=3)
plt.xlabel("t [s]")

plt.plot(t[20::100],y[20::100],"o", markersize=8)
plt.show()

plt.hist(y[20::100])
plt.show()



Translate

このブログを検索

  • ()
  • ()
もっと見る

QooQ