[2019-2] 디지털신호처리 첩 신호 만들기 과제

첩 신호: 주파수 성분이 정현 함수처럼 진동하는 신호

 

clear all
close all

f_sample = 8000;
amp1 = 1;
t = 0:(1/f_sample):amp1;
y1 = amp1*cos(700*pi*t.*t + 440*pi*t);

Fc = 4;
y1 = vco( sin(2*pi*Fc*t), [400 700], f_sample);

Fmax = f_sample/4;
window_length = 128;  
Nfft = 128;   Nover = 100;

nn = 1:4:fix(length(y1));

[B,F,T] = specgram(y1(nn),Nfft,Fmax,[],Nover);

B1 = B; F1 = F; T1 = T; 
subplot(211)
plot(t(1:2000),y1(1:2000))
axis([0 0.05 -3  3 ])
subplot(221),   
imagesc(T,F,db(B,40));

y1 = y1/max(abs(y1));
sound(y1,f_sample)