Фриланс-проєкти

Фриланс-проєкти

Переробити код Матлаб

Десктопні додатки — неправильно зазначені категорії?

В мене є код для ітераційного процесу синхроної апроксимаціїї, потрібно змінити його для асинхронної

ось це код для цієї формули

function rez = cursova()

format longG;

clear;

e = 1e-6;

%----------------------------------------------------------------------------------

x_k = [0;0.5];

disp('---------Метод з синхронною апроксимацією--------')

count = 0;

F = [f1dx1(x_k(1),x_k(2),1) f1dx2(x_k(1),x_k(2),1);

f2dx1(x_k(1),x_k(2),1) f2dx2(x_k(1),x_k(2),1)];

A=inv(F);

x_next=x_k-(A*[f1(x_k(1),x_k(2), 1); f2(x_k(1),x_k(2), 1)]);

A_next=A*([2*[1 0;0 1]-[f1dx1(x_k(1),x_k(2),1) f1dx2(x_k(1),x_k(2),1);

f2dx1(x_k(1),x_k(2),1) f2dx2(x_k(1),x_k(2),1)]*A]);

count = count+1;

disp(x_k)

disp(abs(f1(x_next(1),x_next(2),1)))

disp(abs(f2(x_next(1),x_next(2),1)))

while (norm(x_next-x_k,"inf") >=e)

A=A_next;

x_k=x_next;

x_next=x_k-(A*[f1(x_k(1),x_k(2), 1); f2(x_k(1),x_k(2), 1)]);

A_next=A*([2*[1 0; 0 1]-[f1dx1(x_k(1),x_k(2),1) f1dx2(x_k(1),x_k(2),1);

f2dx1(x_k(1),x_k(2),1) f2dx2(x_k(1),x_k(2),1)]*A]);

count = count+1;

disp(x_k)

disp(abs(f1(x_next(1),x_next(2),1)))

disp(abs(f2(x_next(1),x_next(2),1)))

end

disp('К-ть ітерацій : ')

disp(count)

disp('Приблизний результат')

disp(x_k)

disp('f1(a)')

f1_x=f1(x_next(1),x_next(2),1);

disp(abs(f1_x))

disp('f2(a)')

f2_x=f2(x_next(1),x_next(2),1);

disp(abs(f2_x))

%-----------------------------------------------------------------------------

% перше р-ня системи

function rez = f1(x, y, fun)

if fun == 1

%rez = cos(1/3*(x-y))-2*y;

rez = 2*x-sin(0.5*x-0.5*y);

else

%rez = cos(1/3*(x-y))-2*y;

rez = 2*x-sin(0.5*x-0.5*y);

end

% друге р-ня системи

function rez = f2(x, y, fun)

if fun == 1

%rez = sin(1/3*(x+y))-2*x;

rez = 2*y-cos(0.5*x+0.5*y);

else

% rez = sin(1/3*(x+y))-2*x;

rez = 2*y-cos(0.5*x+0.5*y);

end

% похідна 1го р-ня по X

function rez =f1dx1(x, y, fun)

if fun == 1

%rez = -1/3*sin((x-y)/3);

rez = 2-(0.5*cos(0.5*x-0.5*y));

else

% rez = -1/3*sin((x-y)/3);

rez = 2-(0.5*cos(0.5*x-0.5*y));

end

% похідна 1го р-ня по Y

function rez =f1dx2(x, y, fun)

if fun == 1

% rez = 1/3*(sin((x-y)/3)-6);

rez = 0.5*cos(0.5*x-0.5*y);

else

% rez = 1/3*(sin((x-y)/3)-6);

rez = 0.5*cos(0.5*x-0.5*y);

end

% похідна 2го р-ня по X

function rez =f2dx1(x, y, fun)

if fun == 1

% rez = 1/3*(cos((x+y)/3)-6);

rez = 0.5*sin(0.5*x+0.5*y);

else

% rez = 1/3*(cos((x+y)/3)-6);

rez = 0.5*sin(0.5*x+0.5*y);

end

% похідна 2го р-ня по Y

function rez =f2dx2(x, y, fun)

if fun == 1

% rez = 1/3*cos((x+y)/3);

rez = 0.5*sin(0.5*x+0.5*y)+2;

else

%rez = 1/3*cos((x+y)/3);

rez = 0.5*sin(0.5*x+0.5*y)+2;

end

ПОТРВБНО ПЕРЕПИСАТИ ДЛЯ ЦІЄЇ

Ставки 1

Khrystyna Stadnyk

5 0

Проєкти 5
Оцінка -
Рейтинг 452

Бюджет: 500 UAH Термін: 1 день

Доброго дня, можу зробити ваше завдання протягом дня. Звертайтесь, буду рада співпраці.

Віка Булишин
Львів, Україна

Проєктів 3
Оцінка -
Рейтинг 57

Переробити код Матлаб

Khrystyna Stadnyk

Ставки поки відсутні

Актуальні фриланс-проєкти в категорії Десктопні додатки

Microsoft Power Platform/Power App

Потрібен спеціаліст з розробки інтернет-магазину на Shopify з інтеграцією Amazon

Потрібен Airtable-архітектор для побудови relational schema та нового clean Airtable base

Невелика правка на сайті