版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、信号时域与频域分析实验报告姓名: 杨班级: 机械学号: 213实验数据中,电机转速为1200r/min,采样频率为1280Hz。Hz3为X位移振幅数据,Hz4为Y位移振幅数据,Hz5为速度振幅数据。Matlab中信号特征对应函数编程ma = max(Hz) %最大值mi = min(Hz) %最小值me = mean(Hz) %平均值pk = ma-mi %峰-峰值va = var(Hz); %方差st = std(Hz); %标准差ku = kurtosis(Hz); %峭度rm = rms(Hz); %均方根一、X轴位移测量分析plot(Fs3,Hz3) 时域图:ma =52.0261mi
2、 =56.7010me =1.8200pk =108.7271va =1.3870e+03st =37.2431ku =1.5462rm =37.2693 频域图:fs=1280;x=Hz3;N=length(Hz3);df=fs/N;f=0:df:N*df-df;y=fft(x);y=abs(y)*2/N;figure(1);plot(f,y);xlabel('频率/Hz')ylabel('幅值')频谱幅值取得最大值51.9847um,频率为20Hz,与电机转速对应频率一致,应为电机轴未动平衡所致;二倍频处有较大振幅,可能为轴承间隙过大所致。二、Y轴位移测量分
3、析plot(Fs4,Hz4) 时域图:ma =61.3987mi =-74.6488me =-1.1948pk =136.0475av =42.6109va =2.2428e+03st =47.3582ku =1.5135rm =47.3501 频域图:fs=1280;x=Hz4;N=length(Hz4);df=fs/N;f=0:df:N*df-df;y=fft(x);y=abs(y)*2/N;figure(1);plot(f,y);xlabel('频率/Hz')ylabel('幅值')频谱幅值取得最大值66.6319um,频率为20Hz,与电机转速对应频率一
4、致,应为电机轴未动平衡所致;二倍频处有较大振幅,可能为轴承间隙过大所致。三、速度测量分析plot(Fs5,Hz5) 时域图:ma =1.0907mi =0.1652me =0.6323pk =0.9255av =0.6323va =0.0397st =0.1993ku =1.9936rm =0.6630 频域图:fs=1280;x=Hz5;N=length(Hz5);df=fs/N;f=0:df:N*df-df;y=fft(x);y=abs(y)*2/N;figure(1);plot(f,y);xlabel('频率/Hz')ylabel('幅值')速度频谱幅值取
5、得最大值1.2647mm/s,频率为20Hz,与电机转速对应频率一致。振幅主要集中在20Hz附近,高频较为稳定。速度自相关函数图:fs=1280;x=Hz;N=length(Hz);df=fs/N;f=0:df:N*df-df;a,b=xcorr(Hz,'unbiased');plot(b*df,a);title('速度自相关函数');xlabel('t');四、轴心轨迹测量分析figure(1);plot(Hz3,Hz4)grid on;轴心轨迹图:五、思考题分析1. (1)A,B的波形图:fs=1024;t=0:1/fs:1;x=10*sin
6、(60*pi*t+pi/4)+5*sin(120*pi*t+pi/3);y=12*sin(60*pi*t+pi/2)+4*sin(120*pi*t+pi/6);subplot(2,1,1);plot(t,x);title('A的波形图')xlabel('时间/s')ylabel('幅值')subplot(2,1,2);plot(t,y);title('B的波形图')xlabel('时间/s') ylabel('幅值') figure(2);plot(x,y);title('轴心轨迹'
7、);(2)B的幅值谱,相位谱,功率谱:fs=1024;df=1/fs;t=0:df:1;N=length(t);y=12*sin(60*pi*t+pi/2)+4*sin(120*pi*t+pi/6);f=fft(y,N);mag=abs(f)*2/N;phase=angle(f);power=mag.2;F=0:df*N:(N-1)/N/df;figure(1);plot(F(1:(N+1)/2),mag(1:(N+1)/2);title('B的幅值谱');figure(2);plot(F(1:(N+1)/2),phase(1:(N+1)/2);title('B的相位谱
8、');figure(3);plot(F(1:(N+1)/2),power(1:(N+1)/2);title('B的功率谱');B的幅值谱:B的相位谱:B的功率谱:(3)A和B的互相关函数,A的自相关函数:fs=1024;df=1/fs;t=0:df:1;N=length(t);x=10*sin(60*pi*t+pi/4)+5*sin(120*pi*t+pi/3);y=12*sin(60*pi*t+pi/2)+4*sin(120*pi*t+pi/6);a,b=xcorr(x,y,'unbiased');figure(1);plot(b,a);title(&
9、#39;A,B的互相关图'); a,b=xcorr(x,'unbiased');figure(2);plot(b,a);title('A的自相关图'); A、B的互相关图:A的自相关图:2. 对下列信号进行短时傅里叶变换及EMD分解建立sa.m文件,建立函数function sd=sa(a,b,f,ang,t) if t>bsd=sin(2*pi*f*t+ang)*exp(a*(t-b);else sd=0;endt=0:0.001:2;x=8*sin(60*pi*t+pi/12)+4*sin(120*pi*t+pi/2)+sin(300*pi*t+pi/6)+20*sa(20,0.2,250,pi/3,t);Spec,Freq=STFT(x,542,64,1);短时傅里叶变换图:t=0:0.
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 广播稿400字左右(35篇)
- 高中技术《技术与设计2》模块测试题 一
- 课外书的心得体会范文
- 幼儿园主题方案简单
- 风险合规部工作总结
- 销售培训心得(35篇)
- 居间代理房屋合同(3篇)
- 《技术的未来》教学设计(两篇)
- 苏教版 高中技术《技术与设计1》教案合集
- 26.1 锐角三角函数 同步练习
- 石材加工工艺及技术标准
- 六年级上册美术课件-第1课 建筑艺术的美 ▏人美版 (共27张PPT)
- 医院引流管护理考核评价标准
- 环境规划与管理全套课件完整版电子教案最新板
- 苏教版二年级数学上册《认识线段》课件(市级赛课一等奖)
- 输电线路初步设计评审要点课件
- (完整word版)小餐饮经营食品安全管理制度
- 产后尿潴留的护理个案课件
- 装配式混凝土结构部件吊装监理细则
- 地铁站装饰施工组织设计(181页)
- 杨柳煤矿“三量”动态变化情况分析报告(3)
评论
0/150
提交评论