MATLAB考试试题题库

MATLAB考试试题(1)

产生一个1x10的随机矩阵，大小位于(-55),并且依据从大到小的挨次排列好！(注：要程序和运行

结果的截屏〕

答案：

e=10*rand(lz10)-5;

b=sort(a/descend

”)

1,请产生一个100*5的矩阵，矩阵的每一行都是口2345]

2.变量：A=,ilovematlab,；B=，matlab',请找出：

(A)B在A中的位置。

(B)把B放在A后面，形成C='ilovematlabmatlabr

3,请修改下面的程序，让他们没有for循环语句！

A=[l23;456;789];

[rc]=size(A);

fori=l:l:r

forj=l:l:c

if(A(i,j)>8|A(iJ)<2)

A(ij)=0;

end

end

end

4,请把变量A=[l23;456;789]写到文件里(output.xls),写完后文件看起来是这样的

123456789

5,试从Yahoo网站上获得微软公司股票的2025年9月的每日收盘价。

6.编写M文件，从Yahoo网站批量读取60000.SH至600005.SH在2025年9月份的每日收盘价

(提示：使用字符串函数)。

7.将金牛股份(000937)2025年12月14日至2025年1月10日的交易记录保存到Excel中，

编写程序将数据读入MATLAB中，土一步将数据读入Access数据库文件。

8.资产每日回报率为0.0025,标准差为0.0208,资产现在价值为0.3亿，求5%水平下资产的

10天在险价值(Var)。

9.a=[l2345],b=a(l)*a(5)+a(2)*a(4)+a(3)*a(3)+a(4)*a(2)+a(5)*a(l).试用MATLAB

中最简洁的方法计算b.留意最简洁哦。

1、求以下联立方程的解

3x+4y-7z-12w=4

5x-7y+4z+2w=-3

x+8z-5w=9

-6x+5y-2z+10w=-8

求系数矩阵的秩；

求出方程组的解。

解：⑴

»a=l34-7-12];

5-742；

108-5;

-65-210];

c=[4;-3;9;-8];

b=rank(a)

b=4

(2)»d=a\c

d=-1.4841,-0.6816,0.5337,-1.2429

即：x=-1.4841;y=-0.6816;z=0.5337;w=-1.2429

2、设y=cos[0.5+((3sinx)/(l+x，2))]把x=0~2冗间分为101点，画出以x为横

坐标，y为纵坐标的曲线；

解：»x=linspace(0,2*pi,101);

y=cos(0.5+3.*sin(x)./(l+x.*x));

plot(x,y)

3、设f(x)=x-5-4x“4+3x〃2-

2x+6

<1)取x-[-2,8]之间函数的值(取100个点)，画出曲线，看它有几个零点。

(提示：用polyval函数)

解：»p=[l-43-26];

x=linspace(-2,8,100);

y=polyval(p,x);

plot(x,y);

axis([-2,8,-200,2300]);

为了便于观看，在y=0处画直线，图如下所示:

与y=0直线交点有两个，有两个实根。

(2)用roots函数求此多项式的根

»a=roots(p)

a=3.0000,1.6956,-0.3478+1.0289i,-0.3478-

1.02891

4、在[-10,10；TO,10]范围内画出函数的三维图形。

解：»[X,Y]=meshgrid(-10:0.5:10);

&=sqrt(X.2+Y.2)+eps;

Z=sin(ci).Ai;

mesh(X,Y,Z);

natlab试卷，求答案

一、选择或填空(每空2分，共20分)

1、标点符号()可以使命令行不显示运算结果，()用来表示该行为注释行。

2、以下变量名中()是合法的。

(A)char_l;(B)x*y;(C)x\y;(D)end

，为〜，步长为的向量，使用命令()创立。

A输入矩阵，使用全下标方式用（）取出元素“”，使用单下标方式用（）取出元素

5、符号表达式中独立的符号变量为（）<,

6、M脚本文件和M函数文件的主要区分是（）和（

）。

7、在循环构造中跳出循环，但连续下次循环的命令为（）o

（A）return;（B）break;（C）continue;（D）keyboad

