Банк рефератов содержит более 364 тысяч рефератов, курсовых и дипломных работ, шпаргалок и докладов по различным дисциплинам: истории, психологии, экономике, менеджменту, философии, праву, экологии. А также изложения, сочинения по литературе, отчеты по практике, топики по английскому.
Полнотекстовый поиск
Всего работ:
364150
Теги названий
Разделы
Авиация и космонавтика (304)
Административное право (123)
Арбитражный процесс (23)
Архитектура (113)
Астрология (4)
Астрономия (4814)
Банковское дело (5227)
Безопасность жизнедеятельности (2616)
Биографии (3423)
Биология (4214)
Биология и химия (1518)
Биржевое дело (68)
Ботаника и сельское хоз-во (2836)
Бухгалтерский учет и аудит (8269)
Валютные отношения (50)
Ветеринария (50)
Военная кафедра (762)
ГДЗ (2)
География (5275)
Геодезия (30)
Геология (1222)
Геополитика (43)
Государство и право (20403)
Гражданское право и процесс (465)
Делопроизводство (19)
Деньги и кредит (108)
ЕГЭ (173)
Естествознание (96)
Журналистика (899)
ЗНО (54)
Зоология (34)
Издательское дело и полиграфия (476)
Инвестиции (106)
Иностранный язык (62792)
Информатика (3562)
Информатика, программирование (6444)
Исторические личности (2165)
История (21320)
История техники (766)
Кибернетика (64)
Коммуникации и связь (3145)
Компьютерные науки (60)
Косметология (17)
Краеведение и этнография (588)
Краткое содержание произведений (1000)
Криминалистика (106)
Криминология (48)
Криптология (3)
Кулинария (1167)
Культура и искусство (8485)
Культурология (537)
Литература : зарубежная (2044)
Литература и русский язык (11657)
Логика (532)
Логистика (21)
Маркетинг (7985)
Математика (3721)
Медицина, здоровье (10549)
Медицинские науки (88)
Международное публичное право (58)
Международное частное право (36)
Международные отношения (2257)
Менеджмент (12491)
Металлургия (91)
Москвоведение (797)
Музыка (1338)
Муниципальное право (24)
Налоги, налогообложение (214)
Наука и техника (1141)
Начертательная геометрия (3)
Оккультизм и уфология (8)
Остальные рефераты (21697)
Педагогика (7850)
Политология (3801)
Право (682)
Право, юриспруденция (2881)
Предпринимательство (475)
Прикладные науки (1)
Промышленность, производство (7100)
Психология (8694)
психология, педагогика (4121)
Радиоэлектроника (443)
Реклама (952)
Религия и мифология (2967)
Риторика (23)
Сексология (748)
Социология (4876)
Статистика (95)
Страхование (107)
Строительные науки (7)
Строительство (2004)
Схемотехника (15)
Таможенная система (663)
Теория государства и права (240)
Теория организации (39)
Теплотехника (25)
Технология (624)
Товароведение (16)
Транспорт (2652)
Трудовое право (136)
Туризм (90)
Уголовное право и процесс (406)
Управление (95)
Управленческие науки (24)
Физика (3463)
Физкультура и спорт (4482)
Философия (7216)
Финансовые науки (4592)
Финансы (5386)
Фотография (3)
Химия (2244)
Хозяйственное право (23)
Цифровые устройства (29)
Экологическое право (35)
Экология (4517)
Экономика (20645)
Экономико-математическое моделирование (666)
Экономическая география (119)
Экономическая теория (2573)
Этика (889)
Юриспруденция (288)
Языковедение (148)
Языкознание, филология (1140)

Контрольная работа: Интерполирование функций

Название: Интерполирование функций
Раздел: Рефераты по математике
Тип: контрольная работа Добавлен 22:42:53 07 января 2011 Похожие работы
Просмотров: 3188 Комментариев: 4 Оценило: 1 человек Средний балл: 3 Оценка: неизвестно     Скачать

Содержание

Введение

1. Формула Лагранжа

2. Интерполирование по схеме Эйткена

3. Интерполяционные формулы Ньютона для равноотстоящих узлов

4. Формула Ньютона с разделенными разностями

5. Интерполяция сплайнами

Заключение

Список литературы

Введение

