Как напечатать амплитудно-модулированный сигнал с помощью Python - PullRequest
0 голосов
/ 26 июня 2019

У меня есть пример matlab для этой программы, но я не могу сделать это на Python. Как это должно выглядеть (пример mathlab) https://imgur.com/a/oam3jXl.

import numpy as np
import matplotlib.pyplot as plt
from scipy.signal import hilbert, chirp
duration = int(input("Input duration of signal = "))
A = int(input("Input amplitude of signal = "))
start = int(input("Input start of modulation time = "))
end = int(input("Input lenght of modulation = "))
fs = 400.0
samples = int(fs*duration)
t = np.arange(samples) / fs
main_t = []
result = []
counter=0
for l in range(0,samples):
    main_t.append(0)
    result.append(0)
start_mod = np.arange(start*fs)
end_mod = np.arange(end*fs)
signal = chirp(t, 20.0, t[-1], 100.0)
signal *= (A + 0.5 * np.sin(2.0*np.pi*3.0*t) )

for i in np.arange(0,start*fs):
    signal.insert(i,0)

fig = plt.figure()
ax0 = fig.add_subplot(211)
ax0.plot(main_t, result)
ax0.set_xlabel("time in seconds")
plt.show()

Я ожидаю ввода некоторых данных из консоли, а затем распечатать график, например, с помощью matplotlib. График должен выглядеть как график, который вы можете видеть на изображении (пример mathlab).

1 Ответ

0 голосов
/ 26 июня 2019

Игнорируя все беспорядки с получением ввода (который будет сильно различаться в зависимости от других ограничений), код для генерации амплитудно-модулированного сигнала должен быть простым.

Я начинаю с натяжения numpy, генерируя набор точек, по которым нужно сэмплировать сигнал, вычисляя амплитуду, затем объединяю их:

import numpy as np

x = np.linspace(0, 10, 501)
ampl = np.exp(-(x - 3.5)**2 / 0.8)
y = np.sin(x * 25) * ampl

затем мы можем построить их используя matplotlib, например:

import matplotlib.pyplot as plt

plt.figure(figsize=(10,5))
plt.plot(x, y, label='signal')
plt.plot(x, ampl, ':', label='amplitude')
plt.xlabel('time')
plt.ylabel('value')
plt.legend()

demo plot

Я вытащил морского рожка и использую их стиль ticks, чтобы сделать его немного красивее.

...