- •ПРЕДИСЛОВИЕ
- •Глава 1. НАЧАЛО РАБОТЫ В MATLAB
- •1.1. РАБОЧЕЕ ОКНО
- •1.3. ВЫЧИСЛЕНИЯ В КОМАНДНОМ ОКНЕ
- •1.5. ПЕРЕМЕННАЯ
- •1.6. ВСТРОЕННЫЕ МАТЕМАТИЧЕСКИЕ ФУНКЦИИ
- •1.7. ОЧИЩЕНИЕ КОМАНДНОГО ОКНА И РАБОЧЕГО ПРОСТРАНСТВА
- •1.8. КОМПЛЕКСНЫЕ ЧИСЛА
- •1.10. M-ФАЙЛЫ
- •КОНТРОЛЬНЫЕ ВОПРОСЫ
- •КОНТРОЛЬНЫЕ ЗАДАНИЯ
- •Глава 2. МАССИВЫ
- •2.1. ОПЕРАЦИИ НАД МАТРИЦАМИ В ЛИНЕЙНОЙ АЛГЕБРЕ
- •2.2. СПОСОБЫ ЗАДАНИЯ ВЕКТОРОВ В MATLAB
- •2.3. ОПЕРАЦИИ НАД ВЕКТОРАМИ В MATLAB
- •2.4. СПОСОБЫ ЗАДАНИЯ МАТРИЦ В MATLAB
- •2.5. ОПЕРАЦИИ НАД МАТРИЦАМИ В MATLAB
- •2.6. ИЗВЛЕЧЕНИЕ И ВСТАВКА ЧАСТЕЙ МАТРИЦЫ В MATLAB
- •КОНТРОЛЬНЫЕ ВОПРОСЫ
- •КОНТРОЛЬНЫЕ ЗАДАНИЯ
- •Глава 3. ЛИНЕЙНАЯ АЛГЕБРА
- •3.1. ОПРЕДЕЛИТЕЛЬ КВАДРАТНОЙ МАТРИЦЫ
- •3.2. ОБРАТНАЯ МАТРИЦА
- •3.3. РАНГ МАТРИЦЫ
- •3.4. СИСТЕМЫ ЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ УРАВНЕНИЙ
- •КОНТРОЛЬНЫЕ ВОПРОСЫ
- •КОНТРОЛЬНЫЕ ЗАДАНИЯ
- •Глава 4. ГРАФИКА В MATLAB
- •4.1. ДВУМЕРНАЯ ГРАФИКА
- •4.4. ПОСТРОЕНИЕ ПОВЕРХНОСТЕЙ И КРИВЫХ В ПРОСТРАНСТВЕ
- •4.5. ОФОРМЛЕНИЕ ГРАФИКОВ
- •КОНТРОЛЬНЫЕ ВОПРОСЫ
- •КОНТРОЛЬНЫЕ ЗАДАНИЯ
- •Глава 5. ВЕКТОРНАЯ АЛГЕБРА
- •5.1. ЛИНЕЙНЫЕ ОПЕРАЦИИ НАД ВЕКТОРАМИ
- •5.2. ПРОЕКЦИЯ ВЕКТОРА НА ОСЬ, ОСНОВНЫЕ СВОЙСТВА
- •5.3. БАЗИС И КООРДИНАТЫ ВЕКТОРА. ОРИЕНТАЦИЯ БАЗИСА
- •5.4. ДЛИНА ВЕКТОРА. НАПРАВЛЯЮЩИЕ КОСИНУСЫ ВЕКТОРА
- •5.5. ДЕЙСТВИЯ НАД ВЕКТОРАМИ В КООРДИНАТНОЙ ФОРМЕ
- •5.6. СКАЛЯРНОЕ ПРОИЗВЕДЕНИЕ ВЕКТОРОВ, ЕГО СВОЙСТВА
- •5.7. ВЕКТОРНОЕ ПРОИЗВЕДЕНИЕ ВЕКТОРОВ, ЕГО СВОЙСТВА
- •5.8. СМЕШАННОЕ ПРОИЗВЕДЕНИЕ ВЕКТОРОВ, ЕГО СВОЙСТВА
- •КОНТРОЛЬНЫЕ ВОПРОСЫ
- •КОНТРОЛЬНЫЕ ЗАДАНИЯ
- •Глава 6. АНАЛИТИЧЕСКАЯ ГЕОМЕТРИЯ НА ПЛОСКОСТИ
- •6.1. УРАВНЕНИЯ ПРЯМОЙ НА ПЛОСКОСТИ
- •6.2. ВЗАИМНОЕ РАСПОЛОЖЕНИЕ ДВУХ ПРЯМЫХ НА ПЛОСКОСТИ
- •6.3. РАССТОЯНИЕ ОТ ТОЧКИ ДО ПРЯМОЙ
- •6.4. КРИВЫЕ ВТОРОГО ПОРЯДКА
- •КОНТРОЛЬНЫЕ ВОПРОСЫ
- •КОНТРОЛЬНЫЕ ЗАДАНИЯ
- •Глава 7. АНАЛИТИЧЕСКАЯ ГЕОМЕТРИЯ В ПРОСТРАНСТВЕ
- •7.1. УРАВНЕНИЯ ПЛОСКОСТИ В ПРОСТРАНСТВЕ
- •7.2. ВЗАИМНОЕ РАСПОЛОЖЕНИЕ ДВУХ ПЛОСКОСТЕЙ
- •7.3. РАССТОЯНИЕ ОТ ТОЧКИ ДО ПЛОСКОСТИ
- •7.4. ПРЯМАЯ В ПРОСТРАНСТВЕ
- •7.5. ВЗАИМНОЕ РАСПОЛОЖЕНИЕ ДВУХ ПРЯМЫХ В ПРОСТРАНСТВЕ
- •7.6. ВЗАИМНОЕ РАСПОЛОЖЕНИЕ ПРЯМОЙ И ПЛОСКОСТИ В ПРОСТРАНСТВЕ
- •7.7. РАССТОЯНИЕ ОТ ТОЧКИ ДО ПРЯМОЙ В ПРОСТРАНСТВЕ
- •7.8. ПОВЕРХНОСТИ ВТОРОГО ПОРЯДКА
- •КОНТРОЛЬНЫЕ ВОПРОСЫ
- •КОНТРОЛЬНЫЕ ЗАДАНИЯ
- •БИБЛИОГРАФИЧЕСКИЙ СПИСОК
Обозначим через ax , ay , az проекции вектора |
a OM на координатные |
||||||
оси Ox, Oy, Oz соответственно: |
|
|
|
||||
|
|
|
|
|
|
||
ax |
OM1 |
, ay |
OM 2 |
, az |
OM 3 |
. |
|
Тогда выражение (5.3) примет следующий вид: |
|
|
|
||||
a ax i ay j az k . |
(5.4) |
||||||
По теореме 5.3 представление вектора a в виде (5.4) единственно. Представление вектора a в виде (5.4) называется разложением вектора a
по ортонормированному базису i , j, k .
ax , ay , az – координаты вектора a в базисе i , j, k , т. е. a ax ; ay ; az .
Использование MATLAB при нахождении разложения вектора по заданному базису
Пример 5.7. Разложить вектор |
̅= {11; −7; 13} по базису векторов |
̅ |
−2}. |
̅ = {2, −1,3)}, = {2,1,3} и ̅= {1, −2, |
|
Решение. Введем команды: |
|
>>d=[11;-3;13];
>>a=[2;-1;3];
>>b=[2;1;3];
>>c=[1;-2;-2];
>>A=[a,b,c];
>>A\d
ans =
3.0000
2.0000
1.0000
Ответ. = 3 + 2 + .
5.4. ДЛИНА ВЕКТОРА. НАПРАВЛЯЮЩИЕ КОСИНУСЫ ВЕКТОРА
Найдем длину вектора a в ДПСК.
Из прямоугольного параллелепипеда (рис. 5.16) имеем
OM 2 OM1 2 OM2 2 OM3 2 ,
т. е.
a 2 ax2 a2y az2 .
Получили
119
|
|
a |
a2 |
a2 |
a2 . |
(5.5) |
|
|
|
x |
y |
z |
|
Пример 5.8. Найти длину вектора a 2i 3 j k . |
|
|||||
Решение. ax 2, ay 3, az |
1. По формуле (5.5) имеем |
|
||||
a
22 32 1 2
14 .
Ответ. a
14 .
Пусть α, β, γ – углы между вектором a 0 и координатными осями Ox, Oy, Oz соответственно. По теореме 5.1 имеем
ax |
a |
cos , |
ay |
a |
cos , |
az |
|
a |
|
cos . |
(5.6) |
||||||||||||||||||||||||
Из выражения (5.6) получим: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
cos |
|
ax |
, |
cos |
ay |
|
|
|
, cos |
|
az |
|
|
. |
|
|
|
(5.7) |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
a |
|
|
|
a |
|
|
|
|
|
|
a |
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
Числа cos , cos , cos |
называются направляющими косинусами век- |
||||||||||||||||||||||||||||||||||
тора a . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Подставив выражения (5.6) в равенство |
|
a |
|
2 a2 |
|
a2 |
a2 |
, получим сле- |
|||||||||||||||||||||||||||
|
|
|
|||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
y |
|
z |
|
||
дующее соотношение: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
cos2 cos2 cos2 1. |
|
|
|
||||||||||||||||||||||||||||||
Таким образом, направление вектора в заданной системе координат опре- |
|||||||||||||||||||||||||||||||||||
деляется его направляющими косинусами. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
Легко видеть, что координатами единичного вектора e являются направ- |
|||||||||||||||||||||||||||||||||||
ляющие косинусы: e (cos ; cos ; cos ) . |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
Пример 5.9. Найти направляющие косинусы вектора a 2i 3 j k . |
|||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
Решение. ax 2, ay 3, az |
1, |
a |
|
14 . По формуле (5.7) имеем |
|||||||||||||||||||||||||||||||
cos |
|
2 |
|
|
, cos |
|
3 |
|
|
, cos |
|
1 |
. |
|
|||||||||||||||||||||
|
14 |
|
|
14 |
|
14 |
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
Использование MATLAB для вычисления длины вектора
инаправляющих косинусов вектора
ВMATLAB для вычисления длины вектора существует встроенная функция
norm(a)
120
Пример 5.10. Вычислить длину вектора = {4; 0; 3}: а) с использованием формулы (5.5);
б) с использованием встроенной функции norm.
Решение: а) Воспользуемся формулой (5.5). Введем команды:
>> a=[4 0 3];
>> dlina=sqrt(sum(a.^2)) % найдем длину вектора dlina =
5
б) Введем команды:
>> norm(a) ans =
5
Пример 5.11. Дан вектор = {4; 0; 3}. Найти: а) орт вектора ;
б) направляющие косинусы по формуле (5.7) (результат перевести в градусы);
в) угол наклона вектора к оси ОY (в градусах). Решение. Введем команды:
а) |
б) |
в) |
>> a=[4 0 3]; |
>> cos_a=a(1)/norm(a) |
>> b=acos(a(2)/norm(a))*180/pi |
>> a0=a/norm(a) |
cos_a = |
b = |
a0 = |
0.8000 |
90 |
0.8000 0 0.6000 |
>> cos_b=a(2)/norm(a) |
|
|
cos_b = |
|
|
0 |
|
|
>> cos_g=a(3)/norm(a) |
|
|
cos_g = |
|
|
0.6000 |
|
|
>> cos_a_grad=cos_a*180/pi |
|
|
cos_a_grad = |
|
|
45.8366 |
|
|
>> cos_b_grad=cos_b*180/pi |
|
|
cos_b_grad = |
|
|
0 |
|
|
>> cos_g_grad=cos_g*180/pi |
|
|
cos_g_grad = |
|
|
34.3775 |
|
|
|
|
5.5. ДЕЙСТВИЯ НАД ВЕКТОРАМИ В КООРДИНАТНОЙ ФОРМЕ
Пусть даны векторы a ax , ay , az , b bx ,by ,bz в координатной форме.
Справедливы следующие свойства:
121
1)a b ax bx , ay by , ay by ,
2)λa λax , λay , λaz .
Доказательство.
1)a b ax i ay j az k bx i by j bz k
ax bx i ay by j az bz k ax bx ,ay by ,ay by .
2)λa λ ax i ay j az k λax i λay j λaz k λax ,λay ,λaz .
Равенство векторов: a b ax bx , ay by , az bz .
Условие коллинеарности двух векторов: a || b |
a |
|
ay |
|
a |
||
x |
|
|
z |
. |
|||
b |
b |
|
b |
||||
|
x |
|
|
y |
|
z |
|
Координаты точки и вектора
Пусть в пространстве заданы ДПСК Охyz и М – произвольная точка. Координатами точки М называются координаты ее радиус-вектора OM . Пусть даны две точки
Найдем координаты вектора AB .
Имеем OA AB OB (рис. 5.17). Откуда получим разложение вектора AB в базисе i, j, k :
AB OB OA x2i y2 j z2 k x1i y1 j z1kx2 x1 i y2 y1 j z2 z1 k .
Если известны начало A x1 , y1 , z1 и конец B x2 , y2 , z2 вектора AB , то его координаты равны: AB x2 x1; y2 y1; z2 z1 .
z |
A |
|
|
|
k |
|
B |
|
|
|
|
|
|
|
O |
|
|
i |
|
j |
y |
|
|
|
|
||
|
|
|||
x
Рис. 5.17. OA AB OB
122
Действия над векторами в координатной форме: использование MATLAB
Пример 5.12. Даны точки A(1,1), B(3,2), C(2,0). Требуется: а) найти координаты векторов и ;
б) вычислить сумму и разность векторов в) изобразить векторы и черным,
обозначить вершины.
Решение. Введем команды:
и ;
С зеленым, фиолетовым,
|
а) |
|
б) |
>> A=[1,1]; |
>>summa=AB+BC |
||
>> B=[3,2]; |
summa = |
||
>> C=[2,0]; |
|
1 -1 |
|
>> AB=B-A |
>>raznost=AB-BC |
||
AB = |
|
raznost = |
|
2 |
1 |
3 |
3 |
>> BC=C-B |
|
|
|
BC = |
|
|
|
-1 |
-2 |
|
|
|
|
|
|
в) Построим векторы:
>>grid on, axis equal, hold on >>quiver(1,1,3-1,2-1,1,'k','LineWidth',1) >>quiver(3,2,2-3,0-2,1,'k','LineWidth',1) >>quiver(1,1,2-1,0-1,1,'g','LineWidth',6) >>quiver(2,0,1-2,1-0,1,'m','LineWidth',3) >>text(0.9,1.1,'\bfA') >>text(3.01,1.9,'\bfB') >>text(2.12,0.06,'\bfC')
>>text(1.7,0.12,'\bfAC','Color','green','rotation',-45) >>text(1.2,1,'\bfCA','Color','magenta','rotation',-45)
Изображение векторов появится в окне «Figure 1» (рис. 5.18, а).
Пример 5.13. Дан параллелограмм ABCD, известны координаты трёх его вершин A(1,3), B(3,7), C(8,7). Требуется:
а) найти координаты четвертой вершины D(x,y); б) построить параллелограмм.
Решение. Введем команды:
а) |
б) |
>>A=[1,3]; |
>>A=[1,3]; |
>>B=[3,7]; |
>>B=[3,7]; |
>>C=[8,7]; |
>>C=[8,7]; |
>>AC=C-A; % координаты вектора АС |
>>D=[6,3]; |
>>AB=B-A; % координаты вектора АВ |
>>grid on, axis equal, hold on |
|
|
123
>>AD=AC-AB; % координаты вектора АD |
>>quiver(1,3,3-1,7-3,1,'k','LineWidth',1) |
>>D=AD+A % координаты вершины D |
>>quiver(3,7,8-3,7-7,1,'k','LineWidth',1) |
D = |
>>quiver(8,7,6-8,3-7,1,'k','LineWidth',1) |
6 3 |
>>quiver(1,3,6-1,3-3,1,'k','LineWidth',1) |
|
>>text(1.1,2.8,'\bfA') |
|
>>text(2.8,7.2,'\bfB') |
|
>>text(7.8,7.2,'\bfC') |
|
>>text(6.1,2.8,'\bfD') |
|
|
Изображение параллелограмма появится в окне «Figure 1» (рис. 5.18, б).
а |
б |
Рис. 5.18. Изображение: а – векторов и , где A(1, 1), B(3, 2), C(2, 0); б – параллелограмма ABCD, где A(1, 3), B(3, 7), C(8, 7)
Пример 5.14. Даны точки A(-2; -3), B(-3; -3). Требуется: а) найти координаты вектора 5 ∙ ; б) построить вектор 5 ∙ .
Решение. а) Введем команды:
>>A=[-2,-3]; >>B=[-4,-3];
>>AB=B-A; % координаты вектора АВ >>proizv=5*AB % координаты вектора 5АВ proizv =
-10 0
б) Построим векторы и = 5 ∙ ̅̅̅̅, где точка С – конец вектора 5 ∙ . Введем команды:
>>C=proizv+A % координаты точки С - конца вектора 5АВ
C =
-12 -3
>>grid on, axis equal, hold on
124
