箔式应变片桥路——全半桥

各类全半桥对比---可知线性度

2四个拟合曲线加斜率---排除异常值

（1）全桥正拟合数据

（2）全桥反拟合数据

（3）半桥正拟合数据

（4）半桥反拟合数据

总结：

对比计算后，发现全桥为半桥的两倍，非线性大大减小。

附录代码：

%  这里数据也可以直接导入用xlsread这个算子可以导入excel表格

%全桥＋半桥正实验

x=0:20:200;%x轴上的数据，第一个值代表数据开始，第二个值代表间隔，第三个值代表终止

 a1=[0.002,0.047,0.091,0.136,0.181,0.225,0.270,0.315,0.359,0.403,0.448]; %全桥正数据y值

 b1=[0.005,0.013,0.033,0.054,0.074,0.095,0.113,0.135,0.156,0.176,0.198]; %半桥正数据y值

  subplot(3,2,1)

 plot(x,a1,'-*b',x,b1,'-or'); %线性，颜色，标记

axis([0,200,0,0.6])  %确定x轴与y轴框图大小

set(gca,'XTick',[0:20:200]) %x轴范围0-200，间隔20

set(gca,'YTick',[0:0.04:0.6]) %y轴范围0-0.6，间隔0.04

legend('全桥正','半桥正');   %右上角标注

xlabel('x')  %x轴坐标描述

ylabel('y') %y轴坐标描述

% subspace(2,2,1)  %

%全桥＋半桥反实验

x=0:20:200;%x轴上的数据，第一个值代表数据开始，第二个值代表间隔，第三个值代表终止

 a2=[0.004,0.048,0.092,0.137,0.182,0.226,0.270,0.315,0.359,0.404,0.448]; %全桥反数据y值

 b2=[0.002,0.006,0.016,0.038,0.060,0.082,0.104,0.127,0.151,0.174,0.198]; %半桥反数据y值

  subplot(3,2,2)

 plot(x,a2,'-*b',x,b2,'-or'); %线性，颜色，标记

axis([0,200,0,0.6])  %确定x轴与y轴框图大小

set(gca,'XTick',[0:20:200]) %x轴范围0-200，间隔20

set(gca,'YTick',[0:0.04:0.6]) %y轴范围0-0.6，间隔0.04

legend('全桥反','半桥反');   %右上角标注

xlabel('x')  %x轴坐标描述

ylabel('y') %y轴坐标描述

% subspace(2,2,2)

%全桥正反实验

x=0:20:200;%x轴上的数据，第一个值代表数据开始，第二个值代表间隔，第三个值代表终止

 a3=[0.002,0.047,0.091,0.136,0.181,0.225,0.270,0.315,0.359,0.403,0.448]; %全桥正数据y值

 b3=[0.004,0.048,0.092,0.137,0.182,0.226,0.270,0.315,0.359,0.404,0.448]; %全桥反数据y值

  subplot(3,2,3)

 plot(x,a3,'-*b',x,b3,'-or'); %线性，颜色，标记

axis([0,200,0,0.6])  %确定x轴与y轴框图大小

set(gca,'XTick',[0:20:200]) %x轴范围0-200，间隔20

set(gca,'YTick',[0:0.04:0.6]) %y轴范围0-0.6，间隔0.04

legend('全桥正','全桥反');   %右上角标注

xlabel('x')  %x轴坐标描述

ylabel('y') %y轴坐标描述

%半桥正反实验

x=0:20:200;%x轴上的数据，第一个值代表数据开始，第二个值代表间隔，第三个值代表终止

 a4=[0.005,0.013,0.033,0.054,0.074,0.095,0.113,0.135,0.156,0.176,0.198]; %半桥正数据y值

 b4=[0.002,0.006,0.016,0.038,0.060,0.082,0.104,0.127,0.151,0.174,0.198]; %半桥反数据y值

 subplot(3,2,4)

 plot(x,a4,'-*b',x,b4,'-or'); %线性，颜色，标记

axis([0,200,0,0.6])  %确定x轴与y轴框图大小

set(gca,'XTick',[0:20:200]) %x轴范围0-200，间隔20

set(gca,'YTick',[0:0.04:0.6]) %y轴范围0-0.6，间隔0.04

legend('半桥正','半桥反');   %右上角标注

xlabel('x')  %x轴坐标描述

ylabel('y') %y轴坐标描述

%全桥＋半桥正反对比

