传感器标定程序使用说明书.doc

传感器标定程序使用说明： 建立实验数据文件 在附件的写字板中建立名称为： XXin.txt Yz1input.txt Yz2input.txt Yz3input.txt Yf1input.txt Yf3input.txt Yf3input.txt，并保存于:“d:\传感器实验\”|子目录下。 其中 XXin.txt 是输入量的文件 Yz1input.txt 是输出量的文件，表示第一轮正行程 Yz2input.txt 是输出量的文件，表示第二轮正行程 Yz3input.txt 是输出量的文件，表示第三轮正行程 Yf1input.txt 是输出量的文件，表示第一轮反行程 Yf3input.txt 是输出量的文件，表示第二轮反行程 Yf3input.txt 是输出量的文件，表示第三轮反行程 建立方法：见图，所有输入和输出的实验数据必须按列排放，一行一个数，回车。注意所有的文件长度相等，也即，x和y的个数必须相等。否则计算程序无法执行。 在matlab命令窗口执行 sensor回车即可 出现以下窗口： 使用十字光标在图上做标注，当停止时，回车，往下执行，得到如下图结果： 最后一次按空格键，回到命令窗口看计算结果：包括拟和方程的斜率和截距、线性度、迟滞、重复性和总精度等 程序清单： %读入实验数据 clear load XXin.txt; load Yz1input.txt; load Yz2input.txt; load Yz3input.txt; load Yf1input.txt; load Yf2input.txt; load Yf3input.txt; %变量赋值 X =XXin ; Yz1= Yz1input; Yz2= Yz2input; Yz3=Yz3input; Yf1=Yf1input; Yf2=Yf2input; Yf3=Yf3input; %显示第一轮实验结果 plot(X,Yz1,g,X,Yz1,g*,X,Yf1,g,X,Yf1,go);hold on title(传感器静态特性实验结果) xlabel(x /mm ,/g) ylabel(y/mv); grid gtext(*——正行程数据); gtext(o——反行程数据); gtext(绿色——第一轮实验数据); pause;%停顿，按任意键继续 %显示第二轮实验结果 plot(X,Yz2,b,X,Yz2,b*,X,Yf2,b,X,Yf2,bo);hold on gtext(兰色——第二轮实验数据); pause; %显示第三轮实验结果 plot(X,Yz3,r,X,Yz3,r*,X,Yf3,r,X,Yf3,ro);hold on gtext(红色——第三轮实验数据); pause; %求正反行程三次实验平均值 Ypjz=(Yz1+Yz2+Yz3)/3; Ypjf=(Yf1+Yf2+Yf3)/3; Ypj=(Ypjz+Ypjf)/2; plot(X,Ypj,k,X,Ypj,k+);hold on gtext(+黑色——三轮实验数据平均); pause; %以实验平均值拟和曲线 [k]=polyfit(X,Ypj,1); k %显示拟和曲线的斜率和截距 Yjs=k(1)*X+k(2) %显示拟和曲线 plot(X,Yjs,b,X,Yjs,b^);hold on gtext(三角——拟和曲线计算数据); pause %Yfs是满量程输出 Yfs=max(Yjs); %求实验数据与计算数据的差值 Ycz1=Yz1-Yjs; Ycf1=Yf1-Yjs; Ycz2=Yz2-Yjs; Ycf2=Yf2-Yjs; Ycz3=Yz3-Yjs; Ycf3=Yf3-Yjs; %求各行程的最大差值 mz1=max(abs(Ycz1)); mf1=max(abs(Ycf1)); mz2=max(abs(Ycz2)); mf2=max(abs(Ycf2)); mz3=max(abs(Ycz3)); mf3=max(abs(Ycf3)); l=[mz1 mf1 mz2 mf2 mz3 mf3]; %求线性度 Lmax=max(l); xxdL=Lmax/Yfs*100%; %求迟滞 hmax=max(abs(Ypjz-Ypjf)); chzhH=hmax/Yfs*100%; %求重复性 Yzp1=Ycz1.^2; n=3; Yzp1=sqrt(sum(Yzp1)/(n-1)); Yzp2=Ycz2.^2; Yzp2=sqrt(sum(Yzp2)/(n-1)); Yzp3=Ycz3.^2; Yzp3=sqrt(sum(Yzp3)/(n-1)); Yfp1=Ycf1.^2; Yfp1=sqrt(sum(Yfp1)/(n-1)); Yfp2=Ycf2.^2