国产精品乱码一区-性开放网站-少妇又紧又爽视频-西西大胆午夜人体视频-国产极品一区-欧美成人tv-四虎av在线-国产无遮挡无码视频免费软件-中文字幕亚洲乱码熟女一区二区-日产精品一区二区三区在线观看-亚洲国产亚综合在线区-五月婷婷综合色-亚洲日本视频在线观看-97精品人人妻人人-久久久久久一区二区三区四区别墅-www.免费av-波多野结衣绝顶大高潮-日本在线a一区视频高清视频-强美女免费网站在线视频-亚洲永久免费

 找回密碼
 注冊會員

QQ登錄

只需一步,快速開始

搜索
查看: 5092|回復: 4

[matlab] matlab處理實驗數(shù)據(jù) 實驗:利用轉子試驗臺測試轉子不平衡以及渦動狀態(tài) 故障診斷

[復制鏈接]
1#
發(fā)表于 2011-5-26 11:35:10 | 只看該作者 |倒序瀏覽 |閱讀模式
我有實驗數(shù)據(jù)是dat格式的 要用matlab分析出頻譜圖跟時域圖要得出的圖像跟故障應該有的圖像相似并且?guī)臀医忉屢幌?轉子不平衡的時域圖 頻譜圖各是什么特征我編的程序如下 你看看用得著不%***************1.bch0****************%clcclfclearbch0=load('F:\MATLAB7\lijiande\bch0.dat');fs=100;%設定采樣頻率N=128;n=0:N-1;t=n/fs;bch0=bch0(n+1);%采樣figure(1);subplot(211);plot(t,bch0);%作正弦信號的時域波形axis([0,N/fs,-0.4,0.4]);xlabel('t');ylabel('y');title('bch0采樣t=[0,1.28]');grid;%進行FFT變換并做頻譜圖y=fft(bch0,N);%進行fft變換mag=abs(y);%求幅值f=(0:length(y)-1)'*fs/length(y);%進行對應的頻率轉換figure(1);subplot(212);plot(f,mag);%做頻譜圖axis([0,100,0,4]);xlabel('頻率(Hz)');ylabel('幅值');title('bch0幅頻譜圖N=128');grid;5 b6 k4 o4 n) b% [, K# ]
謝謝大家了~!!!

本帖子中包含更多資源

您需要 登錄 才可以下載或查看,沒有賬號?注冊會員

×
回復

使用道具 舉報

2#
 樓主| 發(fā)表于 2011-5-26 11:47:56 | 只看該作者
我有實驗數(shù)據(jù)是dat格式的 要用matlab分析出頻譜圖跟時域圖6 U( D: e" g) z: L* [
要得出的圖像跟故障應該有的圖像相似
1 {/ l- T9 o+ H并且?guī)臀医忉屢幌?轉子不平衡的時域圖 頻譜圖各是什么特征, ?2 B% Q- n, @# s5 i
我編的程序如下 你看看用得著不/ i# U- @8 k' q7 W
+ ~: F, z: a! \, V: n! m6 @1 g
%***************1.bch0****************%  ?+ z9 a8 k3 Z& L8 }
clc* s4 T# x' B" ~4 n7 T, h% A. r. |
clf* w# F' e  y- g+ l' _) B
clear$ b: o9 D) E2 X( d& C) B
bch0=load('F:\MATLAB7\lijiande\bch0.dat');! N( m6 {! |- i+ D( p2 z
fs=100;%設定采樣頻率
0 @; K* a, b1 o. n$ ^N=128;- h! _9 y: l( k" ]
n=0:N-1;
" `5 K1 k6 J4 B2 p5 Tt=n/fs;
. {' O. Z- ]3 ]- lbch0=bch0(n+1);" d4 G5 `* K( x/ @8 o& e9 P4 f
%采樣; {$ _  M, I, n7 t5 T  d6 m( R
figure(1);0 `( ^  d, P8 v6 F
subplot(211);! S; p* r  I! C& B
plot(t,bch0);%作正弦信號的時域波形
; m* y: T" V" R2 k& X/ Z9 G2 o" C/ @axis([0,N/fs,-0.4,0.4]);
* `1 i( S* X) T  N$ exlabel('t');- i; p' a7 Z7 h: @4 g
ylabel('y');* d' K4 c3 ]9 s
title('bch0采樣t=[0,1.28]');6 o2 F7 {# l/ V2 n5 `
grid;
0 F5 S; D" j( A3 ^. e) l%進行FFT變換并做頻譜圖$ z3 u3 C8 M, k. U, E& d
y=fft(bch0,N);%進行fft變換& z- g, e" X" A! |% T4 y  d
mag=abs(y);%求幅值4 N4 J. S# E1 v- ?3 |$ z1 l
f=(0:length(y)-1)'*fs/length(y);%進行對應的頻率轉換9 z- B* ^0 e, x! e
figure(1);4 H5 b9 B7 }; s2 b0 g
subplot(212);. F! b: b4 ?1 N! `
plot(f,mag);%做頻譜圖# M8 c' z6 r5 j7 j, O3 l: `  V
axis([0,100,0,4]);
0 d. |7 [+ O( J" H" Xxlabel('頻率(Hz)');* ?9 J" n  ~) t) t/ F
ylabel('幅值');# C0 ~2 f5 o% {# ]
title('bch0幅頻譜圖N=128');' G7 t1 ~7 b) g% L
grid;
# d, ~  c) w0 R2 J
; z4 ~' R. t2 u不好意思我寫在寫字板粘貼過來的不知道為什么標點沒了 現(xiàn)在改過來~
3#
發(fā)表于 2011-5-26 12:02:38 | 只看該作者
應該和轉子模型以及試驗工裝模型有關,跨距等。可傳一試驗工裝模型圖,以便分析。
4#
發(fā)表于 2012-5-29 16:45:25 | 只看該作者
你的轉子工作頻率多少
5#
發(fā)表于 2012-5-31 14:35:19 | 只看該作者
把數(shù)據(jù)下下來看看再說吧
您需要登錄后才可以回帖 登錄 | 注冊會員

本版積分規(guī)則

Archiver|手機版|小黑屋|機械社區(qū) ( 京ICP備10217105號-1,京ICP證050210號,浙公網(wǎng)安備33038202004372號 )

GMT+8, 2025-9-3 10:26 , Processed in 0.059985 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.5 Licensed

© 2001-2025 Discuz! Team.

快速回復 返回頂部 返回列表