重采样电阻 非线性后的信号一定和原始数据线性相关吗

君,已阅读到文档的结尾了呢~~
音频重采样算法设计与实现和,音频,算法,和算法,音频重采样,实现音频,和实现,重采样,设计与实现,音频算法
扫扫二维码,随身浏览文档
手机或平板扫扫即可继续访问
音频重采样算法设计与实现
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由:
将文档分享至:
分享完整地址
文档地址:
粘贴到BBS或博客
flash地址:
支持嵌入FLASH地址的网站使用
html代码:
&embed src='/DocinViewer-4.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布,请您等待!
3秒自动关闭窗口查看: 3412|回复: 40
计算阶比跟踪与谐波成分提取整理贴
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
才可以下载或查看,没有帐号?
本帖最后由 impulse 于
11:26 编辑
模拟信号由matlab模拟产生,为一个频率随时间线性增加的信号,基频频率从5Hz变化到30Hz,同时含有0.5倍频和2倍频。
代码如下:
%模拟产生变转速振动信号和键相信号
Ls=;& &%样本数
fs=4096*8;& && &%采样频率
t0 =0:1/fs:(Ls-1)/
t1=t0(length(t0));
f0=5;& &&&%起始频率
f1=30;& & %终止频率
x_1x = 1*chirp(t0,f0,t1,f1,[],280);& && && && & %一倍频信号
x_5x = 0.65*chirp(t0,0.5*f0,t1,0.5*f1,[],115);&&%0.5倍频信号
x_2x = 0.75*chirp(t0,2*f0,t1,2*f1,[],20);& && & %2倍频信号
x_noise=0.15*rand(1,length(x_1x));& && && && &&&%噪声
x=x_1x+x_2x+x_5x;%+x_& && && && && && && &%合成信号
[maxtab, mintab] = peakdet(x_1x, 0.5);&&%通过基频信号(相当于转轴信号)检波,模拟后面的键相信号,函数附在附件中
t=maxtab(:,1)';
tacho=zeros(1,Ls);
tacho(t)=1;& && && && && && && && && &&&%模拟键相信号,假定1为键相脉冲位置
tacho(t+1)=1;
subplot(211)
plot(t0,x);
title('振动信号')
subplot(212)
plot(t0,tacho)
title('转速(键相)信号')
复制代码1、原始信号如下
rawwave.PNG (22.55 KB, 下载次数: 4)
00:30 上传
2、阶比跟踪重采样后的信号
evenwave.PNG (29.7 KB, 下载次数: 5)
00:31 上传
3、原始信号频谱与重采样后信号阶比谱比较
spect.PNG (26.29 KB, 下载次数: 2)
00:32 上传
4、Gabor变换提取谐波波形
gabor.PNG (40.2 KB, 下载次数: 3)
00:33 上传
5、恒带宽零相位滤波器提取谐波波形
cbf.PNG (40.17 KB, 下载次数: 3)
00:34 上传
6、Vold-Kalman滤波提取谐波波形
vkf.PNG (34.73 KB, 下载次数: 1)
00:35 上传
7、上述方法比较,VKF方法堪称完美。
(844 Bytes, 下载次数: 138)
11:25 上传
点击文件名下载附件
赞成: 5 perfect!&
<p id="rate_146" onmouseover="showTip(this)" tip="我来凑热闹,努力向前辈学习&威望 + 2 点
体能 + 5 点
" class="mtn mbn">
本帖被以下淘专辑推荐:
& |主题: 26, 订阅: 2
本帖最后由 ufo158 于
11:43 编辑
顶贴...VKF在计算速度、在临近 交叠阶次的提取上, 效果也很好.
端点上提取结果, 从仿真上看, 优于其他两种, 并且相位想没有偏移。lz VKF带宽怎么取的??常带宽还是百分比?
计算速度上,若用PCG解方程,要慢一些的哈。
小生问题总是很多:
在文献里面看到有一个叫“Constant Percentage Bandwidth”,麻烦解释一下,谢谢。
mxlzhenzhu 发表于
小生问题总是很多:
在文献里面看到有一个叫“Constant Percentage Bandwidth”,麻烦解释一下,谢谢。
恒百分比带宽,滤波器的带宽与中心频率百分比一定。
impulse 发表于
恒百分比带宽,滤波器的带宽与中心频率百分比一定。
非常谢谢你的解释,但是我还是不懂;
一个老师告诉我们滤波器的截止频率就是分析分析频率的最大值。&滤波器的带宽&是不是就是0~fs/2呢,现在的滤波器主要都是数字滤波器了吧?这个带宽就是由那奎斯特频率决定的吧,那中心频率又是啥子呢?
最好提供相关书籍,我自己看去。
mxlzhenzhu 发表于
非常谢谢你的解释,但是我还是不懂;
一个老师告诉我们滤波器的截止频率就是分析分析频率的最大值。&滤波 ...
这是个带通滤波器, 中心频率也就是需要滤出的那部分分量的频率, 滤波带宽即在此中心频率附近也将被滤出的分量的频率范围.
LZ你好,小弟最近在也在学习cot,但是对怎样实现不是很了解。能不能请你分享一下MATLAB程序,让我学习一下阶比跟踪重采样的方法
我也想要matlab程序,学习一下阶比跟踪重采样的方法邮箱
我没有COT的Matlab程序,我的是自己开发的C++商业版本,matlab版本论坛上有一个,虽然效果不是很好,但是基本过程都有
还需要不断学习啊~~~
楼主V5,跟楼主学习太多知识
精华帖。。。学习了
我又来了,我想问下Vold-Kalman filtering的vkf(y,fs,f,p,bw,multiorder)f bw参数应该怎么取?你能给我推荐些这方面的文章吗?谢谢了。
你的这个数据我用VKF是没有提取到你文中的那样效果,而且差距很远.
提取一阶的参数
f0=5;& &&&%起始频率
f1=30;& & %终止频率
inc=(f1-f0)/Ls;
fdn=[f0+inc:inc:f1];
vkf(y,fs,fdn,p,bw);
能给指点下不?
主任可否谈一下:阶比分析及阶比跟踪的具体意义?本人很迷惑,3Q
Powered by>>>>>xa90-BZ72-2
振动信号采集、处理和分析系统
xa90-BZ72-2
振动信号采集、处理和分析系统
产品报价:¥面议
更新时间:访问人数: 356次
产地:国产
厂商性质:生产商品牌:
公司名称: 北京中西远大科技有限公司
型号:xa90-BZ72-
联系人:王璐瑶
(联系我时,请说明是在中国智能制造网上看到的,谢谢!)
该厂商其他产品
振动信号采集、处理和分析系统DASP可自动跟踪信号基频,进行带通滤波。 失真度分析 波形失真度精确分析,谐波阶数可变,提供按国家标准的频域法失真度分析(仅考虑主频的各阶倍频谐波),以及考虑所有非主频倍频的谐波失真的时频法失真度分析。 波形重构分析 对采样后的离散数字信号进行任意的重采样,重采样的采样频率可大于或小于原采样频率,并且两个采样频率可以不是整数倍关系。此重构可以很好地还
振动信号采集、处理和分析系统(16通道,采集频率200KHZ) 型号:xa90-BZ72-2 &库号:M194852 &&硬件USB接口数据采集仪,并口数据采集仪系统功能通道数:16通道、32通道分辨率:16位、12位采样频率:100kHz、200kHz、300kHz程序设计语言主程序及绘图程序:C++语言谱分析及数学预算:FORTRAN数据采集程序:汇编语言软件内容                频谱分析 自功率谱、互功率谱、正富氏变换、逆富氏变换、相干函数、传递函数(频响函数)、自相关、互相关、计算反应谱(含工业规范)数字滤波 富氏带通滤波(低通、带通和高通)、富氏带阻滤波、递归滤波(低通、带通和高通)、褶积滤波(低通、带通 和高通)文件处理 文件加常数、文件乘常数、文件相加、文件相减、文件相乘、文件相除、文件消上限、文件清零、文件叠加、文件取倒数、实虚=&幅相、幅相=&实虚、减均值数学运算 正弦、余弦、正切、双曲正弦、双曲正切、反正弦、反余弦、反正切、常用对数、自然对数、指数、平方根、绝对值、幂、积分、均值积分、微分绘图 时域绘图、频域绘图、时域阵列图、频域阵列图、时域三维绘图、频域周期图、X-Y轨迹图、时域三坐标图、频域三坐标图信号生成 频率扫描函数、幅值扫描函数、渐进渐出方式、生成方波、生成锯齿波、生成正弦波、生成梯形波、生成三角波、生成拍波、生成白噪声、生成复合波形、生工模拟地震波、根据自动频谱生成随机波模态参数识别 频域平滑处理、频域实模态参数识别、振型图形、振型动画积分和微分 频域积分(由加速度、位移)、频域微分(由位移计算速度、加速度)基本处理 归一化、列文件表、文件定义、显示内容、计算均值、计算最大值、文件拷贝、文件移动、删除文件、建立文件、取文件内容、编辑数据内容软件截图工程版测试软件 - 高级型倒频谱分析 冲击响应谱峰值计数分析-波形计数1 统计信号中不同大小的波峰和波谷的数目。常常用于机器寿命设计和试验中,研究载荷的幅值分布特性,即载荷幅值大小及其出现的频次。变程计数分析-波形计数2 统计信号中不同大小的波峰和波谷的数目。在机器寿命设计和试验中,侧重研究载荷时间历程中峰谷之间的距离-即振幅。雨流法计数分析-波形计数3主要是根据研究材料的应力-应变过程进行计数,可以最科学地统计出信号中的所有应力循环和半循环,非常适合于应力疲劳统计和损伤度计算,于是疲劳寿命可以利用累计损伤规律从等幅数据算出。雨流法是当前在机器寿命估计方面应用最为有效的方法之一。小波和小波包分析 典型的时频分析方法,分析结果既包含时间信息,又包含频率信息,适合于非平稳信号分析,尤其对突变信号,可用于机械故障信号(如通过振动)的检测分析、信号识别和降噪、信号滤波检测等。拟小波分析 利用高性能带通卷积滤波实现类似小波方式的分解和重构,效果更好。最大熵和精熵谱分析 信号中任取超短信号(如十几点到几十点)进行分析,用于分析短暂瞬态信号的频率特性,精熵谱还可以消除MEM法幅值不准确以及谱峰分叉的缺点。冲击响应谱分析冲击作用于系统产生的效果,它不是冲击本身的频谱,而是冲击作用下系统响应的频域特性,用于估计系统在冲击作用下的损伤势。倒频谱分析 频谱的再次谱分析,用于在频谱中检测谱的周期性,常用于振动、噪声源识别、机器故障诊断、语音分析、回声剔除等方面的研究,DASP中提供四种倒频谱形式:倒富富、倒熵富、倒富熵和倒熵熵谱。包络谱-共振解调分析用于滚动轴承、齿轮等旋转机械的故障诊断,尤其适合于元件具有表面局部损伤形式的故障。数字跟踪滤波 实现对变频信号进行的带通滤波,对于变频信号(例如升降速过程的转子振动),DASP可自动跟踪信号基频,进行带通滤波。失真度分析 波形失真度精确分析,谐波阶数可变,提供按国家标准的频域法失真度分析(仅考虑主频的各阶倍频谐波),以及考虑所有非主频倍频的谐波失真的时频法失真度分析。波形重构分析对采样后的离散数字信号进行任意的重采样,重采样的采样频率可大于或小于原采样频率,并且两个采样频率可以不是整数倍关系。此重构可以很好地还原信号的特性,而不是简单的插值。非线性灵敏度分析某些传感器具有非线性灵敏度(如热电偶),对于不同大小的物理量具有不同的灵敏度,需要根据其非线性灵敏度性表,将传感器输出的电压量准确地转换为实际物理量。幅频相频曲线测试分析 测试传感器、放大器、滤波器等仪器的频率响应特性,同步测量输入和输出信号,记录各频率点的幅值和相位结果,实时绘制幅度-频率特性曲线和相位-频率特性曲线。&振动信号采集、处理和分析系统(16通道,采集频率200KHZ)&型号:xa90-BZ72-2 &库号:M194852 &&振动信号采集、处理和分析系统(16通道,采集频率200KHZ)&型号:xa90-BZ72-2 &库号:M194852 &&
联系人:王璐瑶
传真:010-
(联系我时,请说明是在上看到的,谢谢!)
产品名称价格地区公司名称更新时间
¥25000北京市
¥25000北京市
面议北京市
面议北京市
面议西安市广东T业大学工程硕士学位论文;视音频格式(wAV)和音频交换文件格式(AIFF;一般音频的常用采样频率包括8kHz,9.6kHz;最后需要考虑的是由有损和可感知压缩算法(如MPE;3.13音频信号的传送环境;实践中,在信号从编码到解码之间有多种可能的传播途;第一种情形是声音文件从一个机器拷贝另一个机器,其;第二种情形是信号仍然保持数字的形式,但采样率发生;第
广东T业大学工程硕士学位论文
视音频格式(wAV)和音频交换文件格式(AIFF)。另一种较低质量音频的量化方式一般采用8bit∥律量化。这些量化方法会使信号产生一些畸变,在8bit∥律量化中显得更为明显。
一般音频的常用采样频率包括8kHz,9.6kHz,10kHz,12kHz,16kHz,22.05kHz和4401kHz。采样频率影响水印数据的隐藏量,因为它给出了可用频谱的上限(如果信号的采样频率为8kHz,则引入的修改分量的频率不会超过4kHz)。对于大多数已有的水印技术而言,可用的数据空间与采样频率的增长至少呈现线性关系。
最后需要考虑的是由有损和可感知压缩算法(如MPEG―AuDIO压缩算法“””71引起的变化。这些变化彻底改变了信号的数据结构。它们仅仅保留了听者能感觉到的特性部分,也就是说,它听起来与原来的相似,但是在最小平方意义上完全不同。数据隐藏的速率依赖于信号的采样率和声音编码的类型和具体的水印算法。
3.13音频信号的传送环境
实践中,在信号从编码到解码之间有多种可能的传播途径。我们只考虑最普通的四种情形p1。
第一种情形是声音文件从一个机器拷贝另一个机器,其中没有任何形式的改变。编码方和解码方的采样率完全一样。
第二种情形是信号仍然保持数字的形式,但采样率发生变化,这…变化保持了大多数信号的幅度和相位,但是改变了信号的时域特征。
第三种情形是信号被转换为模拟形式,通过模拟线路进行传播,在终端被重新采样。在此过程中信号的幅度、量化方式和时域采样率都得不到保持。通常,这种情形下信号的相位值可以得到保持。
第四种情形是信号在空气中传播,经过:麦克风重新采样。这时信号受到未知的非线性改变,会导致相位的改变、幅度改变、不同频率成分的漂移和产生回声等。
在选择水印嵌入算法时,需要考虑信号的表述和传输途径。如果音频信号在传输中没有改变(比如第一种情形),则对水印算法的约束最小。如果音频信
第三章音频数字水印技术
号在传输中发生很大改变(第四种情形),则对水印算法的约束最大,要求算法具有很强的稳健性。
3.1.4对音频数字水印的要求
要想成功地在数字音频媒体中隐藏水印,必须注意以下几方面的要求。1.对数据变换处理操作的稳健性
要求水印本身应能经受得住各种有意无意的攻击。典型的攻击有添加噪声、数据压缩、滤波、重采样、A/D―D/A转换、统计攻击等。
2.听觉相似性
数字水印是在音频载体对象中嵌入一定数量的掩蔽信息,为使得第三方不易觉察这种嵌入信息,需谨慎选择嵌入方法,使嵌入信息前后不产生听觉可感知的变化。这种相似性在理论上可以用“知觉相似性函数”sim(x,y)来描述,实用中利用数字音频中两个信号的相关性函数来估计相似度。
3.是否需要原始数据进行信息提取
根据数据嵌入和提取方案的不同设计,有些方案可以不需要借助于原始数据进行信息提取,这一性能将影响方案的用途和性能。
4.数据提取误码率
数据提取误码率是数据提取结果中错误码位数与嵌入总的位数之比,一般以百分比表示。显然低数据提取误码率也是音频水印方案的一个重要技术指标。因为一方面存在来自物理空间的干扰,另一方面信道中传输的信号会发生衰减和畸变,冉加上人为的数据变换和攻击,都会使数据提取的误码率增加。
5.嵌入数据量指标
根据用途的不同,在有些应用场合中必须保证‘。‘定的嵌入数据量,如利用音频载体进行隐蔽通信。
与图像和视频相比,音频信号在相同的时间间隔内采样的点数少。这使得音频信号可嵌入的信息量要比可视媒体也要少。并且由于人耳听觉系统(}IAs)要比人眼视觉系统(HVs)敏感得多,因此听觉上的不可感知性实现起来要比视觉上困难得多。
广东工业大学工程硕士学位论文
3.1.5数字音频水印的典型应用
随着音频素材在互联网上的增加,数字音频水印技术有着广泛的应用前景:
1)为了便于对音频素材进行查找和检索,可以用水印技术实现元数据的传输,就是用兼容的隐藏的带内方式传送描述性信息。
2)在广播领域中,可以用水印技术执行自动的任务,比如广播节目类型的标识、广告效果的统计分析、广播覆盖范围的分析研究等。其优点是不依赖于特定的频段。
3)用水印技术实现知识产权的保护,包括所有权的证明、访问控制、追踪非法拷贝等。这也是水印技术最初的出发点。
3.2经典的音频水印技术
音频水印技术之间的区别主要体现在数据嵌入/提取方案的不同,早期的方法主要有四种:最低有效位方法、回声隐藏算法、扩频方法、相位编码方法,都是基于时域的方法。分别介绍如下:
3.2.1最低有效位方法
最低有效位(LsB)方法[3],是一种最简单的嵌入方法。任何形式的秘密数据都可以看作是’一串二进制位流,而音频文件的每一个采样数据也是用二进制数来表示,这样可以将每个采样值的最不重要的二进制位(多数情况下为最低位)用代表秘密数据的二进制位代替,以达到在音频信号中嵌入秘密数据的目的。如果将音频信号看作传输的信道,而秘密数据(水印)看作在信道中传输的信号。理想情况下,这种信道容量为lkbps每lkHz,即在无噪声信道中,对8kHz采样率的信号,比特率是8kbps:对44.1kHz采样率的信号,比特率是44.1kbps。伴随这种大信道容量的是可闻噪声的引入,这种噪声的影响效果和音频信号的内容有关。
为了加大对水印攻击的难度,可以使用一段伪随机序列来控制水印的嵌入位置。伪随机信号可以由伪随机序列发生器来产生。当伪随机序列发生器的结构固定时,不同的初始值会产生不同的伪随机序列。这样在收发双方只需要传20
第_三章音频数字水日J技术
送一个初始值(作为密钥),而不需要传送整个伪随机序列值。只要能保证合法用户才能得到该密钥,则根据kerchoff法则可知系统是安全的。任何一个企图提取秘密数据的第三方在不知道密钥的情况下,则不可能达到目的。
这种方法优点是:简单易实现;音频信号罩可编码的数据量大;采用流加密方式分别对数据本身和嵌入过程进行加密,其安全性依赖于密钥,如果选择伪随机性能好的密钥产生机制,则可以做到“~次一密”;信息嵌入和提取算法简单,速度快。缺点:对信道于扰及数据操作的抵抗力很差。如果水印嵌入时不采用冗余技术,噪声、压缩、滤波、重采样都会破坏水印数据;采用冗余技术会使信号速率降低一至二个数量级。在实用中,LsB方法只用于封闭的数字到数字的环境下。
为了提高稳健性,可将水印嵌入到音频信号的较高位,对应于音频信号中的低频分量。但这样带来的结果是大大降低了数据的隐蔽性。为了改善这一点,可以在嵌入过程中根据音频信号的能量进行数据嵌入位的选择,这种方法对平均能量比较高的音频样本更有效。
3.2.2基于回声的水印算法
回声隐藏(EchoHiding)是利用引入回声来将秘密数据嵌入到载体中。它利用人类听觉系统的特性:音频信号在时域的后屏蔽作用,即弱音在强音消失后变得无法听见。弱信号可以在强信号消失后50一200ms的作用面不被人耳觉察。音频信号(载体数据)和经过回声隐藏的水印数据对于人耳来说,前者就像是从耳机里听见的声音,没有回声,而后者就象从扬声器里听到的声音由所处的空间产生的回声。因此,回声算法与其他方法不同,是将秘密数据作为载体数据的环境条件,而不是当作随机噪声嵌入到载体数据中,所以这种方法具有很好的透明性和稳健性。
设音频序列s={s(,z),o≤,z<Ⅳ},按下式即可得到含有回声的音频序列Y:
m,=k‰”聊,,。篡Ⅳ沪,,
其中,m是信号和回声间的延时,一般取m<<Ⅳ,丑为衰减系数。在回声编码中通过修改聊来嵌入水印信息,水印嵌入流程图如图3―1所示。21
广东T业大学工程硕士学位论文
图3―1回声编码水印嵌入流程
FigLlre3-1nowchanofwatermarkembeddingbased帆echohiding
具体方法是:对一个音频数据文件,先将其分成若干包相同样点数的片段,每个片段时间约为几个到几十个毫秒,样点数记为Ⅳ。每段用来嵌入1比特的水印信息。在水印嵌入阶段,对每段信号使用(3―1)式。选择m=m。,则在信号中嵌入水印比特“O”;选择m=聊.,则在信号中嵌入水印比特“l”。延时%或m。是以人耳听不到回声为准则进行选取的。最后将所有含回声的信号段串联成连续信号。
惮夕频Hj吨卫甘
图32回声编码水印提取流程图
Figure3-2nowchan比较萝(蜥o)和歹∽1)ofwate加arkemactingbasedonechohiding
回声编码中水印提取流程图如图3―2所示。对一个音频回声信号,水印的提取关键是确定回声的延时。利用复倒谱可将回声从原始信号中分离出来,但是代表回声的脉冲幅值与载体信号相关度很小。因此它们很难被检测到。使用复倒谱的自相关可以解决这个问题。设回声信号yQ)的复倒谱自相关为萝@)。可以证明,多(”)在回声延时出现极大值,由于引入回声的延时只有坍。和m.两
包含各类专业文献、外语学习资料、中学教育、专业论文、各类资格考试、应用写作文书、幼儿教育、小学教育、高等教育、31基于小波分析的音频数字水印技术研究等内容。 
 41 I II 基于小波变换的数字水印算法研究摘 要数字水印技术是目前信息安全技术...MusiCode 系统提供了音频信号的广播监测,VEIL-II 和 MediaTrax 提供了视频信号...  概括 了数字水印的研究现状和基本原理 , 分析了基于小波变换的数字水印算法的...音频、视频、文本等形式的数字作品极易 得到非法复制和拷贝, 使得盗版侵权日益...  这些行为严重损害了数字作品创 1 基于小波分析的图像水印算法研究 作者和使用者...到目前为止,数 字水印从研究对象上看主要涉及图像水印、视频水印、音频水印、...  基本原理、典型算法几数字水印攻击方法,研究了一 种基于离散小波变换的数字水印...音频和视 频这些信号非常方便,尤其是图像矩阵运算更是方便快捷.另外 MATLAB 内置...  基于小波变换的数字水印技术_信息与通信_工程科技_专业资料。届别学号 2011 届
毕业设计(论文) 基于小波变换的数字水印技术的研究 姓系 名 刘 敏 别...  本文针对基于小波变换的数字水印技术, 提出了一种基于小波域的二值图像 水印算法...由于小波变换的优点, 使小波变换域研究水印处理技术是目前的热点,并且在该 第 ...  计算机毕业设计论文_基于小波变换的数字水印技术研究_工学_高等教育_教育专区 暂无评价|0人阅读|0次下载 计算机毕业设计论文_基于小波变换的数字水印技术研究_工学_...  基于小波变换的静态图像数字水印技术的研究_工学_高等教育_教育专区。I 基于小波...这里的载 体可以是图像、视频、音频,也可以是信道,甚至是某套编制体质和整个...  基于小波变换的数字水印技术研究摘要:算法基于Haar小波变换,把小波系数分块,并计算每个块的平 均值。在一系列信号处理之后,这些块(尤其是大的块)的平均值不会 有...}

我要回帖

更多关于 重新采样双线性三线性 的文章

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信