matlab 三分之一倍频程,三分之一倍频程谱
三分之一倍频程谱是一种频率分析方法,它具有谱线少频带宽的特点。
倍频程实际上是频域分析中频率的一种相对尺度。倍频程谱是由一系列频率点以及对应这些频率点附近的频带内信号的平均幅值(有效值)所构成。这些频率点称为中心频率fc,中心频率附近的频带出于下限频率fl和上限频率fu之间。
三分之一倍频程谱是按逐级式频率进行分析的,它是由多个带通滤波器并联组成,为的是使这些带通滤波器的带宽覆盖整个分析频带。根据国籍电工委员会(IEC)的推荐,三分之一倍频程的中心频率为:
fc = 1000*10^(3n/30)Hz
(n = ...-3,-2,-1,0,1,2,3...)
但在实际应用中,通常采用中心频率是其近似值。按照我国现行标准规定,中心频率为1Hz,1.25Hz,1.6Hz,2Hz,2.5Hz,3.15Hz,4Hz,5Hz,6.3Hz,8Hz,10Hz,...。可以看出,每隔三个中心频率,频率值增加一倍。三分之一倍频程的上下限频率以及中心频率之间的关系为:
fu/fl = 2^(1/3);fc/fl =
2^(1/6);fu/fc = 2^(1/6)
三分之一倍频程带宽为
△f = fu - fl
对于三分之一倍频程谱可以通过两种处理方法得到。一种方法是在整个分析频率从范围,按照不同的中心频率从定义对采样信号进行带通滤波,然后,计算出滤波后数据的均方值或均方根值(有效值),这样,变得到对应每个中心频率的功率谱值或幅值谱值。由于但分之一倍频程谱的滤波带与中心频率的比值是不变的,这种处理方法称为恒定百分比带宽滤波法。另一种方法首先对采样信号进行快速傅立叶变换,计算出功率谱或幅值谱,然后用功率谱或幅值谱的数据,计算每一个中心频率带宽内数矩的平均值,这样便处理得到三分之一倍频程谱值。这种方法显然要比第一种方法处理效率要高的多。
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 三分之一倍频程处理
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear % Remove items from workspace, freeing up system memory
clc % Clear Command Window
clf % Clear current figure window
close all hidden % removal
draws only those lines that are not obscured
% by other objects in the field of view.
format long
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% 加载数据
fun = @(t) sin(50*2*pi*t) +sin(20*2*pi*t) + randn(size(t));
sf =
300; % 采样频率
n =
1000; % 数据长度
t =
0:1/sf:(n-1)/sf; % 建立离散时间列向量
x =
fun(t); % 产生时间序列数据
% 定义三分之一倍频程的中心频率
f = [1.00 1.25 1.60 2.00 2.50 3.15 4.00 5.00 6.300 8.00];
fc = [f,10*f,100*f,1000*f,10000*f];
oc6 =
2^(1/6); % 中心频率与下限频率的比值
nc =
length(fc); %
取中心频率总的长度 nfft = 2^nextpow2(n); % 大于并最接近n的2的幂次方长度
a =
fft(x,nfft); % FFT变换
yc = zeros(1,nc);
for j =
1:nc fl = fc(j)/oc6; %
下限频率 fu =
fc(j)*oc6; % 上限频率
nl = round(fl*nfft/sf+1); % 下限频率对应的序号
nu = round(fu*nfft/sf+1); % 上限频率对应的序号
if fu >
sf/2 % 如果上限频率大于折叠频率则循环中断
m = j-1; break
end
%
以每个中心频率段为通带进行带通频域滤波 b = zeros(1,nfft);
b(nl:nu) = a(nl:nu);
b(nfft-nu+1:nfft-nl+1) = a(nfft-nu+1:nfft-nl+1);
c = ifft(b,nfft);
%计算对应每个中心频率段的有效值
yc(j) = sqrt(var(real(b(1:n))));
end
%绘制输入时程曲线图形
subplot(2,1,1);plot(t,x);
xlabel('时间 (s)');ylabel('加速度 (g)');grid
on; %绘制三分之一倍频程有效值图形
subplot(2,1,2); plot(fc(1:m),yc(1:m));
xlabel('频率 (Hz)');ylabel('有效值'); grid on;
matlab 三分之一倍频程,三分之一倍频程谱相关推荐
- 2倍频程、1倍频程、1/3倍频程介绍
什么是倍频程?有何用? 倍频程,又称倍波程,指在滤波特性曲线上,频率或波长之比为2或1/2的两个频率或波长之间的间隔.常用于音响设备的音频信号分析. 声音信号的频率范围为20Hz~20kHz.分析音频 ...
- Matlab应变片仿真,一种基于Matlab/Adams联合仿真的真实路谱再现系统和方法与流程...
本发明属于汽车系统动力学仿真技术领域,特别是一种基于Matlab/Adams联合仿真的真实路谱再现系统和方法. 背景技术: 汽车系统动力学仿真技术是汽车设计制造中一项不可或缺的技术,尤其是在汽车操纵稳 ...
- MATLAB处理信号得到频谱、相谱、功率谱
MATLAB处理信号得到频谱.相谱.功率谱 第一:频谱 一.调用方法 X=FFT(x): X=FFT(x,N): x=IFFT(X); x=IFFT(X,N) 用MATLAB进行谱分析时注意: (1) ...
- matlab实现三分之一倍频程,三分之一倍频程程序
方法一:%A计权声压级频谱分析 clc; clear; close all; y=wavread('abc.wav'); fs=51200;%采样频率 p0=2e-5;%参考声压 f=[1.00 1. ...
- 倍频程频谱_倍频程科学编程语言速成课程
倍频程频谱 GNU Octave is a free, scientific programming language. It offers a rich mathematical apparatus ...
- matlab怎么选清浊音做短时谱,语音信号处理实验三
实验三语音信号的同态处理 一.实验目的 1.熟悉语音信号求其倒谱与复倒谱的方法: 2.计算一段语音信号的频谱与倒谱: 3.比较由加窗语音信号求得的对数幅度谱与倒谱域滤波后得到的对数幅度谱的各自特点. ...
- matlab怎么选清浊音做短时谱,matlab语音信号处理如何判别清浊音?
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 这是我一个学长给的程序,他说里面要算得基本上已经实现了改一下就可以用了...但是i本人不是主攻matlab的说白了就是不怎么会,,想问下该怎么改来实现判断 ...
- 用MATLAB画出离散幅度和相位谱
Eg:画出x(t)=cos(πt/8),|t|<=4 信号的x(t)的离散幅度和相位谱 main.m echo on fnct='Fcos'; a=-4; b=4; n=24; tol=0.1; ...
- 【光学】基于matlab模拟光纤布拉格光栅FBG反射谱和透射谱仿真
✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,matlab项目合作可私信.
最新文章
- 华为服务器故障灯不开机_总有故障灯亮却不知道是怎么回事?详解这些你不认识的故障灯...
- 郁金香2013网授汇编逆向与外挂 [ 91课,全 ]
- [css] 手写一个满屏品字布局的方案
- js iframe 出现跨越问题
- 21-win10下ElasticSearch.6.1.0安装SQL插件
- View4.5测试参考文档1--测试环境规划
- html音乐静音代码,HTML Audio muted用法及代码示例
- php安装包进行安装吗,php的一键安装包有哪些php环境搭建
- Web — 调色盘打开+div
- Java并发编程实战_盖兹
- getJSON()的使用
- html+css制作小米商城官网初尝试
- html倒计时免费代码,JS倒计时代码汇总
- 360 os3.0 android7.1,360手机N5s发布:升级基于安卓7.1的360 OS 2.0系统
- 如果令 A B C D E F G H I J K L M N O P Q R S T U V W X Y Z 分别等于
- 嵌入式arm板linux路由表设置
- XPO学习(6)----实现一个简单的留言薄(ASP.NET)
- Linux的进程优先级-邹立巍
- 深圳商报 | 跨境电商大会聚焦跨境收款难题
- 微医的第七年:依托互联网医院联合体,如何推动人工智能的战略升级?
热门文章
- LintCode 488.快乐数
- 学习编程语言的第一步,认识什么是计算机!!!
- 成都绿源电动车店寻“新欢”,只这一试,我就认定了
- U8采购入库单API接口示例(参照采购到货)
- QA是干什么的?(职责所在????)
- 魔鬼面试官:用户在电商网站中购买成功了,那么它在微服务中经历了什么?...
- Quick BI企业报表制作——个人笔记
- 微信小程序开发一个多少钱
- 【CSS3】 CSS3实现“图片阴影”效果
- android 画图