《数学实验试题(共24页).doc》由会员分享，可在线阅读，更多相关《数学实验试题(共24页).doc（24页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、精选优质文档-倾情为你奉上电子科技大学成都学院二零零 九 至二零一零 学年第 二 学期期 末 考试数学实验 课程考试题 A卷（120分钟） 考试形式： 闭 考试日期 年 月 日课程成绩构成：平时 10 分， 期中 0 分， 实验 30 分，期末 60 分（本试卷100分）一二三四合计一、单项选择题（共30分）1符号计算与一般数值计算有很大区别，它将得到准确的符号表达式。在MATLAB命令窗口中键入命令syms x，y1=x2；y2=sqrt(x);；int(y1-y2,x,0,1)，屏幕显示的结果是（ ）（A）y1 =x(1/2) （B）ans= 1/3； （C）y2 =x2； （D）ans=
2、 -1/32在MATLAB命令窗口中键入命令B=8,1,6;3,5,7;4,9,2；B*B(:,1)。结果是（ ）（A）ans= （B）ans= （C）ans= （D）ans=91 67 67 6767 91 67 6767 67 91 673在MATLAB命令窗口中，键入命令syms x,F=1/(2+cos(2x);ezplot(diff(F)，结果是（ ）（A）绘出函数F在0，p 的图形； （B）绘出函数F在2p，2p 的图形；（C）绘函数F的导函数在0，p 的图形；（D）绘函数F的导函数在2p，2p 的图形4在MATLAB命令窗口中，键入命令binopdf(3,6,1/2),结果是（A
3、）ans= (B) ans= (C) ans= (D) ans=15/64 5/16 3/32 15用赋值语句给定x数据，计算 对应的MATLAB表达式是（ ）（A）sqrt(7sin(3+2x)+exp(3)log(3) （B）sqrt(7*sin(3+2*x)+exp(3)*log(3)（C）sqr(7*sin(3+2*x)+e3*log(3) （D）sqr(7sin(3+2x)+ e3log(3)6MATLAB语句strcat(int2str(2010),年是, s,年)的功能是（ ）（A）将数据2010转换为符号； （B）将数据2010与符号变量合并；（C）将几个符号变量合并为一个；
4、（D）将符号变量转换为数值变量；7语句L=sqrt(pi); x=fix(1000*L)/1000的功能是（ ）（A）将无理数p 取三位近似； （B）将无理数 取四位近似数（C）将无理数 取三位近似数； （D）将有理数 取四位近似数8MATLAB绘三维曲面需要构建网格数据，语句x,y=meshgrid(-2:2)返回数据中（ ）（A）x是行向量，y是列向量； （B）x是列向量，y是行向量；（C）x是行元素相同的矩阵； （D）x是行向量相同的矩阵9下面有关MATLAB函数的说法，哪一个是错误的（ ）（A）函数文件的第一行必须由function开始，并有返回参数，函数名和输入参数；（B）MATLA
5、B的函数可以有多个返回参数和多个输入参数；（C）在函数中可以用nargout检测用户调用函数时的输入参数个数；（D）如果函数文件内有多个函数，则只有第一个函数可以供外部调用；10设a,b,c表示三角形的三条边，表达式a+bc|a+cb|b+ca， （ ）（A）是三条边构成三角形的条件； （B）是三条边不构成三角形的条件；（C）构成三角形时逻辑值为真； （D）不构成三角形时逻辑值为假二、程序阅读理解（24分）1用plot命令绘制平面曲线v0=515;alpha=45*pi/180;T=2*v0*sin(alpha)/9.8; %第二行t= T*(0:16) /16;x=v0 *t *cos(al
6、pha);y=v0 *t *sin(alpha) -0.5*9.8*t.2;plot(x,y,x,y,r*)Xmax=x(17)（1）对下面有关程序的功能的说法确切的是（ ）（A）以515为初速度的抛射曲线的绘制；（B）以515为初速度的抛射曲线的绘制以及计算射程；（C）以515为初速度以45度为发射角的抛射曲线的绘制以及计算射程；（D）以515为初速度以45度为发射角的抛射曲线的绘制。（2）第二行语句的功能是（ ）（A）计算射程；（B）计算飞行时间；（C）计算航点的横坐标； （D）计算航点的纵坐标。2关于“牟合方盖”的实验程序如下h=2*pi/100;t=0:h:2*pi;r=0:0.05:
7、1;x=r*cos(t);y=r*sin(t);z=sqrt(1-x.2); %第三行meshz(x,y,z),axis offcolormap(0 0 1)view(-47,56),hold onx1=cos(t);y1=sin(t);z1=abs(sin(t);plot3(x1,y1,z1,ro)；（1）下面有关程序的功能的说法确切的是（ ）（A）绘圆柱面x2 + y2 = 1, x2 + z2 = 1的交线；（B）绘圆柱面x2 + y2 = 1, x2 + z2 = 1所围区域的边界曲面；（C）绘圆柱面x2 + y2 = 1, x2 + z2 = 1的交线及所围区域的边界曲面；（D）绘圆
8、柱面x2 + y2 = 1, x2 + z2 = 1的交线及所围区域的边界曲面的上半部分。（2）关于第三行语句错误的解释是（ ）（A）z是矩形域上曲顶柱面高度值； （B）z是与y同型的矩阵；（C）z是圆域上曲顶柱面高度值； （D）z是与x同型的矩阵3中国农历年由天干(10干)和地支(12支)相配而成，计算农历年的MATLAB程序如下year=input(input year:=);S1=辛壬癸甲乙丙丁戊己庚；S2=酉戍亥子丑寅卯辰巳午未申;k1=mod(year,10); k2=mod(year,12);if k1=0,k1=10;end %第五行if k2=0,k2=12;ends1=S1(
9、k1); s2=S2(k2)；strcat(int2str(year),年是, s1,s2,年)（1）输入2010，实验程序的结果将给出（ ）（A）2010年是丁亥年；（B）2010年是庚寅年；（C）2010年是戊子年；（D）2010年是丙戍年（2）第五行语句的功能是（ ）（A）当年份是12的倍数时定位为地支12 （B）当年份是12的倍数时定位为天干12;（C）当年份是10的倍数时定位为地支10；（D）当年份是10的倍数时定位为天干104 非负函数 y =f(x)在有限区间上的图形为上半平面的一条曲线，曲线绕x轴旋转时，产生以x为对称轴的旋转曲面，其体积syms a b xf=exp(a*x)
10、*sin(b*x);f1=subs(f,a,-0.2);f2=subs(f1,b,0.5);V=pi*int(f2*f2,x,0,2*pi)double(V)t=(0:20)*pi/10;theta=t;r=f2(t);x=t*ones(size(t); %第九行y=r*cos(theta); %第十行z=r*sin(theta); %第十一行mesh(z,y,x)colormap(0 0 0)axis offview(-17,54)（1）关于程序的功能确切的说法（ ）(A)计算曲线段 f(x)=exp(a x)sin(b x), 绕X轴旋转的旋转曲面体积(B)计算曲线段 f(x)=exp(-
11、0.2x)sin(0.5x), 绕X轴旋转的旋转曲面体积(C) 计算曲线段 f(x)=exp(a x)sin(b x), 绕X轴旋转的旋转曲面体积并绘图（D）计算曲线段 f(x)=exp(-0.2x)sin(0.5x), 绕X轴旋转的旋转曲面体积并绘图（2）由第九行至第十一行语句可得旋转曲面的方程（ ）(A) （B）（C） (D)三、程序填空（32分）1维维安尼（Viviani）体是圆柱体( x R/2)2 + y2 R2/4被球面x2 + y2 + z2 = R2所割下的立体。下面的实验程序功能是取R=2求体积上半部分，先利用符号计算处理重积分并转换为数值数据，再用蒙特卡罗方法计算体积做对比
12、。完成下面程序填空syms x y;f=sqrt(4-x2-y2);y1=-sqrt(2*x-x2);y2= S1=int(f,y,y1,y2);S2=int(S1,x,0,2)V= P=rand(10000,3);X=2*P(:,1);Y=2*P(:,2)-1;Z=2*P(:,3);II=find(X-1).2+Y.2=1&Z=sqrt(4-X.2-Y.2);V1=8*length(II)/100002对于二重积分 ，D是直线y= x 2和抛物线 y2 = x 所围区域。下面数学实验程序的功能是利用MATLAB符号计算方法将二重积分处理为二次积分进行计算（选择先对x积分后对y积分的秩序），然
13、后将计算结果的符号表达式转换为双精度实数，syms x y;f=x*y*y;x1=y*y;x2= S1=int(f,x,x1,x2);S2=int(S1,y,-1,2);S= 3某年A、B两城镇人各有1000人，以后每年A镇人口10%迁往B镇; B镇人口15%迁往A镇。则有迁移矩阵L，L有两个互异特征值l1，l2，属于两个特征值的特征向量分别为：a 1，a 2，记初始人口分布X0=1000，1000T，则存在不全为零的数c1，c2使得 X0 = c1a 1 + c2a 2（解此方程组求出c1，c2）。下面实验程序主要功能是利用LnX0 = c1l1na 1+ c2l2na 2计算人口变化。完成
14、程序填空n=input(input n:=);A=0.9,0.15;0.1,0.85;X0=1000;1000;P,D=eig(A);C= c1=C(1);c2=C(2);alfa1=P(:,1);alfa2=P(:,2);lamda1=D(1,1);lamda2=D(2,2);Xn= 4判断闰年条件有二能被4整除，但不能被100整除；能被4整除，又能被400整除。n=input(input n:=);if mod(n,4)= &mod(n,100)=0 disp(是闰年)elseif mod(n,400)=0disp(是闰年)elsedisp end四、（14分）线性规划问题建筑公司承建办公
15、楼和住宅楼。建办公楼将获利润500元/平方米，建住宅楼获利润600元/平方米。总建筑面积不少于5000m2，办公楼的面积不能大于5000 m2，住宅楼不能大于3000m2。建立最优化模型，编写使用linprog（）求解问题的简单程序。 (电子科技大学成都学院二零零九至二零一零学年第 二 学期期末考试数学实验 课程考试题 B卷（120分钟） 考试形式： 闭 考试日期 年 月 日课程成绩构成：平时 10 分， 期中 0 分， 实验 30 分，期末 60 分（本试卷满分100分）一二三四合计一、单项选择题（共30分）1在MATLAB命令窗口中，键入命令syms x； int(x*sin(x)。结果是
16、（ ）（A）ans= sin(x)-x*cos(x)； （B）ans= cos(x)+x*sin(x)；（C）ans= sin(x)-cos(x)； （D）ans= -1/2*cos(x)*sin(x)+1/2*x2在MATLAB命令窗口中键入命令A=1 4 2;3 1 2;6 1 5;det(A(1:2,2:3).*A(2:3,2:3)。结果是（ ）（A）ans= -143 （B）ans= 60 （C）ans= 36 （D）ans= -193MATLAB命令x = 3: 2: 100 将创建等差数列，该数列是（ ）（A）以3为初值的98个数,； （B）以 100为终值的98的个数；（C）以9
17、9为终值的49个数； （D）以3为初值的97个数。4设n阶方阵A的特征值为： (i=1，2，n)，称为矩阵A的谱半径, 则下列MATLAB求谱半径命令是（ ）（A）max(abs(eig(A)； （B）abs(max(eig(A)；（C）max(norm(eig(A)； （D）norm(max(eig(A)5在MATLAB命令窗口中输入命令data=4 1 2 3 1 3 1 4 2 4;y=hist(data,4)，结果是（ ）（A）y= 4 1 2 3；（B）y=3 2 3 2；（C）y= 3 2 2 3 ；（D）y= 4 2 1 16数学表达式对应的MATLAB表达式是。（ ）（A）sq
18、rt(7*sin(3+2*x)+exp(3)*log(3) （B）sqrt(7sin(3+2x)+exp(2)log(3) （C）sqrt(7*sin(3+2*x)+e3*log(3) （D）sqrt(7sin(3+2x)+ e2 log(3)7语句L=sqrt(pi); x=fix(100*L)/100的功能是（ ）（A）将无理数p 取三位近似； （B）将取两位近似数（C）将取三位近似数； （D）将无理数p 取两位近似8MATLAB语句x,y=meshgrid(-2:2) 的数据结果中（ ）（A）x是行向量，y是列向量； （B）x是五行五列的矩阵；（C）x是行元素相同的矩阵； （D）x是列向
19、量相同的矩阵9MATLAB的语句colormap(0 0 1) （ ）（A）将三维网面图确定为红色； （B）将三维网面图确定为绿色；（C）将三维网面图确定为蓝色； （D）语句使用格式错误10将带小数的实数处理为整数称为取整，常用四种取整法则是：向正无穷大方向取整、向负无穷大方向取整、向零方向取整和四舍五入取整。MATLAB提供了如下四个取整函数，若a = -2.3，对a 取整的结果是 -2，则不应该选用下面哪个函数。( )（A）round； （B）floor； （C）ceil； （D）fix二、程序阅读理解（24分）1数学实验程序如下L=3/4,1/8,1/8;1/6,2/3,1/6;1/4,
20、1/4,1/2;X2=100;80;120;X=X2;x2=X(2);for k=1:4 X=L*X x2=x2;X(2);endbar(x2) %第八行colormap(1 1 1)（1）实验程序中的循环语句将显示（ ）（A）三阶矩阵L的特征值； （B）方程组X=LX的解；（C）LX的第一分量数据； （D）向量LnX变化规律（2）第八行语句的功能是（ ）（A）绘X的变化曲线； （B）绘图表示方程组X=LX的解；（C）绘LnX第二分量条形图； （D）绘LnX第一分量曲线2十二属相的生肖问题的MATLAB程序如下year=input(input year:=);S=鸡狗猪鼠牛虎兔龙蛇马羊猴; k
21、=mod(year,12); if k=0,k=12;end %第四行s=S(k); s=strcat(int2str(year),年是, s,年)（1）输入2009，实验程序的结果将给出（ ）（A）2009年是龙年；（B）2009年是蛇年；（C）2009年是牛年；（D）2009年是羊年（2）第四行语句的功能是（ ）（A）当年份是12的倍数时定位为猪年；（B）当年份是12的倍数时定位为第12属相；（C）当年份是12的倍数时定位为猴年；（D）当年份是12的倍数时定位为鼠年3一个古典概率问题叙述如下：甲乙丙丁四人按逆时针方向围坐玩麻将甲将两枚均匀骰子同时掷一次，根据骰子点数之和确定先摸谁家面前的麻
22、将例如点数之和为3，7，11时均确定为先摸丙前面的麻将。实验程序如下function Fn=majiang(k)if nargin=0,k=4;endif k4,error(请输入正确编号1到4);endk1=k+4;k2=k+8;N=2000;x=1+fix(6*rand(1,N);y=1+fix(6*rand(1,N);w=find(x+y=k|x+y=k1|x+y=k2); %第七行n=length(w);Fn=n/N %（1）没有输入数据时调用该函数，则程序运行后，将显示（ ）（A）先摸甲面前的麻将的频率； （B）先摸乙面前的麻将的频率；（C）先摸丙面前的麻将的频率； （D）先摸丁面前
23、的麻将的频率（2）第七行语句的功能是（ ）（A）统计2000次随机试验中先摸编号为k的人面前的麻将的频数；（B）统计2000次随机试验中先摸编号为k的人面前的麻将的索引值；（C）统计2000次随机试验中先摸编号为k的人面前的麻将的次数；（D）计算2000次随机试验中先摸编号为k的人面前的麻将的频率。4数学实验程序如下h=439;H=2384;R=6400; a=(h+H+2*R)/2;c=(H-h)/2; e1=c/a; b=sqrt(a*a-c*c); syms e2 t f=sqrt(1-e2*cos(t)2); ft=subs(f,e2,e1*e1); S=int(ft,0,pi/2);
24、 L=4*a*double(S)； V=L/(114*60)； %第九行s1=pi*a*b/(114*60); Vmax=2*s1/(h+R) Vmin=2*s1/(H+R) （1）实验程序的运行后，将显示的数据是（ ）（A）卫星轨道的周长数据； （B）卫星运行的近地速度和远地速度；（C）卫星运行时向径每秒扫过的面积；（D）卫星运行的平均速度数据（2）第九行语句的功能是（ ）（A）计算卫星运行的最小速度； （B）计算卫星运行时向径每秒扫过的面积;（C）计算卫星运行的最大速度； （D）计算卫星运行轨道的平均速度 三、程序填空（32分）1下面实验程序的功能是输入三角形边长数据用海伦公式计算出三角形
25、面积。仔细阅读程序开始部分符号“%”后的注记，根据注记提示的功能完成程序填空function S,C=triangle_area(a,b,c)%海伦公式计算三角形面积（2007-06-20）%triangle_area(a,b,c)：得到边长为a,b,c三角形的面积%如果输入数据不满足三角形两边之和大于第三边则返回出错信息%triangle_area(a,b)：a是等腰三角形的腰长，b是另外一边，返回面积%triangle_area(a)：a是等边三角形的边长，返回面积%S,C=triangle_area(a)：S是等边三角形面积，C是等边三角形周长if nargin=1 b=a;c=a;en
26、dif nargin=2 ; endif error(请输入正确的三角形边长);endC=a+b+c;p=C/2;S=sqrt(p*(p-a)*(p-b)*(p-c);2对于任意正整数n，如果n 只能被1和它自身整除，则称这个数为素数(或质数)。判素数程序的算法思想是试商法，即用2，3,，（n-1）去除n，如果能被这些数中一个整除，则n不是素数，否则是素数。完成下面填空。n=input(input n:=); for k=2:n-1 if mod(n,k)= ,break,end endif k n-1 disp(不是素数)else disp end3汽车紧急刹车问题数据拟合实验 V 20 2
27、5 30 35 40 45 50 55 60 65 70T 20 28 41 53 72 93 118 149 182 221 266V表示刹车时汽车行驶速度(英里/小时),T表示刹车后汽车滑行距离(英尺) 分别做二次多项式和三次多项式数据拟合实验，并绘出数据拟合曲线的图，计算出残差平方和，完成如下实验程序填空v=20 25 30 35 40 45 50 55 60 65 70*1.609;T=20 28 41 53 72 93 118 149 182 221 266*.3048;P2= T2=polyval(P2,v);R2=sum(T-T2).2)figure(2),plot(v,T,*,
28、v,T2)P3=polyfit(v,T,3);T3=polyval(P3,v);R3= figure(3),plot(v,T,*,v,T3)4二阶正交矩阵作用于向量a 时，其效果是将向量a 旋转，旋转角为（逆时针旋转为正）。把一个以原点为中心的正方形旋转pi/24，并做适当缩小，迭代30次形成下图。完成如下程序填空xy=-4 -4;4 -4;4 4;-4 4;-4 -4; A=cos(pi/24) -sin(pi/24);sin(pi/24) cos(pi/24); x=xy(:,1);y=xy(:,2); axis off line(x,y)for k=1:30 xy=.89*xy*A; x
29、= ; y= ; line(x,y),end四、 （14分）线性规划问题投资者拥有1000（万元）用于投资，共有4种投资方式，下表给出了预期收益率：投资方式A1A2A3A4收益率(%)3.51036要求满足如下条件：（1）总投资额不超过现有资金的80%；（2）投资A2不超过投资A1和A4的3倍；（3）投资A1不低于100万元；（4）投资A3不超过300万元；（5）投资A4在50万800万元之间。建立最优化模型，编写使用linprog（）求解问题的简单程序。 电子科技大学二零零八到二零零九学年第二学期期末考试数学实验课程考试题 A卷 (120分钟) 考试形式：闭卷 考试日期：2098年7月8日课
30、程成绩构成：平时10分，期中0分，实验30分，期末60分(本试卷满分100分)所有答案一律写在答题纸上，写在试卷上无效。一、单项选择题(20分)1、三阶幻方又称为九宫图，提取三阶幻方矩阵对角元并构造对角阵用( ) (A) diag(magic(3); (B) diag(magic); (C) diag(diag(magic(3); (D) diag(diag(magic)。2、MATLAB命令P=pascal(3)将创建三阶帕斯卡矩阵，max(P)的计算结果是( ) (A) 1 2 3 (B) 1 2 1 (C) 3 6 10 (D) 1 3 63、命令J=1;1;1*1,2,3;A=j+j-
31、1将创建矩阵( ) (A) ; (B) (C) (D) 4、data=rand(1000,2);x=data(:,1);y=data(:,2);II=find(yx.2);的功能是( ) (A) 统计2000个随机点中落入特殊区域的点的索引值； (B) 统计1000个随机点落入特殊区域的点的索引值； (C) 模拟2000个随机点落入特殊区域的过程； (D) 模拟1000个随机点落入特殊区域的过程。5、MATLAB计算二项分布随机变量分布律的方法是( ) (A) binocdf(x,n,p); (B) normpdf(x,mu,s); (C)binopdf(x,n,p); (D) binornd
32、(x,n,p)。6、MATLAB命令syms e2;f=sqrt(1-e2*cos(t)2);S=int(f,t,0,pi/2)功能是( ) (A) 计算f(x)d 0, 上的积分； (B) 计算f(t)的不定积分符号结果； (C) 计算f(x)的积分的数年结果； (D) 计算f(t)定积分的符号结果。7、y=dsolve(Dy=1/(1+x2)-2*y2,y(0)=0,x);ezplot(y)的功能是( ) (A) 求微分言和特解并绘图； (B) 解代数方程 (C) 求定积分； (D)求微分方程通解。8、X=10000 ;0.5*asin(9.8*X/(5152)的功能是计算关于抛射体问题的
33、( ) (A) 十公里发射角； (B) 十公里飞行时间； (C)最大飞行时间； (D)最大射程。9、theta=linspace(0,2*pi,100) ;r=cos(4*theta) ;polar(theta,r,k)功能是( ) (A) 绘四叶玫瑰线； (B)绘三叶玫瑰线； (C)绘心脏线； (D) 绘八叶玫瑰线。10、北京和纽约的经度分别是：东经118和西经76，根据经度差计算时差用( ) (A) fai1=118;fai2=-76;Dfai=(fai1+fai2)/24; (B) fai1=118;fai2=-76;Dfai=(fai1+fai2)/15; (C) fai1=118;f
34、ai2=-76;Dfai=(fai1-fai2)/24; (D) fai1=118;fai2=-76;Dfai=(fai1-fai2)/15。二、程序阅读题 (40分)1、直方图功能是将数据分为n个类，统计各个类的数据量并绘图。借用现有的直方图命令hist，编写新直方图程序如下。function m=myhist(data,n)if nargin=1,n=7;endXmin=min(data);Xmax=max(data);h=(Xmax-Xmin)/n;m=hist(data,n)/length(data)/h;t=linspace(Xmin,Xmax,n+1);II=1:4:4*n-3;J
35、J=1:n;x(II)=t(JJ);y(II)=zeros(1,n);x(II+1)=t(JJ);y(II+1)=m;x(II+2)=t(JJ+1);y(II+2)=m;x(II+3)=t(JJ+1);y(II+3)=zeros(1,n); plot(x,y,k)(1) 变量data存放了1000个数据，在命令窗口调用myhist(data)的结果是( c) (A) 只绘数据的直方图而不显示被分类后各类的数据量； (B) 只显示被分类后各类的数据量而不绘数据的直方图； (C) 既绘数据直方图也显示被分类后各类的数据量； (D) 根据默认值在数据范围内插入七等分点绘直方图。(2) 关于新直方图绘
36、图程序下面说法不正确的是( D ) (A)h是n等分直方图中小区间长度； (B) 修改程序最后一行可绘红色直方图； (C) 直方图中所有小矩形面积之和为1； (D) 直方图中所有小矩形的高度和为1。 2、3n+1问题反映一个数学猜想：对任一自然数n，按如下法则进行运算：若n为偶数，则将n除2，若n为奇数，则将n乘3加1。重复这种操作，结果终会为1。实验程序如下。function k,N=threeN(n)if nargin=0,n=5;endk=1;N=n;while n=1 r=rem(n,2); if r=0 n=n/2; else n=3*n+1; end N=N,n;k=k+1;end
37、(1)在MATLAB命令窗口中直接调用threeN运行结果为( ) (A)只显示k的最后数值为6； (B) 只显示k的最后数值5； (C) 同时显示k和N的数据； (D) 仅显示N的所有数据。(2)实验程序运行过程中( )(A) 输入变量n不发生改变； (B)N是记录数据变化的一维数组； (C) N记录每次数据变化的单个数据； (D)n是记录数据变化的一维数组。3、将半径为r的球体(密度 )置入水中，球体将浮出水面一定高度。程序如下：function h,Rou=highNu(r)if nargin=0,r=10;endRou=0.3:0.1:1;N=length(Rou);for k=1:N
38、 rouk=Rou(k); P=1,-3*r,0,4*r3*rouk; x=roots(P); II=find(x0); h(k)=2*r-x(II);end(1)在MATLAB命令窗口省略输入调用函数highNu将显示( ) (A) 球体浮出水面的高度数据； (B) 球体的8个不同的密度数据； (C) 球体沉入水下的深度数据； (D) 深度数据和密度数据。(2) 程序中变量x存入如下方程的根( ) (A) ； (B) ； (C) ； (D) 4、一阶常微分方程确定一个平面向量场，初值条件确定了向量场中一条曲线。程序如下：x,y=meshgrid(0:.25:6,0:.05:2);k=y.*(
39、1-y);d=sqrt(1+k.2);px=1./d;py=k./d;quiver(x,y,px,py),hold onu=dsolve(Du=u*(1-u),u(0)=.2); v=dsolve(Dv=v*(1-v),v(0)=1.8);ezplot(u,0,6)ezplot(v,0,6)(1) 程序中所绘向量场对应的一阶常微分方程是( ) (A) ； (B) ； (C) ； (D) 。(2) 关于实验程序下面说法错误的是( ) (A) 程序中第一个初值条件所对应的解曲线在图1中上方； (B) 程序中第二个初值条件所对应的解曲线在图1中上方； (C) 程序绘图原理是根据每一点处曲线切线的单位向量绘图； (D)