3 лаба / src / График статической характеристики преобразования ЦПИ

function res = showPlot()

q = 0.001

lineX=[];
lineY=[];

for i=0:10
   lineX(i+1) = q * i;
   lineY(i+1) = q * i;
end

izmerenia =[
0,
1.16*0.001,
2.18*0.001,
3.19*0.001,
4.12*0.001,
5.27*0.001,
6.27*0.001,
7.29*0.001,
8.3*0.001,
9.33*0.001
];

j = 1;
polY(1) = 0;
for i=1:9
    polX(j) = izmerenia(i);
    polY(j+1) = lineY(i);
    j+=1;
    polX(j) = izmerenia(i+1);
    polY(j+1) = lineY(i+1);
    j+=1;
end

polX(j) = izmerenia(length(izmerenia));
#polY = polY(1: length(polY)-1);
polX;
polY;

errors=[];
for i=2:10
    errors(i-1) = izmerenia(i) - q*(i-1);
end

errorsY=[];
errorsX=[];
i = 1;
for j=1:length(errors)
    errorsX(i) = polX(i+1);
    errorsY(i) = 0 - errors(j) - q;
    i=i+1;
    errorsX(i) = polX(i+1);
    errorsY(i) = 0 - errors(j);
    i=i+1;
end

plot(polX, polY, 'b');
hold on;
plot(lineX, lineY, 'k--');
plot(errorsX, errorsY, 'r');
hold off;
grid on;

xticks(lineX);
yticks(lineY);