MATLAP 操作实例.DOCVIP

数 学 实 验 作 业 数学0802 袁乃坤 20080903059 一 、作图 1、笛卡儿叶形线 令a=1 区间 [-5,5,-5,5] 指令 ezplot(x^3+y^3-3*x*y=0,[-5,5,-5,5]); 三叶玫瑰线 令a=5 x=5*sin(3*h)*cos(h) y=5*sin(3*h)*sin(h) 区间[-2*pi,2*pi]) 指令 ezplot(5*sin(3*h)*cos(h),5*sin(3*h)*sin(h),[-2*pi,2*pi]); 2、椭圆抛物面 (x^2)/16＋(y^2)/25=z 网格图形 指令x=-10:0.1:10; y=-10:0.1:10; [X,Y]=meshgrid(x,y); Z= (X.^2)/16+(Y.^2)/25; surf(X,Y,Z) 平滑图形 指令 x=-10:0.1:10; y=-10:0.1:10; [X,Y]=meshgrid(x,y); Z= (X.^2)/16+(Y.^2)/25; surf(X,Y,Z) shading flat 二 计算 1、(1) 指令 syms x; f=(tan(x)-sin(x))/(x^3); limit(f,x,0) 运行结果 (2) 指令 syms x; f=x*(sqrt(x^2+1)-x); limit(f,x,inf) 运行结果 2、求下列函数的导数或者偏导数： (1) 指令syms x; f=(cos(x))^2*logx; diff(f) 运行结果 指令syms x; f=(cos(x))^2*logx; diff(f,x,4) 运行结果 (2) 指令syms x; z=x^3*y^2-3*x*y^3-x*y+1; diff(z,x,1) 运行结果 指令syms x; z=x^3*y^2-3*x*y^3-x*y+1; diff(z,x,2) 运行结果 指令syms x y; z=x^3*y^2-3*x*y^3-x*y+1; diff(diff(z,x,1),y,1) 运行结果 3、计算下列积分： (1) 指令syms x; f=(1+logx)/x; int(f,x,1,exp(1)) 运行结果 (2) 指令syms x y; f=x*sin(x); int(int(f,x,y,sqrt(y)),y,0,1) 运行结果 (3) 指令syms x; f=(1+sin(x))/sin(x)*(1+cos(x)); int(f,x) 运行结果 (4) 其中为柱面及平面所围成的空间闭区域的整个边界曲面的外侧 指令syms x y z; f1=int((x*sin(y)-z)*x,z,0,3); f2=int(f1,x,0,1); int(f2,y,0,2*pi) 运行结果 4 解非齐次线性方程组 指令： [x1,x2,x3,x4]=solve(4*x1+3*x2-x3+2*x4=6,8*x1+3*x2-3*x3+4*x4=12,3*x1+3*x2-2*x3-2*x4=6,x1,x2,x3,x4) x1 =12/7-15/7*x2 x2 =x2 x3 =-3*x2 x4 =-3/7+9/7*x2 运行结果 5 求矩阵 的行列式、逆、特征值和特征向量 指令a=[1 2 -1;3 4 -2;5 -4 1]; b=det(a) c=inv(a) [v,d]=eig(a) 运行结果 6 解微分方程 指令 y=dsolve(D2y-y=x-2,y(0)=0,y(1)=1,x) 运行结果 三 应用题. 一个圆柱形的容器，内装350升的均匀混合的盐水溶液。如果纯水以每秒14升的速度从容器顶部流入，同时，容器内的混合的盐水以每秒10.5升的速度从容器底部流出。开始时，容器内盐的含量为7千克。求经过时间t后容器内盐的含量。解：设时刻容器内盐的含量为，时刻容器内混合液的体积为： 则时刻溶液的浓度为： 认为在很小的一段时间dt