МІНІСТЕРСТВО ОСВІТИ І НАУКИ УКРАЇНИ
НАЦІОНАЛЬНИЙ УНІВЕРСИТЕТ «ЛЬВІВСЬКА ПОЛІТЕХНІКА»
ІКТА
Кафедра БІT
/
Звіт
до лабораторної роботи №1
з курсу: «Інформаційні технології»
на тему: «Виконання обчислювальних операцій в середовищі Matlab»
Варіант-24
МЕТА РОБОТИ – вивчити засоби для моделювання та виконання обчислювальних операцій в середовищі МАТЛАБ при використанні інформаційних технологій.
Завдання до виконання лабораторної роботи
У всіх завданнях k – номер варіанту (наданий викладачем або остання цифра номера залікової книжки студента).
Завдання 1. Знайти , якщо , ,
Завдання 2. Обчислити скалярний добуток з використанням вхідних даних завдання 1. Перевірити вектори і на ортогональність.
Завдання 3. Знайти детермінант матриці за правилом трикутників.
Завдання 4. Перевірити правильність знаходження , обчисливши визначник як суму добутків елементів:
а) другого стовпця на відповідні алгебраїчні доповнення;
б) третього рядка на алгебраїчні доповнення.
Завдання 5. Визначити ранг матриці із завдання 3.
Завдання 6. Нарисувати графік функції , обчислити першу та другу похідну ( та ), а також неозначений і означений інтеграли і , якщо , , .
Хід роботи
Код програми
%Task_1_
fprintf('Task1\n\n');
k=24;
a=[3.8 (-5+k) (k^2)+8 6-k];
b=[-0.5 -4.5 (-9*k) (k^3)];
Lam=7*k+3;
Mu=10-k^(-4);
c=Lam*a+Mu*b;
fprintf('b=%.lf; %.lf; %.lf; %.lf;)\n',b);
fprintf('Lam=%.lf\n',Lam);
fprintf('Mu=%3f\n',Mu);
fprintf('c=(%3f; %3f; %3f; %3f)\n',c);
%Task2
fprintf('Task2\n\n');
s=sum((Mu * a).* ((k+8)* b));
fprintf('s=%.3f\n',s);
if s==0
fprintf('Ortogonal \n');
else fprintf('Not ortogonal\n');
end;
fprintf('\n \n \n \n \n \n');
%Task3
fprintf( 'Task3\n\n');
A = [8 -k 15; (k - 5) 1.2 (k - 9.5); k -3.1 (k + 1.5)];
DetA=det(A)
DetA = A(1,1)* A(2,2)* A(3,3) + A(1,2) * A(2,3) * A(3,1) + A(2,1) * A(3,2)* A(1,3) - A(3,1) * A(2,2) * A(1,3) - A(2,1) * A(1,2) * A(3,3) - A(1,1) *A(3,2) * A(2,3);
%Task4
fprintf( '4) Perevirka\n\n');
DetA2 = -1^(1+2) * A(1,2) * (A(2,1) * A(3,3) - A(3,1) * A(2,3)) + (-1)^(2+2) * A(2,2) * (A(1,1) * A(3,3) - A(3,1) * A(1,3)) + (-1)^(3+2) * A(3,2) * (A(1,1) * A(2,3) - A(2,1) * A(1,3));
DetA2
DetA3= A(3,1)*(A(1,2)*A(2,3)-A(2,2)*A(1,3))*(-1)^(3+1)+A(3,2)*(A(1,1)*A(2,3)-A(2,1)*A(1,3))*(-1)^(3+2)+A(3,3)*(A(1,1)*A(2,2)-A(2,1)*A(1,2))*(-1)^(3+3);
DetA3
fprintf( '\n\n');
%Task5
fprintf( 'Task5:Rang matrici\n\n');
rang=rank(A)
%Task6
fprintf( 'Task6\n\n');
syms x; %simvolni vichisleniya
f=k*x.^2+(-1)^k;
fprintf('First derivative f(x)=%s \n',char(diff(f)));
fprintf('Second derivative f(x)=% s \n',char(diff(diff(f))));
%Graphic
x=-4*k:0.1:4*k;
fx=k*x.^2+(-1).^k;
y=abs(fx+k-5);
plot(x,y,'r-')
%Integral
syms x;
a=0;b=4*k;
fx=k*x.^2+(-1).^k;
I=int(fx)
Iab=int(fx,a,b)
Результат виконання програми
Task1
b=Lam=Mu=9.999997
c=(644.800002; 3204.000014; 97704.000651; 135161.958333)
Task2
s=-120020251.825
Not ortogonal
Task3
DetA = 2.5649e+003
4) Perevirka
DetA2 = 2.5649e+003
DetA3 = 2.5649e+003
Task5:Rang matrici
rang = 3
Task6
First derivative f(x)=48*x
Second derivative f(x)=48
I = 8*x^3 + x
Iab = 7077984
/
Рис. Графік функції у=