2.3.2 Вид зависимости выхода аммиака от давления

%Зависимость выхода аммиака от давления.

clc, clear

P=[25 25.6 26.26 27.57 28.06 29.55 30.8];

w=[14 13.9 13.79 13.65 13.6 13.61 13.6];

n=length(P);

k=corrcoef(P, w);

X=[ones(n,1) P P.^2];

A=(X*X)^(-1)*X*w

x=24:0.1:31;

f=A(1)+A(2)*x+A(3)*x.^2;

plot(P,w,*,x,f,g);

%Проверка значимости по критерию Стьюдента.

l=n-1;

p=0.135

tk=tinv(p,l)

dw=cov(w);

c=diag((X*X)^(-1));

k=length(A);

for i=1:k

tr(i)=abs(A(i))/sqrt(dw*c(i));

end

tr

%Проверка адекватности по критерию Фишера.

f1=n-1;

f2=n-k;

Fk=finv(p,f1,f2)

Y=X*A;

dost=((Y-w)*(Y-w))/f2;

Fr=dw/dost

%Модель адекватна, т.к. Fr>Fk!

Результаты расчётов:

A =31.9385

-1.2462

0.0211

p = 0.1350

tk = -1.2150

tr =2.1861 1.1836 1.1193

Fk =0.3730

Fr = 61.7702

Графически зависимость отображена на рисунке 2.3.

Рисунок 2.3 - Выход аммиака в зависимости от давления

Математическая модель выхода аммиака от давления имеет следующий вид:

W=31.9385-1.2462x+0.0211x2 (2.20)