Метод квадратичной интерполяции

@buritoz · Регистрация: 13.12.2010

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

Здравствуйте, если есть у кого нибуть данный метод, поделитесь пожалуйста.

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

.

C

x1=1;
dx=0.1;
e1=0.01;
e2=0.01;
while 1
x2=x1+dx;
f1=f(x1);
f2=f(x2);
if (f1>f2)
    x3=x1+(2*dx)
else
    x3=x1-dx
end
f3=f(x3);
fmin=min(f1,f2,f3);
xmin=x1;
fmin=f(x1);
xm=1/2*(((x2^2-x3^2)*f1+(x3^2+x1^2)*f2+(x1^2-x2^2)*f3) /((x2-x3)*f1+(x3-x1)*f2+(x1-x2)*f3));
fm=f(xm);
 
if abs(((fmin-fm)/fm)) < e1 && abs(((xmin-xm)/xm)<e2)
    xkon=xm;
    break;
else
    if x > x1 && x < x3
        x=xmin;
    else
        x1=xm;
    end
end
end

Выручайте

@Mary Fry · 28.12.2014, 21:20

Matlab M

function[f] = laba4 (A1,h,E)
i=0;
q=0;
A2=A1+h;
if fun1(A1)<fun1(A2)
A3=A1-h;
else
A3=A1+2*h;
end
while(fun1(A1)-fun1(q))>E
    A =(A1-A2)*(A2-A3)*(A3-A1);
a=((A3-A2)*fun1(A1)+(A1-A3)*fun1(A2)+(A2-A1)*fun1(A3))/A; 
b=((A2^2-A3^2)*fun1(A1)+(A3^2-A1^2)*fun1(A2)+(A1^2-A2^2)*fun1(A3))/A ;
q=(-b/2*a);
if fun1(A1)<fun1(A2) && (A3<q)&&(q<A1) 
    A2=A1; A1=q;
elseif fun1(A1)<fun1(A2) && (A1<q)&&(q<A2)
    A3=A1; A1=q;
elseif fun1(A2)<fun1(A1) && (A2<q)&&(q<A3)
    A3=A2;A2=q;
else 
    A1=A2; A2=q;
end
end 
f=q;
disp(f);
disp(i);
end

@Nick07 · 29.12.2014, 11:15

Напишите, какие ошибки выдает MATLAB и приведите Ваше толкование этих сообщений

@Nick07 533 / 438 / 47 Регистрация: 17.07.2013 Сообщений: 2,236
	29.12.2014, 11:15	3
	Напишите, какие ошибки выдает MATLAB и приведите Ваше толкование этих сообщений 0

Опции темы