版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、现 代 信 号 处理学院: 仪器与电子学院 专业: 仪器仪表工程 班级:姓名:学号:流水号:时间: 2023年 6月27日 1、设采样周期T=250s采样频率fs =4kHz,用脉冲响应不变法和双线性变换法设计一个三阶巴特沃兹滤波器,其3dB边界频率为fc =1kHz。设计步骤如下:确定所需类型数字滤波器的技术指标。将所需类型数字滤波器的边界频率转换成相应的模拟滤波器的边界频率,转换公式为=2/T tan(0.5)将相应类型的模拟滤波器技术指标转换成模拟低通滤波器的技术指标。设计模拟低通滤波器。通过频率变换将模拟低通转换成相应类型的过渡模拟滤波器。采用脉冲响应不变法和双线性变换法将相应类型的过
2、渡模拟滤波器转换成所需类型的数字滤波器。程序设计如下:B,A=butter(3,2*pi*1000,s);%第一个参数为阶数;第二个参数为角频率;第三个参数为系数B、A按s的正降幂排列。num1,den1=impinvar(B,A,4000); %用冲激响应不变法将模拟滤波器变换成数字滤波器h1,w=freqz(num1,den1);%获取滤波器的频率特性B,A=butter(3,2/0.00025,s);num2,den2=bilinear(B,A,4000); %实现双线性变换,即由模拟滤波器 得到数字滤波器h2,w=freqz(num2,den2);f=w/pi*2000;plot(f,
3、abs(h1),-.,f,abs(h2),-);%x为图形上之x坐标向量,y为其对应的y坐标向量grid;xlabel(频率/Hz )ylabel(幅值/dB)程序中第一个butter 的边界频率21000,为脉冲响应不变法原型低通滤波器的边界频率;第二个butter 的边界频率2/T=2/0.00025,为双线性变换法原型低通滤波器的边界频率.图1 给出了这两种设计方法所得到的频响,虚线为脉冲响应不变法的结果;实线为双线性变换法的结果。脉冲响应不变法由于混叠效应,使得过渡带和阻带的衰减特性变差,并且不存在传输零点。同时,也看到双线性变换法,在z=-1 即=或f=2000Hz 处有一个三阶传输
4、零点,这个三阶零点正是模拟滤波器在=处的三阶传输零点通过映射形成的。MATLAB计算结果如图1所示:图12、设计一数字高通滤波器,它的通带为400500Hz,通带内容许有0.5dB的波动,阻带内衰减在小于317Hz的频带内至少为19dB,采样频率为1000Hz。wc=2*1000*tan(2*pi*400/(2*1000);wt=2*1000*tan(2*pi*317/(2*1000);N,wn=cheb1ord(wc,wt,0.5,19,s);B,A=cheby1(N,0.5,wn,high,s);num,den=bilinear(B,A,1000);h,w=freqz(num,den);f
5、=w/pi*500;plot(f,20*log10(abs(h);axis(0,500,-80,10);grid;xlabel()ylabel(幅度/dB)图2给出了MATLAB计算的结果,可以看到模拟滤波器在=处的三阶零点通过高通变换后出现在=0z=1处,这正是高通滤波器所希望得到的。图23、设计一巴特沃兹带通滤波器,其3dB边界频率分别为f2=110kHz和f1=90kHz,在阻带f3 = 120kHz处的最小衰减大于10dB,采样频率fs=400kHz。w1=2*400*tan(2*pi*90/(2*400);w2=2*400*tan(2*pi*110/(2*400);wr=2*400*
6、tan(2*pi*120/(2*400);N,wn=buttord(w1 w2,1 wr,3,10,s);B,A=butter(N,wn,s);num,den=bilinear(B,A,400);h,w=freqz(num,den);f=w/pi*200;plot(f,20*log10(abs(h);axis(40,160,-30,10);grid;xlabel(频率/kHz)ylabel(幅度/dB)MATLAB计算结果如图3所示:图34、一数字滤波器采样频率fs=1kHz,要求滤除100Hz的干扰,其3dB的边界频率为95Hz和105Hz。低通滤波器传递函数如下:H1a(s)=1/(1+s
7、)w1=95/500;w2=105/500;B,A=butter(1,w1, w2,stop);h,w=freqz(B,A);f=w/pi*500;plot(f,20*log10(abs(h);axis(50,150,-30,10);grid;xlabel(频率/Hz)ylabel(幅度/dB)图4为MATLAB的计算结果图45、用计算机麦克录自己的语音信号,语音信号采样频率为22050,用matlab完成以下分析:1播放语音信号;对信号做1024点FFT变换;做原始语音信号的时域图形;绘制原始语音信号的频率响应图;做原始语音信号的FFT频谱图。fs=22050; %语音信号采样频率为2205
8、0 x1=wavread(D:program-toolMATLABR2023bwork clock001.wav);%读取语音信号的数据,赋给变量x1sound(x1,32768); %播放语音信号y1=fft(x1,1024); %对信号做1024点FFT变换f=fs*(0:511)/1024;figure(1);plot(x1); %做原始语音信号的时域图形title(原始语音信号);xlabel(time n);ylabel(fuzhi n);figure(2);freqz(x1) ;%绘制原始语音信号的频率响应图title(频率响应图);figure(3);subplot(2,1,1)
9、;plot(abs(y1(1:512) ;%做原始语音信号的FFT频谱图title(原始语音信号FFT频谱);subplot(2,1,2);plot(f,abs(y1(1:512);title(原始语音信号频谱);xlabel(Hz);ylabel(fuzhi);MATLAB计算结果如图5、6、7所示:图5图6图72在语音信号中参加随机噪声;播放加噪声后的语音信号,绘制加噪后的语音信号;首先画出语音信号的时域波形,再对语音信号进行快速傅里叶变换,得到信号的频谱特性。给原始的语音信号加上一个高频余弦噪声,频率为5kHz。画出加噪后的语音信号时域和频谱图,与原始信号比照,可以很明显的看出区别程序如
10、下:fs=22050;x1=wavread(qq.wav);f=fs*(0:511)/1024;t=0:1/22050:(size(x1)-1)/22050; %将所加噪声信号的点数调整到与原始信号相同Au=0.3;d=Au*cos(2*pi*5000*t); %噪声为5kHz的余弦信号x2=x1+d;sound(x2,22050); %播放加噪声后的语音信号y2=fft(x2,1024);figure(1)plot(t,x2);title(加噪后的信号);xlabel(time n);ylabel(fuzhi n);figure(2)subplot(2,1,1);plot(f,abs(y1(
11、1:512);title(原始语音信号频谱);xlabel(Hz);ylabel(fuzhi);subplot(2,1,2);plot(f,abs(y2(1:512);title(加噪后的信号频谱);xlabel(Hz);ylabel(fuzhi);得到原始语音信号时域图形,原始语音信号幅度谱,参加噪声之后的语音信号时域图形,参加噪声之后的语音信号幅度谱如图8、9所示:图8图93设计适宜的数字滤波器,将上述加噪声滤掉;播放滤波后的信号;绘制滤波前和滤波后的语音信号及频谱图。双线性变换法设计Butterworth滤波器fs=22050;x1=wavread(qq.wav);t=0:1/22050
12、:(size(x1)-1)/ 22050;Au=0.03;d=Au*cos(2*pi*10000*t);x2=x1+d;wp=0.25*pi;ws=0.3*pi;Rp=1;Rs=15;Fs=32768;Ts=1/Fs;wp1=2/Ts*tan(wp/2); %将模拟指标转换成数字指标ws1=2/Ts*tan(ws/2); N,Wn=buttord(wp1,ws1,Rp,Rs,s);%选择滤波器的最小阶数Z,P,K=buttap(N); %创立butterworth模拟滤波器Bap,Aap=zp2tf(Z,P,K);b,a=lp2lp(Bap,Aap,Wn); bz,az=bilinear(b,
13、a,Fs); %用双线性变换法实现模拟滤波器到数字滤波器的转换H,W=freqz(bz,az); %绘制频率响应曲线figure(1)plot(W*Fs/(2*pi),abs(H)gridxlabel(频率Hz)ylabel(频率响应幅度)title(Butterworth)f1=filter(bz,az,x2);figure(2)subplot(2,1,1)plot(t,x2) %画出滤波前的时域图title(滤波前的时域波形);subplot(2,1,2)plot(t,f1); %画出滤波后的时域图title(滤波后的时域波形);sound(f1, 22050); %播放滤波后的信号F0=
14、fft(f1,1024);f=fs*(0:511)/1024;figure(3)y2=fft(x2,1024);subplot(2,1,1);plot(f,abs(y2(1:512); %画出滤波前的频谱图title(滤波前的频谱)xlabel(Hz);ylabel(fuzhi);subplot(2,1,2)F1=plot(f,abs(F0(1:512); %画出滤波后的频谱图title(滤波后的频谱)xlabel(Hz);ylabel(fuzhi);计算结果如图10、11、12图10图11图126、信号,信号带宽,采样频率,采样点数。用STFT和WD分析其特性。设计过程:用短时傅里叶变换分析
15、其在时域、频域特性,经过STFT和WD运算后的时频分布结果特点,从而将它们识别。设计的源程序如下:k=4;t=5;fc=k*T=20;fs=4*fc;Ts=1/fs;N=T/Ts=400;X=zeros(1,N);T=0:N-1;X=exp(j*k*pi*(t*Ts).2);Subplot(2,2,2);Plot(t-N/2)*fs/N,abs(X);Nw=20;L=Nw/2;Tn=(N-Nw)/L+1;Nfft=32;TF=zeros(Tn,nfft);for i=1:Tnxw=x(i-1)*10+1:i*10+10);temp=fft(xw,nfft);temp=fftshift(temp);TF(i,:)=temp;endsubplot(2,2,3);fnew=(1:nfft)-nfft/2)*fs/nfft;tnew=(1:Tn)*L*Ts;F,T=meshgrid(fnew,tnew);mesh(F,T,abs(TF);subplot(2,2,4);计算结果如图13所示图137、数字信号在125,250和500,625分别有两个频率不同的正弦信号,用STFT和WD分析其特性。k=4;T=5;fc=20;fs=4*fc;Ts=1/fs;N=T/Ts;x=zeros(1
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 临沂市兰山区部分医疗卫生事业单位招聘卫生类岗位人员笔试真题2023
- 2022年吕梁交城县招聘大学生到村工作考试试卷及答案解析
- 2022年湖南湘南学院招聘事业编制人员考试试卷及答案解析
- 2023年淮安生态文化旅游区城市管理工作合同制人员招聘考试真题及答案
- 学校交通安全主题教育倡议书范文
- 预防学生欺凌致家长告知书范文
- 儿童中后期的认知发展-皮亚杰具体运算阶段
- 个案工作的基本原则之非评判
- 第四章第二节第三单元
- 妇女婚姻家庭权益
- 礼貌用语精品课件
- 高教版中职语文上册口语交际《电话交谈》word教案
- 成都的房屋租赁合同3篇(成都市房屋租赁合同 房管局)
- 自动涂胶机机械系统设计和实现 机械制造自动化专业
- 路肩挡土墙施工安全技术交底
- 汽修厂汛期和“五一”期间安全生产工作方案
- DB32∕T 4304-2022 高速公路经营管理安全生产隐患排查治理规程
- 预制管桩、方桩、钻孔灌注桩单桩承载力计算(自动版)
- 输血规范以及流程
- 高二年级日语科答题卡
- 2022年数字化工厂规划报告(三年实施计划)
评论
0/150
提交评论