№4 лабораторная МП
.docxФедеральное агентство связи Ордена
Трудового Красного Знамени
федеральное бюджетное образовательное учреждение высшего образования
«Московский Технический Университет Связи и Информатики»
(МТУСИ)
Кафедра Информатики
Предмет: Введение в математические пакеты прикладных программ
Лабораторная работа по ВвМППП № 4
Тема:
«Решение нелинейных уравнений»
Вариант 00
Выполнил:
Великий преподаватель
__________________________
Москва 0000
Общее задание
Изучите материал Темы 2.1 (п. 2.1.1).
Выберите индивидуальный вариант задания из табл. 2.1.2-1.
Отделите корень нелинейного уравнения f(x)=0 с использованием средств пакета Scilab, для чего:
Построить графики функции f(x) и ее первой производной;
на выбранном отрезке пересечения графика с осью ОХ получить таблицы значений аргумента, функции f(x)и ее первой производной
проверить условие существования единственного корня на выбранном отрезке.
Решите 1-е нелинейное уравнение с использованием функций fsolve(), получив значение корня и значение функции в точке корня.
Задайте вектор коэффициентов для 2-го уравнения.
Сформируйте с использованием функции poly() полином с коэффициентами, хранящимися в векторе.
Вычислите корни полученного полинома, используя функцию roots().
Предъявите результаты выполнение задания на ПК преподавателю.
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