误差理论实验报告3.docx

《误差理论与数据处理》实验报告 实验名称：动态测试数据处理初步 实验目的 动态数据是动态测试研究的重要内容。通过本实验要求学生掌握有关动态数据分析。评价的基本方法，为后续课程做好准备。 实验原理 实验内容和结果 程序及流程 认识确定性信号及其傅立叶频谱之间的关系 用matlab编程画出周期方波信号及其傅立叶频谱，并说明其傅立叶频谱的特点。 fs=30; T=1/fs; t=0:T:2*pi; A=2;P=4; y=A*square(P*t); subplot(2,1,1),plot(t,y) title(方波信号) Fy=abs(fft(y,512)); f2=fs*(0:256)/512; subplot(2,1,2),plot(f2,Fy(1:257)) title(频谱图); set(gcf,unit,normalized,position,[0 0 1 1]); set(gca,xtick,0:0.6:8); axis([0,8,0 300]); 用matlab边城画出矩形窗信号的宽度分别为T=1和T=5两种情况下的时域波形图及其频谱，并分析时域与频域的变化关系。 wlp = 0.35*pi;?whp = 0.65*pi;?wc = [wlp/pi,whp/pi];window1= boxcar(1);? window2=boxcar(5);[h1,w]=freqz(window1,1);? [h2,w]=freqz(window2,5);subplot(411);??stem(window1); axis([0 60 0 1.2]);? title(矩形窗函数(T=1));? subplot(413);??stem(window2); axis([0 60 0 1.2]);? grid;?xlabel(n);??title(矩形窗函数(T=5));?subplot(412);????plot(w/pi,20*log(abs(h1)/abs(h1(1))));? xlabel(w/pi);?ylabel(幅度(dB));??title(矩形窗函数的频谱(T=1));? subplot(414);????plot(w/pi,20*log(abs(h2)/abs(h2(5))));? axis([0 1 -350 0]);?grid;?xlabel(w/pi);?ylabel(幅度(dB));??title(矩形窗函数的频谱(T=5));? 2.认识平稳随机过程自相关函数及其功率谱之间的关系 已知某随机过程x(t)的相关函数为：Rx(t)=e-α 取α=1， 取α=5， 程序： t=0:0.01:1; y1=1.71828.^(-t).*cos(20.*pi.*t); subplot(221) plot(t,y1); title((1)自相关函数) hold on y2=1.71828.^(-5*t).*cos(20.*pi.*t); subplot(222) plot(t,y2); title((2)自相关函数) subplot(2,2,3); pwelch(y1,33,32,[],500); title((1)概率密度函数) subplot(2,2,4); pwelch(y2,33,32,[],500); title((2)概率密度函数) 3.求随机过程的均值、方差和自相关函数 在线纹比长仪上对0-1000mm线纹尺测量六次，所的各段长度对公称值偏差?如下表（个尺寸段单位：mm，表中偏差值单位：um）： 序号 尺寸段 0-100 0-200 0-300 0-400 0-500 0-600 0-700 0-800 0-900 0-1000 1 2 3 4 5 6 0.18 0.30 0.30 0.25 0.30 0.33 0.34 0.38 0.42 0.34 0.38 0.44 0.63 0.70 0.76 0.69 0.73 0.76 1.20 1.26 1.22 1.22 1.30 1.28 1.51 1.55 1.52 1.54 1.58 1.60 2.02 2.10 2.01 1.96 2.03 2.08 2.22 2.26 2.16 2.22 2.28 2.31 2.62 2.66 2.69 2.72 2.71 2.78 2.54 2.56 2.60 2.64 2.69 2.70 2.64 2.66 2.67 2.66 2.71 2.81 编程画出6此实验曲线（散点图或折线图）； X1=[0.18 0.34 0.63 1.20 1.51 2.02 2.22 2.62 2.54 2.64]; X2=[0.30 0.38