介绍
波乘法(调制)用于声波和通信领域。这里我们定义了两个不同频率的波，然后我们用python计算并绘制了相乘的波。
定义 2 个不同频率的波
wave1(波1), f1 = 1,000Hz, 在 1ms 区间内画出波1
wave2(波2), f2 = 5,000Hz, 在 1ms 区间内画出波2
wave(调制) = 波1 x 波2, 在 1ms 区间内画出相乘波
导入函示库 numpy 和 matplotlib
## Import library import numpy as np import matplotlib.pyplot as plt # library for plotting
定义画图函式
def drawGraph(xp, yp, titleLabel = 'Draw Graph', xLabel = 'X Axis', yLabel = 'y Axis'):
'''
Draw graph
Parameters:
xp : x axis data
yp : y axis data
'''
# Draw graph data
fig, ax = plt.subplots(nrows=1, ncols=1) # create figure handle
ax.plot(xp, yp)
ax.set_title(titleLabel) # window title
ax.set_xlabel(xLabel) # x-axis label
ax.set_ylabel(yLabel) # y-axis label
#fig.show()
Wave1(波1), f1 = 1kHz, 在 1ms 区间内划出波1
样品数 = F2(waves) ∗ 10(samples/wave) ∗ 0.001(s)
x 轴是一个时间阵列, [t0,t1,...,tn]
y 轴是时间 t 的正弦波函式频率为 f1,[sin(2πf1t0),sin(2πf1t1),...,sin(2πf1tn)]
y=f(t)=sin(2πf1t)
f1 = 1000 f2 = 5000 nBlock = 10 # 10 points/wave dt = 0.001 # 1ms samples = f2 * nBlock * dt # samples in 5ms # generate t array on x axis xTime = np.arange(0, dt+(1/(f2*nBlock)), 1/(f2*nBlock)) # time base # Draw point p(xt,yt) @ 1kHz yValue1 = np.sin(2*np.pi*xTime*f1) # sine wave of f1 drawGraph(xTime, yValue1, "Draw wave at frequency 1kHz", "Time(s)", "Value")
Wave2(波2), f2 = 5,000Hz, 在 1ms 区间内划出波2
样品数 = F2(waves) ∗ 10(samples/wave) ∗ 0.001(s)
x 轴是一个时间阵列, [t0,t1,...,tn]
y 轴是时间 t 的正弦波函式频率为 f2, [sin(2πf2t0),sin(2πf2t1),...,sin(2πf2tn)]
y=f(t)=sin(2πf2t)
# Draw point p(xt,yt) @ 5kHz yValue2 = np.sin(2*np.pi*xTime*f2) # sine wave of f2 drawGraph(xTime, yValue2, "Draw wave at frequency 5kHz", "Time(s)", "Value")
wave(调制) = wave1(波1) x wave2(波2) in 1ms
wave = wave1 x wave2
y 轴是时间 t 的相乘波函式, [sin(2πf1t0)∗sin(2πf2t0),sin(2πf1t1)∗sin(2πf2t1),...,sin(2πf1tn)∗sin(2πf2tn)]
y=f(t)=sin(2πf1t)∗sin(2πf2t)
# Draw modulation point p(xt,yt) yValue = yValue1 * yValue2 # modulation wave = wave1 x wave2 drawGraph(xTime, yValue, "Draw modulation wave", "Time(s)", "Value")
Start writing here...