Скачиваний:
59
Добавлен:
01.05.2014
Размер:
419.33 Кб
Скачать

Частный случай теоремы 23

X1

X2

Xn

y1

y2

yn

xiyj

Si(x) : Si(x)=yi i=1,2…n

S(x) – многочлен проходящий через все узлы интерполяции – (xi,yi)

<S(x)>X-Xi=yi - остаток от деления S(x) на (x-xi) (Схема Горнера)

S(x)=(x-xi)+c

C(x)=(x-x1) (x-x2)… (x-xn)

Ci(x)=( ) (x-xi-1) (x-xi+1)…

<Di(x)Ci(x)>X-Xi =1

<Di(x) (x-x1) (x-x2)… (x-xi-1) (x-xi+1)… (x-xn)>X-Xi =1

Di(x) (x-x1) (x-x2)… (x-xi-1) (x-xi+1)… (x-xn)=1

S(x)=i=1nDi(x) Ci(x) Ui(x)=i=1nyi - Многочлен Лагранжа

Через n точек можно провести многочлен deg<=n-1

Недостаток: При добавлении данных надо пересчитывать весь многочлен.

  1. Интерпо­ляционная формула Ньютона.

Интерполяционный метод Ньютона.

Подход:

X1

X2

Xn

y1

y2

yn

S1(x) S2(x) … Sn(x) (x-x1)

S1(x)=y1 (x-x2)

…. …

Sk(x)-Sk-1(x) (x-xk-1)

x1,x2….xk-1

* a||b1 и a||b2 (b1;b2)=1 => a||b1b2

по аналогии с *

Sk(x)-Sk-1(x)||(x-x1)(x-x2) … (x-xk-1)

Sk(x)-Sk-1(x)=Ak(x-x1)(x-x2) … (x-xk-1)

Sn(x)=A1+A2(x-x1)+A3(x-x1)(x-x2)+…+An(x-x1)(x-x2) … (x-xn-1)

y1=A1 | x=x1

y2=A1+A2(x2-x1) | x=x2

Пример:

x1 x2 x3

-1 0 1

y1 y2 y3 S3(x) =A1+A2(x-x1)+A3(x-x1)(x-x2)

3 -1 2 3=A1 | x=-1

-1=A1+A2(1) | x=0 | A2=-4

2=A1+A2(2)+A3(2)(1) | A3=7/2

S3(x)=3-4(x+1)+7/2(x+1)x

S4(x)=S3(x)+A4(x-x1)(x-x2)(x-x3)

  1. Приближенная интерполяция (Метод наименьших квадратов).

Приближенная интерполяция

F(x1,x2…xn)

F(x,y)=sin x.cos x или F(x,y,z)=x2z+xy+xy2z3

;

;

Метод наименьших квадратов

S(xk)-yk=k

k=1nk min

Sm(x)=a0+a1x+a2x2+…+amxm

k=1nk2 =k=1n(a0+a1x+a2x2+…+amxm)2

Ф(a0,a1,…,am) min {a0,…,am}

(a,b)ax+b

k=17(axk+b-yk)2=Ф(a,b)

;

Ф(a,b)=(-3a+b+2)2+(-2a+b-4)2+(-a+b)2+(b-1)2+(a+b)2+(2a+b-3)2+(3a+b-2)2

(-3a+b+2)(-3)+(-2a+b-4)(-2)+(-a+b)(-1)+(b-1)(0)+(a+b)(1)+(2a+b-3)(2)+(3a+b-2)(3)=0

  1. Размещения и сочетания без повторений. Биномиальные ко­эффициенты и их свойства. Треугольник Паскаля.

Определение: способ размещения в определенном порядке некоторого числа элементов. Суммы S называют размещения.

S={a,b,c}

{ab,ac,ba,ca,bc,cb}

{ab,ac,bc}

Akn=n(n-1)(n-2)…(n-k+1)

Ann=Pn-n!

