2012年数学实验matlaB复习题
2012年数学实验复习题
一、选择题
1、三阶幻方又称为九宫图，提取三阶幻方矩阵对角元并构造对角阵用( C ) (A) diag(magic(3)); (B) diag(magic);
(C) diag(diag(magic(3))); (D) diag(diag(magic))。
2、data=rand(1000,2);x=data(:,1);y=data(:,2);II=find(y
3、MATLAB 计算二项分布随机变量分布律的方法是( C )
(A) binocdf(x,n,p); (B) normpdf(x,mu,s); (C)binopdf(x,n,p); (D) binornd(x,n,p)。 4、MATLAB 命令syms e2;f=sqrt(1-e2*cos(t)^2);S=int(f,t,0,pi/2)功能是(D) (A) 计算f(x)在[0,pi/2]上的积分； (B) 计算f(t)不定积分符号结果； (C) 计算f(x)积分的数值结果； (D) 计算f(t)定积分的符号结果。 4、y=dsolve(‘Dy=1/(1+x^2)-2*y^2’,’y(0)=0’,’x’);ezplot(y)的功能是( A )
(A) 求微分方程特解并绘图； (B) 解代数方程 (C) 求定积分； (D)求微分方程通解。 6、X=10000 ;0.5*asin(9.8*X/(515^2))的功能是计算关于抛射体问题的(A)
(A) 十公里发射角； (B) 十公里飞行时间； (C)最大飞行时间； (D)最大射程。 7、theta=linspace(0,2*pi,100) ;r=cos(4*theta) ;polar(theta,r,’k’)功能是(D)
(A) 绘四叶玫瑰线； (B)绘三叶玫瑰线； (C)绘心脏线； (D) 绘八叶玫瑰线。 8、MA TLAB 命令A=rand(5,5)；创建55()ij A a ′=，求5
1
max
||ij j
i a =?
用( A )
(A) max(sum(abs(A))); (B) max(sum(abs(A ’))); (C) max(sum(A))); (D) sum(max(A)); 9、MA TLAB 命令x=[1,2,4,5,9];mean(x)，的计算结果是( B ) (A) 4 (B) 4.2 (B) 4.5 (D) 21
10、MATLAB 命令x=rand(10,1)生成10个随机数，将它们从大到小排序，使用( C ) (A) y=sort(x);z=y(10:1); (B) [y,II]=sort(x);z=y(II); (C) y=sort(x);z=y(10:-1;1); (D) [y,II]=sort(x);z=x(II); 11、MATLAB 命令roots([1,0,0,-1])的功能是( D )
(A) 产生向量[1,0,0,1]； (B) 求方程310x +=的根； (C) 求多项式31x -的值 (D) 求方程310x -=的根。
12、MATLAB 命令A=magic(3)创建3阶幻方矩阵，求A 的特征值绝对值最小用( A ) (A) min(abs(eig(A))); (B) min(eig(abs(A))); (C)min(eig(A)); (D) min(abs(A));
13、命令factor()用于分解因式，syms x; f=4*x^3+9*x^2-30*x; factor(diff(f))的结果是( B ) (A) (x-1)*(2*x-5) (B) 6*(x-1)*(2*x+5) (C) 6*(x+1)*(2*x+5) (D) (x+1)*(2*x-5) 14、MATLAB 命令syms x; f=sin(x); V=pi*int(f*f,x,0,2*pi)功能是( C )
(A) 绘出函数f 在[0,2p ]图形； (B) 计算函数f 在[0,2p ]的积分； (C) 计算旋转曲面所围的体积； (D) 计算旋转曲面的表面积。
15、十二属相为“鼠牛虎兔龙蛇马羊猴鸡狗猪”，命令k=mod(2008,12)+1的结果是( D )
(A) k指向第二动物牛；(B) k指向第三动物虎；
(C) k指向第四动物兔；(D) k指向第五动物龙。
16．下面有关MATLAB变量名和函数名的说法中，错误的说法是( D )
(A) 变量名的第一个字符必须是一个英文字母
(B) 变量名可由英文字母、数字和下划线混合组成
(C) 变量名不得包含空格和标点，但可以有下连字符
(D) 变量名和函数名对于英文的大小使用没有区别
17、在MATLAB命令窗口中，键入命令syms x；int(x*sin(x))。结果是（ A）
（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*x
18、在MATLAB命令窗口中键入命令A=[1 4 2;3 1 2;6 1 5];det(A(1:2,2:3).*A(2:3,2:3))。结果是（C ）
（A）ans= -143 （B）ans= 60 （C）ans= 36 （D）ans= -19
19、MATLAB命令x = 3: 2: 100 将创建等差数列，该数列是（ C）
（A）以3为初值的98个数,；（B）以100为终值的98的个数；
（C）以99为终值的49个数；（D）以3为初值的97个数。
20、在MATLAB命令窗口中输入命令data=[4 1 2 3 1 3 1 4 2 4];y=hist(data,4)，结果是（ B）（A）y= 4 1 2 3；（B）y=3 2 3 2；（C）y= 3 2 2 3 ；（D）y= 4 2 1 1
21、MATLAB语句[x,y]=meshgrid(-2:2) 的数据结果中（D ）
（A）x是行向量，y是列向量；（B）x是五行五列的矩阵；
（C）x是行元素相同的矩阵；（D）x是列向量相同的矩阵
22、MATLAB的语句colormap(0 0 1) （D ）
（A）将三维网面图确定为红色；（B）将三维网面图确定为绿色；
（C）将三维网面图确定为蓝色；（D）语句使用格式错误
23、火炮发射炮弹的初始速度和发射角为已知，由此可估算出炮弹在空中的飞行时间Tfly，使用语句Tspan=Tfly*(0:20)/20，将获得一些数据，下面不正确的说法是（D ）A）Tspan为包括发射时刻在内的炮弹在空间飞行的21个不同的飞行时刻；
B）Tspan中任意两个相邻数据之差的绝对值相等；
C）Tspan包含了21个数据，第一个数据为0，最后一个数据为Tfly；
D）Tspan是一个等差数列，公差为Tfly/21
24、北京和纽约的经度分别是：东经118和西经76，根据经度差计算时差用( D)
(A) fai1=118;fai2=-76;Dfai=(fai1+fai2)/24; (B) fai1=118;fai2=-76;Dfai=(fai1+fai2)/15;
(C) fai1=118;fai2=-76;Dfai=(fai1-fai2)/24; (D) fai1=118;fai2=-76;Dfai=(fai1-fai2)/15。
25、用MA TLAB随机产生60个1到365之间的正整数，应该使用下面的哪一条命令(B)
A）fix(365*rand(1,60))；B）1+fix(366*rand(1,60))；
C）1+fix(364*rand(1,60))；D）1+fix(365*rand(1,60))
二、程序阅读
1、3n+1问题反映一个数学猜想：对任一自然数n，按如下法则进行运算：若n为偶数，则将n除2，若n为奇数，则将n乘3加1。重复这种操作，结果终会为1。实验程序如下。function [k,N]=threeN(n)
if nargin==0,n=5;end
k=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
(1)在MA TLAB命令窗口中直接调用threeN运行结果为( )
(A)只显示k的最后数值为6；(B) 只显示k的最后数值5；
(C) 同时显示k和N的数据；(D) 仅显示N的所有数据。
(2)实验程序运行过程中( )
(A) 输入变量n不发生改变；(B)N是记录数据变化的一维数组；
(C) N记录每次数据变化的单个数据；(D)n是记录数据变化的一维数组。
2、关于“牟合方盖”的实验程序如下
h=2*pi/100;t=0:h:2*pi;
r=0:0.05:1;x=r'*cos(t);y=r'*sin(t);
z=sqrt(1-x.^2); %第三行
meshz(x,y,z),axis off
colormap([0 0 1])
view(-47,56),hold on
x1=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）绘圆柱面x2 + y2 = 1, x2 + z2= 1的交线及所围区域的边界曲面的上半部分。
（2）关于第三行语句错误的解释是（）
（A）z是矩形域上曲顶柱面高度值；（B）z是与y同型的矩阵；
（C）z是圆域上曲顶柱面高度值；（D）z是与x同型的矩阵
3、非负函数y =f(x)在有限区间上的图形为上半平面的一条曲线，曲线绕x轴旋转时，产
生以x为对称轴的旋转曲面，其体积
22
[()]
V f x dx
π
π
=?.
syms a b x
f=exp(a*x)*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 off
view(-17,54)
（1） 关于程序的功能确切的说法（ ） (A)计算曲线段 f(x)=exp(a x)sin(b x), 绕X 轴旋转的旋转曲面体积 (B)计算曲线段 f(x)=exp(-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) 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);
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 ）计算卫星运行轨道的平均速度 5、Viviani 体是圆柱体4/)2/(222R y R x ≤+-被球面2222R z y x =++所割立体。下面的数学实验程序功能是取R=2求体积上半部分，先利用符号处理重积分并转换为数值数据，再用蒙特卡罗方法计算体做对比。数学实验程序如下：
syms x y;
f=sqrt(4-x^2-y^2);
y1=sqrt(2*x-x^2);y2=sqrt(2*x-x^2); S1=int(f,y,y1,y2);S2=int(S1,x,0,2) V=double(S2) P=rand(10000,3);
X=2*P(:,1);Y=2*P(:,2);Z=2*P(:,3);
II=find((X-1).^2+Y .^2<=1&Z<=sqrt(4-X.^2-Y .^2));
图1 Vivinai 问题
02x π≤≤02x π≤≤02x π≤≤02x π≤≤()f x =()f x y f =y f =
V1=8*length(II)/10000
(1) 符号计算所用的积分公式是( )
(A) (B) ?
?-----=
2
2
22222
4x x x x dy y x dx
V
(C) ?
-+----=
2
2
222
2
4x x x
x dy y x V (D) ?
-+----=
2
2
22224x x x
x dx y x V
(2) 蒙特卡罗方法选用的随机点变化范围的立方体区域是( B ) (A) )}2,0(),2,0(),2,0(|),,{(∈∈∈=Ωz y x z y x ； (B) )}2,0(),1,1(),2,0(|),,{(∈-∈∈=Ωz y x z y x (C) }20,10),20|),,{(<<<<<<=Ωz y x z y x (D) }20,20),20|),,{(<<<<<<=Ωz y x z y x
6、某厂生产两种产品，产一吨甲产品用A 资源3吨、B 资源4m 3；产一吨乙产品用A 资源2吨，B 资源6m 3，C 资源7个单位。一吨甲产品和乙产品分别价值7万元和5万元，三种资源限制分别为90吨、200m 3和210个单位。生产两种产品使总价值最高的生产方案可用数学实验程序计算。
C=[-7,-5];A=[3 2;4 6;0 7];b=[90;200;210]; Aeq=[];Beq=[];
e0=[0,0];e1=[inf,inf];
[x,fval]=linprog(C,A,b,Aeq,beq,e0,e1); (1) 程序中变量C 表示( )
(A) 目标函数系数； (B) 等式约束系数； (C) 不等式约束系数； (D) 等式约束常向量 (2) 程序中变量A 表示( )
(A) 等式约束矩阵； (B) 不等式约束矩阵； (C) 决策变量的值； (D) 目标函数的最大值
三、程序填空
1、中国农历60年一大轮回，按天干“甲乙丙丁戊已庚辛壬癸”和地支“子丑寅卯辰巳午未申酉戍亥”循环排列而成。已知2009年是农历已丑年，通过简单计算可以找出年份与天干/地支对应的规律。下面数学实验程序对输入年份，计算并输出字符串农历纪年。填空完善程序。
input('input year:=');
if nargin==0, year=2009;end S1=’ 甲乙丙丁戊已庚辛壬癸’; S2=’子丑寅卯辰巳午未申酉戍亥’;
k1= mod(S1-4,10)+1 ①; %定位天干序数%2010庚 s1=S1(k1); k2= mod(S2-4,12) +1 ②; %定位地支序数 s2=S2(k2);
calendar=strcat(int2str(year),’年是’,s1,s2,’年’)
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 ①
end
if k disp('不是素数') else disp ② end 3．汽车紧急刹车问题数据拟合实验 V 20 25 30 35 40 45 50 55 60 65 70 T 20 28 41 53 72 93 118 149 182 221 266 V 表示刹车时汽车行驶速度(英里/小时),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,'*',v,T2) P3=polyfit(v,T,3);T3=polyval(P3,v); R3= ② figure(3),plot(v,T,'*',v,T3) 4、红、绿两队从相距100公里的地点同时出发相向行军。红队速度为10(公里/小时)，绿队速度为8(公里/小时)。开始时，通讯员骑摩托从红队出发为行进中的两队传递消息。摩托车的速度为60(公里/小时)往返于两队之间。每遇一队，立即回驶向另一队。当两队距离小于0.2公里时，摩托车停止，下面数学实验程序模拟计算摩托车跑了多少趟。请填空完善程序。 function k=moto(A,B) if nargin==0,A=0;B=100;end va=10;vb=8;vc=60; f=1;k=0; while (B-A)>0.2 if f==1 tk=(B-A)/(vb+vc); else tk= (B-A)/(vc+va); ①; %计算A 与C 相遇时间 end A= A+va*tk ②; %计算A 点位置 B= B-vb*tk ③; %计算B 点位置 f=-f; k=k+1; end 5、二阶正交矩阵作用于某一向量时，其效果是将该向量旋转，旋转解为α(逆时针旋转为正)。把一个以原点为中心的正三角形旋转50/π，并缩小90%，迭代33次创建图3。完成程序填空： bata=[1/2;7/6;11/6;15/6]*pi; x=cos(bata);y=sin(bata); line(x,y) xy=[x,y]; alfa=pi/50; A=[cos(alfa),-sin(alfa);sin(alfa),cos(alfa)]; for k=1:33 xy= 0.9*xy*A' ①; x=xy(:,1); y= xy(:,2) ②; line(x,y) end 6、长征三号甲运载火箭提供给探月卫星的初始速度不足以将卫星送往月球轨道。为提高到奔月速度，中国航天工程师使用了卫星变轨技术。数学实验程序根据变轨中轨道周期和近地点距离数据，利用开普列第二定律模拟计算计算卫星飞行的最大速度。填空完善下面实验程序。 R=6378;Time=[16,15.63,23.3,50.5,225]*3600; h=[200,600,600,600,600];H=[51000,51000,71000,128000,370000]; a=(h+H+2*R)/2; c= (H-h)/2; ① b=sqrt(a.*a-c.*c); S= pi*a.*b./Time; ② Vmax=2*S./(R+h) 8、为了进入地月转移轨道，嫦娥一号卫星进行了四次变轨调速度。第一次变轨从16小时初始轨道进入16小时轨道，第二次卫星进入24小时轨道，第三次卫星进入48小时轨道，第四次卫星进入116小时地月转移轨道。上面小时数并不是准确轨道周期，变轨目的是将速度从10.3(km/s)逐渐提高到约10.9(km/s)。下面数学实验程序是在区间[10.3,10.9]上插入线性等分点，即每个轨道的最大速度以等差数列出现，然后近似计算出每个轨道的周期参数。填空完善程序。 function satel1() R=6378; h=[200,600,600,600,600]; H=[51000,51000,71000,128000,370000]; a=(h+H+2*R)/2; c=(H-h)/2; b= sqrt(a.^2-c.^2) ①; %计算短半轴数据 E2=(c./a).^2; L=2*pi*a.*(1-E2/4-3*E2.^2/64) format bank Vmax=linspace(10.3,10.9,5) S= 0.5*Vmax.*(R+h) ②; %根据最大速度计算每秒钟扫过的面积 Times=a.*b.*pi./S; myTimes=Times/3600 图3 旋转三角形 MATLAB数学实验练习题 ————————————————————————————————作者：————————————————————————————————日期： “MATLAB”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> solve('exp(x)-3*x^2',0) > ezplot('exp(x)-3*x^2') >> grid on ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 2、求下列方程的根。 1) 5510x x ++= a=solve('x^5+5*x+1',0);a=vpa(a,6) a = 1.10447+1.05983*i -1.00450+1.06095*i -.199936 -1.00450-1.06095*i 1.10447-1.05983*i 2）1sin 02 x x - =至 少三个根 >> fzero('x*sin(x)-1/2', 3) ans = 2.9726 >> fzero('x*sin(x)-1/2',-3) ans = -2.9726 >> fzero('x*sin(x)-1/2',0) ans = -0.7408 3）2sin cos 0x x x -= 所有根 >> fzero('sin(x)*cos(x)-x^2',0) ans = >> fzero('sin(x)*cos(x)-x^2',0.6) ans = 0.7022 3、求解下列各题： 1）30sin lim x x x x ->- >> sym x; >> limit((x-sin(x))/x^3,x,0) ans = 1/6 2） (10)cos ,x y e x y =求 >> sym x >> diff(exp(x)*cos(x),10) ans = (-32)*exp(x)*sin(x) 第一次练习题 1． 求 32 =-x e x 的所有根。（先画图后求解） 2． 求下列方程的根。 1) 0155 =++x x 2） 至少三个根）(0 2 1s i n =- x x 3） 所有根0 c o s s i n 2 =-x x x 3． 求解下列各题： 1） 3 sin lim x x x x ->- 2） ) 10(, cos y x e y x 求= 3） ?+dx x x 2 4 425 4） ）（最高次幂为 展开在将801=+x x 5） )2() 3(1sin y e y x 求 = 4． 求矩阵 ???? ? ? ?--=31 4020 112 A 的逆矩阵1 -A 及特征值和特征向量。 5． 已知,21)(2 2 2)(σ μσ π-- = x e x f 分别在下列条件下画出)(x f 的图形： ）； （在同一坐标系上作图 ，，＝时＝、）；（在同一坐标系上作图，－，＝时、421,0)2(110,1)1(σμμσ=、 6． 画 （1）202004 cos sin ≤≤≤≤???? ?? ? ===u t t z t u y t u x (2) 30,30)sin(≤≤≤≤=y x xy z （3）π π2020sin ) cos 3()cos()cos 3()sin(≤≤≤≤?? ? ??=+=+=u t u z u t y u t x 的图（第6题只要写出程序）. 7绘制曲线x x x sa )sin()(=，其中]10,10[ππ-∈x 。（注意：0=x 处需要特别处理。） 8．作出函数x e x f x cos )(-=的图形；求出方程0=)(x f 在],[020-的所有根；令 n x 为从0向左依次排列的方程的根，输出n n x x --1 ，并指出?)(lim =--∞ >-n n n x x 1 9. 把x cos 展开到2，4，6项，并作出的x cos 和各展开式的图形；并指出用展开式逼 近x cos 的情形。 10. 请分别写出用for 和while 循环语句计算63 263 2 2212+++== ∑ = i i K 的程序。此外， 还请写出一种避免循环的计算程序。 11. 对于0>x ，求1 20 11122 +∞ =∑ ? ? ? ??+-+k k x x k 。（提示：理论结果为x ln ） 第二次练习题 1、 设????? =+=+32/)7(1 1 x x x x n n n ，数列}{n x 是否收敛？若收敛，其值为多少？精确到6位 有效数字。 用两种方法 2、设 ,13 12 11p p p n n x + ++ += }{n x 是否收敛？若收敛，其值为多少？精确到17 位有效数字。 注：学号为单号的取7=p ，学号为双号的取.8=p 3、38P 问 题2 4、编程找出 5,1000+=≤b c c 的所有勾股数，并问：能否利用通项表示 },,{c b a ? 5、编程找出不定方程 )35000(122 2 <-=-y y x 的所有正整数解。（学号为单号 复习题 1、写出3 2、i ｎv(A)表示Ａ的逆矩阵； 3、在命令窗口健入 clc,４、在命令窗口健入clea 5、在命令窗口健入６、x=-1:０.2:17、det(A)表示计算Ａ的行列式的值；8、三种插值方法：拉格朗日多项式插值,分段线性插值，三次样条插值。 9、若A=123456789?? ????????,则fliplr （Ａ）＝321654987?????????? Ａ-３=210123456--??????????A ．＾２=149162536496481?????????? tril(A)=100450789?????????? tri ｕ（Ａ，－1)=123456089??????????ｄiag(A ）＝100050009?????????? Ａ(:,2),＝2 58A(3,:)=369 10、nor ｍcd ｆ（1，1，2)=0.5%正态分布mu=1,s ｉgm ａ＝２,x ＝1处的概率 ｅ45(@ｆ,［ａ,b ］,x0),中参数的涵义是@fun 是求解方程的函数M 文 件,[a,b ］是输入向量即自变量的范围a 为初值,ｘ0为函数的初值，t 为输出指定的[ａ,b],x 为函数值 １5、写出下列命令的功能：te ｘt （1,2,‘ｙ=s ｉｎ(x)’ hold on 16fun ｃｔion 开头; 17 ,4) ３，４） 21、设x 是一向量，则）的功能是作出将Ｘ十等分的直方图 ２2、interp １(［1,2,3]，[3,４,5]，2.5) Ans=4.５ 23、建立一阶微分方程组? ??+='-='y x t y y x t x 34)(3)(2 的函数M 文件。(做不出来） 二、写出运行结果: １、>＞ｅy ｅ（3，4)=1000 01000010 2、>>s ｉze([1,2,3])=1;3 3、设b=ro ｕnd （unifrnd(-５，5,１，4)),则=3 ５ 2 -5 >>[x,m]=ｍin(b)；x ＝－5；m=4 ，［x,n ］=sort(b ） －５ 2 3 ５ 4 3 １ 2 ｍea ｎ(ｂ）=１.2５，m ｅdian(b)＝2.5,range(ｂ）=10 ４、向量ｂ如上题,则 ＞>ａn ｙ(b),all(ｂ<2),ａｌl(b<6） Ans ＝１ 0 １ ５、＞>[5 6;7 8]>[7 8；5 6]=00 11 6、若1234B ??=???? ,则 7、>>ｄｉag(d ｉag （B ）)=10 04 8、＞>[4:-2:１].＊[-1,６]=-4 12 9、>>acos(0.5),a ｔan(1） ans ＝ 1.６5９8 ans= 北京交通大学海滨学院考试试题 课程名称：数学实验2010-2011第一学期出题教师：数学组适用专业: 09机械, 物流, 土木, 自动化 班级：学号：姓名： 选做题目序号： 1.一对刚出生的幼兔经过一个月可以长成成兔, 成兔再经过一个月后可以 繁殖出一对幼兔. 如果不计算兔子的死亡数, 请用Matlab程序给出在未来24个月中每个月的兔子对数。 解: 由题意每月的成兔与幼兔的数量如下表所示： 1 2 3 4 5 6 ··· 成兔0 1 1 2 3 5··· 幼兔 1 0 1 1 2 3··· 运用Matlab程序： x=zeros(1,24); x(1)=1;x(2)=1； for i=2:24 x(i+1)=x(i)+x(i-1); end x 结果为x = 1 1 2 3 5 8 13 21 3 4 5 5 89 144 233 377 610 987 1597 2584 4181 6765 1094 6 7711 2865 7 46368 2.定积分的过程可以分为分割、求和、取极限三部分, 以1 x e dx 为例, 利用 已学过的Matlab 命令, 通过作图演示计算积分的过程, 并与使用命令int() 直接积分的结果进行比较. 解：根据求积分的过程，我们先对区间[0,1]进行n 等分， 然后针对函数x e 取和，取和的形式为10 1 i n x i e e dx n ξ=≈ ∑ ? ，其中1[ ,]i i i n n ξ-?。这里取i ξ为区间的右端点，则当10n =时，1 x e dx ?可用10 101 1.805610 i i e ==∑ 来近似计算， 当10n =0时，100 100 1 01 =1.7269100 i x i e e dx =≈ ∑?，当10n =000时，10000 10000 1 1 =1.718410000 i x i e e dx =≈ ∑ ?. 示意图如下图，Matlab 命令如下： x=linspace (0,1,21); y=exp(x); y1=y(1:20); s1=sum(y1)/20 y2=y(2:21); s2=sum(y2)/20 plot(x,y); hold on for i=1:20 fill([x(i),x(i+1),x(i+1),x(i),x(i)],[0,0,y(i),y(i),0],'b') end syms k;symsum(exp(k/10)/10,k,1,10)；%n=10 symsum(exp(k/100)/100,k,1,100)；%n=100 symsum(exp(k/10000)/10000,k,1,10000)；%n=10000 “”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> ('(x)-3*x^2',0) = -2*(-1/6*3^(1/2)) -2*(-11/6*3^(1/2)) -2*(1/6*3^(1/2)) 3、求解下列各题： 1）30 sin lim x x x x ->- >> x; >> (((x))^3) = 1/6 2） (10)cos ,x y e x y =求 >> x; >> ((x)*(x),10) = (-32)*(x)*(x) 3）2 1/2 0(17x e dx ?精确到位有效数字） >> x; >> ((((x^2),0,1/2)),17) = 0.54498710418362222 4）4 2 254x dx x +? >> x; >> (x^4/(25^2)) = 125*(5) - 25*x + x^3/3 5）求由参数方程arctan x y t ??=? =??dy dx 与二阶导 数22 d y dx 。 >> t; >> ((1^2))(t); >> ()() = 1 6）设函数(x)由方程e所确定，求y′(x)。>> x y; *(y)(1); >> ()() = (x + (y)) 7） sin2 x e xdx +∞- ? >> x; >> ()*(2*x); >> (y,0) = 2/5 8） 08x =展开（最高次幂为） >> x (1); taylor(f,0,9) = - (429*x^8)/32768 + (33*x^7)/2048 - (21*x^6)/1024 + (7*x^5)/256 - (5*x^4)/128 + x^3/16 - x^2/8 + 2 + 1 9） 1sin (3)(2)x y e y =求 >> x y; >> ((1)); >> ((y,3),2) = 清华大学2002至2003学年第二学期数学实验期末考试试题A 数学实验试题 2003.6.22 上午 (A卷；90分钟) 一. 某两个地区上半年6个月的降雨量数据如下（单位：mm）： 月份123456 地区A259946337054 地区B105030204530 在90%的置信水平下，给出A地区的月降雨量的置信区 间： 在90%的置信水平下，A地区的月降雨量是否不小于70（mm）？ 在90%的置信水平下，A、B地区的月降雨量是否相同？ A地区某条河流上半年6个月对应的径流量数据如下（单位：m3）：110，184，145，122，165，143。该河流的径流量y与当地的降雨量x的线性回归方程为；若当地降雨量为55mm，该河流的径流量的预测区间为（置信水平取90%）。 答案：（程序略） (1) [32.35，76.65] (2) 是 (3) 否 (4) y=91.12+0.9857x (5) [130.9,159.7] 二．（10分） （1）（每空1分）给定矩阵，如果在可行域上考虑线性函数，其中，那么的最小值是，最小点为；最大值是，最大点为。 （每空2分）给定矩阵，，考虑二次规划问题，其最优解为，（2） 最优值为，在最优点处起作用约束 为 。 答案：（1）最小值为11/5，最大值为7/2，最小点为（0，2/5，9/5），最大点为（1/2，0，3/2）。 （2）最优解为（2.5556，1.4444），最优值为–1.0778e+001，其作用约束为。 三．（10分）对线性方程组：，其中Ａ＝，b= （3分）当时，用高斯—赛德尔迭代法求解。取初值为，写出迭代第4步的结果=____________________。 （4分）当时，用Jacobi 迭代法求解是否收敛？__________ , 理由是_________________________________________________ 。 （3分）求最大的c, 使得对任意的，用高斯—赛德尔迭代法求解一定收敛，则c应为__________。 答案：（1）x = [ -1.0566 1.0771 2.9897] 注意:在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上）. 第一次练习题 1.求解下列各题： 1)30sin lim x mx mx x ->- 2)(4)cos ,1000.0=x mx y e y 求 3）21/2 0mx e dx ?(求近似值，可以先用inline 定义被积函数，然后用quad 命令) 4）4 224x dx m x +? 5 0x =展开(最高次幂为8). 2.对矩阵21102041A m -?? ?= ? ?-?? ，分别求逆矩阵，特征值，特征向量，行列式，并求矩阵,P D （D 是对角矩阵），使得1A PDP -=。 3. 已知2 1(),()2f x e x μσ=--分别在下列条件下画出)(x f 的图形： (1)/600m σ=,μ分别为0,1,1-(在同一坐标系上作图); (2)0μ=,σ分别为1,2,4,/100m (在同一坐标系上作图). 4.画 （1）sin 020cos 02100x u t t y u t u t z m ??=≤≤?=?≤≤??=? (2) sin()03,03z mxy x y =≤≤≤≤ （3）sin()(/100cos )02cos()(/100cos )02sin x t m u t y t m u u z u π π=+?≤≤?=+?≤≤?=? 的图（第4题只要写出程序）. 5．对于方程50.10200 m x x --=，先画出左边的函数在合适的区间上的图形，借助于软件中的方程求根的命令求出所有的实根,找出函数的单调区间,结合高等数学的知识说明函数为什么在这些区间上是单调的,以及该方程确实只有你求出的这些实根。最后写出你做此题的体会. 第二次练习题 判断迭代收敛速度的程序 x0=1;stopc=1;eps=10^(-8);a=1;c=1;b=2*c;d=a;k=0; f=inline('(a*x+b)/(c*x+d)'); kmax=100; while stopc>eps&k 2011年数学实验复习题 一、选择题 1、三阶幻方又称为九宫图，提取三阶幻方矩阵对角元并构造对角阵用( C ) (A) diag(magic(3)); (B) diag(magic); (C) diag(diag(magic(3))); (D) diag(diag(magic))。 2、data=rand(1000,2);x=data(:,1);y=data(:,2);II=find(y 一元函数微分学 实验1 一元函数的图形（基础实验） 实验目的 通过图形加深对函数及其性质的认识与理解, 掌握运用函数的图形来观察和分析 函数的有关特性与变化趋势的方法,建立数形结合的思想; 掌握用Matlab 作平面曲线图性的方法与技巧. 初等函数的图形 2 作出函数x y tan =和x y cot =的图形观察其周期性和变化趋势. 解：程序代码： >> x=linspace(0,2*pi,600); t=sin(x)./(cos(x)+eps); plot(x,t);title('tan(x)');axis ([0,2*pi,-50,50]); 图象： 程序代码： >> x=linspace(0,2*pi,100); ct=cos(x)./(sin(x)+eps); plot(x,ct);title('cot(x)');axis ([0,2*pi,-50,50]); 图象： cot(x) 4在区间]1,1[-画出函数x y 1 sin =的图形. 解：程序代码： >> x=linspace(-1,1,10000); y=sin(1./x); plot(x,y); axis([-1,1,-2,2]) 图象： 二维参数方程作图 6画出参数方程? ??==t t t y t t t x 3cos sin )(5cos cos )(的图形： 解：程序代码： >> t=linspace(0,2*pi,100); plot(cos(t).*cos(5*t),sin(t).*cos(3*t)); 图象： 极坐标方程作图 8 作出极坐标方程为10/t e r =的对数螺线的图形. 解：程序代码： >> t=0:0.01:2*pi; r=exp(t/10); polar(log(t+eps),log(r+eps)); 图象： 90270 分段函数作图 10 作出符号函数x y sgn =的图形. 解： 程序代码： >> x=linspace(-100,100,10000); y=sign(x); plot(x,y); axis([-100 100 -2 2]); Matlab 数学实验 实验一 插值与拟合 实验内容： 预备知识：编制计算拉格朗日插值的M 文件。 1. 选择一些函数，在n 个节点上（n 不要太大，如5 ~ 11）用拉格朗日、分段线性、三次样条三种插值方法，计算m 个插值点的函数值（m 要适中，如50~100）。通过数值和图形输出，将三种插值结果与精确值进行比较。适当增加n ，再做比较，由此作初步分析。下列函数任选一种。 （1）、 ;20,sin π≤≤=x x y （2）、;11,)1(2/12≤≤--=x x y （3）、;22,c o s 10 ≤≤-=x x y （4）、22),exp(2≤≤--=x x y 2.用电压V=10伏的电池给电容器充电，电容器上t 时刻的电压为 ) (0)()(t e V V V t v ---=，其中0V 是电容器的初始电压，τ是充电常数。试由下面 一组t ，V 数据确定0V 和τ。 实验二 常微分方程数值解试验 实验目的： 1. 用MATLAB 软件求解微分方程，掌握Euler 方法和龙格-库塔方法； 2. 掌握用微分方程模型解决简化的实际问题。 实验内容： 实验三地图问题 1.下图是一个国家的地图，为了计算出它的国土面积，首先对地图作如下测量： 以由西向东方向为x轴，由南到北方向为y轴，选择方便的原点，并将从最西边界点到最东边界点在x轴上的区间适当地划分为若干段，在每个分点的y方向测出南边界点和北边界点的y坐标y1和y2，这样就得到了表中的数据（单位mm）。 根据地图的比例我们知道18mm相当于40km，试由测量数据计算该国土 的近似面积，并与它的精确值41288km2比较。 191 《数学实验》模拟试题一 一、单项选择题 1．符号计算与一般数值计算有很大区别，它得到准确的符号表达式。在MA TLAB 命令窗口中键入命令syms x ，y1=sqrt(x);y2=x^2；int(y1-y2,x,0,1)，屏幕显示的结果是 （A ）y1 =x^(1/2) （B ）ans= 2/3； （C ）y2 =x^2； （D ）ans= 1/3 2．在MA TLAB 命令窗口中键入命令A=[1 4 2;3 1 2;6 1 5];det(A(1:2,2:3).*A(1:2,2:3))。结果是 （A ）ans= -143 （B ）ans= 60 （C ）ans= -16 （D ）ans= -19 3．设n 阶方阵A 的特征值为：i λ (i=1，2，…，n )，称||max )(i i A λρ=为矩阵A 的谱半径, 则下列MA TLAB 求谱半径命令是 （A ）max(abs(eig(A)))； （B ）abs(max(eig(A)))； （C ）max(norm(eig(A)))； （D ）norm(max(eig(A))) 4．MA TLAB 系统运行时，内存中有包括X 和Y 在内的多个变量（数据），要删除所有变量（数据），应该使用的命令是 （A ）clear ； （B ）clc ； （C ）home ； （D ）clear X Y 5．用赋值语句给定x 数据，计算3ln +)2+3sin(72e x 对应的MA TLAB 表达式是 （A ）sqrt(7*sin(3+2*x)+exp(2)*log(3)) （B ）sqrt(7sin(3+2x)+exp(2)log(3)) （C ）sqr(7*sin(3+2*x)+e^2*log(3)) （D ）sqr(7sin(3+2x)+ e^2 log(3)) 6．在MA TLAB 命令窗口中输入命令data=[4 1 2 3 1 3 1 3 2 4];y=hist (data,4)，结果是 （A ） y= 4 1 2 3； （B ）y=3 2 3 2； （C ）y= 1 3 2 4 ； （D ）y= 4 2 1 1 7．在MA TLAB 命令窗口中键入A=magic(6); B=A(2:5,1:2:5) 将得到矩阵B ，B 是 （A ）2行5列矩阵；（B ）4行两列矩阵；（C ）4行3列矩阵；（D ）4行5列矩阵 8．MA TLAB 绘三维曲面需要构建网格数据，语句[x,y]=meshgrid(-2:2)返回数据中 （A ）x 是行向量，y 是列向量； （B ）x 是列向量，y 是行向量； （C ）x 是行元素相同的矩阵； （D ）x 是列元素相同矩阵 9．下面有关MA TLAB 函数的说法，哪一个是错误的 （A ）函数文件的第一行必须由function 开始，并有返回参数，函数名和输入参数； （B ）MA TLAB 的函数可以有多个返回参数和多个输入参数； （C ）如果函数文件内有多个函数，则只有第一个函数可以供外部调用； （D ）在函数中可以用nargin 检测用户调用函数时的输出参数个数 10．将带小数的实数处理为整数称为取整，常用四种取整法则是：向正无穷大方向取 整、向负无穷大方向取整、向零方向取整和四舍五入取整。MA TLAB 提供了如下四个取整函数，若a = -1.4，对a 取整的结果是 -1，则不应该选用下面哪个函数。 （A ）floor ； （B ）round ； （C ）ceil ； （D ）fix ； 二、程序阅读理解 1．如果存在一条曲线L 与曲线簇中每一条曲线相切，则称L 为曲线簇的包络。 简单直线簇的实验程序如下 N=input('input N:='); x=[0:N]/N;y=1-x; Matlab 数学实验 实验一 插值与拟合 实验内容： 预备知识：编制计算拉格朗日插值的M 文件。 1. 选择一些函数，在n 个节点上（n 不要太大，如5 ~ 11）用拉格朗日、分段线性、三次样条三种插值方法，计算m 个插值点的函数值（m 要适中，如50~100）。通过数值和图形输出，将三种插值结果与精确值进行比较。适当增加n ，再做比较，由此作初步分析。下列函数任选一种。 （1）、 ;20,sin π≤≤=x x y （2）、;11,)1(2/12≤≤--=x x y （3）、;22,cos 10≤≤-=x x y （4）、22),ex p(2≤≤--=x x y 2.用电压V=10伏的电池给电容器充电，电容器上t 时刻的电压为 ) (0)()(τt e V V V t v ---=，其中0V 是电容器的初始电压，τ是充电常数。试由下面 一组t ，V 数据确定0V 和τ。 实验二 常微分方程数值解试验 实验目的： 1. 用MATLAB 软件求解微分方程，掌握Euler 方法和龙格-库塔方法； 2. 掌握用微分方程模型解决简化的实际问题。 实验内容： 实验三 地图问题 1. 下图是一个国家的地图，为了计算出它的国土面积，首先对地图作如下测量：以由西向东方向为x 轴，由南到北方向为y 轴，选择方便的原点，并将从最西边界点到最东边界点在x 轴上的区间适当地划分为若干段，在每个分点的y 方向测出南边界点和北边界点的y 坐标y1和y2，这样就得到了表中的数据（单位mm ）。 根据地图的比例我们知道18mm相当于40km，试由测量数据计算该国土 2 实验四狼追兔问题 狼猎兔问题是欧洲文艺复兴时代的著名人物达.芬奇提出的一个数学问题。当一个兔子正在它的洞穴南面60码处觅食时，一只恶狼出现在兔子正东的100码处。当两只动物同时发现对方以后，兔子奔向自己的洞穴，狼以快于兔子一倍的速度紧追兔子不放。狼在追赶过程中所形成的轨迹就是追击曲线。狼是否会在兔子跑回洞穴之前追赶上兔子？ 为了研究狼是否能够追上兔子，可以先考虑求出狼追兔子形成的追击曲线，然后根据曲线来确定狼是否能够追上兔子。 试验五：开放式基金的投资问题 某开放式基金现有总额为15亿元的资金可用于投资，目前共有8个项目可供投资者选择。每个项目可以重复投资，根据专家经验，对每个项目投资总额不能太高，且有个上限。这些项目所需要的投资额已经知道，在一般情况下，投资一年后各项目所得利润也可估计出来（见表一）， 表一: 投资项目所需资金及预计一年后所得利润(单位：万元) “MATLAB”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> solve('exp(x)-3*x^2',0) ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 2、求下列方程的根。 1) 5510x x ++= a=solve('x^5+5*x+1',0);a=vpa(a,6) 1.10447+1.05983*i -1.00450+1.06095*i -.199936 -1.00450-1.06095*i 1.10447-1.05983*i 2） 1 sin0 2 x x-=至少三个根 >> fzero('x*sin(x)-1/2', 3) ans = 2.9726 >> fzero('x*sin(x)-1/2',-3) ans = -2.9726 >> fzero('x*sin(x)-1/2',0) ans = -0.7408 3）2sin cos 0x x x -= 所有根 >> fzero('sin(x)*cos(x)-x^2',0) ans = >> fzero('sin(x)*cos(x)-x^2',0.6) ans = 0.7022 3、求解下列各题： 1）30sin lim x x x x ->- >> sym x; >> limit((x-sin(x))/x^3) ans = 1/6 2） (10)cos ,x y e x y =求 >> sym x; >> diff(exp(x)*cos(x),10) ans = 数学实验 练习2.1 画出下列常见曲线的图形。（其中a=1,b=2,c=3）1、立方抛物线3x y= 解：x=-5:0.1:0;y=(-x).^(1/3); y=-y; x=0:0.1:5; y=[y,x.^(1/3)]; x=[-5:0.1:0,0:0.1:5]; plot(x,y) 2、高斯曲线2x e = y- 解：fplot('exp(-x.^2)',[-5,5]) 3、笛卡儿曲线)3(13,13332 2 2 axy y x t at y t at x =++=+= 解：ezplot('x.^3+y.^3-3*x*y',[-5,5]) x y x.3+y.3-3 x y = 0 或t=-5:0.1:5; x=3*t./(1+t.^2); y=3*t.^2./(1+t.^2); plot(x,y) 4、蔓叶线)(1,13 2 23 2 2x a x y t at y t at x -=+=+= 解：ezplot('y.^2-x.^3/(1-x)',[-5,5]) x y y.2-x.3/(1-x) = 0 或t=-5:0.1:5; x=t.^2./(1+t.^2); y=t.^3./(1+t.^2); plot(x,y) 5、摆线)cos 1(),sin (t b y t t a x -=-= 解：t=0:0.1:2*pi; x=t-sin(t); y=2*(1-cos(t)); plot(x,y) 6、星形线)(sin ,cos 3 23 23 233a y x t a y t a x =+== 解：t=0:0.1:2*pi; x=cos(t).^3; y=sin(t).^3; 2015-2016数学实验练习题 一、选择题 1.清除Matlab工作空间（wordspace）变量的命令是（B ） A. clc B. clear C. clf D.delete 2. 清除当前屏幕上显示的所有内容，但不清除工作空间中的数据的命令是（ A ） A. clc B. clear C. clf D.delete 3. 用来清除图形的命令（ C ） A. clc B. clear C. clf D.delete 4. 在MATLAB程序中，使命令行不显示运算结果的符号是（ A ） A. ; B. % C. # D. & 5. 在MATLAB程序中，可以将某行表示为注释行的符号是（ B ） A. ; B. % C. # D. & 6.在循环结构中跳出循环，执行循环后面代码的命令为 ( B ) A. return B. break C. continue D. Keyboard 7.在循环结构中跳出循环，但继续下次循环的命令为（ C ） A. return B. break C. continue D. Keyboard 8. MATLAB中用于声明全局变量的关键字是( C ) A. inf B. syms C. global D. function 9. 用户可以通过下面哪项获得指令的使用说明（ A ） A. help B. load C. demo D. lookfor 10．在MATLAB命令窗口中键入命令S=zoros(3)；可生成一个三行三列的零矩阵，如果省略了变量名S，MATLAB表现计算结果将用下面的哪一变量名做缺省变量名（ A ） A. ans; B. pi; C. NaN; D. Eps. 11. 9/0的结果是（ B ） A. NAN； B. Inf； C. eps； D. 0 12．在MATLAB中程序或语句的执行结果都可以用不同格式显示，将数据结果显示为分数形式，用下面哪一条命令语句（ D ） A. format long； B. format long e； C. format bank； D. fromat rat 13. 下列MATLAB命令中是构造1行3列的(-1,1)均匀分布随机矩阵的命令的是（D） 第一次练习 教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。 补充命令 vpa(x,n) 显示x 的n 位有效数字，教材102页 fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin lim x mx mx x →-与3 sin lim x mx mx x →∞- syms x limit((902*x-sin(902*x))/x^3) ans = 366935404/3 limit((902*x-sin(902*x))/x^3,inf) ans = 0 1.2 cos 1000 x mx y e =，求''y syms x diff(exp(x)*cos(902*x/1000),2) ans = (46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算 22 11 00 x y e dxdy +?? dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) ans = 2.1394 1.4 计算4 2 2 4x dx m x +? syms x int(x^4/(902^2+4*x^2)) ans = (91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求 syms x diff(exp(x)*cos(902*x),10) ans = -356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x) 1.6 0x =的泰勒展式(最高次幂为4). 电子科技大学二零零六至二零零七学年第二学期期末考试 《数学实验》课程考试题A 卷(120分钟) 考试形式：闭卷 考试日期：2007年7月11日 课程成绩构成：平时10分，期中0分，实验30分，期末60分 (本卷面成绩100) 一、单项选择题(共30分，每小题3分) 1．符号计算与一般数值计算有很大区别，它将得到准确的符号表达式。在MATLAB 命令窗口中键入命令syms x, y1=sqrt(x); y2=x^2; int(y1-y2,x,0,1)，屏幕显示的结果是( D ) (A) y1=x^(1/2) (B)ans=2/3; (C) y2=x^2 (D) ans=1/3 2．在MA TLAB 命令窗口中键入命令A=[1 4 2;3 1 2;6 1 5];det(A(1:2,2:3).*A(1:2,2:3))。结果是( B ) (A) ans= -143 (B) ans=60 (C)ans= -16 (D)ans= -19 3．设n 阶方阵A 的特征值为：),,2,1(n i i =λ，称||max )(i i A λρ=为矩阵A 的谱半径， 则下列MA TLAB 求谱半径 命令是( A ) (A) max(abs(eig(A))); (B) abs(max(eig(A))); (C)max(norm(eig(A))) (D) norm(max(eig(A))) 4．MA TLAB 系统运行时，内存中有包括X 和Y 在内的多个变量(数据)，要删除所有变量(数据)，该使用的命令是( A ) (A) clear (B) clc (C) home (D) clear X Y 5．用赋值语句给定x 数据，计算3ln )23sin(72 e x ++对应的MATLAB 表达式是( A ) (A) sqrt(7*sin(3+2*x)+exp(2)*log(3)) (B) sqrt(7sin(3+2x)+exp(2)log(3)) (C) sqr(7*sin(3+2*x)+exp(2)*log(3)) (D) sqr(7sin(3+2x)+exp(2)log(3)) 6．在MA TLAB 窗口中输入命令data=[4 1 2 3 1 3 1 3 2 4];y=hist(data,4)，结果是( B ) (A) y=4 1 2 3 (B) y=3 2 3 2 (C) y=1 3 2 4 (D) y=4 2 1 1 7．在MA TLAB 命令窗口中键入A=magic(6); B=A(2:5,1:2:5)将得到矩阵B ，B 是( C ) (A) 2行5列矩阵 (B)4行两列矩阵 (C) 4行3列矩阵 (D)4行5列矩阵 8．MA TLAB 绘三维曲面需要构建网格数据，语句[x,y]=meshgrid(-2:2)返回数据中( D ) (A) x 是行向量，y 是列向量 (B) x 是列向量，y 是行向量 (C) x 是行元素相同的矩阵 (D) x 是列元素相同矩阵 9．下面有关MATLAB 函数的说法，哪一个是错误的( D ) (A) 函数文件的第一行必须由function 开始，并有返回参数，函数名和输入参数 (B) MA TLAB 的函数可以有多个返回参数和多个输入参数 (C) 如果函数文件内有多个函数，则只有第一个函数可以供外部调用 (D) 在函数中可以用nargin 检测用户调用函数时的输出参数个数 10．将带小数的实数处理为整数称为取整，常用四种取整法则是：向正无穷大方向取整、向负无穷大方向取整、向零方向取整和四舍五入取整。MA TLAB 提供了如下四个取整函数，若a=-1.4，对a 取整的结果是-1，则不应该选用下面哪个函数。( A ) (A) floor (B) round (C) ceil (D) fix 二、程序阅读理解(共24分，每小题3分) 1．如果存在一条曲线L 与曲线簇中每一条曲线相切，则称L 为曲线簇的包络。简单直线簇的实验程序如下： N=input('input N:='); x=[0:N]/N;y=1-x;MATLAB数学实验练习题
数学实验练习题2012
matlab数学实验复习题(有标准答案)
《数学实验》试题答案
MATLAB实验练习题(计算机)-南邮-MATLAB-数学实验大作业答案
清华大学2002至2003学年第二学期数学实验期末考试试题A
数学实验练习题(MATLAB)
电子科大成都学院2011年数学实验复习题
MATLAB数学实验100例题解
matlab数学实验练习题
数学实验模拟试题
matlab数学实验练习题
MATLAB实验练习题(计算机) 南邮 MATLAB 数学实验大作业答案
数学实验(MATLAB)课后习题答案
数学实验matlab练习题
南邮MATLAB数学实验答案(全)
2007年《数学实验》试卷