전체 글195 [MATLAB 실습 08] ECG 잡음제거 (22.11.16) 1. 300Hz 샘플링 주파수로 심전도를 기록하는 도중에 60Hz 잡음이 섞어 들어갔다. 잡음이 섞인 심전도 데이터(ecg.dat)를 읽어와(load('ecg.dat'):) 60Hz 잡음 신호를 제거를 제거하려고 한다. 필터링 과정을 DFT를 이용하여 제거하려고 한다. (1) 잡음이 낀 심전도 데이터를 600pt DFT를 계산하시오. (2) 60Hz에 해당되는 DFT인덱스를 계산하시오. (3) 60Hz에 해당되는 DFT값을 0으로 값으로 채워 넣으시오. (4) filtered ECG를 IDFT 취하여 원 ECG 신호화 filtered ECG신호를 그리시오. (5) 주파수 영역에서 filtered ECG와 원래의 ECG 신호를 비교하시오, (6) 60Hz 간섭이 제거되었는지 확인하시오. 간섭이 제거되었다... 2023. 6. 14. 09 Responses of Digital Filters (22.11.14) LTI System Characterization ● 임펄스 응답 h[n]을 LTI 필터의 임펄스 응답이라고 하면 컨볼루션을 통해 시스템을 특성화 할 수 있다. h[n]은 LTI 필터의 시간 영역 응답이다. ● 주파수 응답 H(e^jw)는 LTI 필터의 주파수 영역 응답이다. ● Difference Equation LTI 시스템은 아래의 차분 방정식을 충족한다. ● Transfer Function Impulse Response of Digital Filters H(z)가 1차 극점만 있는 z^(-1)의 유리 함수일 때 아래와 같다. 여기서 첫 번째 성분은 M ≥ N인 경우에만 존재한다. 시스템이 causal이면 ROC는 |z| > |max| 이어야 한다. 여기서 pmax는 가장 큰 크기의 극점이다. 이 .. 2023. 6. 14. [MATLAB 실습 07] DFT / FFT (22.11.09) DSP 이론 학습을 마친 후, MATLAB 사용법을 추가적으로 학습하여 DSP 실습을 진행해보았다. 실습은 《MATLAB을 이용한 디지털 신호처리 (저자 Viny K. Ingle, 출판사 Cengage Learning)》 책의 예제와 연습문제로 진행해보았다. [N = 4] [N = 8] [N = 16] [N = 32] N=32일 때 그래프가 부드럽게 이어진다. [N = 40] [N = 50] [N = 60] N이 클수록 분해 능력이 좋아져 입력함수의 스펙트럼을 정확하게 추정할 수 있다. [N = 4] [N = 7] [N = 8] N=7일 때 선형 컨볼루션과 순환 컨볼루션이 같아진다. 2023. 6. 14. 08 Discrete Fourier Series (22.11.07) Discrete Fourier Transform DFT는 주파수 영역에서 이산 시간 유한 기간 신호를 분석하는 데 사용된다. x[n]을 외부 0 ≤ n ≤ N-1이 되는 길이 N의 유한 기간 수열이라고 합니다. DFT 쌍은 아래와 같다. x[n]을 주기 시퀀스 x̄[n]으로 확장하면, X[k] = x̄[k] (0 ≤ k ≤ N-1)이다. 결과적으로 DFT와 DFS는 [0, N-1] 구간 내에서 같다. Properties of DFT DFT 쌍은 [0, N - 1] 내에서 DFS 쌍과 동일하므로 인덱스가 간격 밖에 있을 때, x[n] 및 X[k]의 값을 관리하면 속성이 같다. 1. Linearity 2. Circular Shift of Sequence 3. Duality 4. Symmetry 5. Circ.. 2023. 6. 14. 이전 1 ··· 17 18 19 20 21 22 23 ··· 49 다음
