误差理论实验报告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::8);
>> axis([0,8,0 300]);
用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([0 60 0 ]); 
title('矩形窗函数(T=1)'); 
subplot(413);  
stem(window2);
axis([0 60 0 ]); 
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)'); 

已知某随机过程x(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)自相关函数'