Matlab信号处理倍频程.docx

一、对信号进行倍频程分析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时的滤波结果