机械工程测试技术实验报告.docx
《机械工程测试技术实验报告.docx》由会员分享,可在线阅读,更多相关《机械工程测试技术实验报告.docx(100页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、机械工程测试技术实验报告机械工程测试技术试验报告 本文关键词:机械工程,试验,测试,报告,技术机械工程测试技术试验报告 本文简介:西华高校试验报告试验报告(理工类)课程名称:机械工程测试技术课程代码:8400271学生所在学院:机械工程与自动化学院年级/专业/班:09机电3班学生姓名:王泽学号:312022080307330试验总成果:任课老师:余愚开课学院:机械工程与自动化学院试验中心名称:机械工程专业试验中心第组.西华高校机械工程测试技术试验报告 本文内容:西华高校试验报告实验报告(理工类)课程名称:机械工程测试技术课程代码:8400271学生所在学院:机械工程与自动化学院年级/专业/班:
2、09机电3班学生姓名:王泽学号:312022080307330试验总成果:任课教师:余愚开课学院:机械工程与自动化学院试验中心名称:机械工程专业试验中心第组.西华高校试验报告(机械类)开课学院及试验室:机械工程与自动化学院计算机机房试验时间:2022年5月11日学生姓名王泽学号312022080307330成绩学生所在学院机械工程与自动化学院年级/专业/班09机电3班课程名称机械工程测试技术课程代码8400271试验项目名称周期信号波形的合成和分解项目代码指导教师余愚项目学分一、试验目的1.学习运用Matlab,学会用Matlab供应的函数对信号进行频谱分析;2.加深了解信号分析手段之一的傅立
3、叶变换的基本思想和物理意义;3.视察和分析由多个频率、幅值和相位成肯定关系的正弦波叠加的合成波形;。4.视察和分析频率、幅值相同,相位角不同的正弦波叠加的合成波形;5.通过本试验熟识信号的合成、分解原理,了解信号频谱的含义。二、试验原理按傅立叶分析的原理,任何周期信号都可以用一组三角函数、的组合表示(n=1,2,3,)也就是说,我们可以用一组正弦波和余弦波来合成随意形态的周期信号。对于典型的方波,其时域表达式为:依据傅立叶变换,其三角函数绽开式为:由此可见,周期方波是由一系列频率成分成谐波关系,幅值成肯定比例的正弦波叠加合成的。三、试验内容1.用Matlab编程,绘出7次谐波叠加合成的方波波形
4、图及幅值谱;2.用Matlab编程,变更上述7次谐波中其中两项谐波的幅值绘出合成波形及幅值谱;3用Matlab编程,变更上述7次谐波中其中一项谐波的相位绘出合成波形及幅值谱。四、试验程序及结果分析1、内容1的程序:t=0:0.001:0.512;y=sin(2*pi*50*t)+sin(3*2*pi*50*t)/3+sin(5*2*pi*50*t)/5+sin(7*2*pi*50*t)/7;Y=fft(y);%时间抽取基2法FFT,计算序列x的傅里叶变换Mx=abs(Y)/400;%求肯定值,获得幅值谱f=1000*(0:255)/512;subplot(211),plot(t(1:60),y
5、(1:60),title(Time-domainsignal)subplot(212),plot(f,Mx(1:256),title(Spectrum)波形图及其幅值谱:结论:所谓方波可以由傅立叶绽开得到多个正弦谐波,频率最低的那个叫基波,其他的高频的叫谐波。利用奇数次谐波可以合成方波。并且可以将方波分解成一个直流重量和多数个正弦波重量之和。2、内容2的程序:t=0:0.001:0.512;y=sin(2*pi*50*t)+sin(3*2*pi*50*t)/3+5*sin(5*2*pi*50*t)/5+7*sin(7*2*pi*50*t)/7;Y=fft(y);Mx=abs(Y)/400;f=
6、1000*(0:255)/512;subplot(211),plot(t(1:60),y(1:60),title(Time-domainsignal)subplot(212),plot(f,Mx(1:256),title(Spectrum)波形图及其幅值谱:结论:变更上述7次谐波中其中两项谐波的幅值会使输出的方波的幅值也随之发生变更,若变更的幅值不合适,将其合成的方波的幅值具有不确定性,方波不精确性增大。3、内容3的程序:t=0:0.001:0.512;y=sin(2*pi*50*t)+sin(3*2*pi*50*t)/3+sin(5*2*pi*50*t+30)/5+sin(7*2*pi*50
7、*t)/7;Y=fft(y);%时间抽取基2法FFT,计算序列x的傅里叶变换Mx=abs(Y)/400;%求肯定值,获得幅值谱f=1000*(0:255)/512;subplot(211),plot(t(1:60),y(1:60),title(Time-domainsignal)subplot(212),plot(f,Mx(1:256),title(Spectrum)波形图及其幅值谱:结论:变更上述7次谐波中其中一项谐波的相位,会使合成的方波的拐点的斜率发生变更,若变更的相位不合适,则可能使合成的方波更接近于谐波即正弦波,使合成的方波不精确。第组.西华高校试验报告(机械类)开课学院及试验室:机
8、械工程与自动化学院计算机机房试验时间:2022年5月12日学生姓名王泽学号312022080307330成绩学生所在学院机械工程与自动化学院年级/专业/班09机电3班课程名称机械工程测试技术课程代码8400271试验项目名称用FFT对信号进行频谱分析项目代码指导教师余愚项目学分一、试验目的1.学习运用Matlab,学会用Matlab供应的函数对信号进行频谱分析;2.驾驭采样定理;3.理解加窗对频谱分析的影响;4.理解量化误差对频谱分析的影响;5.驾驭采样点数N、采样频率、数据长度对频谱分析的作用。二、试验原理原理:机械工程测试技术与信号分析第2章,特殊是2.4离散傅立叶变换的内容。设备:PC机
9、;软件:Matlab三、试验内容1.画出x(t)=3sin(2ft)+7sin(10ft)+12sin(15ft)的幅值谱图(f=50Hz)。2.用Mablab设计一程序,能形象地验证离散傅里叶变换中的4个重要问题:(1)采样定理a),其频谱不失真,其频谱失真;b)(工程中常用),可从频域中不失真复原原时域信号;(2)加窗、截断a)信号截断后,其频谱会产生泄漏,出现“假频”;b)信号截断后,降低了频率辨别率;c)采纳适当的窗函数后,可以削减泄漏和提高频率辨别率。(3)量化误差a)对信号进行采样,Hz,采集N64点。用3、8位量化器量化信号每点的幅值,画出原始波形和量化后的信号波形,得出结论。(
10、4)栅栏效应如何才能提高频率辨别率?采样点数N、采样频率起何作用?用例子说明。四、试验程序及结果分析1、内容1的程序:t=0:0.001:0.512;y=3*sin(2*pi*50*t)+7*sin(10*pi*50*t)+12*sin(15*pi*50*t);Y=fft(y);Mx=abs(Y)/400;f=1000*(0:255)/512;plot(f,Mx(1:256),title(Spectrum)Mx=abs(Y)/400;f=1000*(0:255)/51幅值谱:五、试验感想2、内容2的程序:(1)采样定理(1)、试验代码fs=100;N=128;%采样频率和数据点数n=0:N-1
11、;t=n/fs;%时间序列x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t);%信号y=fft(x,N);%对信号进行快速Fourier变换mag=abs(y);%求得Fourier变换后的振幅f=n*fs/N;%频率序列subplot(2,2,1),plot(f,mag);%绘出随频率改变的振幅xlabel(频率/Hz);ylabel(振幅);title(N=128);gridon;subplot(2,2,2),plot(f(1:N/2),mag(1:N/2);%绘出Nyquist频率之前随频率改变的振幅xlabel(频率/Hz);ylabel(振幅);title(N
12、=128);gridon;%对信号采样数据为1024点的处理fs=100;N=1024;n=0:N-1;t=n/fs;x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t);%信号y=fft(x,N);%对信号进行快速Fourier变换mag=abs(y);%求取Fourier变换的振幅f=n*fs/N;subplot(2,2,3),plot(f,mag);%绘出随频率改变的振幅xlabel(频率/Hz);ylabel(振幅);title(N=1024);gridon;subplot(2,2,4)plot(f(1:N/2),mag(1:N/2);%绘出Nyquist频率之前
13、随频率改变的振幅xlabel(频率/Hz);ylabel(振幅);title(N=1024);gridon;(2)Matlab运行仿真图像:其次:用100Hz对信号进行采样源信号为f(t)=0.5*sin(2*pi*15*t1)+2*sin(2*pi*40*t1),用100Hz的频率对f(t)进行采样,其采样图如图1所示,程序如下:(1)、试验代码:fs1=100;t1=-0.1:1/fs1:0.1;fa=0.5*sin(2*pi*15*t1)+2*sin(2*pi*40*t1);%x=0.5*sin(2*pi*15*t)+2*sin(2*pi*40*t)figure(1);plot(t1,f
14、a),xlabel(fs1=100Hz时,fa采样时域图)对信号进行快速离散傅立叶变换将采样信号进行快速离散傅立叶变换(FFT),用300Hz的频率对f(t)进行采样,其采样后快速傅立叶变换频谱图如图f=40;fs=100;N=100;k=0:N-1;t=-0.1:1/fs:0.1;w1=500*k/N;fa=0.5*sin(2*pi*f*t)+2*sin(2*pi*f*t);xfa=fft(fa,N);xf1=(xfa);figure(1);plot(w1,xf1),xlabel(fs=100Hz时,fa经过fft后频谱图.单位:Hz)信号的重建我们可以通过利用内插法把原信号从采样信号中复原
15、出来,视察信号在满意怎样的采样条件下能够复原原信号,下图为复原后的信号.程序如下Wm=180*pi;Wc=Wm;fs=100;Ws=2*pi*fs;n=-800:800;nTs=n/fs;fa=0.5*sin(2*pi*15*nTs)+2*sin(2*pi*40*nTs);Dt=1/fs;t1=-0.1:Dt:0.1;fa1=fa/fs*Wc/pi*sinc(Wc/pi)*(ones(length(nTs),1)*t1-nTsones(1,length(t1);figure(1);plot(t1,fa1);axis(-0.10.1-88);xlabel(fs=100Hz,fa);由抽样定理可知
16、,抽样后的信号频谱是原信号频谱以抽样频率为周期进行周期延拓形成的,周期性在上面两个图中都有很好的体现。但是从10点和50点采样后的结果以及与员连续信号频谱对比可以看出,10点对应的频谱出现了频谱混叠而并非原信号频谱的周期延拓。这是因为N取值过小导致采样角频率,因此经周期延拓出现了频谱混叠。而N取50时,其采样角频率,从而可以实现原信号频谱以抽样频率为周期进行周期延拓,并不产生混叠,从而为下一步通过低通滤波器滤出其中的一个周期(即不失真的原连续信号)打下了基础。(2)加窗、截断(1)、试验代码:n=128;%每个周期抽样的点数N=10;%取10个周期的波形Ts=20/128000;%取样间隔fs
17、=1/Ts;T=n*Ts*N;t=0:Ts:T;fre=50;y=5+220*sqrt(2)*sin(2*pi*fre*t+pi/6)-22*sqrt(2)*sin(2*pi*1.12*fre*t+pi/2)%构造的含有间谐波的波形%*没有加窗函数时分析的构造波形的间谐波*figure(1)subplot(211)plot(t(1:128*N),y(1:128*N)xlabel(t单位:s);ylabel(y单位:V);title(没有加窗函数的时域波形和频率的幅值波形)gridholdonffdat1=fft(y(1:n*N),N*n)/(N*n/2);%对N个周期的波形进行傅立叶变换ftd
18、ata1=abs(ffdat1);ftdata1=ftdata1(1)/2ftdata1(2:end)/sqrt(2);%傅立叶变换后的幅值ftangle1=180*angle(ffdat1)/pi;%傅立叶变换后的相位角tab1=0:1/T:N*n*1/T-1;ftdata1(1:N*n);ftangle1(1:N*n);%统计傅立叶变换%的N*n个点subplot(212)bar(0:fs/(n*N):fs-1),ftdata1(1:N*n);axis(1(fs+1)/250250)xlabel(f(Hz);ylabel(ftdata1单位:V);boxoffgrid%*加hanning窗
19、时分析的构造波形的间谐波*figure(2)subplot(211)plot(t(1:128*N),hanning(n*N).*y(1:128*N)xlabel(t单位:s);ylabel(y单位:V);title(加hanning窗后的时域波形和频率的幅值波形)gridholdonffdat2=fft(hanning(n*N).*y(1:n*N),N*n)/(N*n/2);%对N个周期的波形进行傅立叶变换ftdata2=abs(ffdat2);ftdata2=ftdata2(1)/2ftdata2(2:end)/sqrt(2);%傅立叶变换后的幅值ftangle2=180*angle(ffd
20、at2)/pi;%傅立叶变换后的相位角tab2=0:1/T:N*n*1/T-1;ftdata2(1:N*n);ftangle2(1:N*n);%统计傅立叶变换的N*n个点subplot(212)bar(0:fs/(n*N):fs-1),ftdata2(1:N*n);axis(1(fs+1)/250250)xlabel(f(Hz);ylabel(ftdata2单位:V);boxoffgrid(2)Matlab运行仿真图像:(3)栅栏效应利用500Hz方波信号的频谱分析说明栅栏效应所造成的频谱误差。设定采样频率:fs=5120,MATLAB中默认的FFT计算点数为512,其离散频率点为:fi=i*
21、fs/N=i*5120/512=10,i=0,1,2,3,N/2。位于505H1位置的真实谱峰被拦住看不见,望见的只是它们在相邻频率500Hz或510Hz处能量泄漏的值。若设fs=2560Hz,则频率间隔为df=5Hz,重复上述分析步骤,这时在505Hz位置有谱线,因此可以得到它们的精确值。从时域看,这个条件相当于对信号进行整周期采样,实际中常用此方法来提高周期信号的频谱分析精度。程序如下:t1=0:1/5120:0.1;m1=0:10:5120;t2=0:1/2560:0.2;m2=0:5:2560;y1=square(2*pi*505*t1);transf1=abs(fft(y1)/256
22、);y2=square(2*pi*505*t2);transf2=abs(fft(y2)/256);subplot(2,2,1);plot(t1,y1);title(square(2*pi*505*t1);subplot(2,2,3);stem(m1(1:256),transf1(1:256);title(fs=5120);subplot(2,2,2);plot(t2,y2);title(square(2*pi*505*t2);subplot(2,2,4);stem(m2(1:256),transf2(1:256);title(fs=2560);(2)Matlab运行仿真图像:将上图中的图3和
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 机械工程 测试 技术 实验 报告
限制150内