Gistrec

MatLab лабораторная 2

Mar 8th, 2017
198
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
MatLab 2.14 KB | None | 0 0
  1. % После этого процента, компилятор не учитывает строрку
  2.  
  3. % Задание 1:
  4. % Создаем симметричную матрицу А
  5. A = [1 2 3 4; 2 3 4 5; 3 4 5 6; 4 5 6 7]
  6. % Находим жорданову форму матрицы А
  7. [V, J] = eig(A)
  8. % Проверка ортогональности матрицы перехода
  9. J * J^(-1) % Должна получиться единичная матрица
  10. [U, S, V] = svd(A) % Сингулярное разложение матрицы
  11. % Чтобы убедиться, что жорданова форма совпадает с сингулярным выражением
  12. % Нужно чтобы U было равно V
  13. %%!!! Как создать М-файл прочитаешь тут: http://life-prog.ru/view_zam.php?id=57
  14. %%!!! Буквально первые два абзаца
  15.  
  16. % Задание 2:
  17. % Создаем полином третьей степени, например p(x)= x^3 + 1
  18. % Откуда начинаем: шаг: где заканчиваем
  19. t = -5:0.1:5
  20. x = t.^3 + 1 % точка нужна, хз зачем. Вроде бы из-за того, что мы не матрицу возводим в степень
  21. % Строим наш график
  22. plot (t, x) % Появляется график
  23. % Чтобы найти корни полинома, представляем его в другом виде
  24. p = [3, 0, 0, 1] % Коэффициенты при иксах
  25. % Находим корни полинома
  26. roots (p)
  27. % Создаем вторую функцию
  28. y = t.^2
  29. % Строим две функции на одном графике
  30. plot (t, x, t, y)
  31. clear % Отчищаем прошлые данные
  32. % Хуярим спираль
  33. a=1;
  34. t=0:0.2:20;
  35. r=a*t; % Радиус
  36. x = r.*cos(t);
  37. y = r.*sin(t);
  38. z = t;
  39. plot(x,y);axis equal; % Чтобы график красивый был, а не растянутый. Не нужно переносить строку
  40. [x1, y1]= cart2pol(x,y); % Переводим в полярную систему координат
  41. polar(x1,y1);
  42.  
  43. % Задание 3:
Add Comment
Please, Sign In to add comment