50628

Метод покоординатного спуска

Лабораторная работа

Математика и математический анализ

Напишем программу минимизации функции методом покоординатного спуска: Минимизируем исходную функцию стандартными средствами MatLab

Русский

2014-01-27

56.5 KB

72 чел.

Лабораторные работы

№8

Тема

Метод покоординатного спуска

Ф.И.О.

Пастухова Светлана Владимировна

Группа

403

Вариант

15

Минимизировать функцию f(x,y)=ax + by + exp(cx2 + dy2 )

 N

a

b

c

d

9

22

-0.2

1.44

0.23

Порядок выполнения работы:

  1.  Построим график функции:

ezsurf('13*x-0.2*y+exp(1.44*x^2+0.23*y^2)')

  1.  Минимизируем исходную функцию стандартными средствами MatLab:

x = [0,-6];

min = fminsearch(@my_fun,x)

где x=[0,-6] – начальное приближение, а @my_fun:

function f = my_fun(x)

f =13*x(1)-0.2*x(2)+exp(1.44*x(1)^2+0.23*x(2)^2); 

Получим:

min =[-1.0167, 0.0974]

f(xmin)=  -8.7965

  1.  Напишем программу минимизации функции методом покоординатного спуска:
  2.  %Минимизируем исходную функцию стандартными средствами MatLab:
  3.  x = [0,-6];
  4.  min = fminsearch(@my_fun,x)
  5.  %где
  6.  x=[0,-6] %– начальное приближение, а @my_fun:
  7.  function f = my_fun(x)
  8.  f =13*x(1)-0.2*x(2)+exp(1.44*x(1)^2+0.23*x(2)^2);
  9.  %Получим:
  10.  min =[-1.0167, 0.0974]
  11.  f(xmin)=  -8.7965
  12.  %Напишем программу минимизации функции методом покоординатного спуска:
  13.  aaa=13; bbb=-0.2; ccc=1.44; ddd=0.23;
  14.  x1=0;
  15.  x2=-6;
  16.  e=0.0001;
  17.  k=2;
  18.  h=0.2;
  19.  e1=e/k;
  20.  l=0;
  21.  d=e+1;
  22.  while(d>e)
  23.      d=abs(h);
  24.      h1=h;
  25.      z=aaa*x1+bbb*x2+exp(ccc*x1^2+ddd*x2^2);
  26.      z1=z+1;
  27.      d1=abs(h1);
  28.      while (d1>e)
  29.          d1=abs(h1);
  30.          x1=x1+h1;
  31.          z1=aaa*x1+bbb*x2+exp(ccc*x1^2+ddd*x2^2);
  32.          if (z1>=z) h1=-h1/k; end
  33.          z=z1;
  34.      end
  35.      h1=h;
  36.      z=aaa*x1+bbb*x2+exp(ccc*x1^2+ddd*x2^2);
  37.      d1=e1+1;
  38.      z1=z+1;
  39.      d1=abs(h1);
  40.      while (d1>e)
  41.          d1=abs(h1);
  42.          x2=x2+h1;
  43.          z1=aaa*x1+bbb*x2+exp(ccc*x1^2+ddd*x2^2);
  44.          if (z1>=z)  h1=-h1/k; end
  45.          z=z1;
  46.      end
  47.      h=h/k;
  48.      l=l+1;
  49.  end
  50.  x1
  51.  x2
  52.  aaa*x1+bbb*x2+exp(ccc*x1^2+ddd*x2^2)
  53.  e
  54.  %В результате получим:
  55.  min =[ -1.0165, 0.0981]
  56.  fmin = -8.7965
  57.  n=12
  58.   
  59.  

В результате получим:

min =[ -1.0165, 0.0981]

fmin = -8.7965

n=12


 

А также другие работы, которые могут Вас заинтересовать

49757. ЦИФРОВЫЕ СИСТЕМЫ ПЕРЕДАЧИ НЕПРЕРЫВНЫХ СООБЩЕНИЙ 374.56 KB
  Расчет ширины спектра сигнала модулированного двоичным кодом 10 1. Расчёт отношений мощностей сигнала и помехи необходимых для обеспечения заданного качества приёма 11 2.1 Формирование информационного сигнала 13 2.2 Формирование сигнала синхронизации 15 2.
49758. Выполнение проекта структурированной кабельной системы 344.5 KB
  Современный мир устроен так, что информационное взаимодействие, обмен информацией являются важнейшими компонентами, которые обеспечивают благополучие и развитие общества. Затраты на развитие и поддержание инфраструктуры такого взаимодействия весьма существенны и с целью снижения таких издержек пришли к пониманию необходимости комплексного решения задач информационного взаимодействия
49759. Разработка программы для имитационного моделирования системы массового обслуживания 815.12 KB
  Основные показатели: коэффициент использования системы средняя задержка в очереди среднее время ожидания среднее по времени число требований в очереди и в системе. Емкость накопителя требований r равна 14 дисциплина обслуживания – циклическая с квантом q = 10 секунд. В системе интервалы времени между поступлениями требований являются независимыми случайными величинами со средним временем = 60 секунд. Время обслуживания является случайной величиной некоррелированной с интервалами поступления требований.
49760. ПРОЕКТИРОВАНИЕ ПРИВОДА 591 KB
  Требуемая мощность кВт электродвигателя привода определяем по формуле: где Рв потребляемая мощность измельчителя Уточнение передаточных чисел привода Определяем общее передаточное отношение привода по формуле: Тогда Находим передаточное число редуктора: Тогда Принимаем Уточняем передаточное отношение открытой передачицепной: тогда SH – коэффициент запаса прочности принимаем в соответствии с рекомендациями с....
49761. Разработка базы данных подразделения учета основных средств предприятия 269.12 KB
  Целью данной курсовой работы является разработка базы данных подразделения учета основных средств предприятия. По итогам выполнения данной работы должна получиться информационная система, удовлетворяющая современным требованиям и упрощающая работу работников бухгалтерии по учету основных средств сельскохозяйственного предприятия.
49763. Автоматизация конструирования и технологической подготовки производства РЭС 263.93 KB
  На рисунках 1 и 2 представлены условнографическое отображение микроконтроллера PIC16C8404 P и чертеж его корпуса соответственно. Рисунок 1 – Условнографическое отображение микроконтроллера PIC16C8404 P. На рисунках 1 и 2 представлены условнографическое отображение микросхемы КМ155ИД11 и чертеж её корпуса соответственно. Рисунок 1 – Условнографическое отображение микросхемы КМ155ИД11.
49765. Экономическое обоснование проекта стратегического развития предприятия в сфере дорожного проектирования на примере предприятия ООО «АСК Драфт» 9.38 MB
  Проектирование дорог – это комплексный подход к решению транспортной проблемы в пределах целого региона или страны. Начинается он с этапа принятия решения о строительстве связи в сети дорог. Обязательной задачей этого и последующих этапов проектирования является анализ и сравнение различных вариантов.