Замена одной системы уравнений в программе на другую

@Kinroom · Регистрация: 13.12.2012

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

Здравствуйте. Есть готовая программа, которая работает с численным методом. Она рабочая, нужно только заменить ту систему что в программе, на мою. Система, которую нужно вставить в программу:
4*x1+11*x2^2=9*x3^3=0
14*x1^3+3*x2+28*x3^2=0
11*x1^2+7*x2^3+3*x3=0
Код:

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
F0 = inline('sin(x+y)-1.6*x');
F1 = inline('x^2+y^2-1');
 
disp(F0); disp(F1);
% Их производны
dF0x = inline('cos(x+y)-1.6');
dF0y = inline('cos(x+y)');
dF1x = inline('2*x+0*y');
dF1y = inline('0*x+2*y');
disp(dF0x); disp(dF0y);
disp(dF1x); disp(dF1y);
% Начальное приближение
x=0; y=-1; e=0.000000001;
% Значение функций
root = newton(x,y,e,F0,F1,dF0x,dF0y,dF1x,dF1y);
disp('Решение системы');
disp(root);
%После выполнения программы получим следующее:
%>>
%Inline function:
%F0(x,y) = sin(x+y)-1.6*x
%Inline function:
%F1(x,y) = x^2+y^2-1
%Inline function:
%dF0x(x,y) = cos(x+y)-1.6
%Inline function:
%dF0y(x,y) = cos(x+y)
%Inline function:
%dF1x(x,y) = 2*x+0*y
%Inline function:
%dF1y(x,y) = 0*x+2*y
%Количество итераций
%5
%Решение системы
%-0.6163
%-0.7875

Заранее благодарен)

Зосима · 18.11.2014, 11:06

У тебя система с тремя переменными! рискну предположить:

Matlab M
F1 = inline('4*x+11*y^2-9*z^3');
F2 = inline('14*x^3+3*y+28*z^2');
F3 = inline('11*x^2+7*y^3+3*z');

и производные:

Matlab M
dF1x = inline('4+0*x+0*y+0*z');
dF1y = inline('0*x+22*y+0*z');
dF1z = inline('0*x+0*y-27*z^2');
 
dF2x = inline('42*x^2+0*y+0*z');
dF2y = inline('3+0*x+0*y+0*z');
dF2z = inline('0*x+0*y+56*z');
 
dF3x = inline('22*x+0*y+0*z');
dF3y = inline('0*x+21*y^2+0*z');
dF3z = inline('3+0*x+0*y+0*z');

аналогично и в начальных приближениях нужно еще задать какое-то z
Но самое неприятное это то, что придется менять функцию newton

@Kinroom · 18.11.2014, 14:04 **[ТС]**

Зосима, спасибо, да там 3 переменные) С функцией я разберусь думаю) С математикой все туго просто)

Новые блоги и статьи Все статьи Все блоги /
делаю науч статью по влиянию грибов на сукцессию anaschu 13.03.2026 прикрепляю статью	SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .
Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .	Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .

Зосима 5245 / 3573 / 379 Регистрация: 02.04.2012 Сообщений: 6,477 Записей в блоге: 18
	18.11.2014, 11:06
	У тебя система с тремя переменными! рискну предположить: Matlab M F1 = inline('4x+11y^2-9z^3'); F2 = inline('14x^3+3y+28z^2'); F3 = inline('11x^2+7y^3+3z'); и производные: Matlab M dF1x = inline('4+0x+0y+0z'); dF1y = inline('0x+22y+0z'); dF1z = inline('0x+0y-27z^2'); dF2x = inline('42x^2+0y+0z'); dF2y = inline('3+0x+0y+0z'); dF2z = inline('0x+0y+56z'); dF3x = inline('22x+0y+0z'); dF3y = inline('0x+21y^2+0z'); dF3z = inline('3+0x+0y+0z'); аналогично и в начальных приближениях нужно еще задать какое-то z Но самое неприятное это то, что придется менять функцию newton 1

@Kinroom 1 / 1 / 0 Регистрация: 13.12.2012 Сообщений: 99
	18.11.2014, 14:04 [ТС]
	Зосима, спасибо, да там 3 переменные) С функцией я разберусь думаю) С математикой все туго просто) 0