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);