№4 лабораторная ВВПП

Федеральное агентство связи Ордена

Трудового Красного Знамени

федеральное бюджетное образовательное учреждение высшего образования

«Московский Технический Университет Связи и Информатики»

(МТУСИ)

Кафедра Информатики

Предмет: Введение в математические пакеты прикладных программ

Лабораторная работа по ВвМППП № 4

Тема:

«Решение нелинейных уравнений»

Вариант 00

Выполнил:

Великий преподаватель

__________________________

Москва 0000

Общее задание

1. Изучите материал Темы 2.1 (п. 2.1.1).

2. Выберите индивидуальный вариант задания из табл. 2.1.2-1.

3. Отделите корень нелинейного уравнения f(x)=0 с использованием средств пакета Scilab, для чего:

• Построить графики функции f(x) и ее первой производной;

• на выбранном отрезке пересечения графика с осью ОХ получить таблицы значений аргумента, функции f(x)и ее первой производной

• проверить условие существования единственного корня на выбранном отрезке.

4. Решите 1-е нелинейное уравнение с использованием функций fsolve(), получив значение корня и значение функции в точке корня.

5. Задайте вектор коэффициентов для 2-го уравнения.

6. Сформируйте с использованием функции poly() полином с коэффициентами, хранящимися в векторе.

7. Вычислите корни полученного полинома, используя функцию roots().

8. Предъявите результаты выполнение задания на ПК преподавателю.

x=2:0.1:4;

deff('[y]=f(x)','y= x.^2-3.4*log(1+x)-1');

deff('[y1]=f1(x)','y1=2*x-3.4./(1+x)');

fx=f(x);f1x=f1(x);

Z=[x;fx;f1x]'

Z =

2. -0.7352818 2.8666667

2.1 -0.4367672 3.1032258

2.2 -0.1147128 3.3375

2.3 0.2306636 3.569697

2.4 0.5991635 3.8

2.5 0.9906059 4.0285714

2.6 1.4048249 4.2555556

2.7 1.8416684 4.4810811

2.8 2.3009964 4.7052632

2.9 2.7826797 4.9282051

3. 3.2865992 5.15

3.1 3.8126443 5.3707317

3.2 4.3607126 5.5904762

3.3 4.9307089 5.8093023

3.4 5.5225446 6.0272727

3.5 6.1361369 6.2444444

3.6 6.7714086 6.4608696

3.7 7.4282875 6.6765957

3.8 8.1067059 6.8916667

3.9 8.8066003 7.1061224

4. 9.5279111 7.32

z=f(x);

plot(x,z,x,0*z)

--> xtitle('Графическая функция f(x)','x','f(x)');

scf(1);

z=f1(x);

plot(x,z);

xtitle('Графическая производной f(x)','x','f1(x)');

t=fsolve(2,f)

t =

2.2339682

f(t)

ans =

8.882D-16

a=[6 0 -18 0 1];

p=poly(a,'x','c')

p =

2 4

6 -18x +x

X=roots(a)

X =

1.7156269

-1.7156269

-0.2379587

0.2379587