9.4Matlab在建筑、通信、经管等方面的应用-福州理工学院教务管理系统.ppt

高 等 数 学;*第四节 Matlab在建筑、通信、经管等方面的应用; *第四节 Matlab在建筑、通信、经管等方面的应用; *第四节 Matlab在建筑、通信、经管等方面的应用;一、Matlab在建筑方面的应用 1．Matlab在空间解析几何中的相关应用 (1) 曲线切线与法平面的作法 设空间曲线Γ的参数方程为： 这里假定 都在[α,β]上可导。 在曲线Γ上取对应于t=t0的一点M0(x0,y0,z0)及对应于t=t0+?t 的邻近一点M0(x0+?x,y0+?y,z0+?z)。作曲线的割线MM0，其方程为： 当点M沿着Γ趋于点M0时割线MM0的极限位置就是曲线在点M0处的切线?。考虑 当M →M0即?t→0时。得曲线在点M0处的切线方程为： ; 曲线的切向量：切线的方向向量称为曲线的切向量．向量 就是曲线 在点M0处的一个切向量。 法平面：通过点M0而与切线垂直的平面称为曲线Γ 在点M0 处的法平面，其法平面方程为 下面利用MATLAB求空间曲线的切线与法平面． 例1 画出曲线x = t-sin(t)；y = 1-cos(t)；z =4sin(t/2) 在 t=π/2 处的切线与法平面． 解：输入命令： syms t x y z x1=t-sin(t);y1=1-cos(t);z1=4*sin(t/2); S1=jacobian([x1,y1,z1],t); t=pi/2; x0=t-sin(t);y0=1-cos(t);z0=4*sin(t/2); s0=subs(S1); ;syms t F=-[x;y;z]+[x0;y0;z0]+s0*t；G=[x-x0,y-y0,z-z0]*s0； t=-2*pi:pi/50:2*pi; x=-3:0.2:3;y=-3:0.2:3;[x,y]=meshgrid(x,y); tt=-6:0.2:6; x1=t-sin(t);y1=1-cos(t);z1=4*sin(t/2); x2=321329765582243/562949953421312+tt; y2=1+tt;z2=2*2^(1/2)+tt*2^(1/2); z=(x-884279719003555/562949953421312+y-4)/-sqrt(2); plot3(x1,y1,z1) hold on plot3(x2,y2,z2) mesh(x,y,z) axis([-3 3 -4 4 -17 22]) view(-37.5+20*(i-1),25) grid on 显示图形： ;(2) 曲面切平面与法线的作法 设曲面Σ的方程为F(x,y, z)=0，M0(x0,y0, z0)是曲面Σ上的一点，并设函数F(x,y,z)的偏导数在该点连续且不同时为零．在曲面Σ上，通过点M0任意引一条曲线Γ．假定曲线Γ的参数方程式为： t→t0对应于点M0(x0,y0, z0)且 不全为零．曲线在点的切向量为： 考虑曲面方程F (x, y,z)=0两端在t→t0的全导数： 引入向量 ;易见T与n是垂直的．因为曲线Γ是曲面Σ上通过点M0的任意一条曲线．它们在点M0的切线都与同一向量n垂直．所以曲面上通过点M0的一切曲线在点M0的切线都在同一个平面上．这个平面称为曲面Σ在点M0的切平面．这切平面的方程式是： 曲面的法线： 通过点M0(x0,y0, z0)而垂直于切平面的直线称为曲面在该点的法线．法线方程为： 曲面的法向量：垂直于曲面上切平面的向量称为曲面的法向量．向量 就是曲面Σ在点M0处的一个法向量． ;下面利用Matlab求空间曲面的切平面与法线 例2 求椭圆抛物面z=x2+3y2在点（1，1，4）处的切平面及法线方程． 解：输入命令： syms t x y z F=x^2+3*y^2-z;x0=1;y0=1;z0=4; w=[x,y,z]; s1=jacobian(F,w); v1=subs(s1,x,x0);z2=subs(v1,y,y0);n=subs(z2,z,z0); F=[x-x0,y-y0,z-z0]*n' G=-[x;y;z]+[x0;y0;z0]+n'.*t F = 2*x-4+6*y-z G = -x+1+2*t -y+1+6*t -z+4-t;u=[0:0.1:2]';v=0:p