Цель работы: изучение и сравнительный анализ методов интерполяции функций; реализация этих методов в виде машинных программ на языке высокого уровня и практическое решение задач интерполяции на ЭВМ.

При разработке математического обеспечения САПР часто приходится иметь дело с функциями f (x ), заданными в виде таблиц, когда известны некоторое конечное множество значений аргумента и соответствующие им значения функции. Аналитическое выражение функции f (x ) при этом неизвестно, что не позволяет определять ее значения в промежуточных точках аргумента, отсутствующих в таблице. В таком случае решается задача интерполирования, которая формулируется следующим образом.

На отрезке [a , b ] заданы n + 1 точки x 0 , x 1 , ..., xn , которые называются узлами интерполяции, и значения некоторой функции f (x ) в этих точках f (x 0 ) = y 0 ,f (x 1 ) = y 1 , ..., f (x n ) = y n . Требуется построить интерполирующую функцию F (x ), принимающую в узлах интерполяции те же значения, что и f (x ), т.е. такую, что F (x 0 ) = y 0 ,F (x 1 ) = y 1 , ...,F (x n ) = yn .

Геометрически это означает, что нужно найти кривую y = F (x ) некоторого определенного типа, проходящую через заданную систему точек Mi (xi , yi ) дляi = . Полученная таким образом интерполяционная формула y = F (x ) обычно используется для вычисления значений исходной функции f (x ) для значений аргумента x , отличных от узлов интерполяции. Такая операция называется интерполированием функции f (x ). При этом различают интерполирование в узком смысле, когда x принадлежит интервалу [x 0 , xn ], и экстраполирование, когда x не принадлежит этому интервалу.

В такой общей постановке задача интерполирования может иметь бесчисленное множество решений. Чтобы получить единственную функцию F (x ), необходимо предположить, что эта функция не произвольная, а удовлетворяет некоторым дополнительным условиям.

В простейшем случае предполагается, что зависимость y = f (x ) на каждом интервале (xi , xi +1 ) является линейной. Тогда для каждого участка (xi , xi +1 ) в качестве интерполяционной формулы y =F (x ) используется уравнение прямой, проходящей через точки Mi (xi , y i ) и Mi +1 (xi +1 , y i +1 ), которое имеет вид

. ( 1)

Припрограммировании процедур линейной интерполяции следует учитывать, что процесс решения задачи интерполирования с использованием формулы (1) включают два этапа: выбор интервала (xi , xi +1 ), которому принадлежит значение аргумента х ; собственно вычисление значения y =F (x ) по формуле (1).

На практике в качестве интерполирующей функции F (x ) обычно используется алгебраический многочлен

Pn (x ) = a 0 + a 1 x + a 2 x 2 + ... +an xn

степени не выше n , такой, что Pn (x 0 ) = y 0 ,Pn (x 1 ) = y 1 ,...,Pn (xn ) = yn . Наиболее известными методами построения интерполяционного многочлена Pn (x )являются метод Лагранжа, итерационные и разностные методы.


1. Формула Лагранжа

Интерполяционная формула Лагранжа обеспечивает построение алгебраического многочлена Pn (x ) для произвольно заданных узлов интерполирования. Для n + 1 различных значений аргумента x 0 , x 1 , ..., xn и соответствующих значений функции f (x 0 ) = y 0 ,f (x 1 ) = y 1 , ..., f (x n ) = y n интерполяционная формула Лагранжа имеет вид

,

где х - значение аргумента функции, расположенного в интервале [x 0 , xn ].

Необходимо отметить, что формула Лагранжа, в отличие от других интерполяционных формул, содержит явно yi (i = ), что бывает иногда важно.

Пример 1. Построить интерполяционный многочлен Лагранжа для функции, заданной следующей таблицей.

x 0 = 0, x 1 = 1, x 2 = 2, x 3 = 5,
y 0 = 2, y 1 = 3, y 2 = 12, y 3 = 147.

Для случая четырех узлов интерполяции (n = 3) многочлен Лагранжа представляется следующим образом:

Заменив переменные xi , yi (i = )их числовыми значениями, получим интерполяционный многочлен


Интерполирование по формуле Лагранжа связано с большим объемом вычислений, значительная часть которых повторяется при получении нескольких значений Pn (x ) для одной функции f (x ). В том случае, когда формула Лагранжа используется для многократного получения значений одной функции при различных значениях аргумента, можно значительно уменьшить объем вычислений. Для этого формула Лагранжа представляется в виде


