搜索
误差理论与数据处理姓名:****班级：研硕控制8班学号：*****日期：2015年12月2日学校:矿大第1页共8页误差理论与数据处理概述：现代科学与技术工作离不开实验，而实验又离不开对具体量值的测量，并且测量总是存在误差的。人们总是希望减少或消除测量中存在的误差来提高测量精度。误差理论正是研究误差的来源、性质、规律，减小误差的措施，评估测量结果精度的办法，如何合理的设计实验去减小甚至消除误差，以及提高实验或测量精度的理论。数据处理泛指对数据，特别是实验观测数据进行的一切加工及运算，包括由测量数据求待测值，估其测量精度，列表或作图表示数据之间的对应关系与数据变化规律，以及通过所得到的数据求取相关未知量、求取相关量之间的未知关系式等等。一、数据处理的几种滤波算法1、高斯滤波高斯函数具有五个重要的性质，这些性质使得它在早期图像处理中特别有用．这些性质表明，高斯平滑滤波器无论在空间域还是在频率域都是十分有效的低通滤波器，且在实际图像处理中得到了工程人员的有效使用。（1）二维高斯函数具有旋转对称性，即滤波器在各个方向上的平滑程度是相同的．一般来说，一幅图像的边缘方向是事先不知道的，因此，在滤波前是无法确定一个方向上比另一方向上需要更多的平滑．旋转对称性意味着高斯平滑滤波器在后续边缘检测中不会偏向任一方向．（2）高斯函数是单值函数．这表明，高斯滤波器用像素邻域的加权均值来代替该点的像素值，而每一邻域像素点权值是随该点与中心点的距离单调增减的．这一性质是很重要的，因为边缘是一种图像局部特征，如果平滑运算对离算子中心很远的像素点仍然有很大作用，则平滑运算会使图像失真．（3）高斯函数的付立叶变换频谱是单瓣的．正如下面所示，这一性质是高斯函数付立叶变换等于高斯函数本身这一事实的直接推论．图像常被不希望的高频信号所污染(噪声和细纹理)．而所希望的图像特征（如边缘），既含有低频分量，又含有高频分量．高斯函数付立叶变换的单瓣意味着平滑图像不会被不需要的高频信号所污染，同时保留了大部分所需信号．第2页共8页（4）高斯滤波器宽度(决定着平滑程度)是由参数σ表征的，而且σ和平滑程度的关系是非常简单的．σ越大，高斯滤波器的频带就越宽，平滑程度就越好．通过调节平滑程度参数σ，可在图像特征过分模糊(过平滑)与平滑图像中由于噪声和细纹理所引起的过多的不希望突变量(欠平滑)之间取得折衷．（5）由于高斯函数的可分离性，大高斯滤波器可以得以有效地实现．二维高斯函数卷积可以分两步来进行，首先将图像与一维高斯函数进行卷积，然后将卷积结果与方向垂直的相同一维高斯函数卷积．因此，二维高斯滤波的计算量随滤波模板宽度成线性增长而不是成平方增长．二维高斯函数：2222222),(ryxAeAeyxG当时，；时，一般用宽度小于的滤波器，即当时，由连续Gaussian分布求离散模板，需采样、量化，并使模板归一化，举例结果如下：图1原图像lenna图2有噪声的lenna第3页共8页图3高斯滤波，σ2=1图4高斯滤波，σ2=3可以看到高斯滤波虽然能够在一定程度上去掉噪声，但也使得图象变得模糊不清，效果并不能令人满意。Matlab程序如下：%%%%%%%%%%%%%Themain.mfile%%%%%%%%%%%%%%%clc;%ParametersoftheGaussianfilter:n1=5;sigma1=3;n2=5;sigma2=3;theta1=0;[w,map]=imread('lenna.gif');x=ind2gray(w,map);filter1=d2gauss(n1,sigma1,n2,sigma2,theta1);y=imnoise(x,'gauss',0.01);f1=conv2(x,filter1,'same');rf1=conv2(y,filter1,'same');figure(1);subplot(2,2,1);imagesc(x);title('lenna');subplot(2,2,2);imagesc(y);title('noisylenna');subplot(2,2,3);imagesc(f1);title('smooth');subplot(2,2,4);imagesc(rf1);title('noisecancel');colormap(gray);%%%%%%%%%%%%%%Endofthemain.mfile%%%%%%%%%%%%%%%%Functiond2gauss.m:%Thisfunctionreturnsa2DGaussianfilterwithsizen1*n2;thetais第4页共8页%theanglethatthefilterrotatedcounterclockwise;andsigma1andsigma2%arethestandarddeviationoftheGaussianfunctions.functionh=d2gauss(n1,std1,n2,std2,theta)r=[cos(theta)-sin(theta);sin(theta)cos(theta)];fori=1:n2forj=1:n1u=r*[j-(n1+1)/2i-(n2+1)/2]';h(i,j)=gauss(u(1),std1)*gauss(u(2),std2);endendh=h/sqrt(sum(sum(h.*h)));%Functiongauss.m:functiony=gauss(x,std)y=exp(-x^2/(2*std^2))/(std*sqrt(2*pi));%%%%%%%%%%%%%%Endofthefunctions%%%%%%%%%%%%%%%%2、最小二乘法滤波算法已知在不同的温度T下，测定铜棒的长度l如表1-1所示。表1-1i12345678iT/℃1015202530354045il/cm2000.362000.502000.722000.802001.072001.252001.482001.60Matlab利用原始数据画折线图如图5，程序如下：clc,clear;T=[1015202530354045];L=[2000.362000.502000.722000.802001.072001.252001.482001.60];plot(T,L,'m');第5页共8页gridon;xlabel('T/℃');ylabel('L/cm');title('T-LLinechart');legend('T-L');图5T-LLineChart由折线图可知，铜棒的长度l随温度T呈线性变化,设laTb=+，用最小二乘法给出参数a和b的最小二乘估计值。Matlab实现最小二乘参数估计LN=[2000.362000.502000.722000.802001.072001.252001.482001.60]';TN=[10,1;15,1;20,1;25,1;30,1;35,1;40,1;45,1]；ab=inv(TN'*TN)*TN'*LN；x=10:1:45；plot(x,y,'b',T,L,'m');gridon;xlabel('T/℃');ylabel('L/cm');title('T-LLinechart');legend('L=aT+b','T-L');a=ab(1)a=第6页共8页0.0368b=ab(2)b=2.0000e+003原始数据折线图与l=aT+b函数图形(图6)对比：图6折线图与直线图对比所以铜棒的长度l与温度T的线性关系式为：l=0.0368T+2000。3、双边滤波算法双边滤波（Bilateralfilter）是一种可以保边去噪的滤波器。之所以可以达到此去噪效果，是因为滤波器是由两个函数构成。一个函数是由几何空间距离决定滤波器系数。另一个由像素差值决定滤波器系数。Matlab程序如下:I=imread('einstein.jpg');I=double(I)/255;w=5;sigma=[30.1];%bilateralfilterstandarddeviationsI1=bfilter2(I,w,sigma);subplot(1,2,1);Imshow(I);subplot(1,2,2);第7页共8页Imshow(I1);仿真结果如图7所示。图7二、几种算法的优缺点高斯滤波(Gaussfilter)实质上是一种信号的滤波器，其用途为信号的平滑处理，数字图像用于后期应用，其噪声是最大的问题，因为误差会累计传递等原因，大多图像处理教材会在很早的时候介绍Gauss滤波器，用于得到信噪比SNR较高的图像（反应真实信号）。高斯平滑滤波器对于抑制服从正态分布的噪声非常有效。但同时使图像变得模糊，效果不是太好。最小二乘法（又称最小平方法）是一种数学优化技术。它通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据，并使得这些求得的数据与实际数据之间误差的平方和为最小。最小二乘法还可用于曲线拟合。其他一些优化问题也可通过最小化能量或最大化熵用最小二乘法来表达。但它依赖先验公式，当数据无法用先验公式模拟和数据的数量巨大时，就不适宜使用。双边滤波Bilateralblur相对于传统的高斯blur来说很重要的一个特性即可可以保持边缘（EdgePerseving），这个特点对于一些图像模糊来说很有用。一般的高斯模糊在进行采样时主要考虑了像素间的空间距离关系，但是却并没有考虑像素值之间的相似程度，因此这样我们得到的模糊结果通常是整张图片一团模糊。Bilateralblur的改进就在于在采样时不仅考虑像素在空间距离上的关系，同时加入了像素间的相似程度考虑，因而可以保持原始图像的大体分块进而保持边缘。在于游戏引擎的postblur算法中，bilateralblur常常被用到。第8页共8页参考文献：[1]费业泰.误差理论与数据处理（第4版）[M].北京:机械工业出版社,2003.[2]冈萨雷斯.数字图像处理[M].电子工业出版社,2003.[3]张志涌等，MATLAB与仿真（2010a）,北京航空航天大学出版社,2010（即将出版）.