最近处理一批线性增量式的数据,即线性回归的数据,用了一些常见的代码,存放之处,以供日后参考。
%欧阳的事情
A = importdata('test.xlsx');
x1 = A(1,:);
x2 = A(2,:);
y=[1 2 3 4 5 6];
figure(1)
plot(y,A(1,:),'-','LineWidth',1);
hold on;
plot(y,A(2,:),'-','LineWidth',1);
hold on;
plot(y,A(3,:),'-','LineWidth',1);
hold on;
plot(y,A(4,:),'-','LineWidth',1);
hold on;
plot(y,A(5,:),'-','LineWidth',1);
hold on;
plot(y,A(6,:),'-','LineWidth',1);
title('累计量/g');
hold on;
%计算差值
y1 = A(:,1);
y2 = A(:,2);
y3 = A(:,3);
y4 = A(:,4);
y5 = A(:,5);
y6 = A(:,6);
dy6 = y6 - y5;
dy5 = y5 - y4;
dy4 = y4 - y3;
dy3 = y3 - y2;
dy2 = y2 - y1;
dy1 = y1;
B = [dy1 dy2 dy3 dy4 dy5 dy6];
figure(2)
% ylim([0,2]);
plot(y,B(1,:),'-','LineWidth',1);
hold on;
plot(y,B(2,:),'-','LineWidth',1);
hold on;
plot(y,B(3,:),'-','LineWidth',1);
hold on;
plot(y,B(4,:),'-','LineWidth',1);
hold on;
plot(y,B(5,:),'-','LineWidth',1);
hold on;
plot(y,B(6,:),'-','LineWidth',1);
title('每格量/g');
hold on;
axis([ -inf inf 0 2]);
mean1 = zeros(1,6);
s1 = zeros(1,6);
for i = 1:6
AA = mean(B(i,:));
mean1(i) = AA;
end
figure(3);
plot(y , mean1,'LineWidth',1);
hold on;
mean0 = mean(mean1);
plot(y ,mean0);
plot([mean0 mean0 mean0 mean0 mean0 mean0],'--k','LineWidth',1);
title('平均值/g');
axis([ -inf inf 0 2]);%限制坐标轴的范围
|
请发表评论