误差理论实验报告3.doc

: .
《误差理论与数据处理》实验报告
实验名称：动态测试数据处理初步
一、 实验目的
动态数据是动态测试研究的重要容。通过木实验要求学生掌握 有关动态数据分析。评价的基木方法，为后续课程做好准备。
二、 实验原理
三、 实验容和结果
1. 程序及流程
1. 认识确定性信号及其傅立叶频谱之间的关系
1. 用mat lab编程画岀周期方波信号及其傅立叶频谱，并说明其 傅立叶频谱的特点。
>> fs=30;
>> T= 1/fs;
>> t=0:T:2*pi;
>> A=2;P=4;
>> y=A*square(P*t);
>> subplot(2z 1,l),plot(t,y)
» title(方波信号')
>> Fy=abs(fft(y,512));
>> f2=fs*(0:256)/512;
>> subplot(2,l,2),plot(f2,Fy(l:257))
» title(濒谱图')；
>> setCgcf/unif/normalized'/positlon^O 0 1 1]);
>> set(gca/xtick',0::8);
>> axis([0,8z0 300]);
lab边城画出矩形窗信号的宽度分别为T=1和T=5两种 情况下的时域波形图及其频谱，并分析时域与频域的变化关系。
wlp = *pi;
whp = *pi;
wc = [wlp/pi,whp/pi];
windowl = boxcar(1);
window2=boxcar(5);
[hi ,w]=freqz(window1,1);
[h2,w]=freqz(window2,5);
subplot(411);
stem(windowl);
axis([0 60 0 ]);
titled矩形窗函数(T=1»；
subplot(413);
stem(window2);
axis([0 60 0 ]);
grid;
xlabelfn1);
titled矩形窗函数(T=5»； subplot(412);
plot(w/pi,20*log(abs(h 1 )/abs(h 1 (1))));
xlabelCw/pi')；
ylabel(幅度(dB));
title。矩形窗函数的频i§(T=1D;
subplot(414);
plot(w/pi,20*log(abs(h2)/abs(h2(5))));
axis([0 1 -350 0]);
grid;
xlabelfw/pi1)；
ylabel(幅度(dB 门；
title(矩形窗函数的频谱(T=5));
2. 认识平稳随机过程自相关函数及其功率谱之间的关系
已知某随机过程x(t)的相关函数为：Rx(t)=e-a|x|cosu)oT,画岀下 列两种情况下的自相关函数和功率谱函数。
1. 取a =1, coo = 2tc♦ 10；
2. 取a = 5, a)o = 2tt♦ 10；
程序:>> t=0:：l;
yl=