二、（此题12分）利用MATLAB数值运算，求解线性方程组（将程序保存为test02.m文件）

三、（此题20分）利用MATALAB符号运算完成（将程序保存为test03.m文件）：

（1）创立符号函数

（2）求该符号函数对的微分；

（3）对趋向于求该符号函数的极限：

（4）求该符号函数在区间上对的定积分：

（5）求符号方程的解。

网、（此题20分）编写MATALAB程序，完成以下任务（将程序保存为test04.m文件）：

（1）在区间上均匀地取20个点构成向量：

（2）分别计算函数与在向量处的函数值：

（3）在同一图形窗口绘制曲线与，要求曲线为黑色点画线，曲线为红色虚线圆圈：并在图中恰

当位置标注两条曲线的图例；给图形加上标题“y1andy2-0

五、（此题15分）编写M函数文件，利用for循环或while循环完成计算函数的任务，并利用该

函数计和时的和（将总程序保存为test05.m文件）。

六、［此题13分）求解线性规划模型：

的MATLAB命令为x=linp「og

（c,A,b,Aeq,beq,VLB,VUB）

试编写MATLAB程序，求解如下线性规划问题（将程序保存为怕stD6.m文件）:

问题补充：

卷子的地址

看不见符号，能做就做了一些.

1、标点符号（；）可以使命令行不显示运算结果，（%）用来表示该行为注释行。

2、以下变量名中（A）是合法的。

（A）char_l;（B）x*y;（C）x\v;（D）end

3、为〜，步长为的向量，使用命令（此题题意不清）创立。

A输入矩阵，使用全下标方式用（此题题意不清）取出元素“”，使用单下标方式用（此

题题意不清）取出元素”工

5、符号表达式中独立的符号变量为（）o

6、M脚本文件和M函数文件的主要区分是（变量生存期和可见性）和

（函数返回值）0

7、在循环构造中跳出循环，但连续下次循环的命令为（C）<

（A）return;（B）break;（C）continue:（D）keyboad

二、（此题12分）利用MATLAB数值运算，求解线性方程组（将程序保存为test02.m文件）

三、（此题20分）利用MATALAB符号运算完成（将程序保存为【est03.m文件）：

（1）创立符号函数symsx

（2）求该符号函数对的微分；

（3）对趋向于求该符号函数的极限；

（4）求该符号函数在区间上对的定积分；

（5）求符号方程的解。

四、（此题20分）编写MATALAB程序，完成以下任务（将程序，呆存为test04.m文件）：

（1）在区间上均匀地取20个点构成向量：

（2）分别计算函数与在向量处的函数值；

（3）在同一图形窗口绘制曲线与，要求曲线为黑色点画线，曲线为红色虚线圆圈：并在图中恰

当位置标注两条曲线的图例：给图形加上标题“ylandy2w

五、（此题15分）编写M函数文件，利用for循环或while循环完成计算函数的任务，并利用该

函数计算时的和（将总程序保存为test05.m文件）。

六、（此题13分）求解线性规划模型：

F勺MATLAB命令为

x=linprog（c,A,b,Aeq,beq,VLB,VUB）

试编写MATLAB程序，求解如下线性规划问题（将程序保存为test06.m文件）：

[例2.1]SIS0系统的状态空间表达式为(2-3)式，求系统的传递函数。A=[0

10;001;-4-3-2];B=[1;3;-6];C=[100];D=0;

[num,den]=ss2tf(a,b,c,d,u)

[num,den]=ss2tf(A,B,C,D,1)

[例2.2]从系统的传递函数(2-4)式求状态空间表达式。

num=[153];

den=[1234];

[A,B,C,D]=tf2ss(num,den)

[例2.3]对上述结果进展验证编程。

中将[例2.2]上述结果赋值给A、B、C、D阵;

R=[-2-3-4;100;010]；B=[1;0;0]：C=[153]：D=0：

[num,den]=ss2tf(A,B,C,D,1)

[例2.4]给定系统G(s)=S3+2S2+S+3,求系统的零极点增益模型和状态空间模型，并求其

