小白提问~~~

Simoncolosi

在numpy和matplotlib.pyplot模块里面这一段是什么意思啊求大神解

y = np.sin(2*np.pi*ff*t)

jinlovelive

把大概的一段代码贴出来啊...现在只能看出是在计算Y轴的值，2π后面没代码不清楚

Simoncolosi

jinlovelive 发表于 2019-9-5 16:26
把大概的一段代码贴出来啊...现在只能看出是在计算Y轴的值，2π后面没代码不清楚

import matplotlib.pyplot as plt

import numpy as np

import seaborn
Fs =120; # sampling rate采样率

Ts = 1.0/Fs; # sampling interval 采样区间
t = np.arange(0,1,Ts) # time vector,这里Ts也是步长
ff1 = 10; # frequency of the signal 采样频率
ff2 = 50
y = 10*np.sin(2*np.pi*ff1*t)+10*np.sin(2*np.pi*ff2*t)
n = len(y) # length of the signal 信号长度
k = np.arange(n)
T = n/Fs
frq = k/T # two sides frequency range 两侧频率范围
frq1 = frq[range(int(n/2))] # one side frequency range单侧频率范围
YY = np.fft.fft(y)
Y = np.fft.fft(y)/n # fft computing and normalization 快速傅立叶变换计算与归一化
Y1 = Y[range(int(n/2))]
fig, ax = plt.subplots(4, 1)
ax[0].plot(t,y)
ax[0].set_xlabel('Time')
ax[0].set_ylabel('Amplitude')
ax[1].plot(frq,abs(YY),'r') # plotting the spectrum 绘制光谱
ax[1].set_xlabel('Freq (Hz)')
ax[1].set_ylabel('|Y(freq)|')
ax[2].plot(frq,abs(Y),'G') # plotting the spectrum
ax[2].set_xlabel('Freq (Hz)')
ax[2].set_ylabel('|Y(freq)|')
ax[3].plot(frq1,abs(Y1),'B') # plotting the spectrum
ax[3].set_xlabel('Freq (Hz)')
ax[3].set_ylabel('|Y(freq)|')
plt.show()
其实我更想问这一段红色的 大佬求讲解。。。

jinlovelive

哟，通信工程的啊
y=10sin(2pi*ff1*t)+10*sin(2pi*ff2*t)
每次采样的值为这两个信号在时间t的值的和
红色那段代码就是画出坐标图，比如：
ax[0].plot(t,y)
ax[0].set_xlabel('Time')
ax[0].set_ylabel('Amplitude')
建立坐标轴ax中的第一条线即画出(t,y)的对应关系，X轴为t，标签是Time,Y轴为y标签为Amplitude，默认黑色
下面以此类推，'G','B','R'分别代表绿色、蓝色、红色

2369520493

这也算小白吗？

Simoncolosi

jinlovelive 发表于 2019-9-5 17:00
哟，通信工程的啊
y=10sin(2pi*ff1*t)+10*sin(2pi*ff2*t)
每次采样的值为这两个信号在时间t的值 ...

通信倒不是这是我用的别人的代码在理解和询问呢~哈哈 我是在做关于故障检测的。还是个本科菜鸟。
谢谢大佬的回答

Simoncolosi

2369520493 发表于 2019-9-5 17:13
这也算小白吗？

这是我在看别人的代码呀，哥哥 我真是小白
我觉得要学习的话也一定要看别人的代码才行

jinlovelive

Simoncolosi 发表于 2019-9-5 17:43
这是我在看别人的代码呀，哥哥 我真是小白
我觉得要学习的话也一定要看别人的代码才行

搞不懂你要学matlab还是python了

Simoncolosi

jinlovelive 发表于 2019-9-5 17:50
搞不懂你要学matlab还是python了

我是想学Python因为都说很好啊。反正Python有这方面的功能的。而且社区环境这么好

jinlovelive

Simoncolosi 发表于 2019-9-5 17:56
我是想学Python因为都说很好啊。反正Python有这方面的功能的。而且社区环境这么好

欢迎加入python

Simoncolosi

jinlovelive 发表于 2019-9-5 18:43
欢迎加入python

我要熬到小甲鱼变成老甲鱼

dragov

Simoncolosi 发表于 2019-9-5 22:57
我要熬到小甲鱼变成老甲鱼

应该熬到甲鱼汤……