matlab选修结课作业,matlab在高等数学中的应用结课作业

matlab在高等数学中的应用结课作业

MatLab软件在作图中的应用 1、描点作图(MatLab及在电子信息课程中的应用131页) a=0:0.01:3; >> y=a； plot(y); 2、显函数作图(fplot)(MatLab及在电子信息课程中的应用184页) fplot(@(x)[tan(x),sin(x),cos(x)], 2*pi*[-1 1 -1 1]) 3、隐函数作图(ezplot)(MatLab及在电子信息课程中的应用185页) ezplot( x^2+y^2-9 );axis equal 4、极坐标作图(MatLab及在电子信息课程中的应用50页) t=0:pi/50:2*pi; r=sin(t).*cos(t); polar(t,r, -* ); 空间曲面作图(三维曲面图)(MatLab及在电子信息课程中的应用61页) t=0:pi/20:2*pi; [x,y,z]= cylinder(2+sin(t),30); subplot(2,2,1); surf(x,y,z); subplot(2,2,2); [x,y,z]=sphere; surf(x,y,z); subplot(2,1,2); [x,y,z]=peaks(30); surf(x,y,z); MatLab软件在高等数学应用 1、极限问题( 例的几何解释: 下面的示意图中共描了40个点，取如图1所示， N可以取2或大于2的正整数，即当n>2时，的值落在直线y=0.5,y=1.5之间；取如图2所示， N可以取20，即当n>20时，的值落在直线y=0.9,y=1.1之间。可见N的取值依赖于的取值，且不唯一。 源程序 clf subplot(1,2,1) hold on grid on n=1:40; %描40个点 m=1+(-1).^(n-1)./n; plot(n,m, . ) fplot( 0.5 ,[0,40]) %取0.5 fplot( 1.5 ,[0,40]) axis([0,40,0,2]) subplot(1,2,2) hold on grid on plot(n,m, . ) fplot( 0.9 ,[0,40]) %取0.1 fplot( 1.1 ,[0,40]) axis([0,40,0,2]) 2、导数问题 函数在x=0处连续,但在x=0处不可导.如图所示. clf; >> subplot(1,2,2) hold on fplot( abs(x) ,[-1,1]) 3、不定积分 用MatLab软件，计算下列不定积分 (1) syms x int( x^3*exp(-x^2) ,x) ans =-1/2*x^2/exp(-x^2)-1/2/exp(-x^2) ans = -1/2*(1+x^2)*exp(-x^2) ans = -1/2*x^2/exp(-x^2)-1/2/exp(-x^2) 定积分 利用MatLab软件找出满足定积分中值定理的点,使得 。 syms x y=sqrt(1-x^2); zhi=int(y,0,1) %计算 z=y-zhi; zf=char(z); fzero(zf,0.5) ; 执行结果：zhi = 1/4*pi 4、求函数的泰勒展开式问题 syms x taylor(exp(x),x,4,1) 执行结果：ans = exp(1)+exp(1)*(x-1)+1/2*exp(1)*(x-1)^2+1/6*exp(1)*(x-1)^3 MatLab软件在线性代数应用 1、行列式问题 计算行列式的值,其中a,b,c,d是参数. syms a b c d A=[a,1,0,0;-1,b,1,0;0,-1,c,1;0,0,-1,d] det(A)； 执行结果： A = [ a, 1, 0, 0] [ -1, b, 1, 0] [ 0, -1, c, 1] [ 0, 0, -1, d] ans = a*b*c*d+a*b+a*d+c*d+1 2、矩阵运算问题(包括加，减，乘，转置，求逆，求秩) 矩阵的加、减、乘(线性代数13页) >> A=[1 2 3 4 5];B=3:7;C=linspace(2,4,3); AT=A ;BT=B ; >> E1=A+B,E2=A-B,F=AT-BT,G1=3*A,G2=B/3 E1 = 4 6 8 10 12 E2 = -2 -2 -2 -2 -2 F = -2 -2 -2 -2 -2 G1 = 3 6 9 12 15 G2 = 1.0000 1.3333 1.6667 2.0000 2.3333 矩阵的转置求逆(线性代数76页) t1=0:6; t2=0:10; t3=1:10; t4=0:6; m=3; y1=sin(t1); y2=cos(t2); y3=zeros(1,10); for k=4:10 y3(k)=y1(k-3); end; y4=fliplr(y1); subplot(2,2,1),stem(t1,y1, r ); subplot(2,2,2),stem(t2,y2, g ); subplot(2,2,3),stem(t3,y3); subplot(2,2,4),stem(t4,y4); 矩阵的求秩 ,求矩阵A的秩，并求A的一个最高阶非零子式。 A=[3,2,0,5,0;3,-2,3,6,-1;2,0,1,5,-3;1,6,-4,-1,4]; rank(A) rref(A) 执行结果： ans = 3 ans = 1.0000 0 0.5000 0 3.5000 0 1.0000 -0.7500 0 -0.2500 0 0 0 1.0000 -2.0000 0 0 0 0 0 3、线性方程组的求解问题 求下列非齐次方程组的通解。 (3) A=[1,-1,-1,1;1,-1,1,-3;1,-1,-2,3]; B=[0;1;-1/2] ; jfch(A,B) 执行结果: ra = 2 the special solution is : ss = 0.5000 0.5000 0 0 the basic solution is : bs =