где - лагранжевы коэффициенты, определяемые как


Вычисление лагранжевых коэффициентов выполняется по следующей схеме, удобной при использовании ЭВМ. Составляется таблица разностей:

Произведение элементовi -й строки обозначается через Ki . Отсюда лагранжевы коэффициенты вычисляются по формуле


где Пn +1 (x ) = (x - x 0 )(x - x 1 )…(x - xn ) - произведение элементов главной диагонали таблицы (эти элементы подчеркнуты). Тогда формула Лагранжапринимает вид:


Использование формулы (2) позволяет сократить значительную часть вычислений по определению лагранжевых коэффициентов Li (n ) (x )при различных значениях аргумента. Для этого произведение элементов i -й строки таблицы разностей представляется как Ki = (xxi )Di , где Di - произведение всех элементов строки, кроме расположенного на главной диагонали. Величина Di (i= )не зависит от значения аргумента x и может быть вычислена для заданной функции только один раз.

2. Интерполирование по схеме Эйткена

Итерационные методы интерполирования основаны на повторном применении некоторой простой интерполяционной схемы. Наиболее известным из итерационных методов является метод Эйткена, в основе которого лежит многократное применение линейной интерполяции.

В соответствии со схемой Эйткена линейная интерполяция по точкам Mi (xi , yi ) и Mi +1 (xi +1 , yi +1 ) сводится к вычислению определителя второго порядка



При интерполировании по трем и более точкам последовательно вычисляются многочлены



В общем случае интерполяционныймногочлен n -й степени, принимающий в точках xi значения yi (i = ), записываются следующим образом:


(3)

Основным достоинством схемы Эйткена является возможность постепенного увеличения числа используемых значений xi до тех пор, пока последовательные значения P 0,1,2,…,n (x ) и P 1,2,…,n -1 (x ) не совпадут в пределах заданной точности. Иначе говоря, вычисления прекращаются при выполнении условия

|P 0,1,2,…, n (x ) - P 1,2,…, n -1 (x )| < e (k £n ).

При использовании ЭВМ вычисления по формуле (3) реализуются в виде рекурсивной подпрограммы - функции РХ(I, J) с формальными параметрами I, J, определяющими индексы крайних узлов интерполирования, которые используются для получения значения соответствующего многочлена Pi ,i +1,…, j (x ).

Для хранения вычисленных значений P (x )используется двумерный массив M размером N*N элементов, где N - максимальное число узлов интерполирования. Каждому возможному значению P (x ) соответствует один из элементов M(I, J), расположенный выше главной диагонали (I < J) и определяемый сочетанием индексов крайних узлов интерполирования.

Например, значению многочлена P 1,2 (x ) соответствует элемент M(1,2), значению P 2,3,4 (x ) - элемент M(2, 4) и т.д. Симметричные элементы M(J, I), расположенные ниже главной диагонали (J > I), показывают, вычислены ли соответствующие значения P (x ) на данный момент, и определяются как


Схема рекурсивной процедуры PX приведена на рис. 1, где Х- массив значений узлов интерполирования, Y- массив значений функциивузлах интерполирования, Z- значение аргумента. Параметры X, Y, Z, M должны быть описаны как общие для главной программы и подпрограммы PX.

3. Интерполяционные формулы Ньютона для равноотстоящих узлов

Узлы интерполирования x 0 , x 1 , ..., xn называются равноотстоящими, если , гдеh - шаг интерполирования. При этом для некоторой функции f (x ) таблично задаются значения yi = f (xi ), где xi = x 0 + ih .



Существуют две формулы Ньютона для случая равноотстоящих узлов интерполирования, которые называются соответственно первой и второй интерполяционными формулами Ньютона и имеют вид:

;

,

В этих формулах Di yj - конечные разности, где i - порядок разности, j - ее порядковый номер, а параметры t и q определяются следующим образом:

t = (x - x 0 ) / h ; q = (x - xn ) / h .

Конечные разности первого порядка вычисляются как Dyj = yj +1yj , где

j = , для более высоких порядков используется известная формула

(i = 2, 3, ...; j = ).

