基于matlab的误差数据处理实验报告.docVIP

《误差理论与数据处理》实验美扎巴·测控三班 实验一：MATLAB软件基础 （一）实验目的：熟悉MATLAB软件的用户环境；了解MATLAB软件的一般目的命令；掌握MATLAB教组操作与运算函数；掌握MATLAB软件的基本绘图命令；掌握MATLAB语言的几种循环、条件和开关选择结构。 （二）实验内容： MATLAB软件的数组处理及运算操作 E=eye(3,3); R=rand(3,2); O=zeros(2,3); S=[2,0;0,4]; A=[E,R;O,S] B=[E,R+R*S;O,S^2] C=A*A 所以B=C，原结论成立。 直接使用MATLAB软件进行作图练习 （1）t=-1:0.01:1; x=sin(2*pi*t); y=cos(2*pi*10*t); plot(t,x,t,y) xlabel(t);ylabel(函数值) legend(正弦函数,余弦函数) （2） 1) x=-10:0.01:10; y=normpdf(x,0,1); plot(x,y) （3）[x,y]=meshgrid(-10:0.05:10); z=sin(pi*sqrt(x.^2+y.^2)); mesh(x,y,z); 3.用MATLAB语言编写命令M-文件和函数M-文件 a=input(请输入a的值) x0=a./2 x1=(x0+a./x0)./2 while(abs((x0-x1)1e-5)) x0=x1; x1=(x0+a./x0)./2; end digits(8) vpa(x1) 实验二：测量数据的统计分析 （一）实验目的：通过对测量数据进行统计分析，学习掌握测量数据统计分析的基本方法。 （二）实验内容： 1. x=normrnd(10,5,500,1); mu=mean(x) mu =9.7672 sigma=std(x) sigma =4.8754 va=var(x) va =23.7697 hist(x) y=normpdf(x,mu,sigma); plot(x,y) 2. x=-15:0.01:15; y1=normpdf(x,0,1); y2=normpdf(x,0,4); y3=normpdf(x,10,1); plot(x,y1,y2,y3); 3. x=randn(500,1); mu=mean(x); va=var(x); cs=skewness(x); ck=kurtosis(x); hist(x); sigma=std(x); y=normpdf(x,mu,sigma); plot(x,y) cs cs =0.1117 ck ck =3.0089 mu mu =0.0730 va va=0.9981 4. x=-5:0.1:5; y1=tpdf(x,5); y2=tpdf(x,10); y3=tpdf(x,20); z=normpdf(x,0,1); plot(x,y1,x,y2,x,y3,x,z) 5. x=0:0.2:30; y1=chi2pdf(x,5); y2=chi2pdf(x,10); y3=chi2pdf(x,20); plot(x,y1,x,y2,x,y3) 6.x=0:0.01:10; y1=fpdf(x,4,5); y2=fpdf(x,10,20); y3=fpdf(x,50,50); plot(x,y1,x,y2,x,y3) 实验三：等精度和非等精度直接测量数据处理 （一）实验目的：通过本实验使学生掌握等精度和非等精度直接测量数据的基本处理方法；学习如何发现和处理测量列中的随机误差、系统误差和粗大误差，如何科学地表达测量结果。 （二）实验内容： l=[24.674,24.675,24.673,24.676,24.671,24.678,24.672,24.674]; x1=mean(l); v=l-x1; a=sum(v); ah=abs(a); bh=ah-(8/2)*0.0001; xt=sum(v(1:4))-sum(v(5:8)); bz=sqrt((sum(v.^2)/7)); p=sort(l) g0=2.03; g1=(x1-p(1))/bz; g8=(p(8)-x1)/bz; sc=bz/(sqrt(8)); t=2.36; jx=t*sc l1=x1+jx; l2=x1-jx 结果：p =24.6710 24.6720 24.6730 24.6740 24.6740 24.6750 24.6760 24.6780 jx =0.0019 l2 =24.6723 实验四：分析发现测量数据中的系统误差 （一）实验目的：通过本