S3+0.5S2+2S+1

单位脉冲响应及单位阶跃响应。

解：

rum=(1213];den=(10.521);

sys=tf(num,den)%系统的传递函数模型

Transferfunction:

sA3+2sA2+s+3

sA3+0.5sA2+2s+1

sysl=tf2zp(num,den)%系统的零极点增益模型sysl=

sys2=tf2ss(sys)%系统的状态空间模型模型；或用［a,b,c,d］=tf2ss(num,den)形式

impulse(sys2)%系统的单位脉冲响应

step(sys2)%系统的单位阶跃响应

［例3.1］对下面系统进展可控性、可观性分析。

解：

a=[-l-22;0-11;10-l];b=[20l]w;c=[l20]

Qc=ctrb(a,b)%生成能控性判别矩阵

rank(Qc)%求矩阵Qc的秩

ans=3先满秩，故系统能控

Qo=obsv(a,c)%生成能观测性判别矩阵

rank(Qo)舟求矩阵Qo的秩

ans=3舟满秩，故系统能观测

[例3.2]系统状态空间方程描述如下：

试判定其稳定性，并绘制出时间响应曲线来验证上述推断。

解：

A=[-10-35-50-24:1000;0100;0010]：

B=[l;0;0;0];C=[l72424];D=[0];

[z,p,k]=ss2zp(A,B,C,D,1);

Flagz=0;

n=length(A);

fori=l:n

ifreal(p(i))>0

Flagz=l;

end

end

disp("系统的零极点模型

为")；z,p,k系统的零极点模型为

ifFlagz==l

dispC系统不稳定”)；

elsedisp("系统是稳定

的“);end

运行结果为：

系统是稳定的

step(A,B,C,D)%系统的阶跃响应

资源与环境工程学院2025级硕士争论生《MatLab及其应用》试题

留意，每题的格式均须包含3个局部

a.程序(含程序名及完整程序)：

b.运行过程：

c.运行结果：

(1)求解线性规划问题：