Ckn=Akn/Pk=n!/(n-k)!k!

Свойства:

1)

  1. k=0nCkn=2n

  2. k=0n(-1)kCkn=0

  3. Ckn=Cn-kn

  4. Ckn=Ck-1n-1+Ckn-1

  5. C00=Cnn=1

Треугольник Паскаля:

С0i, С1i, … Сii (x+y)i

1

1 1

1 2 1

1 3 3 1

1 4 6 4 1

1 5 10 10 5 1

  1. Тождество Коши:

Сkm+n=s=0kCsmCk-sn

m=n=k

Ck2k=s=0k(Csn)2

Доказательство:

M,n,k – надо выбирать

s=0kCsmCk-sn

  1. Размещения и сочетания с повторениями. Кодирование с ис­правлением ошибок. Граница Хемминга.

Размещение и сочетание с повторениями.

n-k

Akn=nk n, n, n … n - k раз

Сочетания

X={a1,a2,…an)

a1 x1 xi>=0

a2 x2 x1+x2+…+xn=k

.. ..

an xn

ai - xi+1 > 0

a1a1a1| a2a2a2|…|ananan|

x1+1 x2+1 xn+1

Cn-1n+k-1

Кодирование с исправлениями (кодирование Хэмминга)

a=(a1,…,an) ai{0,1}

P<2km r – не более, чем в r разрядах могут существовать ошибки.

P – число слов.

a=(a1,…,am) | d(a,b)={ai+bi} – количество разрядов, которые не совпадают друг с b=(b1,…,bm) | другом. – Расстояние от a до b.

Пример: a=(0,1,1,0) |

b=(1,1,0,1) | => d(a,b)=3. = функция расстояния между объектами

(A,B)M (A,B) M*M  R+ (декартово произведение).

  1. (A,B)>0

(A,B)=0 <=> A=B

  1. (A,B)=(B,A)

  2. (A,B)<=(A,C)+(C,B)

Докажем, что d(a,b) – расстояние т.е. обладает 1,2,3.

  1. очевидно

  2. очевидно

  3. d(a,b)<=d(a,c)+d(c,b)

ai=bi; 0<= …

aibi; 1<= Пусть нет вклада ai=ci |

ci=bi | => ai=bi Значит наше предположение не верно aibi;

Теорема: 26

S={a(1), a(2), … , a(p)} –

D(a(i), a(j))>=2r+1

a(i) – передали i разрядное слово, не более, чем в r разрядах ошибка.

b – переданное слово.

Известно d(a(i),b)<=r

Ai={b|d(a(i),b)<=r - окрестность с центром в точке b.

b<i

Окрестности не пересекаются:

Пусть это не так, тогда:

d(a(i),w)<=r | => d(a(i), b(i))<=d(a(i),w)+d(a(j),w)<=2+2=22 Значит наше пред-

d(a(j),w)<=r | положение неверно.

. .

. . . слова

Окрестность, любое слово на

расстоянии <r стало P<rn слов.

2m, r  P?

Ai ai Ai={a, d(a,ai)<=r}

|Ai|=C0m+ C1m+ C2m+…+Crm=rCkm – число слов в любой окрестности.

PrCkm<=2m => p<=2m/rCkm

Пример: Кодирования по Хэммингу.

m=3, r=1, p<=8/C03+C13=2

a(0)=(0,0,0) |

a(1)=(1,1,1) | d(a(0), a(1))=3=2k+1

A0={(0,0,0), (1,0,0), (0,1,0), (0,0,1)}

A1={(1,1,1), (1,1,0), (0,1,1), (1,0,1)}

(0,1,1) (1,1,1)

(0,1,0) (1,1,0)

(0,0,1) (1,0,1)

(0,0,0) (1,0,0)

  1. Полиномиальное кодирование. Код Шеннона-Фэно и алгоритм Хаффмена.

Соседние файлы в папке Шпаргалка по дискретной математике