mavzu. Matlabda interpolyatsiya va approksimatsiya masalasini yechish. Reja

Misol: ko‘phad ildizlarini topamiz.
Ko‘phad ildizlarini topishga teskari protsedura, ya’ni ko‘phadlarni tiklash, р=poly(c) funksiyasi asosida amalga oshiriladi, bu yerda с – ko‘phad ildizlari vector ustun; р – ko‘phad koeffitsentlari.
Ko‘phadning hosilasi dp=polyval(р) funksiyasi yordamida topiladi, bu yerda р – berilgan ko‘phad koeffitsentlari vektori; dp – ko‘phadning hosilasi koeffitsentlari vektori.
8
>> P=[1 -5 6]
P =
1 -5 6
>> C=roots(P)
C =
3.0000
2.0000
>>
Ko‘pgina amaliy masalalarini yechishda funksiyalarni approksimatsiya yoki interpolyatsiya yordamida yaqinlashtirish masalalarini yechish talab etiladi.
Approksimatsiya deganda bir funksiya (approksimatsiya-lanuvchi) ni berilgan qiymatlari va ma’lum kriteriy asosida boshqa eng yaxshi yaqinlashuvchi funksiyaga almashtirish tushuniladi.
Injenerlik amaliyotida odatda tekis va o‘rta kvadratik yaqinlashish kriteriysi qo‘llaniladi.
9
APPROKSIMATSIYA NIMA?
10
Approksimatsiya – bu lotincha approximatus so’zidan olingan bo’lib, (ingliz tilida approximation) juda yaqin degan ma’noni bildiradi.
Fanda yoki tadqiqotlar jarayonida “Approksimatsiya” tushunchasi mavjud modeldan (funksiyadan) foydalanish qiyin bo’lgan hollarda, model yoki jarayon sifatida qaraladi.
Approksimatsiya hisoblash jarayonini osonlashtiradi.
Shunday ma’lumotlar borki, ularni aniq funksiya bilan ifodalash imkoniyati yo’q, bunday hollarda ham approksimatsiyalash masalasidan foydalaniladi;
f(x) -> g(x)
Matlabда аппроксимацияловчи функция сифатида n – тартибли кўпҳад, аппроксимация критерийси сифатида ўрта квадратик четланиш ишлатилади.
Аппроксимациялаш функцияси қуйидаги кўринишга эга:
р=polyfit(x,y,n),
бу ерда:
x,y – бир хил ёки турли қадамдаги тугун нуқталар ва шу нуқтадаги берилган қийматлар;
n – аппроксимацияловчи полином тартиби;
р – аппроксимацияловчи полином коэффициентлари вектори.
11
12
13
14
15
Мисол. функциянинг бир хил қадамдаги тугун нуқталардаги қийматлари асосида 5-тартибли кўпҳад билан аппроксимация қилиш.
>> x=pi/8:pi/8:4*pi; y=sin(x)./x;
>> p=polyfit(x,y,5); fa=polyval(p,x);
>> subplot(3,1,1:2),plot(x,y,'-o',x,fa,':*'), grid, hold on;
>> error=abs(fa-y);subplot(3,1,3),plot(x,error,'--p')
16
Интерполяция деганда бир функциянинг кам сонли тугун нуқталари (интерполяция тугунлари)да берилган қийматлардан фойдаланиб, қийматлари берилган функциянинг тугун нуқталардаги қийматлари билан устма-уст тушувчи ва тугун нуқталар орасидаги ихтиёрий нуқтада функциянинг қийматларини ҳисоблашга имкон берувчи яқинлашувчи полином билан алмаштириш тушунилади.
17