误差理论实验报告3.docx

误差理论实验报告3《误差理论与数据处理》实验报告实验名称:动态测试数据处理初步实验目的动态数据是动态测试研究的重要内容。通过本实验要求学生掌握有关动态数据分析。评价的基本方法,为后续课程做好准备。实验原理实验内容和结果程序及流程认识确定性信号及其傅立叶频谱之间的关系用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',[0011]);>>set(gca,'xtick',0::8);>>axis([0,8,0300]);用matlab边城画出矩形窗信号的宽度分别为T=1和T=5两种情况下的时域波形图及其频谱,并分析时域与频域的变化关系。wlp=*pi; whp=*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([]); title('矩形窗函数(T=1)'); subplot(413);  stem(window2);axis([]); 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([01-3500]); grid; xlabel('w/pi'); ylabel('幅度(dB)');  title('矩形窗函数的频谱(T=5)'); (t)的相关函数为:Rx(t)=e-ατcosω0τ,画出下列两种情况下的自相关函数和功率谱函数。取α=1,ω0=2π*10;取α=5,ω0=2π*10;程序:>>t=0::1;y1=.^(-t).*cos(20.*pi.*t);subplot(221)plot(t,y1);title('(1)自相关函数')holdony2=.^(-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,