Програма знаходження статичної моделі електричної печі у вигляді ряду Пайка і Сільверберга
clc
clear
P=[20 40 60 80 100 ];
U= [23.5 28.2 32.9];
T=[73 60 50 41 36; 89 71 61 55 47; 112 82 70 63 58],
d=input ('Acceptable error of approximation d=');
n=length(P);%number of columns
m=length(U);%number of rows
N=n*m;
Tser=mean(T(:));%mean value of temperature
y=T-Tser;
g0=mean(y),%averaged columns
h0=mean(y'),%averaged rows
%Calculation of functions g0 and h0
for j=1:3
for i=1:5
y1(j,i)=y(j,i)-g0(i)-h0(j);
end
end
A1=polyfit(P,g0,3);
g0r=polyval(A1,P);
figure(1); plot(P,g0,'*',P,g0r); grid;
xlabel('P , kPa'); ylabel('g0');
A2=polyfit(U,h0,1);
h0r=polyval(A2,U);
figure(2); plot(U,h0,'*',U,h0r); grid;
xlabel('U , V'); ylabel('h0');
for j=1:3
for i=1:5
Tr(j,i)=g0r(i)+h0r(j);
end
end
sigma=std(y1(:));
sigma0=sigma,
y2=y1;
AG=[]; AH=[]; SIGMA=[];
%Calculation of functions gs and hs
hs=y1(:,1);
while sigma>d,
for k=1:10
for i=1:5
gs(i)=y2(:,i)'*hs/(hs'*hs);
end
for j=1:3
hs(j)=y2(j,:)*gs'/(gs*gs');
end
end
g1=gs,
h1=hs',
Ag=polyfit(P,g1,3);
g1r=polyval(Ag,P);
figure(3); plot(P,g1,'*',P,g1r); grid;
xlabel('P , kPa'); ylabel('gs')
AG=[AG; Ag ];
Ah=polyfit(U,h1,2);
h1r=polyval(Ah,U);
figure(4); plot(U,h1,'*',U,h1r); grid
xlabel('U , V'); ylabel('hs');
AH=[AH; Ah];
y2=y2-h1'*g1;
for j=1:3
for i=1:5
Tr(j,i)=Tr(j,i)+h1r(j)*g1r(i);
end
end
sigma=std(y2(:));
SIGMA=[SIGMA; sigma];
hs=y2(:,1);
end
disp('Polynomial coefficients of function g0' )
A1,
disp('Polynomial coefficients of function h0')
A2
disp('Polynomial coefficients of functions gk')
AG,
disp('Polynomial coefficients of functions hk')
AH,
SIGMA
Ta=Tr+Tser,
delta=Ta-T,
maxdelta=max(abs(delta(:))),
Tmax=max(T(:));Tmin=min(T(:));
delzv=maxdelta/(Tmax-Tmin)*100,
k=0;
for j=1:m
k=k+1;
for i=1:n
U1(j,i)=U(k);
P1(j,i)=P(i);
end
|
|
|
end
figure(7), mesh(U1,P1,Ta),
ylabel('P , kPa'); zlabel('Ta , oC'); xlabel('U , V');
title('Relation Ta = f (P,U)');
ЗМІСТ
| Основні теоретичні відомості…………………………………………….. | 3 |
| 1. Мeтoд пoслiдoвнoгo вилучeння склaдoвих функцій з вихiднoї змінної ………………………………………………………………. | 4 |
| 2. Апроксимація функції двох змінних рядом Пайка і Сільверберга …………………………………………………………. | 6 |
| 3. Послідовна апроксимація за кожною вхідною величиною | 9 |
| 4. Побудова лінійної статичної моделі багатовимірного об’єкта методом найменших квадратів……………………………………. | 9 |
| Опис лабораторної установки……………………………………….…… | 10 |
| Послідовність виконання роботи …………………………………….… | 12 |
| Приклад знаходження статичної моделі…………………………….…… | 14 |
| Протокол лабораторної роботи…………………………………………… | 18 |
| Контрольні запитання……………………………………………………... | 19 |
| Список літератури…………………………………………………………. | 19 |
| Додаток…………………………………………………………………… | 20 |
Навчальне видання
Побудова статичної моделі багатовимірного ОБ’ЄКТА
Регулювання за результатами експериментальних досліджень
|
|
|
Методичні вказівки
для самостійної підготовки та інструкція
до лабораторної роботи № 2
з дисципліни “Ідентифікація та моделювання технологічних об’єктів”
для студентів базового напряму 6.050202
“Автоматизація та комп’ютерно-інтегровані технології”
Укладачі: Крих Ганна Бориславівна,
Кореньков Валентин Павлович,
Матіко Галина Федорівна.
Дата добавления: 2018-04-04; просмотров: 235; Мы поможем в написании вашей работы! |
Мы поможем в написании ваших работ!