传感器试验程序MATLAB

霍尔传感器特性试验数据处理程序

x=[Ov1：0.1:1]；%位移

v11=[1.42,2o85,4o32,5.71,7.10,8。52,9。91,11。06,11。19/1。19]：*揄出电压(1次正行程)

v12=-[—1.41,—2.84,—4O44,-5.91,—7.42,-8.81,-10.06,—10o13,一10。14,一10.14]；%输出电压(1次正行程)

v21=[1.43,2o85,4o27,5.63,7.02,8.40,9。73,11。09,11.20,11。20]；*输出电压(2次正行程)

v22=-[—1.44,—2.89,-4。42,—5.88,—7.37,—8。80.-10.09,-10.13,-10.14,-10o14]；%输出电压(2次反行程)

v31=[1.40,2.81,4o27,5。65,7。02,8。46,9.76,11.15,11。20,11。20];%输出电压(3次正行程)

v32=—[―1o33,—2。76,—4.18,-5.83,-7。02,—8。35,—9。84,—10。12,-10.14,—10.14]比输出电压(3次反行

程)

部％结论：k1为一次行程拟合的曲线的斜率(第一个值)和截距(第二个值),

%%%Imd1为1次行在灵锹度，fxxwul为1次行程非线性误差，uuwul为1次行程迟滞误差，ufx为这个实验的重复性误差

$%%其余的k2,Imd2ooo类推。。。

k1=polyfit(x,v11,1)%拟合直线的斜率和极距(1次行程)

y=polyval(k1,x)；

hold

figure(1)

plot(x,v11,',x,y)

gridon

xlabel('位移值(mm)')；

ylabel输出电压(mV),);

title。霍尔传感器输出特性图(1次正行程)')：

Imd1=k1(1)%灵敏度

y1=v11-y;

I1=max(y1);%最大非线性绝对误差

mlc=max(v11)—min(v11);

fxxwc1=l1/mlc*100$非线性误差(算出来的值写的时候记得带正负号和百分号)

figure(2)

