python 波形图_Python实现正弦信号的时域波形和频谱图示例【基于matplotlib】

本文实例讲述了Python实现正弦信号的时域波形和频谱图。分享给大家供大家参考,具体如下:

# -*- coding: utf-8 -*-

# 正弦信号的时域波形与频谱图

import numpy as np

import matplotlib.pyplot as pl

import matplotlib

import math

import random

row = 4

col = 4

N = 500

fs = 5

n = [2*math.pi*fs*t/N for t in range(N)] # 生成了500个介于0.0-31.35之间的点

# print n

axis_x = np.linspace(0,3,num=N)

#频率为5Hz的正弦信号

x = [math.sin(i) for i in n]

pl.subplot(221)

pl.plot(axis_x,x)

pl.title(u'5Hz的正弦信号',fontproperties='SimHei')

pl.axis('tight')

#频率为5Hz、幅值为3的正弦+噪声

x1 = [random.gauss(0,0.5) for i in range(N)]

xx = []

#有没有直接两个列表对应项相加的方式??

for i in range(len(x)):

xx.append(x[i]*3 + x1[i])

pl.subplot(222)

pl.plot(axis_x,xx)

pl.title(u'频率为5Hz、幅值为3的正弦+噪声',fontproperties='SimHei')

pl.axis('tight')

#频谱绘制

xf = np.fft.fft(x)

xf_abs = np.fft.fftshift(abs(xf))

axis_xf = np.linspace(-N/2,N/2-1,num=N)

pl.subplot(223)

pl.title(u'频率为5Hz的正弦频谱图',fontproperties='SimHei')

pl.plot(axis_xf,xf_abs)

pl.axis('tight')

#频谱绘制

xf = np.fft.fft(xx)

xf_abs = np.fft.fftshift(abs(xf))

pl.subplot(224)

pl.title(u'频率为5Hz的正弦频谱图',fontproperties='SimHei')

pl.plot(axis_xf,xf_abs)

pl.axis('tight')

pl.show()

运行效果:

希望本文所述对大家Python程序设计有所帮助。

  • 0
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值