示波器 信号处理.docVIP

机电工程学院 信号处理课程设计 学 号：专 业学生姓名：教师：2013年月  y1=1.5*sin(2*pi*fc*t); %正弦波信号 y2=randn(1,N); %均匀白噪声信号 n=0:N-1; y=y1+y2; %正弦波信号叠加白噪声信号 subplot(311); %创建子图 plot(t,y,r); title(均匀白噪声与正弦波信号叠加图); ylabel(y); xlabel(t); %定义各坐标轴:横坐标为t,纵坐标为y grid on; %划分网格 rt=xcorr(y,125,unbiased); %自相关函数 X=fft(t); %傅立叶变换 mag=abs(X); f=n*fs/N; subplot(312); %创建子图 plot(rt); title(自相关函数图形); ylabel(rt); xlabel(t); %定义各坐标轴:横坐标为t,纵坐标为rt grid on; %划分网格 subplot(313); %创建子图 plot(f,mag); title(叠加后信号的FFt变换图形); ylabel(F(jw));xlabel(w); %定义各坐标轴:横坐标为w,纵坐标为F(jw) grid on; %划分网格 原始波形、自相关函数图形和FFT变换后的图形 3. 信号中主要混有频率低于0.5w和高于4.5w的信号。 4.（1）在信号处理中，自相关检测淹没在随机噪声干扰信号在信号处理中， 题2: IIR 滤波器设计 得分 Butterworth低通IIR数字滤波器设计 （1）源程序 fp=50;rp=3;fs=80;Fs=1000;rs=20 %初始化相应参数，给定需要的衰减，单位为dB wp=2*pi*fp/Fs; ws=2*pi*fs/Fs;Fs=Fs/Fs; %令Fs=1 wap=tan(wp/2); was=tan(ws/2); %实现由数字滤波器的频率到模拟滤波器的频率的转换[n,wn]=buttord(wap,was,rp,rs,s); [n,wn]=buttord(wap,was,rp,rs,s) %求模拟滤波器的阶次 [z,p,k]=buttap(n); %设计模拟低通Butterworth滤波器，得到极零点 [bp,ap]=zp2tf(z,p,k); %由极零点得到滤波器分子分母多项式的系数[bs,as]=lp2lp(bp,ap,wap); [bs,as]=lp2lp(bp,ap,wap); %归一化低通到实际低通的转换 [bz,az]=bilinear(bs,as,Fs/2); %实现s到z的转换 [h,w]=freqz(bz,az,256,Fs*1000); %求出滤波器的频率响应 subplot(311) plot(w,abs(h)); title(Butterworth低通滤波器) xlabel(频率(Hz));ylabel(幅度) grid on; fs=5;fc=29/100;t=(0:0.2:49.99);N=250; %初始化 y1=1.5*sin(2*pi*fc*t); %正弦波信号 y2=randn(1,N); %均匀白噪声信号 n=0:N-1; y=y1+y2;