数字信号处理基础据
大理学院实验报告
课程名称生物医学信号处理
实验名称数字信号处理基础据
2015—2016学年度第
三、实验内容
在做频谱分析的时候，为了提高频域分辨率，常常要增加采样点数。如果数据量略有不足，传统的方法是在数据尾部补0来解决，称为高密度频谱（The High Density Spectrum ）。但是补0并不能够提高频域分辨率，我们认为填入适当的现有数据会更好，称为高分辨率频谱（The High Resolution Spectrum）。Matlab 中做频谱分析分函数为fft 。
理想的数字滤波器有：低通、高通、带通和带阻数字滤波器。为了抑制高频噪声，常常还需要平滑数据，这种滤波器也称为平滑滤波器或者中值滤波器。Matlab 中常使用函数medfiltl 函数。
四、实验结果与分析
（包括实验原理、数据的准备、运行过程分析、源程序（代码）、图形图象界面等）
例2-9：
%【例2-9】
N1=10;%定义采集数据的长度为10 N2=10;%定义采集数据的长度为10 N3=100;%定义采集数据的长度为100 N4=50;%定义采集数据的长度为50 N5=80;%定义采集数据的长度为80 % 做10点的DFT for i=1:N1
x_1(i)=cos(0.48*pi*i)+cos(0.52*pi*i);%用该信号采集10点 end
subplot(3,2,1),stem(0:N1-1,x_1);
title('x(n),0
subplot(3,2,2),stem(0:(2/N1):1,abs(f_1(1:N1/2+1))); title('DFT Magnitude')%输出进行dft 后图像 xlabel('frequency in π units') % 做100点的DFT ，其中90个点补零 for j=1:N2
x_2(j)=cos(0.48*pi*j)+cos(0.52*pi*j);%用该信号采集10点 end
for j=N2+1:N3
x_2(j)=0;%其他点补0 end
subplot(3,2,3),stem(0:N3-1,x_2);
title('x(n),0
subplot(3,2,4),stem(0:(2/N3):1,abs(f_2(1:N3/2+1)));