clc;
clear;
%generate signal
fs = 10; %采样频率
N = 10000; %数据总长度
n = 0:N-1;
t = n/fs; %时间序列
for i=1:N
if t(i)<500
x(i)=sin(10*t(i));
else
x(i)=sin(20*t(i));
end
end
x;
subplot(411)
plot(t,x);
axis([0,1000,-1,1]);
%fourier Transform
y_FFT = fft(x,N); %对信号快速Fourier变换
mag1 = abs(y_FFT); %Fourier变换后的振幅
f = n*fs/N; %频率序列
subplot(412)
plot(f,mag1);
axis([0,10,0,3000]);
xlabel('频率/HZ');
ylabel('振幅');
title('Fourier Transform');
%STFT, Hamming windows
h=hanning(N);
[y_STFT,ff,tt,pp] = spectrogram(x,h,100,N,fs);
mag2 = abs(y_STFT);
subplot(413);
n1 = 1:50
clear;
%generate signal
fs = 10; %采样频率
N = 10000; %数据总长度
n = 0:N-1;
t = n/fs; %时间序列
for i=1:N
if t(i)<500
x(i)=sin(10*t(i));
else
x(i)=sin(20*t(i));
end
end
x;
subplot(411)
plot(t,x);
axis([0,1000,-1,1]);
%fourier Transform
y_FFT = fft(x,N); %对信号快速Fourier变换
mag1 = abs(y_FFT); %Fourier变换后的振幅
f = n*fs/N; %频率序列
subplot(412)
plot(f,mag1);
axis([0,10,0,3000]);
xlabel('频率/HZ');
ylabel('振幅');
title('Fourier Transform');
%STFT, Hamming windows
h=hanning(N);
[y_STFT,ff,tt,pp] = spectrogram(x,h,100,N,fs);
mag2 = abs(y_STFT);
subplot(413);
n1 = 1:50