搜索
一、对信号进行倍频程分析1、流程图2、程序代码%对信号进行倍频程分析，时间2012-6-14，编程人员韩宝安clc;%清空clearall;%清除所有变量closeall;%关闭所有窗口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);%将读取的数据存定义输入信号xn显示信号xn滤波器设计显示信号滤波后结果初始化结束于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:');%确定倍频程数nN=5;%滤波器阶数k=0;%循环次数初始化w2=22.5;%初始化w2while(w2=Fs/2)%k循环加1，直到w2Fs/2w1=w2;w2=w1*2^n;k=k+1;endw2=22.5;form=1:1:k-1%m从1每次加1，直到m=k-1w1=w2;%确定带通下截止频率w1w2=w1*2^n;%确定带通上截止频率w2centerf(m)=(w1+w2)/2;%计算中心频率centerfwn=[w1w2]/(Fs/2);%得到归一化角频率[b,a]=butter(N,wn);%设计带通巴特沃斯滤波器y=filter(b,a,xn);%得到滤波结果ysm=std(y);%计算y的标准偏差e(m)=sm/length(y);%能量值endfigure;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。图2n=1时的滤波结果②当倍频程数n=1/3时，滤波结果如图3。图3n=1/3时的滤波结果③当倍频程数n=1/6时，滤波结果如图4。图4n=1/6时的滤波结果（3）切比雪夫滤波器结果说明：由于切比雪夫滤波器的实现代码基本和巴特沃斯滤波器的一致，这里就不在这里赘述，只需将[b,a]=butter(N,wn)改为[b,a]=cheby1（N,r,wn）,其中r为通带内的波纹大小，一般取0.5。①当倍频程数n=1时，滤波结果如图5。图5n=1时的滤波结果②当倍频程数n=1/3时，滤波结果如图6。图6n=1/3时的滤波结果③当倍频程数n=1/6时，滤波结果如图7。图7n=1/6时的滤波结果