北理工_数据分析_实验6_信号频谱分析与信号滤波

北京理工大学
现代数据分析
实验6实验报告
主讲：李明
学生：李经
2012/11/14
实用文档
实验6 数据拟合 (1)
6.1 实验目的 (1)
6.2 实验内容 (1)
6.2.1 信号的频谱分析 (1)
6.2.2 滤波 (1)
6.3 实验代码及结果 (1)
6.3.1 信号的频谱分析&& 滤波 (1)
(图表页)
图1. 噪声信号 (3)
图2. 功率谱 (3)
图3. 显示信号 (4)
图4. 滤波结果 (4)
实用文档
实验6数据拟合
6.1实验目的
熟悉信号的频谱分析与信号滤波。

6.2实验内容
6.2.1信号的频谱分析
建立一个含10Hz和100Hz幅值为1的正弦信号（sin），然后叠加一个幅值为1的随机信号，利用Matlab分析其频谱。

6.2.2滤波
接上题，滤除噪声信号和100Hz正弦信号。

(启动Demo，在地址栏输入Signal Processing Demo: filtdem)
6.3实验代码及结果
6.3.1信号的频谱分析&& 滤波
代码：
fs=1000;
t = 0:1/fs:1;
x = sin(2*pi*10*t) + sin(2*pi*100*t);
y = x + randn(size(t));
实用文档
plot(y(1:50));
Y = fft(y,256);
Pyy = Y.*conj(Y)/256;
f = fs/256*(0:127);
figure(2);
plot(f,Pyy(1:128));
xlabel('Frequency (Hz)')
figure(3);
plot(t,y);
[b,a] = ellip(4,0.1,40,[5 15]*2/fs); sf = filter(b,a,y);
figure(4);
plot(t,sf);
axis([0 1 -1 1]);
结果：
实用文档
图1.噪声信号
图2.功率谱实用文档
图3.显示信号
图4.滤波结果实用文档。