Grupė.
Darbo tikslas – susipažinti su interpoliavimo polinomų galimybėmis aprašant objektų statines charakteristikas. Įgyti praktinius įgūdžius darbui su modernia programine įranga MATLAB, kuri procesų valdymo specialistų tarpe labiausiai paplitusi pasaulyje.
% Ivedame ekspermentinius duomenius
xe=[0.0 9.4 20.6 49.4 80.4 89.2 94.9 100.0];
ye=[0.0 0.3 2.5 7.7 11.4 13.7 14.2 14.7];
% Ekspermentinius taskus atvaizduojame grafiskai
figure(1); plot(xe,ye,’*’);grid;
% Parenkame 4 interpoliavimo mazgus, 3 eiles Lagranzo polinomui
xi=[0.0 20.6 49.4 100.0]; yi=[0.0 2.5 7.7 14.7];
% Apskaiciuojame polinomo koeficientus
a0=yi(1)
a1=(yi(2)-a0)/(xi(2)-xi(1))
a2=(yi(3)-a0-a1*(xi(3)-xi(1)))/((xi(3)-xi(1))*(xi(3)-xi(2)))
a3=(yi(4)-a0-a1*(xi(4)-xi(1))-a2*(xi(4)-xi(1))*(xi(4)-xi(2)))/...
((xi(4)-xi(1))*(xi(4)-xi(2))*(xi(4)-xi(3)))
% Grafiskai atvaizduojame gauto poliomo ir ekspermento vertes
i=0;
for int=0.0:2:100
i=i+1; x(i)=int;
p3(i)=a0+...
a1*(int-xi(1))+...
a2*(int-xi(1))*(int-xi(2))+...
a3*(int-xi(1))*(int-xi(2))*(int-xi(3));
end
figure(2);
plot(x,p3,xe,ye,’o’,xi,yi,’r*’); grid;
legend(‘3-io laipsnio polinomas’,’ekspermentiniai taskai’,’interpoliavimo taskai’)
title(‘Voztuvo charakteristika – 3 laipsnio polinomas’,’FontName’,’Times NewRoman baltic’);
xlabel(‘Voztuvo atidarymas, %’,’FontName’,’Times NewRoman baltic’);
ylabel(‘Srautas, m3/h’,’FontName’,’Times NewRoman baltic’);
% Modeliavimo tikslumui padidinti ivedame papildoma interpoliavimo taska
xi(5)=89.2; yi(5)=13.7;
% Apskaiciuojame nauja koeficienta a4
a4=(yi(5)-a0-a1*(xi(5)-xi(1))-a2*(xi(5)-xi(1))*(xi(5)-xi(2))-...
a3*(xi(5)-xi(1))*(xi(5)-xi(2))*(xi(5)-xi(3)))/...
((xi(5)-xi(1))*(xi(5)-xi(2))*(xi(5)-xi(3))*(xi(5)-xi(4)))
% Grafiskai atvaizduojame naujas polinomo vertes
i=0; clear x; clear p;
for int=0.0:2:100
i=i+1; x(i)=int;
p4(i)=a0+...
a1*(int-xi(1))+...
a2*(int-xi(1))*(int-xi(2))+...
a3*(int-xi(1))*(int-xi(2))*(int-xi(3))+...
a4*(int-xi(1))*(int-xi(2))*(int-xi(3))*(int-xi(4));
end
figure(3);
plot(x,p4,’b’,xe,ye,’o’,xi,yi,’r*’); grid;
legend(‘4-io laipsnio polinomas’,’ekspermentiniai taskai’,’interpoliavimo taskai’)
title(‘Voztuvo charakteristika – 4 laipsnio polinomas’,’FontName’,’Times NewRoman baltic’);
xlabel(‘Voztuvo atidarymas, %’,’FontName’,’Times NewRoman baltic’);
ylabel(‘Srautas, m3/h’,’FontName’,’Times NewRoman baltic’);
% Nuokrypu nuo pagrindinio ir papildomo ekspermento tasku skaiciavimas
xe=[];
ye=[];
xe=[0.0 9.4 20.6 49.4 80.4 89.2 94.9 100.0 4.0 17.0 47.0 72.0 85.0 97.0];
ye=[0.0 0.3 2.5 7.7 11.4 13.7 14.2 14.7 0.2 1.4 6.2 10.0 12.6 14.5];
% Pagrindinio ir papildomo ekspermentinius taskus atvaizduojame grafiskai
figure(4); plot(xe,ye,’*’); grid;
%Nuokrypu kvadratu sumos skaiciavimas trecio laipsnio polinomui;
IIIkvsum=0;
for i=1:1:14;
nuokrIII(i)=(ye(i)-(a0+a1*(xe(i)-xi(1))+a2*(xe(i)-xi(1))*(xe(i)-xi(2))+a3*(xe(i)-xi(1))*(xe(i)-xi(2))*(xe(i)-xi(3))))^2;
IIIkvsum=IIIkvsum+nuokrIII(i)
end
%Nuokrypu kvadratu sumos skaiciavimas ketvirto laipsnio polinomui;
Ivkvsum=0;
for i=1:1:14;
nuokrIV(i)=(ye(i)-(a0+a1*(xe(i)-xi(1))+a2*(xe(i)-xi(1))*(xe(i)-xi(2))+a3*(xe(i)-xi(1))*(xe(i)-xi(2))*(xe(i)-xi(3))+a4*(xe(i)-xi(1))*(xe(i)-xi(2))*(xe(i)-xi(3))*(xe(i)-xi(4))))^2;
Ivkvsum=Ivkvsum+nuokrIV(i)
end
2) Interpoliavimo polinomų koeficientai (3 ir 4 – os eilės):
(panaudojus pagrindinio eksperimento duomenis)
Polinomo koeficientai:
a0...
Šį darbą sudaro 609 žodžiai, tikrai rasi tai, ko ieškai!
★ Klientai rekomenduoja
Šį rašto darbą rekomenduoja mūsų klientai. Ką tai reiškia?
Mūsų svetainėje pateikiama dešimtys tūkstančių skirtingų rašto darbų, kuriuos įkėlė daugybė moksleivių ir studentų su skirtingais gabumais. Būtent šis rašto darbas yra patikrintas specialistų ir rekomenduojamas kitų klientų, kurie po atsisiuntimo įvertino šį mokslo darbą teigiamai. Todėl galite būti tikri, kad šis pasirinkimas geriausias!
Norint atsisiųsti šį darbą spausk ☞ Peržiūrėti darbą mygtuką!
Mūsų mokslo darbų bazėje yra daugybė įvairių mokslo darbų, todėl tikrai atrasi sau tinkamą!
Panašūs darbai
Atsisiuntei rašto darbą ir neradai jame reikalingos informacijos? Pakeisime jį kitu nemokamai.
Pirkdamas daugiau nei vieną darbą, nuo sekančių darbų gausi 25% nuolaidą.
Išsirink norimus rašto darbus ir gauk juos akimirksniu po sėkmingo apmokėjimo!