Исследование модели фрактального броуновского движения
Статистическое моделирование системы с возмущениями в виде фрактального броуновского движения. Интерполяция, экстраполяция и прогнозирование процесса по наблюдениям в двух точках. Анализ дифференциальной системы на основе фильтрации Калмана-Бьюси.
Рубрика | Математика |
Вид | дипломная работа |
Язык | русский |
Дата добавления | 01.05.2015 |
Размер файла | 586,1 K |
Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже
Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.
if (scale > 249 ) {
scale = 249;
}
Image1->Canvas->Pen->Color = clSilver;
Image1->Canvas->MoveTo(5,250);
Image1->Canvas->LineTo(405,250);
Image1->Canvas->MoveTo(405,250);
Image1->Canvas->LineTo(405,255);
Image1->Canvas->MoveTo(5,250-scale);
Image1->Canvas->LineTo(5,250+scale);
Image1->Canvas->MoveTo(5,250-scale);
Image1->Canvas->LineTo(10,250-scale);
Image1->Canvas->MoveTo(5,250+scale);
Image1->Canvas->LineTo(10,250+scale);
Image1->Canvas->MoveTo(5,250-Y[0]*scale);
Image1->Canvas->Pen->Color = clBlue;
for (int i=1; i<=n; i++) {
Image1->Canvas->LineTo(5+i*2,250-Y[i]*scale);
}
Image1->Canvas->MoveTo(5,250-X[0]*scale);
Image1->Canvas->Pen->Color = clBlack;
for (int i=1; i<=n; i++) {
Image1->Canvas->LineTo(5+i*2,250-X[i]*scale);
}
Image1->Canvas->MoveTo(5,250-M[0]*scale);
Image1->Canvas->Pen->Color = clRed;
for (int i=1; i<=n; i++) {
Image1->Canvas->LineTo(5+i*2,250-M[i]*scale);
}
}
//---------------------------------------------------------------------------
void __fastcall TForm1::btnShowClick(TObject *Sender)
{
CanvasClear();
// Вычисление масштаба шкалы значений ФБД
double maxb = 0;
for (int i=0; i<=n; i++) {
if (maxb < fabs(B[i])) {
maxb = fabs(B[i]);
}
}
int scale = floor(250/maxb);
if (scale > 249 ) {
scale = 249;
}
// Отрисовка осей координат
Image1->Canvas->Pen->Color = clSilver;
Image1->Canvas->MoveTo(5,250);
Image1->Canvas->LineTo(405,250);
Image1->Canvas->MoveTo(405,250);
Image1->Canvas->LineTo(405,255);
Image1->Canvas->MoveTo(5,250-scale);
Image1->Canvas->LineTo(5,250+scale);
Image1->Canvas->MoveTo(5,250-scale);
Image1->Canvas->LineTo(10,250-scale);
Image1->Canvas->MoveTo(5,250+scale);
Image1->Canvas->LineTo(10,250+scale);
Image1->Canvas->MoveTo(5,250);
Image1->Canvas->Pen->Color = clBlack;
// Построение графика ФБД
for (int i=0; i<=n; i++) {
Image1->Canvas->LineTo(5+i*2,250-B[i]*scale);
}
}
//---------------------------------------------------------------------------
void __fastcall TForm1::btnFBM2Click(TObject *Sender)
// Моделирование ФБД (альтернативный вариант)
{
Memo1->Clear();
CanvasClear();
double mu[N/2+1];
double Ak[N/2+1];
double Bk[N/2+1];
double beta[n+1];
AnsiString buf;
// Считывание значений Mu[k] для данного N и параметра Харста H
for (int i=0; i<=N/2; i++) {
buf = Memo2->Lines->Strings[i];
mu[i] = StrToFloat(buf);
}
// Генерирование случайных величин V[k]
randomize();
for (int i=0; i<=N/2; i++) {
Ak[i] = sqrt(2.0)*RandG(0,1)/sqrt(N);
Bk[i] = sqrt(2.0)*RandG(0,1)/sqrt(N);
}
// Вычисление реализаций ФГШ
for (int i=0; i<n; i++) {
beta[i]=0;
for (int j=1; j<=N/2; j++) {
beta[i] += cos(mu[j]*i)*Ak[j] + sin(mu[j]*i)*Bk[j];
}
}
// Вычисление значений ФБД
B[0] = 0;
Memo1->Lines->Add(AnsiString(B[0]));
for (int i=1; i<=n; i++) {
B[i] = B[i-1] + beta[i-1];
Memo1->Lines->Add(AnsiString(i)+" : "+AnsiString(B[i]));
}
//ScaleMotion();
// Вычисление масштаба шкалы значений ФБД
double maxb = 0;
for (int i=0; i<=n; i++) {
if (maxb < fabs(B[i])) {
maxb = fabs(B[i]);
}
}
int scale = floor(250/maxb);
if (scale > 249 ) {
scale = 249;
}
// Отрисовка осей координат
Image1->Canvas->Pen->Color = clSilver;
Image1->Canvas->MoveTo(5,250);
Image1->Canvas->LineTo(405,250);
Image1->Canvas->MoveTo(405,250);
Image1->Canvas->LineTo(405,255);
Image1->Canvas->MoveTo(5,250-scale);
Image1->Canvas->LineTo(5,250+scale);
Image1->Canvas->MoveTo(5,250-scale);
Image1->Canvas->LineTo(10,250-scale);
Image1->Canvas->MoveTo(5,250+scale);
Image1->Canvas->LineTo(10,250+scale);
Image1->Canvas->MoveTo(5,250);
Image1->Canvas->Pen->Color = clBlack;
// Построение графика ФБД
for (int i=0; i<=n; i++) {
Image1->Canvas->LineTo(5+i*2,250-B[i]*scale);
}
}
//---------------------------------------------------------------------------
void __fastcall TForm1::btnErrClick(TObject *Sender)
{
double maxb;
int scale;
Memo1->Clear();
CanvasClear();
if (ind2 == 0) {
double eps[n+1];
double maxe=0;
for (int i=0; i<=n; i++) {
eps[i] = X[i] - M[i];
if (i>10) {
if (fabs(eps[i]) > maxe) {
maxe = fabs(eps[i]);
}
}
}
// Вычисление масштаба шкалы значений функций
maxb = 0;
for (int i=0; i<=n; i++) {
if (maxb < fabs(eps[i])) {
maxb = fabs(eps[i]);
}
}
scale = floor(250/maxb);
if (scale > 249 ) {
scale = 249;
}
Image1->Canvas->Pen->Color = clSilver;
Image1->Canvas->MoveTo(5,250);
Image1->Canvas->LineTo(405,250);
Image1->Canvas->MoveTo(405,250);
Image1->Canvas->LineTo(405,255);
Image1->Canvas->MoveTo(5,250-scale);
Image1->Canvas->LineTo(5,250+scale);
Image1->Canvas->MoveTo(5,250-scale);
Image1->Canvas->LineTo(10,250-scale);
Image1->Canvas->MoveTo(5,250+scale);
Image1->Canvas->LineTo(10,250+scale);
Image1->Canvas->Pen->Color = clRed;
Image1->Canvas->MoveTo(5,250-maxe*scale);
Image1->Canvas->LineTo(405,250-maxe*scale);
Image1->Canvas->MoveTo(5,250+maxe*scale);
Image1->Canvas->LineTo(405,250+maxe*scale);
Image1->Canvas->MoveTo(5,250-eps[0]*scale);
Image1->Canvas->Pen->Color = clBlack;
for (int i=1; i<=n; i++) {
Image1->Canvas->LineTo(5+i*2,250-eps[i]*scale);
}
ind2 = 1;
ShowMessage("maxe = "+AnsiString(maxe));
}
else {
maxb = 0;
for (int i=0; i<=n; i++) {
if (maxb < fabs(X[i])) {
maxb = fabs(X[i]);
}
if (maxb < fabs(Y[i])) {
maxb = fabs(Y[i]);
}
if (maxb < fabs(M[i])) {
maxb = fabs(M[i]);
}
}
scale = floor(250/maxb);
if (scale > 249 ) {
scale = 249;
}
Image1->Canvas->Pen->Color = clSilver;
Image1->Canvas->MoveTo(5,250);
Image1->Canvas->LineTo(405,250);
Image1->Canvas->MoveTo(405,250);
Image1->Canvas->LineTo(405,255);
Image1->Canvas->MoveTo(5,250-scale);
Image1->Canvas->LineTo(5,250+scale);
Image1->Canvas->MoveTo(5,250-scale);
Image1->Canvas->LineTo(10,250-scale);
Image1->Canvas->MoveTo(5,250+scale);
Image1->Canvas->LineTo(10,250+scale);
Image1->Canvas->MoveTo(5,250-Y[0]*scale);
Image1->Canvas->Pen->Color = clBlue;
for (int i=1; i<=n; i++) {
Image1->Canvas->LineTo(5+i*2,250-Y[i]*scale);
}
Image1->Canvas->MoveTo(5,250-X[0]*scale);
Image1->Canvas->Pen->Color = clBlack;
for (int i=1; i<=n; i++) {
Image1->Canvas->LineTo(5+i*2,250-X[i]*scale);
}
Image1->Canvas->MoveTo(5,250-M[0]*scale);
Image1->Canvas->Pen->Color = clRed;
for (int i=1; i<=n; i++) {
Image1->Canvas->LineTo(5+i*2,250-M[i]*scale);
}
ind2=0;
}
}
//
Размещено на Allbest.ru
...Подобные документы
Байесовские алгоритмы оценивания (фильтр Калмана). Постановка задачи оценивания для линейных моделей динамической системы и измерений. Запись модели эволюции и модели измерения в матричном виде. Составление системы уравнений, описывающей эволюцию системы.
курсовая работа [3,0 M], добавлен 14.06.2011Оценка вероятности простоя цеха в виде схемы движения заявок или в виде соответствия "состояния системы"-"события". Выбор единицы моделирования и погрешности измеряемых параметров. Создание блок-схемы и листинга программы, отладка модели на языке GPSS.
лабораторная работа [213,6 K], добавлен 15.04.2012Принципы и этапы построения математической модели движения неуправляемого двухколесного велосипеда. Условия устойчивого движения. Вопрос гироскопической стабилизации движения. Модель движения велосипеда с гиростабилизатором в системе Matlab (simulink).
статья [924,5 K], добавлен 30.10.2015Анализ динамических процессов в системе на основе использования построенной аналитической модели. Моделирование с использованием пакета расширения Symbolic Math Tolbox. Построение модели в виде системы дифференциальных уравнений, записанных в форме Коши.
курсовая работа [863,4 K], добавлен 21.06.2015Моделирование непрерывной системы контроля на основе матричной модели объекта наблюдения. Нахождение передаточной функции формирующего фильтра входного процесса. Построение графика зависимости координаты и скорости от времени, фазовой траектории системы.
курсовая работа [1,5 M], добавлен 25.12.2013Вывод уравнения движения маятника. Кинетическая и потенциальная сила энергии. Определение всех положений равновесия. Исследование на устойчивость. Аналитический и численный расчет траектории системы. Изображение траектории системы разными способами.
контрольная работа [344,2 K], добавлен 12.04.2016Изучение актуальной задачи математического моделирования в биологии. Исследование модифицированной модели Лотки-Вольтерра типа конкуренция хищника за жертву. Проведение линеаризации исходной системы. Решение системы нелинейных дифференциальных уравнений.
контрольная работа [239,6 K], добавлен 20.04.2016Назначение и принципы действия корреляционно-экстремальной навигационной системы, особенности ее программно-аппаратной реализации, целесообразность статистического моделирования. Описание технологического процесса разработки и отладки программы.
магистерская работа [1,5 M], добавлен 06.12.2013Использование системы MathCAD как средства описания алгоритмов решения основных математических задач. Рассмотрение законов Кеплера и понятия о всемирном тяготении. Аналитические и численные решения задачи трех тел (материальных точек), вывод уравнений.
курсовая работа [287,2 K], добавлен 04.06.2013Построение квадратичной двумерной стационарной системы, нахождение состояний равновесия, исследование бесконечно-удаленной части плоскости. Необходимые и достаточные условия существования у системы двух частных интегралов. Построение траектории в круге.
дипломная работа [118,3 K], добавлен 07.09.2009Исследование зависимости погрешности решения от погрешностей правой части системы. Определение корня уравнения с заданной точностью. Вычисление точностных оценок методов по координатам. Сплайн интерполяция и решение дифференциального уравнения.
контрольная работа [323,4 K], добавлен 26.04.2011Построение многофакторной корреляционно-регрессионной модели доходности предприятия: оценка параметров функции регрессии, анализ факторов на управляемость, экономическая интерпретация модели. Прогнозирование доходности на основе временных рядов.
дипломная работа [5,1 M], добавлен 28.06.2011Теоретические основы моделирования: понятие модели и моделирования. Моделирование в решении текстовых задач. Задачи на встречное движение двух тел. Задачи на движение двух тел в одном направлении и в противоположных направлениях. Графические изображения.
курсовая работа [98,9 K], добавлен 03.07.2008Вычислительные методы линейной алгебры. Интерполяция функций. Интерполяционный многочлен Ньютона. Узлы интерполяции. Интерполяционный многочлен Лагранжа. Интерполяция сплайнами. Коэффициенты кубических сплайнов.
лабораторная работа [70,5 K], добавлен 06.02.2004Изучение основных принципов функционирования системы оптимального слежения. Моделирование привода антенны на основе экспериментальных данных, полученных при проведении исследований динамических характеристик и параметров привода РЛС в НПО "Горизонт".
дипломная работа [1,5 M], добавлен 24.11.2010Схема блоков модели Карааслана, система дифференциальных уравнений, методы решения. Блоки и биохимические законы системы Солодянникова, переход между фазами. Моделирование патологий, графики экспериментов. Построение комплексной модели гемодинамики.
дипломная работа [4,1 M], добавлен 24.09.2012Спектральная теория графов. Теоремы теории матриц и их применение к исследованию спектров графов. Определение и спектр предфрактального фрактального графов с затравкой регулярной степени. Связи между спектральными и структурными свойствами графов.
дипломная работа [272,5 K], добавлен 05.06.2014Задачи о пифагоровых треугольниках с целочисленными значениями сторон. Практическое использование задач в геодезии, в атомных и молекулярных структурах и в астрономических расчетах. Число вариантов представления исходного числа в виде двух сомножителей.
статья [29,9 K], добавлен 26.08.2013Математические модели технических объектов и методы для их реализации. Анализ электрических процессов в цепи второго порядка с использованием систем компьютерной математики MathCAD и Scilab. Математические модели и моделирование технического объекта.
курсовая работа [565,7 K], добавлен 08.03.2016Изучение понятия, классификации, свойств математических моделей. Особенности работы с функциями, переменными, графикой, программированием (интерполяция, регрессия) в системе MathCad. Проведение алгоритмического анализа задачи и аппроксимация результатов.
курсовая работа [4,5 M], добавлен 15.02.2010