Добавил:
Вуз:
Предмет:
Файл:
function res = showPlot()
q = 0.01
lineX=[];
lineY=[];
lineY(1) = q/5;
for i=0:10
lineX(i+1) = q * i;
lineY(i+2) = -(q/5) * i;
end
polX = [
13.5*0.001,
24.5*0.001,
34.3*0.001,
44.4*0.001,
54.5*0.001,
64.7*0.001,
74.4*0.001,
84.9*0.001,
95.5*0.001
];
polY = [
-8.5*0.001,
-9.5*0.001,
-9.3*0.001,
-9.4*0.001,
-9.5*0.001,
-9.7*0.001,
-9.4*0.001,
-9.9*0.001,
-10.5*0.001
];
xi = 0:0.01:1;
p = polyfit(polX, polY, 1);
plot(xi,polyval(p,xi), "b--");
hold on;
plot(polX, polY, 'r');
hold off;
grid on;
[maxx, nmaxx] = max(polX);
[miny, nminy] = min(polY);
maxx = maxx + 0.05*maxx;
miny = miny + 0.05*miny;
xlim([0 maxx]);
ylim([miny q/5]);
xticks(lineX);
yticks(lineY);
a = p(2)
b = p(1)
q = 0.01
lineX=[];
lineY=[];
lineY(1) = q/5;
for i=0:10
lineX(i+1) = q * i;
lineY(i+2) = -(q/5) * i;
end
polX = [
13.5*0.001,
24.5*0.001,
34.3*0.001,
44.4*0.001,
54.5*0.001,
64.7*0.001,
74.4*0.001,
84.9*0.001,
95.5*0.001
];
polY = [
-8.5*0.001,
-9.5*0.001,
-9.3*0.001,
-9.4*0.001,
-9.5*0.001,
-9.7*0.001,
-9.4*0.001,
-9.9*0.001,
-10.5*0.001
];
xi = 0:0.01:1;
p = polyfit(polX, polY, 1);
plot(xi,polyval(p,xi), "b--");
hold on;
plot(polX, polY, 'r');
hold off;
grid on;
[maxx, nmaxx] = max(polX);
[miny, nminy] = min(polY);
maxx = maxx + 0.05*maxx;
miny = miny + 0.05*miny;
xlim([0 maxx]);
ylim([miny q/5]);
xticks(lineX);
yticks(lineY);
a = p(2)
b = p(1)
Соседние файлы в папке src