滤波器设计Word文件下载.docx
- 文档编号:18617691
- 上传时间:2022-12-29
- 格式:DOCX
- 页数:12
- 大小:243.22KB
滤波器设计Word文件下载.docx
《滤波器设计Word文件下载.docx》由会员分享,可在线阅读,更多相关《滤波器设计Word文件下载.docx(12页珍藏版)》请在冰豆网上搜索。
常用的窗函数有矩形窗、三角形(Bartlett)窗、汉宁(Hanning)窗、汉明(Hamming)窗、布莱克曼(Blackman)窗、凯泽(Kaiser)窗等,前五种窗函数都是各以一定的主瓣加宽为代价来换取某种程度的旁瓣抑制,而凯泽窗全面的反应了这种主瓣宽度与旁瓣衰减之间的交换关系,可以在这两者之间选择它们的比重,因此本文选取凯泽窗进行滤波器的设计。
利用凯泽窗进行滤波器设计时,截止频率、过渡带宽度、阻带最小衰减及滤波器阶数是设计中的重要指标参数,他们共同影响滤波器的性能。
对于带通数字滤波器,若通带过窄会滤除部分目标信号分量,过宽则会引入过多的噪声,导致信号失真。
对于船体结构应力信号的滤波器设计,本文主要通过考虑结构应力响应和船舶可能遭遇的海况特点来确定滤波器的截止频率;
过渡带宽度、阻带最小衰减及滤波器阶数的关系为
(2-19)
式中:
为过渡带宽度,
为阻带最小衰减,N为阶数。
过渡带宽度越窄滤波器的频域特性越好,越接近理想滤波器,但由式(3-44)可以看出减小过渡带宽度会增加滤波器阶数,阶数越大滤波时计算量越大,效率越低;
增大阻带最小衰减能够更好的抑制阻带噪声信号,但同样会增加滤波器阶数。
本文利用综合评判法根据不同的情况综合考虑过渡带宽度、阻带最小衰减及滤波器阶数这三个因素来确定最佳的滤波器设计参数。
2.2.2 FIR滤波器的参数确定
低海况下船舶在波浪上的载荷响应为线性响应系统,因此可以根据船体结构的应力响应函数确定结构响应明显的频带,并将此频带范围作为滤波器的通带。
下面以某船为例对该方法进行说明。
某船总长241m,船宽32m,设计吃水10.8m,排水量48200t,航速18kn,沿船长重量分布如图2.2所示。
采用三维线性势流理论计算结构波浪载荷响应。
图2.3为船中剖面在不同浪向beta与频率ω下的弯矩响应,可以看出船中剖面弯矩响应主要集中在自然频率为0.3~1.1rad/s的范围内;
图2.4为浪向角为0°
时沿船长均匀分布的21个剖面的弯矩频率响应,可以看出各剖面弯矩响应趋势一致,只是幅值有所不同,其他浪向的各剖面弯矩频率响应规律与0°
浪向时一致,这说明不仅在船中范围内,在船体的任意部位其剖面弯矩响应主要频率范围均在0.3~1.1rad/s。
考虑到实际监测得到的结构应力响应是在遭遇频率下测得的,可按式(2-20)将自然频率转化为遭遇频率
(2-20)
为航向角;
U为航速,m/s;
g为重力加速度,9.81m/s2。
将上述浪向、航速和自然频率带入式(2-20)可得弯矩响应的主要频率范围,可知其最大值约为3.4rad/s,即0.5Hz。
由于随浪工况下当船舶的航速与波速接近时船体与波浪处于相对静止状态,不会产生波浪载荷响应。
因此,滤波器截止频率的下限需通过其他途径确定。
图2.2全船重量分布
(a)满载工况
(b)压载工况
图2.3船中剖面弯矩频响曲线
图3.4各剖面弯矩频响曲线
该船航行区域为北大西洋海域,根据海况资料[8]可知,北大西洋海区波浪周期主要分布在3.5~18.5s,其分布如图3.5所示,由此可知船舶在该海域驻泊时所能遭遇的最大波浪周期不超过20s。
对于随浪航行导致的遭遇周期大于20s的波浪,由于其周期较长,由该船的频响函数可知,船体结构响应很小,因此滤波器截止频率的下限可取为1/20s=0.05Hz。
图2.5北大西洋海况周期分布
当船舶遭遇中高级海况时,由于船体和波浪之间的剧烈相对运动,船体结构容易发生砰击及甲板上浪等强非线性波浪载荷作用,砰击产生的瞬间水动力会引起船体梁的颤振效应,严重的颤振响应会在船中区域的强力甲板和船底板上产生瞬间高应力[9],对于船体结构应力的长期监测而言这种由砰击引起的高应力是不可避免的,同时对这种砰击应力的记录有助于砰击载荷的理论与试验研究。
本文提及的船体结构应力监测系统的硬件设备不仅能够感知砰击载荷引起的瞬间高应力,同时在确定砰击现象发生后能够自动将系统的采样频率跳跃至更高的频率,以详细记录砰击引起的结构应力变化。
上文提及的低海况带通滤波器仅保留了由于波浪循环作用引起的结构应力,却滤除了由于砰击效应引起的结构振动响应。
由于高海况下船体波浪载荷的非线性响应明显,在总载荷响应中由砰击引起的结构振动响应占有较高的比例[10],因此上述滤波器在该工况下已不适用。
考虑到砰击载荷可能引起的高频振动,高海况下的滤波器不仅需要保留由于波浪循环作用引起的结构应力,还需增加结构的振动应力分量。
由于船体高阶振动对总载荷的贡献不大,因此仅保留船体的低阶振动分量即可。
本文在进行滤波器设计时保留船体结构的一阶和二阶振动响应。
船体结构的振动模态可以通过迁移矩阵法、有限元法或试验方法获得,由于前两种方法或多或少总会引入误差,因此本文通过实船试验的方法获得船体结构的一阶和二阶固有频率。
由试验测得的应力时域及其频谱如图2.6所示,根据图2.6(b)可知结构的一阶和二阶固有频率分别为2.257Hz、7.537Hz,则其高频通带范围可取为2.0~2.5Hz和7.25~7.75Hz。
(a)时域
(b)频谱
图2.6实测结果
2.2.3 FIR滤波器设计
本文选取窗函数方法中的凯泽窗设计FIR滤波器。
通常情况下,可以先设计理想滤波器,不论是低通、带通还是高通滤波器,都可以
对于船体结构应力实时监测系统,低频工况下滤波器通带为0.05~0.5Hz,过渡带为0.03Hz,阻带最小衰减为50dB;
高频工况下滤波器通带为2.0~2.5Hz,7.25~7.75Hz,过渡带宽为0.05Hz,阻带最小衰减为40dB。
在低频采样工况下(采样频率取10Hz),采用加凯泽窗的窗函数设计,可以由公式:
(2-21)
和公式:
(2-22)
可以手动计算参数
和
,也可以采用MATLAB自带的函数kaiserord()来估计滤波器的阶数等参数。
[N,wn,bta,ftype]=kaiserord([0.030.050.50.53],[010],[10^-2.50.0510^-2.5],Fs);
该段代码就可以估算一个凯泽窗的
(对应代码中bta参数),再由以下代码:
h1=fir1(n,wn,ftype,kaiser(n+1,bta),'
noscale'
);
就能得到相应的带通滤波器。
滤波器的脉冲响应如图2.7所示,滤波器幅频特性和相频特性如图2.8所示。
图2.7低频采样工况下FIR滤波器脉冲响应
图2.8低频采样工况下FIR滤波器幅频相频特性
由图2.8可知,滤波器通带范围内,该FIR滤波器的相位具有线性的特性,与上述关于FIR滤波器的线性相位特性相吻合。
MATLAB仿真使用多个正弦信号来摸你合成信号,假设有用的信号频率为0.08Hz,噪声信号分别为0.03Hz和0.55Hz。
上述带通滤波器的通带是0.05Hz~0.5Hz,理论上可以将噪声信号滤除,得到0.08Hz的有用信号。
MATLAB核心代码如下:
Fs=10;
[n,wn,bta,ftype]=kaiserord([0.030.050.50.53],[010],[10^(-2.5)0.05
10^(-2.5)],Fs);
[hh1,f]=freqz(h1,1,1024,Fs);
freqz(h1,1,1024,Fs);
滤波前不加噪声和加噪声的信号时域以及频域特性,如图2.9所示,滤波之后的信号时域和频域特性,如图2.10所示。
由图可以看出,滤波后的信号频域分量为0.08Hz,达到了滤波的效果。
图2.9低频采样工况下原始信号
图2.10低频采样工况下滤波后信号
在高频采样工况下(采样频率去1000Hz),滤波器脉冲响应如图2.11所示,幅频相频特性如图2.12所示,可以看出在三个带通范围内,相位呈现线性。
图2.11高频采样工况下FIR滤波器脉冲响应
图2.12高频采样工况下FIR滤波器幅频相频特性
如图2.13和2.14所示,分别为原是信号和滤波后信号,可以看出滤波后信号能够较好的保留有用的信号,将噪声信号滤除。
图2.13高频采样工况下原始信号
图2.14高频采样工况下滤波后信号
Welcome!
!
欢迎您的下载,
资料仅供参考!
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 滤波器 设计