车轮村庄 发表于 2011-10-28 11:24

旋转机械波德图的做法,或者是源代码

请问有没有做过旋转机械伯德图的朋友 ,如果有的话,请教下做法!

impulse 发表于 2011-10-29 11:08

回复 1 # 车轮村庄 的帖子

知道bode图纵横坐标的意义不就会画了嘛。

车轮村庄 发表于 2011-10-29 12:19

回复 2 # impulse 的帖子

具体做的时候出现了问题?幅值能求出来 ,但是做相位的时候有问题 ,等不到准确的值,我用了内插法校正FFT任然不能满足要求。

impulse 发表于 2011-10-29 14:13

回复 3 # 车轮村庄 的帖子

这是算法本身问题,不是画图问题

车轮村庄 发表于 2011-10-30 09:03

回复 4 # impulse 的帖子

对,就是这个意思 ,现在的问题是吗,我想知道大概用什么办法去做 。我做相位用的fft内插法校正,在人家的论文里面看人家也是这么做的 ,但得不到想要的结果 。我就想问下看大家还知不知道其他的方法,或者说我的方法哪里错了,我好修改。

impulse 发表于 2011-10-30 22:57

本帖最后由 impulse 于 2011-10-30 23:03 编辑

回复 5 # 车轮村庄 的帖子

%计算振动某阶特征量(幅值、相位)
%wave-振动波形
%sample_per_cycle-每转样本点数
%order-阶比
function =getWaveAttr_Order(wave,sample_per_cycle,order)
a=0;
b=0;
deg0=2.0*pi/sample_per_cycle;
deg=0;
sample_len=length(wave);
for i=1:sample_len
    b=b+wave(i)*sin(deg*order);
    a=a+wave(i)*cos(deg*order);
    deg=deg+deg0;
end
amp=2*2*sqrt(a*a+b*b)/sample_len;
if(b>1e-10 | b<-1e-10)
    ph=atan(a/b)*180/pi;
    if(b>0)
      ph=ph+270;
    end
    if(b<0)
      ph=ph+90;
    end
else
    if(a>0)
      ph=0;
    else
      ph=180;
    end
end
phase=360-ph;给你一段代码,直接从波形计算某阶振动幅值、相位。可以减小采用FFT方法时非整周期采样带来的误差。

车轮村庄 发表于 2011-11-1 16:49

回复 6 # impulse 的帖子

谢谢!我先试试

车轮村庄 发表于 2011-11-1 17:33

回复 6 # impulse 的帖子

仔细看了下你给的程序 ,基本上可以算是稳定转速的情况,或者说是等角度采样的情况。我想请教下你 ,旋转机械启停过程中 ,转频变化的情况下怎样得到转频相关的幅值和相位!我现在是不知道该怎样去实现这个了,用FFT分小段每段1024点(信号是恒时间间隔采样的)做频谱,能够取得幅值,与东方所的软件对比临界转速只相差3转/Min。现在是不知道相位该怎么弄

impulse 发表于 2011-11-1 17:46

本帖最后由 impulse 于 2011-11-1 17:48 编辑

回复 8 # 车轮村庄 的帖子

旋转机械振动为了减少泄露,基本都是采用等角度采样,硬件或者软件方式,目前国外主流是软件重采样。你要是采用非同步采样,首先要得到转频,其次要采用校正算法。

车轮村庄 发表于 2011-11-1 23:24

回复 9 # impulse 的帖子

恩 ,我也是这样做的 ,首先计算到转频,然后对振动(相位)信号做FFT,找到计算得到的转频附近的最大幅值点,做内插法FFT校正等到相位和幅值,但是结果不对。我的校正法用模拟的正弦信号验证了,是正确的。我是合肥工业大学毕业的,我看你现在也在合肥,呵呵,是不是知道合肥工大啊!我的QQ是244299541,希望能在QQ上请教你

impulse 发表于 2011-11-1 23:53

回复 10 # 车轮村庄 的帖子

你的转频怎么计算的?如果同时采集了转速键相信号,完全可以采用重采样方法实现等角度采样。

车轮村庄 发表于 2011-11-2 11:26

回复 11 # impulse 的帖子

转频是通过简单的数一小段时间的脉冲个数得到的,比较粗糙。

impulse 发表于 2011-11-2 13:00

回复 12 # 车轮村庄 的帖子

既然这样,你可以采用软件重采样来做,看过重大的一篇硕士论文做过这个。

车轮村庄 发表于 2011-11-2 17:01

回复 13 # impulse 的帖子

恩 ,找找看

车轮村庄 发表于 2011-11-2 22:01

回复 14 # 车轮村庄 的帖子



上面两图是我自己做的和对比软件做的 ,幅值还可以相位完全不对 。如果用你说的论文的那种方法貌似太复杂了点 ,参考软件好像是用的FFT做的 ,因为他的界面上这个    ,纠结就在这里。我用FFT做的就是不对 ,不知道问题的根源 ,急需指导
页: [1] 2
查看完整版本: 旋转机械波德图的做法,或者是源代码