Минимизация методом парабол. (перевести из c++ в матлаб)

@Mixirist · Регистрация: 27.06.2014

Author24 — интернет-сервис помощи студентам

Вот программа.

Кликните здесь для просмотра всего текста

C++

#include<iostream.h>
#include<math.h>
#include<conio.h>
#include<stdio.h>
#include "stdafx.h"
double f(double x)
{double f;
f=pow(x, 3)-x+exp(-x);
return f;
}
void main()
{double a=0, b=1, c, ya, yb, yc, s, t, yt, x, y;
const double e=0.00001;
int N;
cout<<"Method parabol "<<endl<<endl;
do
{cout<<"Vv. c ";
cin>>c;
}while ((a>=c) || (b<=c) || (f(a)<=f(c)) || f(b)<=f(c));
ya=f(a);
yb=f(b);
yc=f(c);
N=3;//kol-vo experimentov
while (b-a>2*e)
  {
   s=c+0.5*( (b-c)*(b-c)*(ya-yc) - (c-a)*(c-a)*(yb-yc) )/( (b-c)*(ya-yc) + (c-a)*(yb-yc) );
   if (s==c) t=(a+c)/2;   //see here
     else t=s;
   yt=f(t); N++;
 
   if (t<c)
     {
     if (yt<yc) {b=c; yb=yc; c=t; yc=yt;}
       else if (yt>yc) {a=t; ya=yt;}
        else {a=t; ya=yt; b=c; yb=yc; c=(a+b)/2; yc=f(c); N++;}
      }
   else if (t>c)
      {
       if (yt<yc) {a=c; ya=yc; c=t; yc=yt;}
          else if(yt>yc) {b=t; yb=yt;}
         else {a=c; ya=yc; b=t; yb=yt; c=(a+b)/2; yc=f(c); N++; }
         }
   }
x=(a+b)/2; y=f(x);
 
printf("\nx = %.5f \n", x);
printf("y = %.5f \n \n", y);
cout<<"Kol-vo tochek: "<<N;
}

Буду признателен если все получится.

Эх, забыли про мою темку)))

Зосима · 30.06.2014, 11:59

держи

Matlab M

clear, clc
 
f = @(x) x.^3 - x + exp(-x);
 
a=0;
b=1;
ezplot(f,[a, b]), hold on
 
e=0.00001;
disp('Метод парабол')
c = (a+b)/2;
% c = input('Введите с: ');
% while (a>=c) || (b<=c) || (f(a)<=f(c)) || (f(b)<=f(c))
%     c = input('Введите с: ');
% end
 
ya=f(a);
yb=f(b);
yc=f(c);
N=3; % kol-vo experimentov
while b-a>2*e
    s=c+0.5*( (b-c)*(b-c)*(ya-yc) - (c-a)*(c-a)*(yb-yc) )/( (b-c)*(ya-yc) + (c-a)*(yb-yc) );
   if (s==c) 
       t=(a+c)/2;   % see here
   else
       t=s;
   end
   yt=f(t); 
   N = N+1;
  if (t<c)
    if (yt<yc) 
        b=c; yb=yc; c=t; yc=yt;
    elseif (yt>yc) 
        a=t; ya=yt;
    else
        a=t; ya=yt; b=c; yb=yc; c=(a+b)/2; yc=f(c); 
        N = N+1;
    end
      
  elseif (t>c)
      if (yt<yc) 
          a=c; ya=yc; c=t; yc=yt;
      elseif (yt>yc) 
          b=t; yb=yt;
      else
          a=c; ya=yc; b=t; yb=yt; c=(a+b)/2; yc=f(c);
              N=N+1; 
      end
  end
end
x=(a+b)/2; y=f(x);
 
fprintf('\nx = %.5f \n', x);
fprintf('y = %.5f \n \n', y);
fprintf('Kol-vo tochek: %i \n', N);
 
plot(x,y,'or'), grid on

@Mixirist · 02.07.2014, 08:33 **[ТС]**

Спасибо)))

@Mixirist 0 / 0 / 0 Регистрация: 27.06.2014 Сообщений: 7
	02.07.2014, 08:33 [ТС]	3
	Спасибо))) 0

Минимизация методом парабол. (перевести из c++ в матлаб)

Решение