专题一 在数字信号处理中的应用 信号处理涉及大量的运算,可以说离开计算机及相应的软件,就不可能解决任何稍微复杂的实际应用问题
借助程序可以很好理解抽象的傅里叶变换的作用及其意义
示例 1 通过傅里叶变换找出被随机信号淹没的信号频率 设采样频率 f=1000Hz,信号频率分布在50~100Hz,然后叠加随机噪声 源程序: fs = 1000; % 采样速率 Ts = 1/fs; t = 0:Ts:0
6; signal = sin(2*pi*50*t)+sin(2*pi*100*t); % 信号频率分布在50~100Hz noise = 2*rand(size(t)); sn = signal+noise; % 染噪信号 NFFT = 512; Ys = fft(signal,NFFT); Yn = fft(sn,NFFT); Ps = abs(Ys); Pn = abs(Yn); f = 1000*(0:256)/NFFT; subplot(2,2,1); plot(t,signal); xlabel('时间(秒)'); title('不含噪声信号'); subplot(2,2,2); plot(t,sn); xlabel('时间(秒)'); title('含噪声信号'); subplot(2,2,3); plot(f,Ps(1:257)); xlabel('频率(Hz)'); title('不含噪声信号的幅值谱'); axis([0,300,0,300]); subplot(2,2,4); plot(f,Pn(1:257)); xlabel('频率(Hz)'); title('含噪声信号的幅值谱'); axis(