/Python/Pseudo-random input signal python

Pseudo-random input signal python

import numpy as np
import matplotlib.pyplot as plt

nstep = 300

# random signal generation

a_range = [0,2]
a = np.random.rand(nstep) * (a_range[1]-a_range[0]) + a_range[0] # range for amplitude

b_range = [2, 10]
b = np.random.rand(nstep) *(b_range[1]-b_range[0]) + b_range[0] # range for frequency
b = np.round(b)
b = b.astype(int)

b[0] = 0

for i in range(1,np.size(b)):
    b[i] = b[i-1]+b[i]

# Random Signal
i=0
random_signal = np.zeros(nstep)
while b[i]<np.size(random_signal):
    k = b[i]
    random_signal[k:] = a[i]
    i=i+1

# PRBS
a = np.zeros(nstep)
j = 0
while j < nstep:
    a[j] = 5
    a[j+1] = -5
    j = j+2

i=0
prbs = np.zeros(nstep)
while b[i]<np.size(prbs):
    k = b[i]
    prbs[k:] = a[i]
    i=i+1

plt.figure(0) 
plt.subplot(2,1,1)
plt.plot(random_signal, drawstyle='steps',label='Random Signal')
plt.legend()
plt.subplot(2,1,2)
plt.plot(prbs, drawstyle='steps', label='PRBS')
plt.legend()
plt.show()

Source: stackoverflow.com

Similar codes

0 vote105 views1 answer
Django auth user model
Python
11 votes132 views4 answers
Invert a dictionary python
Python
8 votes109 views3 answers
Python check if string is int
Python
71 votes68 views10 answers
Sus sus sus
Python
0 vote100 views1 answer
Print convert hex as string ascii
Python
1 vote117 views1 answer
X= [10] def List_ex(): x.append(20) def add_list(): x=[30,40] x.append(50) print (x) List_ex() print (x) add_list() print (x)
Python
0 vote131 views1 answer
DJANGO model instance get by variable
Python Django
0 vote73 views1 answer
Calendly python
Python
0 vote111 views1 answer
Tf.io path copy
Python
1 vote106 views1 answer
Python tipi array
Python