MIT-多变量微积分-1

网友笔记-不全配套     视频   习题视频 [mathjax] 多变量微积分(推荐)    MIT官方资料

MIT-多变量微积分-2
MIT-多变量微积分-3
MIT-多变量微积分-4

向量和点乘

向量的表示

A=a1i^+a2j^+a3k^=a1,a2,a3向量的方向dir(A)

向量的加法:几何意义和数值计算。

向量乘以一个标量:直接放大或者缩小。

 

点乘定义和几何意义


余弦定理和点乘的等价性

 Law of cosines |C|2=|A|2+|B|22|A||B|cosθ|C|2=CC=(AB)(AB)=AAABBA+BB=|A|2+|B|22AB因此AB=|A||B|cosθ即通过余弦定理推导出了矩阵叉乘的定义;同样可以由矩阵叉乘的定义反向推导余弦定理。

 

点乘的应用

(1)  求解空间中任意两个向量之间的夹角(判断垂直关系);求解向量长度。
 (2)  说明x+2y+3z=0是一个平面?
答:可以看成是与向量<1,2,3>垂直并且过原点的一个平面。
(3)  求解一个向量在某个方向上的投影分量。
(4)  求解平面图形的面积
首先看一种比较容易想到的方法:
接下来看如何利用点乘来计算(实现正余弦转换)
先求出A,它和A垂直,但是大小一样。那么对AB做点乘的绝对值,就是AB围城的平行四边形的面积,如果乘以二分之一,那么就是对应三角形的面积。

求解A
注意我们是利用AB=a2,a1b1,b2=a1b2a2b1而这个结果正好对应行列式|a1a2b1b2|的值。当然,如果值出现负数,那么就要取绝对值符号,小问题。
面积求解总结一下:

 

补充(六面体体积)

三个三维向量组成的行列式计算det(A,B,C)=|a1a2a3b1b2b3c1c2c3|=a1|b2b3c2c3|a2|b1b3c1c3|+a3|b1b2c1c2|三个向量组成的平行六面体的体积计算

 

叉乘

定义

A×B=|ijka1a2a3b1b2b3|=|a2a3b2b3|i|a1a3b1b3|j+|a1a2b1b2|k其实上面不是真正意义上的行列式,因为计算出的结果是一个向量,但是我们借助了行列式计算的规则。从定义上来说,应该直接写成|a2a3b2b3|i|a1a3b1b3|j+|a1a2b1b2|k,但是|ijka1a2a3b1b2b3|这样写有助于记忆而已。

定理

 

右手定则计算

i×j=k|ijk100010|=0i0j+1k

再探平行六面体体积

上面右侧,其实是按照点乘的定义(aibi=a1b1+a2b2+a3b3)进行计算的;三个行列式对应着三个维度的值,是我们从“叉乘定义”那里得到的。

应用

(1)  已知三个点,判断任意一个点是否在这三个点组成的平面内。
方法1:三个向量围成的平行六面体的体积为零。
方法2:原来的三个点组成的平面的法向量N垂直于P1P
法向量N=P1P2×P1P3
注:
(1)A×A=0
(2)  前面的点乘和叉乘的混合叫作“混合积”,这是和行列式一样的东西。

 

矩阵

基本知识

1. 改变坐标系,其实就是线性变换。
2. 两个矩阵相乘,我们可以这样思考运算规则:
3. 线性变换AB的先后顺序

4. 向量逆时针旋转90度

横轴线性变换之后,变成了原来纵轴的方向;纵轴线性变换,变成了原来横轴的负方向。总的效果就是整个坐标轴逆时针旋转了90度。R2表示逆时针旋转了两次90度,最终效果就是xy前面都加一个负号。

 

计算逆矩阵

第一步:计算余子式矩阵;
第二步:计算代数余子式矩阵;(特定元素变号)
第三步:将代数余子式的转置;
第四步:将第三步的矩阵除以原矩阵的行列式的值。

 

解特定法向量的平面

例:求解法向量为N=1,5,10而且过原点的平面?
解:
如果不是要求过原点,而是过特定点,比如
 Plane through P0=(2,1,1) and N=1,5,10
解:这里利用的是对该平面上任意一点P,它和平面上已知点组成的向量一定在该平面上,所以垂直于给定的法向量。注:这里的3除以法向量的模长就是原来过零点的该平面沿着法向量方向移动的距离。

 

