nihe.asv

峰值检测

clear all;close all;
Y=[1.0000	1.0486	1.0765	1.1070	1.1136	1.1071	1.1173	1.1179];
X=[0 10 20 50 100 200 500 1000];
% X=log10(X)
x1=X(1:4);
y1=Y(1:4);
x2=X(:6);
y2=Y(3:6);
p1=polyfit(x1,y1,1);
p2=polyfit(x2,y2,1);
scatter(X,Y)
hold on;
plot(x1,p1(1)*x1+p1(2),'k-');
plot(x2,p2(1)*x2+p2(2),'k-');
legend('最小二乘拟合')
xlabel('浓度');
ylabel('光强');
y01=p1(1)*X(1:3)+p1(2);
y02=p2(1)*X(4:6)+p2(2);
Y1=[y01 y02];
sum=0;
A=[];
tmp=0;
for i=1:6
    A(i)=(Y1(i)-Y(i))^2;
    tmp=tmp+A(i);
end
  B=sqrt(mean(A)*length(A))
  B2=sqrt(tmp)