版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
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. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- DB11T 1087-2014 公共建筑装饰工程质量验收标准
- 河南省周口市淮阳实验小学等校2024-2025学年二年级上学期月考语文试卷(10月份)
- 湖北省黄冈市部分学校2024-2025学年七年级上学期期中地理试卷(含答案)
- 家用电动蒸汽拖把产业深度调研及未来发展现状趋势
- 叶酸膳食补充剂市场发展预测和趋势分析
- 提供婚礼司仪服务行业市场调研分析报告
- 加热板市场发展预测和趋势分析
- 婴儿绒线鞋产业链招商引资的调研报告
- 化妆用杏仁油产业规划专项研究报告
- 医用洗液产业深度调研及未来发展现状趋势
- (完整word版)水利工程设计概(估)算编制规定
- (完整word版)小餐饮经营食品安全管理制度
- 产后尿潴留的护理个案课件
- 中医护理发展史课件(PPT 35页)
- 装配式混凝土结构部件吊装监理细则
- 地铁站装饰施工组织设计(181页)
- 色彩的基础知识课件.PPT
- 杨柳煤矿“三量”动态变化情况分析报告(3)
- 动火作业及动火工作票管理规定
- 变电站综合自动化电子教案
- 桥梁伸缩缝施工及质量保证要点
评论
0/150
提交评论