- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
matlab频域积分,matlab数值积分实现 (时域频域积分)学步
园
最近做有关加速度的数据处理,需要把加速度积分成 移,⽹上找了找相关资料,发现做这个并不多,把最近做的总结⼀下吧 !
积分操作主要有两种⽅法 :时域积分和频域积分,积分中常见的问题就是会产⽣⼆次趋势。关于积分的⽅法,在国外⼀个论坛上有⼈提出了
如下说法,供参考。
Double integration of raw acceleration data is a pretty poor estimate for displacement. The reason is that at each
integration, you are compounding the noise in the data.
If you are dead set on working in the time-domain, the best results come from the following steps.
1. Remove the mean from your sample (now have zero-mean sample)
2. Integrate once to get velocity using some rule (trapezoidal, etc.)
3. Remove the mean from the velocity
4. Integrate again to get displacement.
5. Remove the mean. Note, if you plot this, you will see drift over time.
6. To eliminate (some to most) of the drift (trend), use a least squares fit (high degree depending on data) to determine
polynomial coefficients.
7. Remove the least squares polynomial function from your data.
A much better way to get displacement from acceleration data is to work in the frequency domain. To do this, follow these
steps...
1. Remove the mean from the accel. data
2. Take the Fourier transform (FFT) of the accel. data.
3. Convert the transformed accel. data to displacement data by dividing each element by -omega^2, where omega is the
frequency band.
4. Now take the inverse FFT to get back to the time-domain and scale your result.
This will give you a much better estimate of displacement.
说到底就是频域积分要⽐时域积分效果更好,实际测试也发现如此。原因可能是时域积分时积分⼀次就要去趋势,去趋势就会降低信号的能
量,所以最后得到的结果常常⽐真实幅值要⼩。下⾯做⼀些测试,对⼀个正弦信号的⼆次微分做两次积分,正弦频率为50Hz,采样频率
1000Hz,恢复效果如下
时域积分
频域积分
可见恢复信号都很好(对于50Hz是这样的效果)。
分析两种⽅法的频率特性曲线如下
时域积分
频域积分
可以看到频域积分得到信号更好,时域积分随着信号频率的升⾼恢复的正弦幅值会降低。
对于包含两个正弦波的信号,频域积分正常恢复信号,时域积分恢复的⾼频信息有误差 ;对于有噪声的正弦信号,噪声会使积分结果产⽣⼤
的趋势项(不是简单的⼆次趋势),如下图
对此可以⽤滤波的⽅法将⼤的趋势项去掉。
测试的代码如下
% 测试积分对正弦信号的作⽤
clc
clear
close all
%% 原始正弦信号
ts = 0.00 1;
fs = 1/ts;
t = 0:ts:1000
您可能关注的文档
最近下载
- 2024新人教版九年级化学上册《第三单元 物质构成的奥秘》大单元整体教学设计.docx
- 基于化学核心素养的初中化学大单元教学设计.docx VIP
- 小学一年级安全教育《校园安全篇》课件.ppt VIP
- 2024考研408计算机基础综合真题及解析.docx VIP
- 小学新课标体育与健康水平三田径大单元18课时教案.docx
- 政治-黑龙江省哈尔滨市第三中学2024-2025学年高三10月月考试题和答案.docx
- 深静脉血栓VTE防治中心年度工作报告总结.pptx
- 【PPT课件】2024中国慢性阻塞性肺疾病基层诊疗与管理指南解读.pptx
- 骨质疏松症的治疗药物考核试题与答案.docx
- 2024年全国新能源汽车关键技术技能大赛考试题库(试题及答案).docx VIP
文档评论(0)