python實現(xiàn)周期方波信號頻譜圖
更新時間：2018年07月21日 14:33:26 作者：落葉_小唱
這篇文章主要介紹了python 周期方波信號頻譜圖，小編覺得挺不錯的，現(xiàn)在分享給大家，也給大家做個參考。一起跟隨小編過來看看吧
在學(xué)習(xí)傅里葉變換的時候遇到了求周期方波信號頻譜圖的例子，在書上和網(wǎng)上查閱了一些資料，發(fā)現(xiàn)大都是討論的都是下圖左邊的周期信號的頻譜，課程老師的PPT中也只列出了另一種周期信號頻譜圖的結(jié)論，沒有在進(jìn)行傅里葉變換，自己便根據(jù)定義推導(dǎo)了一遍，貼在這里作記錄和分享之用。
關(guān)于傅立葉級數(shù)展開的另一討論在我的另一篇文章http://m.fzitv.net/article/144194.htm
2016年11月21號更新
在第二個周期方波信號的傅里葉變換里，注意是
轉(zhuǎn)換為sin函數(shù)下為
之前寫錯了，今天更正。
對于這兩種方波信號，我們也可以編程驗證一下，就是用正弦函數(shù)去逼近方波信號，Python的實現(xiàn)代碼如下：
# 分析傅里葉級數(shù)分解之后cos和sin的和項的圖像輸出
from numpy import mgrid,sin,cos,array,pi
from matplotlib.pyplot import plot,show,title,legend,xlabel,ylabel
x = mgrid[0:10:0.02] # 這里類似于MATLAB用冒號產(chǎn)生步長為0.02的序列，但是語法和MATLAB不同
# 下面的這段循環(huán)實現(xiàn)y=sin(x)+sin(3x)+...+sin(19x)
def cos_square():
y1 = 0;
for i in range(0,20,1):
b = (-1)**(i)*cos((2*i+1)*x)/(2*i+1)
y1=b+y1 # 這種求和的方法是從C語言移植過來的
plot(x,y1,'orange',linewidth=0.6)
title('cos_square')
xlabel('Time')
ylabel('Amplitude')
show()
def sin_square():
y2 = 0
for i in range(0,20,1):
b = sin((2*i+1)*x)/(2*i+1)
y2=b+y2 # 這種求和的方法是從C語言移植過來的
plot(x,y2,'g',linewidth=0.6)
title('sin_square')
xlabel('Time')
ylabel('Amplitude')
show()
cos_square()
sin_square()
輸出結(jié)果如下：
由結(jié)果可以看兩種不同的方波級數(shù)展開是正確的
以上就是本文的全部內(nèi)容，希望對大家的學(xué)習(xí)有所幫助，也希望大家多多支持腳本之家。
相關(guān)文章
在本篇文章里小編給大家整理的是關(guān)于Python3爬蟲中Splash的知識總結(jié)內(nèi)容，需要的朋友們可以學(xué)習(xí)參考下。2020-07-07
這篇文章主要介紹了基于python監(jiān)控程序是否關(guān)閉,文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值,需要的朋友可以參考下2020-01-01
這篇文章主要為大家詳細(xì)介紹了PyQt5+requests實現(xiàn)車票查詢工具，具有一定的參考價值，感興趣的小伙伴們可以參考一下2019-01-01
python模塊shutil函數(shù)應(yīng)用示例詳解教程
這篇文章主要為大家介紹了python模塊中shutil函數(shù)的應(yīng)用示例詳解，有需要的朋友可以借鑒參考下，希望能夠有所幫助，祝大家多多進(jìn)步，早日升職加薪2021-10-10
這篇文章主要給大家介紹了關(guān)于如何利用 python 對目錄下的文件進(jìn)行過濾刪除的相關(guān)資料，文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值，需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧。2017-12-12
在本篇文章里小編給大家整理了一篇關(guān)于python集合的新增元素方法整理內(nèi)容，有需要的朋友們可以學(xué)習(xí)下。2020-12-12
?蒙特卡洛就是產(chǎn)生隨機(jī)變量，帶入模型算的結(jié)果，尋優(yōu)方面，本文主要介紹了python 蒙特卡羅模擬法實現(xiàn)，文中通過示例代碼介紹的非常詳細(xì)，具有一定的參考價值，感興趣的小伙伴們可以參考一下2022-03-03
Python實現(xiàn)yaml與json文件批量互轉(zhuǎn)
這篇文章主要為大家詳細(xì)介紹了如何利用Python語言實現(xiàn)yaml與json文件的批量互轉(zhuǎn)，文中的示例代碼講解詳細(xì)，感興趣的小伙伴可以動手嘗試一下2022-07-07
這篇文章主要介紹了如何通過Python實現(xiàn)標(biāo)簽云算法,文中通過示例代碼介紹的非常詳細(xì)，對大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價值,需要的朋友可以參考下2019-07-07
弱引用在很多語言中都存在,最常用來解決循環(huán)引用問題,本文就來和大家一起探索一下python中的弱引用,感興趣的小伙伴可以跟隨小編一起學(xué)習(xí)一下2023-09-09