Решение нелинейного уравнения методом бисекции

[Pavel.V. 2]

 

 Доброго времени суток всем участникам форума. За неимением опыта работы со скриптами на apdl в Ansys обращаюсь за помощью к знатокам. Задача в следующем: необходимо написать макрос решения нелинейного уравнения с помощью метода бисекции. Суть метода известна, для меня вопрос каким образом это записать в скрипте? Само уравнение: X^3-cos(x)+1=0 взято из примера решения в mathcad для сравнения результатов. Заданная точность  Е=10^-10,  интервал поиска [-0,6 ; -0,4].  Заранее благодарен за помощь!

[soklakov 1 475]

1 час назад, Pavel.V. сказал:

Суть метода известна, для меня вопрос каким образом это записать в скрипте?

APDL - это язык из последовательных команд примерно одинакового синтаксиса:

NAME,ARG1,ARG2,ARG3

где NAME  - имя команды, а дальше идут ее аргументы, которые могут быть условно произвольными, в зависимости от команды.

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

*IF

*ELSE

*DO

*ENDDO

*WHILE

звездочка (*) здесь - это часть имени команды.

Подробности по командам есть в полной мере в справке.

Изменено 19 декабря 2017 пользователем soklakov

[Борман 2 375]

@soklakov

Успехи налицо.

 

Файл с функцией FNC.txt

!=================================

FNC=ARG**3-cos(ARG)+1

!=================================

 

Файл с солвером. Результат содержится в переменной SOL если не достигнуто ограничение по итерациям.

!=================================

*afun,rad

x1=-0.6
x2=-0.4

MaxIter=100

*do,i,1,MaxIter,1

     *if,abs(x1-x2),lt,1e-10,then
         SOL=(x1+x2)/2
         *exit
     *endif

     dx=(x2-x1)/2
     xi=x1+dx

     ARG=x1
     /INPUT,FNC.txt
     RES1=FNC

     ARG=xi
     /INPUT,FNC.txt
     RESi=FNC

     *if,RES1*RESi,lt,0,then
         x2=xi
     *else
         x1=xi
     *endif

*enddo

ARG=
DX=
FNC=
MAXITER=
RES1=
RESi=
Xi=

!=================================

 

[Pavel.V. 2]

 

@soklakov Спасибо за помощь! 

[soklakov 1 475]

1 час назад, Борман сказал:

Успехи налицо.

Успехи? В чем?

[Pavel.V. 2]

Борман, за предоставленный пример скрипта отдельное спасибо! 

[Борман 2 375]

2 часа назад, soklakov сказал:

Успехи? В чем?

Твои в механикл "мать его" апдле. А ты про что подумал?

[piden 2 874]

- Как выучить основные команды на апдл?

- Попить пиваса с Борманом!

 

 

- А что будет, если с Боманом дунуть? 

 

[soklakov 1 475]

15 часов назад, Борман сказал:

Твои в механикл "мать его" апдле. А ты про что подумал?

Про то и подумал. Только думаю... успехи мои в апдле за последние десять лет мало изменились, как и сам апдл. А циклы и условные переходы писать еще в школе учат.

[Radon 14]

Только следует знать, что для расчётов apdl оооочень медленный, медленнее того же визуал бейсика на порядки (даже не в разы), так что решение примерно 100000 уравнений в нём занимает часы. Так что может придётся выносить код при больших вычислениях наружу, на c++ или каком другом быстром языке писать и организовывать обмен информацией

[Fedor 1 597]

Скриптовые языки не для того придумывают. Еще бы на змеюке скорость искали :)