minZ=-Ax[+X2+7X3

stx+x—x=5

l23

3x-x+x<4

I23一

X|+x2-4犬34-7

x,x>0

I2一

问各外分别取何值时，z有何微小值。(10分)

答：fprintf("线性规划问题求解\n");

f=[-4;l;7];

A=[3/-l/l;l,l,-4;];

b=[4,・7]”；

Aeq=[14,-1];

beq=[5]”；

lb=[0,0,];

ub=[];

x=linprog(f,A,b,Aeq,beq,lb,ub);

x

z=f"*x;

fprintf("MINz=%f\n”,z);

运行结果：线性规划问题求解

Optimizationterminatedsuccessfully,

x=

2.2500

6.7500

4.0000

MINz=25.750000

0.5x,x<2

⑵编写一个函数，使其能够产生如下的分段函数：/(X)=J15-0.25X,2<X<6,

0.5>6<x

并调用此函数，绘制在x=[0,+2]范围的/6)・/(»2)曲线。(10分)

答：functiony=f(x)

ifx<=2

y=0.5*x;

elseifx>6

V=0.5;

elsey=1.5-0.25*x;

end

end

运行结果x=2

f(x)=l

x=0:0.05:2;

y=diag(A2(x)"*A2(x+2));

P>ot(x,y);

xlabel(M\bfx°);

ylabel(M\bfy");

⑶将一个屏幕分4幅，选择适宜的步长在右上幅与左下幅绘制出以下函数的图形。

（10分）

①历而J；（曲线图）；②=+Zl；(-2<x<2,-4<y<4)(曲面

42

图)。

答：»subplot(2,2,2)

»ezplot(w(cos(x))"(1/2)w,[-pi/2pi/2])

»ylabel(My")

»subplot(2,2,3)

»x=-2:0.5:2;

»y=-4:1:4;

»czsurfc(Mx〃2/2-2+y"2/4"2")

(4)A是一个雉度mxn的矩阵.写一段程序，算出A中有多少个零元素(10分)

答：»A=input("请输入一个矩阵")

[m,n]=size(A);

sig=O;

fori=l:m

forj=l:n

ifA(i,j)==O

sig-sig+1;

end

end

end

请输入一个矩阵[012;102;000]

A=

012

102

000

»sig

sig=

5

⑸向量A%鸟，写一段程序，找出人中的最小元素〔io分)

答：A=input(w请输入一个向量”)

[m,n]=sizeA

min=A(l,n);

fori=l:n

ifA(lJ)<min

min=A(l,i)

end

end

请输入一个向量[123-520]

A=

123-520

min=

-5

B.应用题(50分)

依据专业方向特色和相关科研工作需求，经过与导师商量后，结合一个课题具

体任务，编写一份Matlab应用工作报告。报告由：a课题任务要求，b技术路

线，c程序，d运行结果，e总结、等局部构成，完成的报告经导师给出简洁评

语并签字后缴来。

a,课题任务：争论了一种生物质，油菜秸秆对水溶液中金属离子铜的吸附

行为，分别从pH,用量，温度几个方面考察秸秆的吸附性，并对分析的最正确

条件进展了探讨。同时从吸附热力学和吸附动力学角度探讨吸附机理。结果表

明，100ml溶液pH=5.30,秸秆用量0.75g时，桔秆对铜的吸附量可到达6mg/g

左右。

b,技术路线：通过试验，获得一系列的数据，然后通过Matlab来做各种关

系图。从图中找到g各种关系式。

C,程序:x=[2.202.723.444.135.38]

y=[2.393.836.076.395.84]

plot(x,y);

xlabel("pH“);

ylabel(“吸附量”)

图1

x=[0.50.751.01.251.5]

y=[6.056.195.334.694.02]

plot（x,y）;

xlabel（"秸秆用量g”）；

ylabel（“吸附量”）

S5

特秆用包］

图2

通过数据图，得到比较抱负的试验条件pH和秸秆用量，接下来做动力

学和等温线。

>x=[0.1670.5123458]

y=[0.0332010.0860590.16977903220610.48C7690.6441220.809061

1.269841]

plot(x,y);

xlabel(“时间t");

ylabelC时间/吸附

图3

x=[0.23630.154960.136190.129060.133730.13315]

y=[0.252180.047070.020250.012670.008810.00706]

Plot(x,y);

图4

x=[0.626540.809770.865850.88920.873770.87564]

y=[0.598291.32731.695891.897372.055032.15149]

plot(x,y);

xlabel(^^Lg吸附

量”);ylabel(°Lg平

d,总结：从图1和图2,分析看可以得到比较抱负的对于本次试验的pH和

秸秆用量。后面试验是在前面的根底上得到的。图3是吸附动力学反响速率

图，从图中可以看到线性拟合程度很好，符合二级反响速率方程。图4和图5

是吸附等温线作图，看以看出图4的线性拟合较图5的好，说明符合Langmuir

吸附等温模型。

[例2.1]SIS0系统的状态空间表达式为(2-3)式，求系统的传递函数。A=[0

10;001;-4-3-2];B=[1;3;-6];C=[100];D=0;

[num,den]=ss2tf(a,b,c,d,u)

[num,den]=ss2tf(A,B,C,D,1)

[例2.2]从系统的传递函数(2-4)式求状态空间表达式.

num=[153];

den=[1234]：

LA,B,C,D]=tf2ss(num,den)

[例2.3]对上述结果进展验证编程。

%将[例2.2]上述结果赋值给人、8、C、D阵；

A=[-2-3-4;100：010]；B=[1;0;0]：C=[153]：D=0；

[num,den]=ss2tf(A,B,C,D.1)

[例2.4]给定系统G(s)=S3+2S2+S+3,求系统的零极点增益模型和状态空间模

S3+0.5S2+2s+1

型，并求其单位脉冲响应及单位阶跃响应。

解：

num=[l213];den=[l0.521];

sys=tf(num,den)%系统的传递函数模型

Transferfunction:

sA3+2sA2+s+3

sA3+0.5sA2+2s+1

sysl=tf2zp(num,den)%系统的零极点增益模型sysl=

sys2=tf2ss(sys)%系统的状态空间模型模型：或用归,6(:,5=1£255(加111,(^11)形式

impulse(sys2)%系统的单位脉冲响应

step(sys2)%系统的单位阶跃响应

［例3.1］对下面系统进展可控性、可观性分析。

解：

a=[-l-22;0-11:10-l];b=[201]”；c=[l20]

Qc=ctrb(a,b)%生成能控性判别矩阵

rank(Qc)先求矩阵Qc的秩

ans=3先满秩，故系统能控

Qo=obsv(a,c)%生成能观测性判别矩阵

rank(Qo)务求矩阵Q。的秩

ans=3时的秩，故系统能观测

［例3.2］系统状态空间方程描述如下:

试判定其稳定性，并绘制出时诃响应曲线来验证上述推断。

解：

A=[-10-35-50-24;1000;0100;0010];

B=[l;0;0;0];C=[l72424];D=[0];

[z,p,k]=ss2zp(A,B,C,D,1)：

Elagz=0;

n=length(A);

fori=l:n

ifreal(p(i))>0Flagz=l;

end

end

disp("系统的零极点模型

为")；z,p,k系统的本极点模型为

ifFlagz==l

disp("系统不稳定”)；

elsedisp（M系统是稳定的”）;

end

运行结果为：

系统是稳定的

step(A,B,C,D)%系统的阶跃响应。

1、使用以下哪一个函数可以产生单位矩阵(B)

A.zerosB.eyeC.randD.diag

2、以下哪一个函数是求模函数(I))

A.remB.signC.fixD.mod

3、使用以下哪一个函数可以交换矩阵左右对称位置上的元素(A)

A.fliplrB.flipdimC.flipudI),find

4、使用以下哪一个函数可以比较字符串，且比较时无视字符的大小写(D)

A.strncmpB.strcmpC.strncmpiD.strcmpi

5、要利用图形方式显示元胞数组，则应当使用以下哪一个函数(B

A.cellfunB.celIplotC.celIdispD.cel12mat

6、以下哪一个函数可以猎取构造字段的数据(B)

A.fieldnamesB.getfieldC.setfieldD.rmfield

7、执行以下哪一条命令后，图形窗体的轴将显示坐标网格线(A)

A.gridonB.holdonC.gridoffD.holdoff

8、进展格式化绘图时，使用哪一个函数可以添加图例(B)

A.titleB.legendC.labelD.text

9、使用以下哪一条指令可以将图形窗体分割成二行三列，并且将第一行其次列

的绘图区域设置为当前的绘图区域(B)

A.subplot(2,3,1)B.subplot(2,3,2)

C.subplot(2,3,4)D.subplot(2,3,5)

10、使用以下哪一个函数可以绘制三维网线图(C)

A.surfB.plotC.meshD.plot3

1、»A=[l23:456];B=[25;83](2分)

B=

25

83

2、假设向量八=[9876543210]

»A([l:3:5])(2分)

ans=

96

»A([l35])(2分)

ans=

975

3、»A=ones(2,2);A(:)=1:4;

»A*A'(2分)

ans=

1014

1420

»B=A.*A*(2分)

B=

16

616

4、使用三元组法，将以下满屿转变为稀疏矩阵

1500220-15

0113000

S=000-600

9100000

0002800

解：

»ir=[l4221351]：

»jc=[l1234446];

»data=[159111322-628-15];

>>s=sparse(ir,jc,data,5,6)

(1,1)15

(4,1)91

(2,2)11

(2,3)3

(1,4)22

(3,4)-6

(5,4)28

(1,6)-15

5、»A=reshape(l:24,4,6);

»A(:,[2•54])=[]

A=

11721

21822

31923

42024

6、使用函数struct创立一个构造。

此构造名为Student；

有三个字段，分别为name、age、grade；

有两条记录，分别为‘Way'、23、3和'Deni'、21、1

解：

>>_________________________________________________________________

Student=struct(Mnamew,[nWayDeniage”,{23,21}grade”,

{3,1})7、绘出下幅图

»plot([l4285])

8

7

6

5