东华大学样题及答案(finished).doc

东华大学高等数学实验试题A 考试时间：90分钟 （附参考解答） 班级 学号 姓名 得分 上机考试说明： 开考前可将准备程序拷到硬盘, 开考后不允许用移动盘，也不允许上网； 领座考生试卷不同，开卷，可利用自己备用的书和其他资料，但不允许讨论，也不允许借用其他考生的书和资料。 解答(指令行，答案等)全部用笔写在考卷上。 计算题（60分） 要求：写出M函数（如果需要的话）、MATLAB指令和计算结果。 解线性方程组 并求系数矩阵的行列式。 指令行：A=[5 1 -1 0;1 0 3 -1;-1 -1 0 5;0 0 2 4];b=[1;2;3;-1];x=A\b,d=det(A) 结果：x1=1.4, x2= -5.9, x3=0.1, x4= -0.3. 行列式=70. 设 f(x,y) = 4 sin (x 3y)，求 。 指令行：syms x y; f=diff(4*sin(x^3*y),x); f=diff(f,y); f=subs(f,x,2); f=subs(f,y,3) 结果：1063.6 求方程 3x4+4x3-20x+5 = 0 的所有解。 指令行：roots([3 4 0 –20 5]) 结果：-1.5003 - 1.5470i, -1.5003 + 1.5470i, 1.4134, 0.2539 matlab求导命令diff调用格式: diff(函数) ， 求的一阶导数； diff(函数， n) ， 求的n阶导数(n是具体整数)； diff(函数，变量名)， 求对的偏导数； diff(函数， 变量名，n) ，求对的n阶偏导数； 在符号表达式或矩阵中进行符号替换 Syntax R = subs(S)R = subs(S, new)R = subs(S, old, new) R = subs(S, old, new) 利用new的值代替符号表达式中old的值，old为符号变量或是字符串变量名。new是一个符号货数值变量或表达式。的近似值。 方法一：指令行：syms x; s=int(1/sqrt(2*pi)*exp(-x^2/2),0,1); vpa(s,5)结果：0.34135 方法二：指令行：x=0:0.01:1; y=1/sqrt(2*pi)*exp(-x.^2/2);trapz(x,y) 结果：0.3413 方法三：M函数ex4fun.m function f=ex4fun(x) f=1/sqrt(2*pi)*exp(-x.^2/2); 指令行：s=quadl(@ex4fun,0,1) 结果：0.3413 int函数用以演算一函数的积分项， 这个函数要找出一符号式 F 使得diff(F)=f。如果积分式的解析式 (analytical form, closed form) 不存在的话或是MATLAB无法找到，则 int 传回原输入的符号式。相关的函数语法有下列 4个：int(f) 传回f对预设独立变数的积分值int(f,t) 传回f对独立变数t的积分值int(f,a,b) 传回f对预设独立变数的积分值，积分区间为[a,b]，a和b为数值式int(f,t,a,b) 传回f对独立变数t的积分值，积分区间为[a,b]，a和b为数值式int(f,m,n) 传回f对预设变数的积分值，积分区间为[m,n]，m和n为符号式matlab控制运算精度用的是digits和vpa这两个函数 digits用于规定运算精度，比如： digits(20); 这个语句就规定了运算精度是20位有效数字。但并不是规定了就可以使用，因为实际编程中，我们可能有些运算需要控制精度，而有些不需要控制。vpa就用于解决这个问题，凡是用需要控制精度的，我们都对运算表达式使用vpa函数。例如： digits(5); a=vpa(sqrt(2)); 这样a的值就是1.4142，而不是准确的1.4142135623730950488016887242097 又如： digits(5); a=vpa(sqrt(2)); b=sqrt(2); 这样a的值是1.4142，b没有用vpa函数，所以b是1.4142135623730950488016887242097...... vpa函数对其中每一个运算都控制精度，并非只控制结果。 digits(11); a=vpa(2/3+4/7+5/9); b=2/3+4/7+5/9; a的结果为1.7936507936，b的结果为1.79