Получаемые конечные разности удобно представлять в табличной форме записи, например, в виде табл. 1, которая называется горизонтальной таблицей конечных разностей.

Таблица 1

x y Dy D2 y D3 y D4 y
x 0 Y 0 Dy 0 D2 y 0 D3 y 0 D4 y 0
x 1 Y 1 Dy 1 D2 y 1 D3 y 1 D4 y 1
x 2 Y 2 Dy 2 D2 y 2 D3 y 2
x 3 Y 3 Dy 3 D2 y 3 -
x 4 Y 4 Dy 4 - -
x 5 Y 5 - - -

Пepвая формула Ньютона применяется для интерполирования вперед и экстраполирования назад, т.е. в начале таблицы разностей, где строки заполнены и имеется достаточное число конечных разностей. При использовании этой формулы для интерполирования значение аргумента x должно лежать в интервале [x 0 , x 1 ]. При этом за x 0 может приниматься любойузел интерполяции xk с индексом , где m - максимальный порядок конечных разностей.

Вторая формула Ньютона применяется для интерполирования назад и экстраполирования вперед, т.е. в конце таблицы конечных разностей. При этом значение аргумента x должно находиться в интервале [xn -1 , xn ], причем за xn может приниматься любой узел интерполирования .

Одно из важнейших свойств конечных разностей заключается в следующем. Если конечные разности i –го порядка (i < n ) постоянны, то функция представляет собой полином i –й степени. Следовательно, формула Ньютона должна быть не выше i -й степени. При использовании ЭВМ вычисление конечных разностей завершается при выполнении условий

где L - число значащих цифр после запятой в представлении значений функции.

Необходимо отметить, что формулы Ньютона являются видоизменениямиформулы Лагранжа. Однако в формуле Лагранжа нельзя пренебречь ни одним из слагаемых, так как все они равноправны и представляют многочлены n -й степени. В формулы Ньютона в качестве слагаемых входят многочлены повышающихся степеней, коэффициентами при которых служат конечные разности, разделенные на факториалы. Конечные разности, как правило, быстро уменьшаются, что позволяет в формулах Ньютона пренебречь слагаемыми, коэффициенты при которых становятся малыми. Это обеспечивает вычисление промежуточных значений функции достаточно точно спомощью простых интерполяционных формул.


4. Формула Ньютона с разделенными разностями

Первая и вторая формулы Ньютона предполагают, что узлы интерполирования являются равноотстоящими. Однако, в общем случае функция f (x ) может быть задана таблицей, в которой узлы находятся на произвольном расстоянии друг от друга , где значения hi (i = ) являются различными.

При таких условиях первая и вторая интерполяционные формулы Ньютона неприменимы. В данном случае, для решения задачи интерполяции применяются не конечные, а разделенные разности.

Разделенная разность первого порядка определяется:

Для вычисления разделенных разностей высших порядков используется формула:

Разделенные разности удобно представлять диагональной таблицей, вид которой для n = 4 соответствует табл. 2.

Таблица 2

Интерполяционный многочлен Ньютона, использующий разделенные разности, имеет вид:

где , Пk (x ) = 1.

Представленная формула позволяет повышать точность вычислений постепенно, добавляя разделенные разности более высоких порядков. Следует отметить, что при этом все полученные результаты сохраняются, т.е. не вычисляются заново, а только наращиваются. Это следует из соотношения

Оценка погрешности интерполирования выполняется по формуле


5. Интерполяция сплайнами

Пусть задана таблица значений функции f (xi ) = yi (), в которой они расположены по возрастанию значений аргумента: x 0 < x 1 < … < xn . Чтобы построить кубический сплайн, требуется определить коэффициенты ai 0 , ai 1 , ai 2 , ai 3 , которые задают интерполяционный кубический многочлен

на каждом интервале интерполирования [xi -1 , xi ], .

Таким образом, необходимо определить 4n коэффициентов aij (, ), для чего требуется 4n уравнений. Необходимые уравнения определяются следующими условиями.

1. Условия непрерывности функции:

2. Условия непрерывности 1-х и 2-х производных функции:

3. Граничные условия:

Часто используются граничные условия видаПолучаемый при этом сплайн называется естественным кубическим сплайном.

