logo

您所在位置网站首页 > 海量文档  > 计算机 > matlab

Matlab信号处理倍频程.docx 7页

本文档一共被下载: ,您可全文免费在线阅读后下载本文档。

  • 支付并下载
  • 收藏该文档
  • 百度一下本文档
  • 修改文档简介
全屏预览

下载提示

1.本站不保证该用户上传的文档完整性,不预览、不比对内容而直接下载产生的反悔问题本站不予受理。
2.该文档所得收入(下载+内容+预览三)归上传者、原创者。
3.登录后可充值,立即自动返金币,充值渠道很便利
一、对信号进行倍频程分析1、流程图定义输入信号xn显示信号xn滤波器设计显示信号滤波后结果2、程序代码%对信号进行倍频程分析,时间2012-6-14,编程人员韩宝安clc; %清空clear all; %清除所有变量close all; %关闭所有窗口pref=2e-5;%构造输入信号bookName = '011b_A_12高度_1-signal时域信号.xls'; % book名sheetName = 'sheet1'; %sheet名range = 'C85:C4180'; %Excel表中的A85:F2132之间的数据tmp = xlsread(bookName, sheetName, range); %将读取的数据存于tmp中Fs=4096; %采样频率xn = tmp'; %输入信号xnt=1/Fs:1/Fs:1; %定义时间tplot(t,xn); %以t1为横坐标,x1为纵坐标绘制图像xlabel('t/ s','fontsize',15); %横坐标轴标记为t/sylabel('xn / (Pa)','fontsize',15); %纵坐标轴标记为xn/(pa)%滤波器设计n=input('请确定倍频程数n: '); %确定倍频程数n N=5; %滤波器阶数 k=0; %循环次数初始化 w2=22.5; %初始化w2while(w2<=Fs/2) %k循环加1,直到w2>Fs/2 w1=w2; w2=w1*2^n; k=k+1;end w2=22.5;for m=1:1:k-1 %m从1每次加1,直到m=k-1 w1=w2; %确定带通下截止频率w1 w2=w1*2^n; %确定带通上截止频率w2 centerf(m)=(w1+w2)/2; %计算中心频率centerf wn=[w1 w2]/(Fs/2); %得到归一化角频率 [b,a]=butter(N,wn); %设计带通巴特沃斯滤波器 y=filter(b,a,xn); %得到滤波结果y sm=std(y); %计算y的标准偏差 e(m)=sm/length(y); %能量值 end figure; stem(centerf,20*log10(e/pref)); xlabel('centerf/ Hz','fontsize',15); %横坐标轴标记为centerf/ Hzylabel('e / (Pa)','fontsize',15); %纵坐标轴标记为e/(pa)3、程序运行结果(1)程序运行时,输入信号xn的波形如图1。图1输入信号波形(2)巴特沃斯滤波器滤波结果①当倍频程数n=1时,滤波结果如图2。图2 n=1时的滤波结果②当倍频程数n=1/3时,滤波结果如图3。图3 n=1/3时的滤波结果③当倍频程数n=1/6时,滤波结果如图4。图4 n=1/6时的滤波结果(3)切比雪夫滤波器结果说明:由于切比雪夫滤波器的实现代码基本和巴特沃斯滤波器的一致,这里就不在这里赘述,只需将[b,a]=butter(N,wn)改为[b,a]=cheby1(N,r,wn),其中r为通带内的波纹大小,一般取0.5。①当倍频程数n=1时,滤波结果如图5。图5 n=1时的滤波结果②当倍频程数n=1/3时,滤波结果如图6。图6 n=1/3时的滤波结果③当倍频程数n=1/6时,滤波结果如图7。图7 n=1/6时的滤波结果

发表评论

请自觉遵守互联网相关的政策法规,严禁发布色情、暴力、反动的言论。
用户名: 验证码: 点击我更换图片

“原创力文档”前称为“文档投稿赚钱网”,本站为“文档C2C交易模式”,即用户上传的文档直接卖给(下载)用户,本站只是中间服务平台,本站所有文档下载所得的收益归上传人(含作者)所有【成交的100%(原创)】。原创力文档是网络服务平台方,若您的权利被侵害,侵权客服QQ:3005833200 电话:19940600175 欢迎举报,上传者QQ群:784321556