x=0:20:200;%x轴上的数据，第一个值代表数据开始，第二个值代表间隔，第三个值代表终止

 a5=[0.002,0.047,0.091,0.136,0.181,0.225,0.270,0.315,0.359,0.403,0.448]; %全桥正数据y值

 b5=[0.004,0.048,0.092,0.137,0.182,0.226,0.270,0.315,0.359,0.404,0.448]; %全桥反数据y值

  subplot(3,2,5)

 plot(x,a5,'-*b',x,b5,'--pc'); %线性，颜色，标记

hold on

 a6=[0.005,0.013,0.033,0.054,0.074,0.095,0.113,0.135,0.156,0.176,0.198]; %半桥正数据y值

 b6=[0.002,0.006,0.016,0.038,0.060,0.082,0.104,0.127,0.151,0.174,0.198]; %半桥反数据y值

%   subplot(3,2,1)

 plot(x,a6,'-or',x,b6,'--og');

 hold on

axis([0,200,0,0.6])  %确定x轴与y轴框图大小

set(gca,'XTick',[0:20:200]) %x轴范围0-200，间隔20

set(gca,'YTick',[0:0.04:0.6]) %y轴范围0-0.6，间隔0.04

legend('全桥正','全桥反','半桥正','半桥反');   %右上角标注

xlabel('重量（g）')  %x轴坐标描述

ylabel('电压（V）') %y轴坐标描述

% 计算灵敏度

% Gx=[0,20,40,60,80,100,120,140,160,180,200];

% %存储全桥数据

% Gx1=[0.002,0.047,0.091,0.136,0.181,0.225,0.270,0.315,0.359,0.403,0.448];

% Gx2=[0.004,0.048,0.092,0.137,0.182,0.226,0.270,0.315,0.359,0.404,0.448];

% %存储半桥数据

% Gx3=[0.005,0.013,0.033,0.054,0.074,0.095,0.113,0.135,0.156,0.176,0.198]; %半桥正数据y值

% Gx4=[0.002,0.006,0.016,0.038,0.060,0.082,0.104,0.127,0.151,0.174,0.198]; %半桥反数据y值

%为防止数据有误差，先进行一次拟合算法处理数据 然后得到结果

% [Gx,Gx1]=size(c1);

% [Gx,Gx2]=size(c2);

% [Gx,Gx3]=size(c3);

% [Gx,Gx4]=size(c4);

%以下为算法方法，可以自行带入，data表格需要表格

% [n,m]=size(data);

% %计算“斜率”和“截距”

% k=((n*sum(data(:,1).*data(:,2)))-sum(data(:,2)*sum(data(:,1))))/((n*sum(data(:,1).^2))-(sum(data(:,1))).^2);

% b=((sum(data(:,1).^2)*sum(data(:,2)))-(sum(data(:,1))*sum(data(:,1).*data(:,2))))/((n*sum(data(:,1).^2))-(sum(data(:,1))).^2);

%

% %隐式方程

% y=@(x)k*x+b;

%

% %拟合值

% nhdata(:,1)=data(:,1);

% nhdata(:,2)=data(:,1).*k+b;

%

% %计算  总体平方和 SST

% SST=sum((data(:,2)-mean(data(:,2))).^2);

%

% %计算  误差平方和 SSE

% SSE=sum((data(:,2)-nhdata(:,2)).^2);

%

% %计算  回归平方和 SSR

% SSR=sum((nhdata(:,2)-mean(data(:,2))).^2);

%

% %拟合优度

% R=SSR/SST;

%

%

%

% %画图

% fplot(y,[2,8]);

% hold on;

% xlabel(['拟合优度：',num2str(R)]);

% plot(data(:,1),data(:,2),'o');

% grid on;

% k1=0.446/200;

% k2=0.444/200;

% k3=0.193/200;

% k4=0.196/200;

% eLq=0.002/0.448*100;%全桥

% sprintf("0.4464%");

% eLb=0.003/0.198*100;%半桥

% sprintf("1.5152%");

% eHq=0.001/0.448*100; %全桥

% sprintf("0.2232%");

% eHb=0.013/0.198*100;%半桥

% sprintf("6.5657%");

%

% sprintf("0.4464%");%quan

% sprintf("1.5152%");%ban

% sprintf("0.2232%");%quan

% sprintf("6.5657%");%ban

线性度 = (最大偏差值 / 测量长度) x 100%
eL=±(ΔLmax/yF.S)×100%
ΔLmax ——最大非线性误差； yF.S ——理论满量程输出。
然后迟滞是Eh=(△Hmax x yF.S)×100%
△Hmax代表正反行程间的最大误差值

附录数据图片

全桥                         半                       桥