Задача определения кубического сплайна существенно упрощается при использовании многочлена Эрмита. Кубический многочлен Эрмита на интервале [xi- 1 , xi ] определяется с помощью значений функции yi -1 , yi и ее производных y ¢i -1 , y ¢i . Так как значения производных в общем случае могут быть неизвестны, обозначим их как y ¢i -1 = Si -1 ; y ¢i = Si . При построении сплайна переменные S i называются наклонами сплайна в соответствующих точках xi .

Запишем многочлен Эрмита для интервала [xi- 1 , xi ], где hi = xi - xi- 1 :

При таком выборе кубического многочлена автоматически выполняются условия непрерывности функции и ее первых производных:

Чтобы определить сплайн, нужно задать условия непрерывности второй производной:

Для записи этих условий в развернутом виде определим кубический многочлен Эрмита на интервале [xi , xi +1 ], где hi +1 = xi +1 - xi :

Определим вторые производные многочленов Qi (x ) и Qi +1 (x ) в точке x = xi :

(4)

(5)

Отсюда условие непрерывности вторых производных имеет вид:

(6)

Это условие порождает систему линейных уравнений относительно наклонов сплайна Si , которая содержит n - 1 уравнение и n + 1 переменную. Чтобы определить два недостающих уравнения используются граничные условия. Например, для естественного кубического сплайна:

Указанные граничные условия могут быть получены из уравнения (5) для i = 0 и из уравнения (4) для i = n соответственно. В развернутом виде:


(7)

Решение системы линейных уравнений, образованной условиями (6) и (7), позволяет вычислить наклоны сплайна Si (i = ) и определить кубический сплайн путем записи многочлена Эрмита для каждого интервала [xi- 1 , xi ],i = .


Заключение

В вычислительной математике существенную роль играет интерполяция функций, т.е. построение по заданной функции другой (как правило, более простой), значения которой совпадают со значениями заданной функции в некотором числе точек. Причем интерполяция имеет как практическое, так и теоретическое значение. На практике часто возникает задача о восстановлении непрерывной функции по ее табличным значениям, например полученным в ходе некоторого эксперимента. Для вычисления многих функций оказывается эффективно приблизить их полиномами или дробно-рациональными функциями. Теория интерполирования используется при построении и исследовании квадратурных формул для численного интегрирования, для получения методов решения дифференциальных и интегральных уравнений.


Список литературы

1. В.В. Иванов. Методы вычислений на ЭВМ. Справочное пособие. Изд-во "Наукова думка". Киев. 1986.

2. Н.С. Бахвалов, Н.П. Жидков, Г.М. Кобельков. Численные методы. Изд-во "Лаборатория базовых знаний". 2003.

3. И.С. Березин, Н.П. Жидков. Методы вычислений. Изд. ФизМатЛит. Москва. 1962.

4. К. Де Бор. Практическое руководство по сплайнам. Изд-во "Радио и связь". Москва. 1985.

5. Дж. Форсайт, М.Мальком, К. Моулер. Машинные методы математических вычислений. Изд-во "Мир". Москва. 1980.

