可视化计算离线作业(6)

太阳 | 地球

??

??

5.（10分）设计一个低通滤波器，从混合信号：

x(t)=sin(2*pi*10*t) + cos(2*pi*100*t) + 0.2*randn(size(t)) 中获取10Hz的信号（10分）。

clear

hh=figure('numbertitle','off','name','低通滤波器<2015年秋 >');

randn('state',1)

ws=1000;%采样频率

t=0:1/ws:0.4

x=sin(2*pi*10*t)+cos(2*pi*100*t)+0.2*randn(size(t));%生成带噪声的多频信号wn=ws/2;%Nyquest 频率

[B,A]=butter(10,30/wn);

y=filter(B,A,x);

plot(t,x,'b-',t,y,'r.','MarkerSize',10);

legend('Input','Output',0);

title(' ')