应用matlab对信号进行频谱分析及滤波.pdfVIP

应用MATLAB对信号进行频谱分析及滤波 fs input(please input the fs:);%设定采样频率 N input(please input the N :);%设定数据长度 t 0:0.001:1; f 100;%设定正弦信号频率 %生成正弦信号 x sin (2*pi*f*t); figure(1); subplot(211); plot(t,x);%作正弦信号的时域波形 axis([0,0.1,-1,1]); title(正弦信号时域波形); z square(50*t); subplot(212) plot(t,z) axis([0,1,-2,2]); title(方波信号时域波形);grid; %进行FFT变换并做频谱图 y fft(x,N);%进行fft变换 mag abs(y);%求幅值 f (0:N-1)*fs/N ;%横坐标频率的表达式为f (0:M-1)*Fs/M ; figure(2); subplot(211); plot(f,mag);%做频谱图 axis([0,1000,0,200]); title(正弦信号幅频谱图); y1 fft(z,N);%进行fft变换 mag abs(y1);%求幅值 f (0:N-1)*fs/N ;%横坐标频率的表达式为f (