久久久国产一区二区_国产精品av电影_日韩精品中文字幕一区二区三区_精品一区二区三区免费毛片爱

 找回密碼
 注冊(cè)會(huì)員

QQ登錄

只需一步,快速開(kāi)始

搜索
查看: 5070|回復(fù): 4

[matlab] matlab處理實(shí)驗(yàn)數(shù)據(jù) 實(shí)驗(yàn):利用轉(zhuǎn)子試驗(yàn)臺(tái)測(cè)試轉(zhuǎn)子不平衡以及渦動(dòng)狀態(tài) 故障診斷

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

本帖子中包含更多資源

您需要 登錄 才可以下載或查看,沒(méi)有賬號(hào)?注冊(cè)會(huì)員

×
回復(fù)

使用道具 舉報(bào)

2#
 樓主| 發(fā)表于 2011-5-26 11:47:56 | 只看該作者
我有實(shí)驗(yàn)數(shù)據(jù)是dat格式的 要用matlab分析出頻譜圖跟時(shí)域圖
5 S; h$ l+ I! m要得出的圖像跟故障應(yīng)該有的圖像相似
1 W* J" A/ l3 P/ q. h) n并且?guī)臀医忉屢幌?轉(zhuǎn)子不平衡的時(shí)域圖 頻譜圖各是什么特征# L4 X% k; p7 y2 O) m
我編的程序如下 你看看用得著不
4 _, i) h& M& \* G' S  m+ U7 \3 x# T$ u) n- ]+ s
%***************1.bch0****************%
0 n5 _. P* r" f' @5 A0 D/ V6 w+ jclc& Q  Z# s3 k$ V
clf
4 A1 ]9 R0 A- G% e3 U' W. O' Iclear. g0 x& a. @8 Q$ n! M6 d
bch0=load('F:\MATLAB7\lijiande\bch0.dat');3 ^' O! e. h' M+ ?  f( w7 E
fs=100;%設(shè)定采樣頻率
# s8 y$ U8 p& h$ l/ U: Q5 gN=128;
2 m4 [0 R9 z& w. P+ on=0:N-1;* ?2 {7 ^% p/ x6 j- B' ]  d
t=n/fs;
' q8 u; J$ M6 e" v7 Q' Xbch0=bch0(n+1);
& ^4 U4 q$ u, Y/ e7 h%采樣( e. o' u4 G3 }# J
figure(1);$ Z( F1 f7 r5 z, `! ?
subplot(211);
+ |- [4 x  x* M+ Q! vplot(t,bch0);%作正弦信號(hào)的時(shí)域波形1 e1 {% w$ a0 X- L/ ~& ]
axis([0,N/fs,-0.4,0.4]);
8 e2 H+ k, Y5 O8 r: P$ Xxlabel('t');
* x) A' M4 @* S* P0 Pylabel('y');
) b! V! r+ j0 E( q0 u' Atitle('bch0采樣t=[0,1.28]');* s5 w5 e9 n' e. ?6 P- t* b+ ~
grid;7 R9 s& `+ V- ]& a( P
%進(jìn)行FFT變換并做頻譜圖% t$ L7 z2 o7 Z1 K% |9 ^/ Q
y=fft(bch0,N);%進(jìn)行fft變換
5 q$ \/ f3 I& J) W/ jmag=abs(y);%求幅值! S* Y1 a! c5 R2 {/ g
f=(0:length(y)-1)'*fs/length(y);%進(jìn)行對(duì)應(yīng)的頻率轉(zhuǎn)換* m8 ^3 f4 k, \
figure(1);
! E2 \% P7 h/ ?; g9 i3 g) csubplot(212);" X' h+ q3 l( x' A3 A: {7 u# J
plot(f,mag);%做頻譜圖& ]; a+ H) _  w3 e) b( u+ z
axis([0,100,0,4]);
" n0 s( _4 c0 W: zxlabel('頻率(Hz)');
  F  |% t  f. l8 s4 ^' ?  @6 zylabel('幅值');2 |( r1 B" b- O; m$ \, Q" _
title('bch0幅頻譜圖N=128');8 f; M) ~. e' p5 N0 o+ J; H* H
grid;) K/ T& S5 {/ ~  e
% s; Z9 _# m- S9 M& w
不好意思我寫在寫字板粘貼過(guò)來(lái)的不知道為什么標(biāo)點(diǎn)沒(méi)了 現(xiàn)在改過(guò)來(lái)~
3#
發(fā)表于 2011-5-26 12:02:38 | 只看該作者
應(yīng)該和轉(zhuǎn)子模型以及試驗(yàn)工裝模型有關(guān),跨距等。可傳一試驗(yàn)工裝模型圖,以便分析。
4#
發(fā)表于 2012-5-29 16:45:25 | 只看該作者
你的轉(zhuǎn)子工作頻率多少
5#
發(fā)表于 2012-5-31 14:35:19 | 只看該作者
把數(shù)據(jù)下下來(lái)看看再說(shuō)吧

本版積分規(guī)則

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

GMT+8, 2025-8-21 22:08 , Processed in 0.082122 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.5 Licensed

© 2001-2025 Discuz! Team.

快速回復(fù) 返回頂部 返回列表