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

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

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

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

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 [amp phase]=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上请教你