MATLAB数学实验答案(全).pdf

第一次练习 教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。 补充命令 vpa(x,n) 显示x 的n 位有效数字，教材102 页 fplot( ‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中 为你的学号的后3 位（1-9 班）或4 位（10 班以上） m mx −sin mx mx −sin mx 1.1 计算lim 与lim 3 3 x →0 x x → x syms x limit((902*x-sin(902*x))/x^3) ans = 366935404/3 limit((902*x-sin(902*x))/x^3,inf)//inf 的意思 ans = 0 mx 1.2 y e x cos ，求y 1000 syms x diff(exp(x)*cos(902*x/1000),2)//diff 及其后的2 的意思 ans = (46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1 1 2 2 1.3 计算 ex +y dxdy 0 0 dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1)//双重积分 ans = 2.1394 4 x 1.4 计算 dx m2 +4x2 syms x int(x^4/(902^2+4*x^2))//不定积分 ans = atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 y e x cosmx , 求y (10) //高阶导数 syms x diff(exp(x)*cos(902*x),10) ans = -356485076957717053044344387763*cos(902*x)*exp(x)-39523230242776424948220058 84*sin(902*x)*exp(x) m 1.6 给出 +x 在x 0 的泰勒展式(最高次幂为4). 1000.0 1 syms x taylor(sqrt(902/1000+x),5,x)//泰勒展式 ans = - (9765625*451^(1/2)*500^(1/2)*x^4)/82743933602 +(15625*451^(1/2)*500^(1/2)*x^3) - (125*451^(1/2)*500^(1/2)*x^2)/406802 + (451^(1/2)*500^(1/2)*x)/902 +(451^(1/2)*500^(1/2))/500 1.7 Fibonacci 数列{xn }的定义是x1 1,x 2 1 ,xn xn−1 =+xn−2 (n 3,4, ) 用循环语句编 程给出该数列的前20 项（要求将结果用向量的形式给出）。//已知数列的递推公式求数 列的项 x=[1,1]; for n=3:20//3 比20 默认步长为多少??写50 时为什么会出错??