python 실현 주기 방 파 신호 스펙트럼 도

푸 리 엽 의 변환 을 공부 할 때 주기 적 인 방 파 신호 스펙트럼 도 를 구 하 는 예 를 만 났 다.책 과 인터넷 에서 일부 자 료 를 찾 아 보 니 대부분이 다음 그림 왼쪽 의 주기 신호 스펙트럼 을 토론 하 는 것 을 발견 했다.과정 선생님 의 PPT 에서 도 다른 주기 신호 스펙트럼 도 에 대한 결론 만 나 왔 고 푸 리 엽 의 변환 을 하지 않 았 기 때문에 자신 은 정의 에 따라 한 번 유도 했다.여기에 붙 여 기록 하고 공유 하 는 용도 로 쓰 인 다.
푸 립 엽 급수 에 관 한 또 다른 토론 은 나의 다른 글 에 있다.

2016 년 11 월 21 일 업데이트
두 번 째 주기 방 파 신호 의 푸 리 엽 변환 에서 주의 하 는 것 은?

sin 함수 로 변환

전에 잘못 썼 으 니 오늘 정정 하 겠 습 니 다.
이 두 가지 방 파 신호 에 대해 우 리 는 프로 그래 밍 을 통 해 검증 할 수 있다.바로 사인 함수 로 방 파 신호 에 접근 하 는 것 이다.Python 의 실현 코드 는 다음 과 같다.

#            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        0.02   ,     MATLAB  



#          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()
출력 결 과 는 다음 과 같 습 니 다.


결 과 를 보면 두 가지 서로 다른 방 파 급수 가 전개 되 는 것 이 정확 하 다.
이상 이 바로 본 고의 모든 내용 입 니 다.여러분 의 학습 에 도움 이 되 고 저 희 를 많이 응원 해 주 셨 으 면 좋 겠 습 니 다.

좋은 웹페이지 즐겨찾기