Оценить/Добавить комментарий
Имя
Оценка
Комментарии:
Где скачать еще рефератов? Здесь: letsdoit777.blogspot.com
Евгений08:14:29 19 марта 2016
Кто еще хочет зарабатывать от 9000 рублей в день "Чистых Денег"? Узнайте как: business1777.blogspot.com ! Cпециально для студентов!
08:57:12 29 ноября 2015
если бы еще программы на Паскале
12:37:05 08 декабря 2011
не все формулы видно =(
21:07:46 22 ноября 2011Оценка: 3 - Средне

Работы, похожие на Контрольная работа: Интерполирование функций
Математические основы теории систем
ОГЛАВЛЕНИЕ Оглавление 1 Введение 3 Объект и устройство 3 Задачи управления 4 Матричный формализм в теории систем 6 Линейные операторы 6 Инвариантное ...
(6) ѭn=M[Xn]= = xn f(x)dx
(14) ѭ0,..,0,1,0,..,0=M[(x1)°...(xi-1)° xi (xi+1)°...(xn)°=M[xi]
Раздел: Рефераты по математике
Тип: реферат Просмотров: 1245 Комментариев: 3 Похожие работы
Оценило: 1 человек Средний балл: 5 Оценка: неизвестно     Скачать
Сравнительный анализ численных методов
Министерство образования и науки Республики Казахстан Карагандинский Государственный Технический Университет Кафедра ____САПР_ ПОЯСНИТЕЛЬНАЯ ЗАПИСКА ...
Допустим, в n+1 точке заданы значения x0,x1,.xn и соответствующие им значения f (x0), f (x1), ., f (xn).
Интерполяция заключается в выборе функции ѭ (х), значения которой в узлах интерполяции совпадают со значениями f (xi).
Раздел: Рефераты по информатике, программированию
Тип: дипломная работа Просмотров: 1573 Комментариев: 2 Похожие работы
Оценило: 0 человек Средний балл: 0 Оценка: неизвестно     Скачать
Численные методы решения типовых математических задач
ФЕДЕРАЛЬНОЕ АГЕНСТВО ПО ОБРАЗОВАНИЮ ГОСУДАРСТВЕННОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ Тульский государственный университет Кафедра ...
Докажем, что многочлен Лагранжа является интерполяционным многочленом, проходящим через все узловые точки, т.е. в узлах интерполирования xi выполняется условие Ln(xi) = yi.
Чем больше узлов интерполирования на отрезке [x0,xn] , тем точнее интерполяционный многочлен приближает заданную табличную функцию, т.е. тем точнее равенство:
Раздел: Рефераты по математике
Тип: курсовая работа Просмотров: 4370 Комментариев: 2 Похожие работы
Оценило: 0 человек Средний балл: 0 Оценка: неизвестно     Скачать
Сравнительный анализ численных методов
Министерство образования и науки Республики Казахстан Карагандинский государственный технический университет Кафедра САПР Курсовая работа по ...
При этом потребуем, чтобы каждый многочлен обращался в нуль во всех узлах интерполяции, за исключением одного (i-го), где он должен равняться единице.
Предположим, что в окрестности точки xi функция F(x)дифференцируема достаточное число раз.
Раздел: Рефераты по математике
Тип: курсовая работа Просмотров: 882 Комментариев: 2 Похожие работы
Оценило: 0 человек Средний балл: 0 Оценка: неизвестно     Скачать
Программа вычисления минимума заданной функции
1. Индивидуальное задание Вычислить минимум функции F(x)=L(x1)x2-2.5 L(x2)x-3 на отрезке [a;b] с точностью ѭ. L(x1), L(x2) значения интерполяционного ...
Этим условием является g(xi,a0,a1,.an)=f(xi) при i=0,n, которое означает, что аппроксимируемая функция f(x) совпадает с g(xi,a0,a1,.an) в т.н. узлах интерполяции x0,x1,.,xn.
интерполяционные формулы Ньютона Выберем для решения задачи интерполяции интерполяционный многочлен Лагранжа, так как его построение просто в алгоритмизации, не требует вычисления ...
Раздел: Рефераты по информатике, программированию
Тип: курсовая работа Просмотров: 251 Комментариев: 2 Похожие работы
Оценило: 0 человек Средний балл: 0 Оценка: неизвестно     Скачать
Выбор и построение интерполирующей функции
Министерство науки и образования Украины Сумской государственный университет кафедра информатики Численные методы Курсовая работа на тему: " Выбор ...
Если и значение близко к , то выгодно использовать первый интерполяционный многочлен Ньютона, тогда и Таким образом, первая интерполяционная формула Ньютона применяется для ...
В данной программе есть два массива: одномерный массив X, в нем хранятся значения узлов интерполирования хi и двумерный массив Р, в нем хранятся значения многочленов степени не ...
Раздел: Рефераты по математике
Тип: курсовая работа Просмотров: 577 Комментариев: 3 Похожие работы
Оценило: 1 человек Средний балл: 5 Оценка: неизвестно     Скачать

Все работы, похожие на Контрольная работа: Интерполирование функций (1299)

Назад
Меню
Главная
Рефераты
Благодарности
Опрос
Станете ли вы заказывать работу за деньги, если не найдете ее в Интернете?

Да, в любом случае.
Да, но только в случае крайней необходимости.
Возможно, в зависимости от цены.
Нет, напишу его сам.
Нет, забью.



Результаты(150629)
Комментарии (1838)
Copyright © 2005-2016 BestReferat.ru bestreferat@mail.ru       реклама на сайте

Рейтинг@Mail.ru