Как решить нелинейную систему уравнений

@Matlabist · Регистрация: 28.11.2010

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

Здравствуйте,
есть система уравнений, x вектор

f1(x)=0
...
fn(x)=0

Якобиан можно посчитать.
Как записать это в матлабе, чтоб можно было создать еще dll и использовать полученную функцию в visual studio c++?

@Галина Борисовн · 16.01.2011, 10:43

Посмотрите здесь

http://matlab.exponenta.ru/opt... 0/1_11.php
http://matlab.exponenta.ru/opt... 10/2_2.php
http://matlab.exponenta.ru/opt... 10/2_1.php

@Matlabist · 16.01.2011, 13:09 **[ТС]**

Взял пример с книги

Matlab M

1
2
3

function у = funsc(х)
y=[x(l)+x(2)-sin(pi*x(l)); х(1)-х(2)-cos(pi*x(2))];
[х,f] = fsolve(@funsc,[0.2;l],optimset('Display','off))

Пишет

??? Error: File: Untitled2.m Line: 1 Column: 10
The input character is not valid in MATLAB statements or
expressions.

Добавлено через 17 минут
Заменил буквы на латинские,

Matlab M

1
2
3

function y = funsc(x)
y=[x(1)+x(2)-sin(pi*x(1)); x(1)-x(2)-cos(pi*x(2))];
[x,f] = fsolve(@funsc,[0.2;l],optimset('Display','off'))

Теперь пишет
??? Input argument "x" is undefined.

Error in ==> Untitled2 at 2
y=[x(1)+x(2)-sin(pi*x(1)); x(1)-x(2)-cos(pi*x(2))];

@Sergik1 · 16.01.2011, 15:15

Сохраните первые 2 строки как отдельный файл. Третью строку запишите в новом чистом файле.
Вы записали в первых строках ФАЙЛ-функцию. В последней строке ссылаетесь на ФАЙЛ.
Оба файла должны лежать в одной папке.

@Галина Борисовн · 16.01.2011, 15:21

Matlab M

1
2
3

% ПП-функция. Отдельный m-файл
function y = funsc(x)
y=[x(1)+x(2)-sin(pi*x(1)); x(1)-x(2)-cos(pi*x(2))];

Основная(вызывающая) программа

Matlab M

clc
 clear
x0 = [-5; -4];
options=optimset('Display','iter');
[x,fval] = fsolve(@funsc,x0,options)

Оба m-файла сохранить в одной папке.

@Matlabist · 21.01.2011, 20:00 **[ТС]**

А как можно создать dll из полученной программы? Так что можно было использовать полученную функцию. Т.е. задаю параметры функции например a,b f(x)=a*x+ln(x)+b. А функция дает решение, вызывается из программы на c++.

@Галина Борисовн · 21.01.2011, 20:07

Посмотрите здесь
http://matlab.exponenta.ru/for... 382&nt=All

@Matlabist · 22.01.2011, 08:34 **[ТС]**

Я все правильно сделал?
f(x)=ln((V-Vn1)/(V*m+Jtc/(ron*exp(-wtc*(pn-p0)))))/wtc-ln(Vn1/(V-V*m))/wn
функция должна находить нули f(x) вблизи точки Vn0.

Matlab M

function [ Vn1_out ] = fun2( wtc,V,m,Jtc,ron,pn,p0,wn,Vn0 )
Vn1=Vn0;
Vn1=fsolve(ln((V-Vn1)/(V*m+Jtc/(ron*exp(-wtc*(pn-p0)))))/wtc-ln(Vn1/(V-V*m))/wn);
Vn1_out=Vn1;
end

Добавлено через 10 часов 36 минут
Пытаюсь сделать одну функцию которая решает нелинейное уравнение, чтобы сделать одну функцию dll.

Matlab M

clc
 clear
x0 = 10;
%options=optimset('Display','iter','PlotFcns',@optimplotstepsize);
%[x,fval] = fsolve(@fun,x0,options)
wtc=1e-9;
V=2e4;
m=0.5;
ron=1e3;
pn=1e6;
p0=1e6;
wn=0.5e-9;
Vn0=m*V;
Jtc=1000;
Vn1=fun2(wtc,V,m,Jtc,ron,pn,p0,wn,Vn0)

Matlab M

1
2
3

function [ Vn1_out ] = fun2( wtc,V,m,Jtc,ron,pn,p0,wn,Vn0 )
Vn1_out=fsolve(log((V-Vn1)/(V*m+Jtc/(ron*exp(-wtc*(pn-p0)))))/wtc-log(Vn1/(V-V*m))/wn,Vn0);
end

Пишет

Код

??? Undefined function or variable 'Vn1'.

Error in ==> fun2 at 2
Vn1_out=fsolve(log((V-Vn1)/(V*m+Jtc/(ron*exp(-wtc*(pn-p0)))))/wtc-log(Vn1/(V-V*m))/wn,Vn0);

Error in ==> main at 15
Vn1=fun2(wtc,V,m,Jtc,ron,pn,p0,wn,Vn0)

@Matlabist 2 / 2 / 1 Регистрация: 28.11.2010 Сообщений: 102
		1
	Как решить нелинейную систему уравнений 16.01.2011, 10:00. Показов 5973. Ответов 7 Метки нет (Все метки) Здравствуйте, есть система уравнений, x вектор f1(x)=0 ... fn(x)=0 Якобиан можно посчитать. Как записать это в матлабе, чтоб можно было создать еще dll и использовать полученную функцию в visual studio c++? 0

@Галина Борисовн 2834 / 2131 / 86 Регистрация: 02.05.2010 Сообщений: 3,195
	16.01.2011, 10:43	2
	Посмотрите здесь http://matlab.exponenta.ru/opt... 0/1_11.php http://matlab.exponenta.ru/opt... 10/2_2.php http://matlab.exponenta.ru/opt... 10/2_1.php 0

@Sergik1 128 / 127 / 10 Регистрация: 09.11.2010 Сообщений: 200
	16.01.2011, 15:15	4
	Сохраните первые 2 строки как отдельный файл. Третью строку запишите в новом чистом файле. Вы записали в первых строках ФАЙЛ-функцию. В последней строке ссылаетесь на ФАЙЛ. Оба файла должны лежать в одной папке. 1

@Matlabist 2 / 2 / 1 Регистрация: 28.11.2010 Сообщений: 102
	21.01.2011, 20:00 [ТС]	6
	А как можно создать dll из полученной программы? Так что можно было использовать полученную функцию. Т.е. задаю параметры функции например a,b f(x)=a*x+ln(x)+b. А функция дает решение, вызывается из программы на c++. 0

@Галина Борисовн 2834 / 2131 / 86 Регистрация: 02.05.2010 Сообщений: 3,195
	21.01.2011, 20:07	7
	Посмотрите здесь http://matlab.exponenta.ru/for... 382&nt=All 1

	22.01.2011, 08:34