MATLAB综合练习题Word文档下载推荐.docx
- 文档编号:22548064
- 上传时间:2023-02-04
- 格式:DOCX
- 页数:23
- 大小:115.94KB
MATLAB综合练习题Word文档下载推荐.docx
《MATLAB综合练习题Word文档下载推荐.docx》由会员分享,可在线阅读,更多相关《MATLAB综合练习题Word文档下载推荐.docx(23页珍藏版)》请在冰豆网上搜索。
6、MATLAB中生成平面网格坐标矩阵的函数为
(A)meshgrid(B)gridmesh(C)subplot(D)sub
7、给定
个点可唯一确定一个
次多项式,需调用的命令
(A)poly(x,y,n)(B)polyfit(x,y,n)(C)polyfit(n)(D)poly(n)
8、MATLAB中‘legend’函数的功能为
(A)给图形加标题(B)在指定位置给图形加注释
(C)给图形变色(D)给图形添加图例
9、已知
,则在MATLAB中调用命令A(3)的结果为
(A)3(B)2(C)无结果(D)7
10、MATLAB中调用命令[L,U,P]=lu(A)后,矩阵A满足
(A)
(B)
(C)
(D)
11、MATLAB中调用命令v=eye(2,1)后的结果为
(A)v=[1;
0](B)v=[1,0](C)v=[0;
1](D)v=[1,0;
0,1]
12、MATLAB中计算积分
,则下列程序中的空格处应填的命令是
symt
;
f=3/((t-1)^2*(t-2));
F=_____________,y=double(F)
(A)(f,t,2,6)(B)qad(f,t,2,6)(C)int(t,2,6)(D)int(f,t,2,6)
13、如果在MATLAB的命令窗口输入命令a=input('
a=?
'
),出来的结果为
(A)a=?
(B)没有任何显示(C)a=(D)’a=?
’
14、’try’语句为MATLAB所特有的试探性语句,下面的程序MATLAB能执行吗,执行的命令为A=[1,2,3;
4,5,6];
B=[7,8,9;
10,11,12];
tryC=A*B;
catchC=A.*B;
end
(A)不执行(B)执行的是’catch’后的命令
(C)执行,但没有结果(D)执行的是’try’后的命令
15、MATLAB中与linspace(3,7,11)等价的命令为
(A)3:
7(B)3:
0.4:
7(C)3:
4/11:
7(D)3:
7/10:
7
16、对于下面的程序,在MATLAB中运行后f的形式为
a=0;
b=3*pi;
n=1000;
h=(b-a)/n;
x=a:
h:
b;
f=exp(-0.5*x).*sin(x+pi/6)
(A)数(B)矩阵(C)向量(D)没有结果
17、已知
,MATLAB中若想由x生成矩阵
可调用命令
(A)A=[x,x.^2,x.^3](B)A=[x,x^2,x^3](C)A=[x;
x^2;
x^3](D)A=[x;
x.^2;
x.^3]
18、已知数据点为x=-6:
6,y=sin(x),x0=-6
:
0.25
6,过
、
做分段三次样条插值,并求该多项式再x0处的函数值,在MATLAB中应使用的命令是
(A)interp1(x,y,x0,’pchip’)(B)interp1(x,y,x0,pchip)
(C)interp1(x,y,x0,’spline’)(D)interp1(x,y,x0,’cubic’)
19、下面命令中不能求解线性方程组
的命令是
(A)A/b(B)A\b(C)inv(A)*b(D)A^(-1)*b
20、MATLAB中命令exp
(1)的结果为
(A)0(B)1(C)e(D)inf
21、MATLAB中命令sum([412;
036])的结果为
(A)16(B)448(C)79(D)程序出错,没有结果
22、MATLAB中生成希尔伯特矩阵的函数为
(A)invhilb(n)(B)company(n)(C)eye(n)(D)hilb(n)
23、MATLAB中给变量w赋空矩阵的语句为
(A)w=zeros()(B)w=[](C)w=zeros(0,0)(D)clearw
24、在MATLAB中,如需查询一个函数的基本用法,命令为
(A)help(B)version(C)demo(D)path
25、A=[121;
223;
500],则MATLAB中求A的逆矩阵的命令为
(A)det(A)(B)uv(A)(C)A*(-1)(D)inv(A)
26、b=[5,23,7,2,3,6,0],要把b向量的元素按从小到大的顺序排列,MATLAB中用的命令为
(A)sort(B)put(C)fprintf(D)mesh
27、已知A=[40],MATLAB中C=[A,ones(size(A))]为
(A)4011(B)4010(C)4000(D)4040
28、在MATLAB中,列出当前工作空间变量名的函数是
(A)path(B)clc(C)who(s)(D)draw
29、计算矩阵、向量范数的命令为
(A)length(B)norm(C)size(D)anorm
30、MATLAB对图形窗口进行灵活分割的命令为
(A)plot(B)subplot(C)aplot(D)ploval
31、MATLAB中生成一个2行4列的单位矩阵的命令为
(A)eye(2,4)(B)eye(2,2)(C)eye
(2)(D)zeros
32、已知x=[4,2,],求以x为根的多项式的语句是
(A)x=[4,2,];
polval(x)(B)x=[4,2,];
poly(x)
(C)p=roots(x);
poly(p)(D)x=[4,2,];
roots(x)
33、在MATLAB中,清除所有变量的命令为
(A)clear(B)clc(C)delete(D)holdoff
34、命令floor(3.4)及ceil(3.4)的值为多少
(A)34(B)44(C)33(D)43
35、在MATLAB中,打开示例程序的命令为
(A)demo(B)help(C)clear(D)plot
36、MATLAB中保持图形窗口不关闭的命令为
(A)alter(B)holdon(C)holdoff(D)cd
37、在MATLAB中,强行终止程序运行的命令为
(A)ctrl+v(B)ctrl+shift(C)ctrl+c(D)shift+c
38、
,则A.*B=
(A)运算不执行,因有错(B)683(C)574(D)663
39、双纵坐标函数为
(A)plot(B)yplotty(C)plot3(D)plotty
40、已知数据点为x=[1,1.4,1.6,1.8,2.0]'
,y=[3.75,3.99,2.07,1.52,4.22]'
,拟合一条过这些点的3次曲线,应使用的命令是
(A)polyfit(x,y)(B)polyfit
(2)(C)polyfit(x,y,1)(D)polyfit(x,y,3)
41、MATLAB中表示无穷的符号是
(B)int(C)inf(D)+in
42、MATLAB中给图形加网格线的命令为
(A)gridoff(B)holdon(C)boxon(D)gridon
43、MATLAB中程序A=[123;
-121;
352];
B=[346;
589;
453];
C=(B>
1)&
(B<
7),运行后结果为
(A)C=[1,1,1;
1,0,0;
1,1,1](B)C=[1,0,0;
1,0,1]
(C)无结果(D)C=[0,1,0;
1,1,1]
44、MATLAB中生成一个从0到20的等差为0.5的数列应使用的命令为
(A)x=linspace(0,20)(B)x=0:
20(C)x=linspace(0,20,41)(D)x=linspace(0,20,0.5)
45、积分区间被等分为奇数个小区间时,合理使用1/3辛普森法和3/8辛普森法的是
(A)只用3/8辛普森法
(B)1/3辛普森法和3/8辛普森法都不能使用,用cotes方法
(C)前三个区间用3/8辛普森法,剩下的区间用1/3辛普森法(D)只用1/3辛普森法
46、MATLAB中c=[1234],要得到c的转置向量所用命令为
(A)b=c’(B)b=[1,2,3,4](C)b=inv(c)(D)b==c
47、MATLAB中计算积分
f=5/((t-1)*(t-2)*(t-3));
(A)(f,t,4,5)(B)qad(f,t,4,5)(C)int(t,4,5)(D)int(f,t,4,5)
48、MATLAB表示圆周率的符号是
(B)int(C)pi(D)没有符号可表示
49、MATLAB中a=[1,3,2],则b=a.^2的结果为
(A)b=[1,6,4](B)b=[1,9,4](C)b=[2,6,4](D)执行出错,无结果
50、求解
在x=1附近的根的错误的语句为
(A)fzero('
x^3-3*x+1=0'
1)(B)fzero('
x^3-3*x+1'
1)
(C)fzero('
[1,2])(D)P=[1,-3,1];
roots(p)
三、1.阅读下述MATLAB脚本程序
A=[3579;
0241;
6448];
B=ones(3,3);
b=A([12],2:
3);
B(4,4)=5;
A([1,2],:
)=[];
C=[A;
B(:
1:
3)'
];
D=reshape(C,2,8);
写出b、A、B、C、D的结果。
2.用MATLAB编写函数,计算
,函数的输入参数为x,输出参数为f,函数名为fun。
functionf=fun(x)
x=input('
请输入x的值:
);
ifx<
1&
x>
=0
y=x;
elseifx<
2&
=1
y=2-x;
else
y=0;
y
3.用如下迭代公式求
,迭代的终止条件为
迭代初值
迭代次数不超过100次。
分别对迭代结果和准确
值进行比较,并统计迭代次数。
迭代公式:
试对此迭代问题编写M函数文件。
function[x,n]=sqrt_a(a)
x=1.0;
fork=1:
100
m=x;
x=x/2+a/(2*x);
ifabs(x-m)<
=10^(-5)
break
x;
n=k;
s=(x-sqrt(a));
ifabs(s)<
disp('
正确'
else
错误'
end
3.阅读一绘制地形图的程序,写出所标出的语句的功能。
fid=fopen('
F:
\work\dengzhixian.txt'
'
r'
__打开数据文件_
data=fscanf(fid,'
%f%f%f'
[3,inf]);
_从文件中读取数据_
data=data'
fclose(fid);
____关闭文件____
x=data(:
1);
y=data(:
2);
z=data(:
3);
xb=min(x);
xe=max(x);
yb=min(y);
ye=max(y);
Xl=linspace(xb,xe,200);
___将区间离散化_
Yl=linspace(yb,ye,200);
[X,Y]=meshgrid(Xl,Yl);
_生成二维网格数据__
method='
cubic'
Z=griddata(x,y,z,X,Y,method);
__曲面拟合插值计算_
figure;
___打开一个窗口__
holdon;
___保持打开的图形窗口_
contour(X,Y,Z);
__绘制等高线图形__
colormapcool;
__设定图形色彩_
boxon;
___对坐标系的右边和上边加边框__
xlabel('
x'
_坐标横轴添加标签为x轴__
ylabel('
y'
axisequal;
__坐标轴的长度单位设成相等__
title('
等值线图:
(单位:
scale)'
__添加标题_
holdoff_关闭图形保持__
4阅读下述MATLAB脚本程序
)=[]
答案:
b=[57;
24]
A=[6448]
B=1110
1110
0005
C=6448
1110
D=61414180
11111100
5.下面的函数完成什么功能?
functionf=factor(n)
ifn<
f=1;
f=factor(n-1)*n;
6.下面的程序完成功能是什么?
t=0:
pi/50:
4*pi;
y0=exp(-t/3);
y=exp(-t/3).*sin(3*t);
plot(t,y,'
-r'
t,y0,'
b'
)
t'
gridon;
答:
绘制两条曲线,第1条是蓝色的虚线;
第2条是红色的实线。
给坐标轴添加标签。
7.下面的函数主要完成什么功能?
a=0;
i=1;
while(a<
100000)
a=a+pow2(i);
i=i+1;
m=i-1
求满足下式的最小m值
8.编写一个函数,计算下面函数的值,给出x的值,调用该函数后,
返回y的值。
functiony=fun(x)
y=sin(x);
=3&
y=-x+6;
9.用如下迭代公式求
迭代次数不超过200次。
function[x,n]=sqrt_a(a)
x=1.5;
200
=10^(-8)
ifabs(s)<
10以下为对一水准网进行间接平差的程序,试阅读并填写所标出的
语句的功能。
functionshuizhunjjpc()__函数定义行_
水准网间接平差示例'
已知高程'
)___提示输入已知高程_
Ha=5.015
Hb=6.016
观测高差'
L=[1.359;
2.009;
0.363;
1.012;
0.657;
-0.357]
fid=fopen(‘E:
\szjj.txt’,‘r’)__打开数据文件_
temp=fscanf(fid,‘%f%f%f%f’,[4,6])__从文件中读取数据_
__关闭数据文件_
temp=temp‘;
B=temp(:
[1,2]);
_读取矩阵的第1、2列_
l=temp(:
S=temp(:
4);
__读取矩阵的第4列__
S=S’;
disp(‘C是单位线路公里数,S是线路长度'
C=ones(1,6)__生成1行6列全1阵_
P=C./S;
_计算2个向量的点除_
P=diag(P)__计算观测值的权阵_
disp(‘参数的解’)
x=inv(B'
*P*B)*B'
*P*l__解算法方程_
误差V(mm),各待定点的高程平差值L1(m)'
V=B*x-l_计算观测值的改正数向量_
LP=L+V/1000_计算观测值的平差值__
精度评定'
n=6;
_输入观测值个数_
t=2;
delta=sqrt(V'
*P*V/(n–t))__计算验后单位权中误差_
四、已知多项式
和
,编写脚本程序计算:
1)多项式a(x)+b(x)
2)多项式a(x)*b(x)。
3)多项式a(x)的根。
五、用MATLAB编写函数,计算
六、编程求
的值,要求:
(1)使用循环语句编程,但不能使用数组运算符和sum;
(2)使用数组运算符和sum编程,不能使用循环结构。
7、已知超定方程组
写出用Matlab中语言求解上述线性方程组的程序。
八、Hilbert矩阵是非常著名的病态矩阵,写出其定义,并编出Matlab程序求
到
的Hilbert矩阵的条件数、行列式。
九、数值求解定积分时,当积分区间被等分为奇数个小区间时,如何合理使用1/3辛普森法和3/8辛普森法。
十、写出将矩阵
进行LU分解,并求出该矩阵的行列式的语句
x=-5:
5;
y=1./(1+x.^2);
t=-5:
0.05:
y1=1./(1+t.^2);
n=length(t);
fori=1:
n
z=t(i);
s=0;
for
Lk=1;
u=x(k);
for
if
Lk=Lk*(z-x(j))/(u-x(j));
;
y2(i)=s;
plot(x,y,'
ko'
t,y1,t,y2,'
十一、写出
个插值节点
的Lagrange插值多项式的公式,并将下面的Lagrange插值多项式的Matlab实现程序补充完整。
(每空3分,共15分)
Lagrange插值多项式
公式:
十二、已知
的表达式
,当
时,求
的值,请补充完整。
n=100;
i=1:
n;
f=;
y=
十三、
是数据点的横坐标数组,
是拟合数据点
的多项式,
是一个矩阵,它的第
行即为
的幂系数,请补充完整。
functionp=shape_pw(x)
np=length(x);
forj=1:
np
y=zeros(1,np);
y(j)=1;
p(j,:
)=;
十四、下面的程序是解线性方程组的一个函数文件,请补充完整。
function[x,y]=line_solution(A,b)
[m,n]=size(A);
y=[];
ifnorm(b)>
0
ifrank(A)==rank([A,b])
ifrank(A)==
x=A\b;
else
disp('
原方程组有有无穷个解,其齐次方程组的基础解系为y,特解为x'
y=null(A,'
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB 综合 练习题