Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- % После этого процента, компилятор не учитывает строрку
- % Задание 1:
- % Создаем симметричную матрицу А
- A = [1 2 3 4; 2 3 4 5; 3 4 5 6; 4 5 6 7]
- % Находим жорданову форму матрицы А
- [V, J] = eig(A)
- % Проверка ортогональности матрицы перехода
- J * J^(-1) % Должна получиться единичная матрица
- [U, S, V] = svd(A) % Сингулярное разложение матрицы
- % Чтобы убедиться, что жорданова форма совпадает с сингулярным выражением
- % Нужно чтобы U было равно V
- %%!!! Как создать М-файл прочитаешь тут: http://life-prog.ru/view_zam.php?id=57
- %%!!! Буквально первые два абзаца
- % Задание 2:
- % Создаем полином третьей степени, например p(x)= x^3 + 1
- % Откуда начинаем: шаг: где заканчиваем
- t = -5:0.1:5
- x = t.^3 + 1 % точка нужна, хз зачем. Вроде бы из-за того, что мы не матрицу возводим в степень
- % Строим наш график
- plot (t, x) % Появляется график
- % Чтобы найти корни полинома, представляем его в другом виде
- p = [3, 0, 0, 1] % Коэффициенты при иксах
- % Находим корни полинома
- roots (p)
- % Создаем вторую функцию
- y = t.^2
- % Строим две функции на одном графике
- plot (t, x, t, y)
- clear % Отчищаем прошлые данные
- % Хуярим спираль
- a=1;
- t=0:0.2:20;
- r=a*t; % Радиус
- x = r.*cos(t);
- y = r.*sin(t);
- z = t;
- plot(x,y);axis equal; % Чтобы график красивый был, а не растянутый. Не нужно переносить строку
- [x1, y1]= cart2pol(x,y); % Переводим в полярную систему координат
- polar(x1,y1);
- % Задание 3:
Add Comment
Please, Sign In to add comment