误差理论实验报告3..docxVIP

《误差理论与数据处理》实验报告实验名称：动态测试数据处理初步实验目的动态数据是动态测试研究的重要内容。通过本实验要求学生掌握有关动态数据分析。评价的基本方法，为后续课程做好准备。实验原理实验内容和结果程序及流程认识确定性信号及其傅立叶频谱之间的关系用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)=，画出下列两种情况下的自相关函数和功率谱函数。取；取；程序： t=0:0.01:1;y1=1.71828.^(-t).*cos(20.*pi.*t);subplot(221)plot(t,y1);title((1)自相关函数)hold ony2=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-1000-2000-3000-4000-5000-6000-7000-8000-9000-10001234560.180.300.300.250.300.330.340.380.420.340.380.440.630.700.760.690.730.761.201.261.221.221.301.281.511.551.521.541.581.602.022.102.011.962.032.082.222.262.162.222.282.312.622.662.692.722.712.782.542.562.602.642.692.702.642.662.672.662.712.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 0.70 1.26 1.55 2.10 2.26 2.66 2.56 2.66];X3=[0.30 0.42 0.67 1.22 1.52 2.01 2.16 2.69 2.60 2.67];X4=[0.25 0.34 0.69 1.22 1.54 1.96 2.22 2.72 2.64 2.66];X5=[0.30