plot(x,v11/X,v12/-o')%迟滞特性图

hold

gridon

xlabel('位移值(mm),);

ylabel('输出电压(mV)');

title('迟滞特性图(1次行程)，)；

y2=v11—v12：

I11=max(y2)然正反行程最大差值

ccwc1=l11/(2*(max([v11,v12])—min([v11,v12])))*100%迟滞误差(算出来的值写的时候记得带百分号)

k2=polyfit(x,v21,1)%拟合直线的斜率和截距(2次行程)

y=polyval(k2,x);

figure(3)

plot(x,v21,'—,x,y)

gridon

xlabel('位移值(mm),)；

ylabel('输出电压(mV),);

title('霍尔传感器输出特性图(2次正行程)')；

Imd2=k1(1)%灵敏度

y1=v11-y；

I1=max(y1);%最大非线性绝对误差

mIc=max(v21)—min(v21);

fxxwc2=11/mlc大100%非线性误差(算出来的值写的时候记得带正负号和百分号)

figure(4)

plot(x,v21,'-*',x,v22,'-o')%迟滞特性图

hold

gridon

xlabel('位移值(mm)')；

ylabel('揄出电压(mV)');

title('迟滞特性图(2次行程)’)；

y2=v21-v22；

I11=max(y2)器正反行程最大差值

ccwc2=l11/(2*(max([v21,v22])-min([v21,v22])))*100%迟滞误差(算出来的值写的时候记得带百分号)

k3=polyfit(x,v31,1)%拟合直线的斜率和截距(1次行程)

y=polyval(k3,x)：

figure(5)

plot(x,v31,，一*',x,y)

gridon

xIabeI('位移值(mm)');

ylabel('输出电压(mV)');

title(・霍尔传感器输出特性图(3次正行程)');

Imd3=k1(1)%灵敏度

y1=v11—y;

I1=max(y1)：%最大非线性绝对误差

mlc=max(v31)-min(v31);

fxxwc3=l1/mlc*100%非线性误差(算出来的值写的时候记得带正负号和百分号)

figure(6)

plot(x,v31/—*,,x,v32,,—o')%迟滞特性图

hold

gridon

xlabel('位移值(mm)')：

ylabel('输出电压(mV),);

title。迟滞特性图(3次行程)');

y2=v31-v32；

11kmax(y2)：%正反行程最大差值

ccwc3=l11/(2*(max([v31,v32])-min([v31,v32])))*100%迟滞误差(算出来的值写的时候记得带百分号)

cfx=max(max([abs(v11—v21),abs(v11—v31),abs(v21-v31)]),irax([abs(v12-v22),abs(v12—v32),abs(v22一

v32)]))/(2*(max([v11,v12])—min([v11,v12])))*100与重复性：算出来的值写的时候记得带百分号)

反射式光纤位移传感器特性试验数据处理程序

x=[0o5:0.5:5]；%位移

v11=[0.48,0o81,0.96,0.99,0。94,0。86,0。77,0.68。6,0。53];%输出电压（1次正行程）

v12=[0o41,0.48,0.78,0。96,0.99,0。95,0。86,0.77,0.68,0。6]：%输出电压（1次正行程）

v21=[Oo53,Oo84,0.98,0.99,0。93,0。85,0。75,0。67,0。59,0。52]；%输出电压(2次正行程)

v22=[Oo07,0.51,0.83,0.97,0.99,0.93,0。85,0.75,0.67,0。59]；%输出电压(2次反行程)

v31=[0.4,0。75,Ov95,1,0.96,0.88,0„78,0v69,0v61,0.53]；%揄出电压(3次正行程)

v32=[0.07,0.39,0.74,0。94,1,0。96,0。88,0。79,0.69,0.61];%输出电压(3次反行程)

嬲%结论：k1q为一次行程拟合的前坡曲线的斜率(第一个值)和截距(第二个值)，

$%%k1h为一次行程拟合的后坡曲线的斜率(第一个值)和截距(第二个值)

%%%Imdlq为1次行程前坡灵敏度，fxxwclq为1次行程前坡非线性误差，ccwcl为1次行程迟滞误差，cfx为这个实脸的重

复性误差

%%%Imdlh为1次行程后坡灵敏度，fxxwclh为1次行程后坡非线性误差，其余的k2,Imd2。。。类推.

k1=polyfit(x(1:5),v11(1:5),1)舟拟合直线的斜率和截距(1次行程前坡)

k2=pulyriL(x(6：10),v11(6:10),1)%拟合直线的斜率和他距(1次行程后坡)

yq=polyvaI(k1,x(1:5)):

yh=poIyvaI(k2,x(6:10))；

hold

figure(1)

plot(x,v11,x(1:5),yq,x(6:10),yh,'black')

gridon

xIabeI('位移值(mm)');

ylabel('输出电压(mV),)：

title('反射式光纤位移传感器输出特性图(1次正行程)‘)；

Imd1q=k1(1)%灵敏度(前坡)

Imd1h=k2(1)%灵敏度(后坡)

y1q=v11(1:5)—yq；

y1h=v11(6:10)-yh;

11q=max(y1q);%最大非线性绝对误差

I1h=max(y1h)方最大非线性绝对误差

mlcq=max(v11(1:5))―min(v11(1:5));

mlch=max(v11(6:10))-min(v11(6:10))：

fxxvvc1q=l1q/mlcq*100%非线性误差(算出来的值写的时候记得带正负号和百分号)

fxxwc1h=l1h/mlch*100*非线性误差(算出来的值写的时候记得带正负号和百分号)

figure(2)

plot(x,v11/v12,'-o')舟迟滞特性图

hold

gridon

xIabeI('位移值(mm)');

ylabel输出电，压(mV)')：

titled迟滞特性图(1次行程)0;

y2=v11-v12；

I11=max(y2)：%正反行程最大差值

ccwc1=l11/(2*(max([v11,v12])—min([v11,v12])))*100%迟滞误差(异出来的值写的时候记得带百分号)

k3=polyfit(x(1:5),v21(1:5),1)%拟合直线的斜率和截距(2次行程前坡)

k4=polyfit(x(6：10),v21(6：10),1)%拟合直线的斜率和截距(2次行程后坡)

yq=polyval(k3,x(1:5));

yh=polyvaI(k4,x(6:10));

hold

figure(3)

pIot(x,v21,'—'A-',x(1：5),yq,x(6:10),yh,'bIack*)

gridon

xlabel(位移值(mm),);

ylabel('输出电压(mV),)；

title('反射式光纤位移传感器输出特性图(2次正行程)')；

Imd2q=k3(1)%灵敏度(前坡)

Imd2h=k4(1)$灵敏度(后坡)

y1q=v21(1:5)一yq；

y1h=v21(6:10)—yh；

I1q=max(y1q)；%最大酢线性绝对误差

I1h=max(y1h)：*最大非线性绝对误差

mlcq=max(v21(1:5))―min(v21(1:5));

mlch=max(v21(6:10))-min(v21(6:10))：

fxxwc2q=l1q/mlcq*100%非线性误差(算出来的值写的时候记得带正负号和百分号)

fxxwc2h=l1h/mlch*100%非线性误差(算出来的值写的时候记得带正负号和百分号)

figure(4)

plot(x,v21,x,v22,o')%迟滞特性图

hold

gridon

xIabeI('位移值(mm)');

ylabel('输出电压(mV),)：

title。迟滞特性图(2次行程)’)；

y2=v21-v22；

I11=max(y2)：%正反行程最大差值

ccwc2=H1/(2*(max([v21,v22])—min([v21,v22])))*100%迟滞误差(算出来的值写的时候记得带百分号)

k5:polyfit(x(1:5),v31(1：5),1)%拟合直线的斜率和截距(3次行程前坡)

k6=polyfit(x(6：10),v31(6:10),1)%拟合直线的斜率和截距(3次行程后坡)

yq=polyval(k5,x(1:5));

yh=poIyvaI(k6,x(6：10))；

hold

figure(5)

plot(x,v31,'—5),yq,x(6:10),yh,'black')

gridon

xlabel位移值(mm)')；

ylabelV揄出电压(mV),)：

title。反射式光纤位移传感器输出特性图(3次正行程)’)；

Imd3q=k5(1)%灵敏度(前坡)

Imd3h=k6(1)%灵敏度(后坡)

y1q=v31(1:5)-yq；

y1h=v31(6:10)-yh;

I1q=max(y1q);*最大非线性绝对误差

I1h=max(y1h);%最大非线性绝对误差

mlcq=max(v31(1:5))-min(v31(1：5));

mIch=max（v31（6：10））~min（v31（6:10））;

fxxwc3q=l1q/mlcq*100%非线性误差（算出来的值写的时候记得带正负号和百分号）

fxxwc3h=lIh/mlch^lOO%非线性误差（算出来的值写的时候记得带正负号和百分号）

figure(6)

plot(x,v31,*',x,v32/-o')%迟滞特性图

hold

gridon

xlabel('位移值(mm),);

ylabel('输出电压(mV)')：

title(,迟滞特性图(3次行程)');

y2=v31—v32;

111二max(y2)；%正反行程最大差值

ccwc3=l11/(2*(max([v31,v32])—min([v31,v32])))*100%迟滞误差(算出来的值写的时候记得带百分号)

cfx=max(max([abs(v11-v21),abs(v11—v31),abs(v21-v31)]),max([abs(v12—v22),abs(v12-v32),abs(v22—v32)]))/

(2*(max([v11,v12])-min([v11,v12])))*100%重复性(算出来的值写的时候记得带百分号)

电涡流传感器特性与位移测量试验数据处理程序

x=[0o25：0.25：1O25];先位移

v11=-[1O98,2.81,3o55,4.28,4。96]：%输出电压（1次正行程）

v12=—[0„01,1.95,2。81,3.54,4。26];%榆出电压（1次正行程）

v21=-[0o01,1O97,2O8,3.56,4O3];先输出电压（2次正行程）

v22=—[0,08,0.01,1o96,2.79,3.55];%输出电压（2次反行程）

v31=—[2,09,2.92,3.69,4。39,5.08]；%输出电压（3次正行程）

v32=-[0o2,2o05,2.9,3。66,4.37]；%输出电压（3次反行程）

%嬲结论:k1为一次行程拟合的曲线的斜率（第一个值）和截距（第二个值），

%%%Imd1为1次行程灵敏度，fxxwcl为1次行程非线性误差，ccwc1为1次行程迟滞误差，cfx为这个实脸的重复性误差

题疑其余的k2,类推。.（＞

k1=polyfit（x,v11,1）%拟合直线的斜率和截距（1次行程）

y=polyval(k1,x)；

hold

figure(1)

plot(x,v11,',x,y)

gridon

xIabeI('位移值(mm)')；

yIabeI(,输出电压(mV),);

title('电涡流传感器输出特性图(1次正行程)’);

Imd1=k1(1)%灵敏度

y1=v11—y;

I1=max(y1);%最大非线性绝对误差

mIu=max(v11)—rniri(v11);

fxxwc1=l1/mlc*100%非线性误差(算出来的值写的时候记得带正负号和百分号)

figure(2)

plot(x,x,v12/—o')%迟滞特性图

hold

gridon

xlabel('位移值(mm),);

ylabel('输出电压(mV),)：

title。迟滞特性图(1次行程)’)；

y2=v11-v12;

11kmax(y2)；%正反行程最大差值

ccwc1=H1/(2*(max([v11,v12])-min([v11,v12])))*10。％迟滞误差(算出来的值写的时候记得带百分号)

k2=polyfit(x,v21,1)%拟合直线的斜率和截距(2次行程)

y=polyval(k2,x);

figure(3)

plot(x,v21,'—*',x,y)

gridon

xlabel('位移值(mm)')；

yIabeI('输出电压(mV),);

title('电涡流传感器输出特性图(2次正行程)’)；

Imd2=k1(1)%灵敏度

y1=v11-y;

I1=max(y1);%最大非线性绝对误差

mlc=max(v21)—min(v21);

fxxwc2=l1/mlc*100%非线性误差(算出耒的值写的时候记得带正负号和百分号)

figure(4)

plot(x,v21,'—x,v22,'-o')%迟滞特性图

hold

gridon

xlabel('位移值(mm),);

ylabel('输出电压(mV),);

title('迟滞特性图(2次行程)’)：

y2=v21—v22;

I11=max(y2)；舟正反行程最大差值

ccwc2=H1/(2*(max([v21,v22])—min(