matlab基本运算点乘,MATLAB基本运算和仿真(1)ppt课件

《MATLAB基本运算和仿真(1)ppt课件》由会员分享，可在线阅读，更多相关《MATLAB基本运算和仿真(1)ppt课件(16页珍藏版)》请在人人文库网上搜索。

1、MATLAB基本运算和仿真,1,一、矢量的点乘dot和叉乘cross运算 用法：点乘dot(a,b)；叉乘cross(a,b)，其中a和b是矢量(向量)。,例1.1 已知 和 ，求(1) ；(2) ；(3)A和B之间的夹角。,解 代码如下 A=1,3,5; %矢量A B=2,4,6; %矢量B c=dot(A,B) %矢量A和B的点积 D=cross(A,B) %矢量A和B的叉积 a=sqrt(dot(A,A); %矢量A的模值 b=sqrt(dot(B,B); %矢量B的模值 theta=acosd(c/(a*b) %矢量A、B间的夹角,运行结果为 c = 44 D = -2 4 -2 th。

2、eta = 6.3532,2,二、标量场的等值线contour命令和矢量场的矢量线quiver命令,1. 对标量函数z=z(x,y)，contour命令的用法 contour(X,Y,Z,n,v) 必须输入的矩阵： X、Yx、y的取值范围矩阵， Z 对应于X和Y的数值矩阵。 选填参量： n绘制等值线的条数， v向量，该向量中元素的个数是n，元素的数值是等值线的值。,例1.2 在 ， 的范围内画出标量场 的9条等值线。,解 代码如下 X,Y=meshgrid(-2:0.2:2,-2:0.2:2); % x,y-2,2，间隔为0.2 Z=X.*exp(-X.2-Y.2); % 点(x，y)处的函数。

3、值z v=-0.4,-0.3, -0.2,-0.1,0,.1,0.2,.3,0.4; % 设置向量，有9个等值线的值 C,h=contour(X,Y,Z,v); % C是等高线矩阵，h是等高线句柄 clabel(C,h) % 标出等值线数值,3,4,2. 对二维矢量函数，quiver 命令的用法 (1)单个矢量： quiver(x,y,u,v) 表示在点(x，y)用箭头图形绘制矢量图，u是x方向分量，v是y方向分量。 (2)多点绘图： quiver(X,Y,U,V,scale) 其中 X、Yx、y的取值范围矩阵， U对应(x,y)点上x方向的分量， V对应(x,y)点上y方向的分量， scal。

4、e用来控制矢量线长度的实数，不填写时默认为1，有时为避免重叠，可设置较小的值。 3. 对三维矢量函数，用 quiver3 命令 quiver3(X,Y,Z,U,V,W),5,三、梯度gradient、散度divergence和旋度curl运算,1. 对二元函数F=f(x,y)，它的梯度gradient调用格式为 Fx,Fy=gradient(F) 其中F是二维矩阵，Fx是F在x方向的偏导数，Fy是F在y方向的偏导数，并假定自变量的间距是1。 对三元函数F=f(x,y,z)，它的梯度gradient调用格式为 Fx,Fy,Fz=gradient(F) 其中F是三维矩阵，Fx是F在x方向的偏导数，。

5、Fy是F在y方向的偏导数，Fz是F在z方向的偏导数，并假定自变量的间距是1。,例1.3 绘出例1.2中标量场的梯度的矢量图。 解 代码如下 X,Y=meshgrid(-2:0.2:2,-2:0.2:2); % x,y-2,2 ，间隔为0.2 Z=X.*exp(-X.2-Y.2); % 计算点(x，y)处的函数值z DX,DY=gradient(Z,.2,.2); % 计算梯度 contour(X,Y,Z) % 绘制等值线 hold on % 图形保持 quiver(X,Y,DX,DY) % 绘制矢量图,6,7,2. 对一个矢量场，计算它的散度用divergence命令，格式为 div=dive。

6、rgence(X,Y,Z,U,V,W),例1.4 画出矢量场 ，及其散度 的等值线。,解 利用坐标变换公式,矢量场A在直角坐标中的表达式为,8,代码如下： %矢量场A的矢量图及其散度的等值线 x=linspace(-2,2,21);y=linspace(-2,2,21); % 产生等间隔的21个点 X,Y=meshgrid(x,y); % 点(x，y)坐标 Ax= exp(-X.2-Y.2).*X./sqrt(X.2+Y.2+eps); % 矢量场A的x分量 Ay= exp(-X.2-Y.2).* Y./sqrt(X.2+Y.2+eps); % 矢量场A的y分量 div=divergence(。

7、X,Y,Ax,Ay); % 计算A的散度 quiver(X,Y,Ax,Ay,.8) % 画A的矢量图 Xlim(-2 2); % 将显示范围限制在网格内Ylim(-2 2); hold on % 图形保持 contour(X,Y,div) % 画A的等值线,9,10,3. 对一个矢量场，计算它的旋度用curl命令，格式为 curlx,curly,curlz,cav = curl(X,Y,Z,U,V,W),例1.4 画出矢量场 ，及其旋度A的矢量线。,解 利用坐标变换公式,矢量场A在直角坐标中的表达式为,11,代码如下： %矢量场A的矢量图 theta=0:pi/5:2*pi;rho=0.4:0。

8、.4:2; % 角和半径的范围和间隔 t,r=meshgrid(theta,rho); % 在圆柱坐标下建立网格 U=- exp(-r.2).* sin(t); % 矢量场A的x分量 V= exp(-r.2).*cos(t); % 矢量场A的y分量 x,y=pol2cart(t,r); % 将网格转化为直角坐标 quiver(x,y,U,V,.2) % 画场矢量A的二维矢量图 Xlim(-2 2); % 将显示范围限制在网格内 Ylim(-2 2);,12,13,%矢量场A的旋度的矢量图 x,y,z=meshgrid(-2:2); % 建立-2,2区间间隔为1的网格 A=exp(-x.2-y.。

9、2); % 计算函数 U=-A.*y./sqrt(x+eps).2+(y+eps).2); % 矢量场A的x分量 V=A.* x./sqrt(x+eps).2+(y+eps).2); % 矢量场A的y分量 W=A.*0; % 矢量场A的z分量 curlx,curly,curlz=curl(x,y,z,U,V,W); % 计算矢量场A的旋度 quiver3(x,y,z,curlx,curly,curlz) % 画旋度的三维矢量图 Xlim(-2 2); % 将显示范围限制在网格内Ylim(-2 2); Zlim(-2,2);,注意：代码中的eps时间算计中的最小数，为10-52，是为了避免出现分母为零的情况而加的。,14,15,练习题,为邻边构成的平行六面体的体积。,试求A、B之间的夹角；,在C上的分量；并用matlab验证。,利用matlab(1)在xy平面上画出 f 的等值线族；(2)画出f 的矢量图。,利用matlab求以矢量,2. 已知矢量,3. 已知二维标量场,16。