三种不同采样频率下的抽样信号及其频谱图

clc

clear all

fs = input('采样频率为：fs=');

N = 1024;

t = (0:N-1)/fs;

y = 0.5*sin(20*pi*t) + 0.6*sin(180*pi*t) + 0.8*cos(300*pi*t); subplot(2,1,1)

plot(t,y);grid on

title(['抽样信号时域波形图，抽样频率fs=', num2str(fs)]);xlabel('t / ms');ylabel('y')

f = (0:N-1)*fs/N;

y1 = fft(y,N); %对信号进行快速Fourier 变换

y2 = abs(y1)/N*2

subplot(2,1,2);

plot(f(1:N/2),y2(1:N/2)); grid on

title('信号的频谱图');xlabel('f / Hz');ylabel('Rex（f）');