平面的线性系统

三个平面交于一点:唯一解(系数矩阵可逆);
交于一条线:无数组解(系数矩阵行列式为零);
不相交:无解(系数矩阵不可逆,行列式也等于零)。

判断到底是无数组解还是无解,要通过消元法。

 

齐次方程组

零解是trivial solution(平凡解)
(a)  det(A)0,那么只有零解;
(b) det(A)=0,那么det(N1,N2,N3)=0N1,N2,N3 coplanar 三个平面的法向量共面,而垂直于这个面而且过原点的直线就是方程组的解。这条直线一方面垂直这些法向量,那么和这些面都平行,另一方面,这条线过原点,因此这条线必然同时处在三个平面之内。

 

参数方程(直线/摆线/弧长/开普勒第二定律)

参数方程表示空间中的直线

已知三维空间中两个点的坐标,求解两点连成的直线上的所有点的位置坐标?
解:将这条直线看作是“移动的点的轨迹”
{x(t)+1=t2y(t)2=tz(t)2=3t{x(t)=1+2ty(t)=2+tz(t)=23tQ0Q(t)=tQ0Q1Q(t)=Q0+tQ0Q1

平面和直线的交点

给定平面x+2y+4z=7,上面谈的直线和这个平面的交点坐标是多少?
解:

摆线的参数方程

摆线为例(Cycloid),晚上自行车车轮上固定的灯的轨迹就是摆线。特点就是无论这个轮子的转速是快是慢,都不影响轨迹线,因此我们不采用时间t作为变量参数,那么可以采用车轮水平移动的距离或者车轮转动的角度。下面我们以车轮转动的角度θ作为参数,构建摆线的参数方程。摆线的演示动画
最终的参数方程为OP=aθasinθx(θ),aacosθy(θ)问题:Question What happens near bottom?
答:在θ=0处进行泰勒展开(简单起见,让a=1),取几个高阶项,那么x(θ)θ(θθ36)θ36y(θ)1(1θ22)θ22显然θ=0处摆线的切线斜率无穷大。

摆线运动的速度和加速度


我们前面求解的是关于角度θ的参数方程,现在我们要求速度和加速度,那么必须将摆线方程转变成和时间有关的方程。假设这个轮子的半径为1,运动速度是1个单位长度每单位时间,轮子和地面接触的痕迹的长度就是轮子水平方向移动的距离。换句话说,轮子转动的角速度,为每单位时间转动1个角度,那么对应的弧长正好就是1

所以摆线时间的参数方程为r(t)=tsint,1cost速度矢量的求解方式为 Velocity vector v=drdt=dxdt,dydt,dzdt因此速度矢量的参数方程为v=1cost,sint( at t=0v=0!)速度的大小为|v|=(1cost)2+sin2t=12cost+cos2t+sin2t=22cost加速度a=dvdt e.g. cycloid a=sint,cost
上图表示上一个周期运动突然结束,然后速度突然变为零,同时在该瞬间获得一个向上的加速度。注意下式不成立|drdt|d|r|dt|drdt|表示的是速度矢量的模长,而后者就是表示位置向量tsint,1cost取模长对时间求导而已,没有很好的物理内涵。

 

弧长(Arc length)

s= distance travelled along trajectory,弧长对时间求导就是速度矢量的模长|dsdt|= speed =|v|
这里引入一个新的概念“单位正切矢量”T^,从轨迹上看,每一个点的单位正切矢量表示的是质点在该处的速度方向,那么点乘在该点移动的微小距离(标量),那么得到的就是移动前后位移矢量的微小变化。


 

开普勒第二定律

开普勒定律(Kepler's laws of planetary motion)是开普勒1609年发表的成果,其中的开普勒第二定律描述的是:
Motion of planets in a plane, and the area is swept out by the line from sun to planet at a constant rate. 就是单位时间内行星和太阳的连线扫过的面积相等。后来,牛顿用万有引力定律科学上解释了这一现象。下面的推导,发现要实现上述定律,只需要位移矢量和加速度矢量的叉乘为零,换句话说二者的方向共线,这显然满足,因为星星运动的加速度只来自于和太阳之间的万有引力,所以这个加速度的方向和位移矢量的方向恰好共线而且反向。矢量运算真强大!

拓展资料:
(1) 行星的运动—上海八中
(2) 行星轨道为什么是椭圆—科普时报

Leave a Reply