MATLAB-考试试题-(1)汇总
MATLAB 考试试题 (1)
产生一个1x10的随机矩阵，大小位于（-5 5），并且按照从大到小的顺序排列好！（注：要程序和运行结果的截屏）
答案：
a=10*rand(1,10)-5;
b=sort(a,'descend')
1.请产生一个100*5的矩阵，矩阵的每一行都是[1 2 3 4 5]
2. 已知变量：A=’ilovematlab’；B=’matlab’, 请找出：
（A） B在A中的位置。
（B）把B放在A后面，形成C=‘ilovematlabmatlab’
3. 请修改下面的程序，让他们没有for循环语句！
A=[1 2 3; 4 5 6; 7 8 9];
[r c]=size(A);
for i=1:1:r
for j=1:1:c
if (A(i,j)>8 | A(i,j)<2)
A(i,j)=0;
end
end
end
4. 请把变量A=[1 2 3; 4 5 6; 7 8 9]写到文件里(output.xls)，写完后文件看起来是这样的
1 2 3 4 5 6 7 8 9
5.试从Yahoo网站上获得微软公司股票的2008年9月的每日收盘价。
6.编写M文件，从Yahoo网站批量读取60000.SH至600005.SH在2008年9月份的每日收盘价（提示：使用字符串函数）。
7. 将金牛股份（000937）2005年12月14日至2006年1月10日的交易记录保存到Excel中，编写程序将数据读入MATLAB中，进一步将数据读入Access数据库文件。
8.已知资产每日回报率为0.0025，标准差为0.0208，资产现在价值为0.8亿，求5％水平下资产的10天在险价值(Var)。
9.a=[1 2 3 4 5],b=a(1)*a(5)+a(2)*a(4)+a(3)*a(3)+a(4)*a(2)+a(5)*a(1).试用MATLAB 中最简单的方法计算b，注意最简单哦。
1、求下列联立方程的解
3x+4y-7z-12w=4
5x-7y+4z+ 2w=-3
x +8z- 5w=9
-6x+5y-2z+10w=-8
求系数矩阵的秩；
求出方程组的解。
解：（1）
>> a=[3 4 -7 -12];
5 -7 4 2 ;
1 0 8 -5;
-6 5 -2 10];
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)/(1+x^2))]把x=0~2π间分为101点，画出以x为横坐标，y为纵坐标的曲线；
解：>> x=linspace(0,2*pi,101);
y=cos(0.5+3.*sin(x)./(1+x.*x));
plot(x,y)
3、设f(x)=x^5-4x^4+3x^2-
2x+6
（1）取x=[-2,8]之间函数的值（取100个点），画出曲线，看它有几个零点。
（提示：用polyval 函数）
解：>>p=[1 -4 3 -2 6];
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.0289i
4、在[-10，10；-10，10]范围内画出函数的三维图形。
解：>>[X,Y]=meshgrid(-10 : 0.5 :10);
a=sqrt(X.^2+Y.^2) +eps;
Z=sin(a)./a;
mesh(X,Y,Z);
matlab试卷，求答案
一、选择或填空（每空2分，共20分）
1、标点符号（）可以使命令行不显示运算结果，（）用来表示该行为注释行。
2、下列变量名中（）是合法的。
(A) char_1 ; (B) x*y ; (C) x\y ; (D) end
3、为～，步长为的向量，使用命令（）创建。
4、输入矩阵，使用全下标方式用（）取出元素“ ”，使用单下标方式用（）取出元素“ ”。
5、符号表达式中独立的符号变量为（）。
6、M脚本文件和M函数文件的主要区别是（）和（
）。
7、在循环结构中跳出循环，但继续下次循环的命令为（）。
(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）在同一图形窗口绘制曲线与，要求曲线为黑色点画线，曲线为红色虚线圆圈；并在图中恰当位置标注两条曲线的图例；给图形加上标题“y1 and y2”。
五、（本题15分）编写M函数文件，利用for循环或while循环完成计算函数的任务，并利用该函数计算时的和（将总程序保存为test05.m文件）。
六、（本题13分）已知求解线性规划模型：
的MATLAB命令为
x=linprog（c,A,b,Aeq,beq,VLB,VUB）
试编写MATLAB程序，求解如下线性规划问题（将程序保存为test06.m文件）：
问题补充：
卷子的地址
看不见符号,能做就做了一些.
1、标点符号（ ; ）可以使命令行不显示运算结果，（ % ）用来表示该行为注释行。
2、下列变量名中（ A ）是合法的。
(A) char_1 ; (B) x*y ; (C) x\y ; (D) end
3、为～，步长为的向量，使用命令（本题题意不清）创建。
4、输入矩阵，使用全下标方式用（本题题意不清）取出元素“ ”，使用单下标方式用（本题题意不清）取出元素“ ”。
5、符号表达式中独立的符号变量为（）。
6、M脚本文件和M函数文件的主要区别是（变量生存期和可见性）和
（函数返回值）。
7、在循环结构中跳出循环，但继续下次循环的命令为（ C ）。
(A) return; (B) break ; (C) continue ; (D) keyboad
二、（本题12分）利用MATLAB 数值运算，求解线性方程组(将程序保存为test02.m 文件)
三、（本题20分）利用MATALAB 符号运算完成（将程序保存为test03.m 文件）： （1）创建符号函数 syms x （2）求该符号函数对 的微分； （3）对 趋向于 求该符号函数的极限； （4）求该符号函数在区间 上对 的定积分； （5）求符号方程 的解。
四、（本题20分）编写MATALAB 程序，完成下列任务（将程序保存为test04.m 文件）： （1）在区间 上均匀地取20个点构成向量 ； （2）分别计算函数 与 在向量 处的函数值；
（3）在同一图形窗口绘制曲线 与 ，要求 曲线为黑色点画线， 曲线为红色虚线圆圈；并在图中恰当位置标注两条曲线的图例；给图形加上标题“y1 and y2”。
五、（本题15分）编写M 函数文件，利用for 循环或while 循环完成计算函数 的任务，并利用该函数计算 时的和（将总程序保存为test05.m 文件）。 六、（本题13分）已知求解线性规划模型：
的MATLAB 命令为
x=linprog （c,A,b,Aeq,beq,VLB,VUB ）
试编写MATLAB 程序，求解如下线性规划问题（将程序保存为test06.m 文件）：
[例2.1] 已知SISO 系统的状态空间表达式为(2-3)式，求系统的传递函数。 A=[0 1 0;0 0 1;-4 -3 -2];B=[1;3;-6];C=[1 0 0];D=0; [num,den]=ss2tf(a,b,c,d,u) [num,den]=ss2tf(A,B,C,D,1)
[例2.2] 从系统的传递函数(2-4)式求状态空间表达式。 num =[1 5 3]; den =[1 2 3 4];
[A,B,C,D]=tf2ss(num,den)
[例2.3] 对上述结果进行验证编程。 %将[例2.2]上述结果赋值给A 、B 、C 、D 阵；
A =[-2 -3 -4;1 0 0; 0 1 0]；
B =[1;0;0]；
C =[1 5 3]；D=0； [num,den]=ss2tf(A ，B ，C ，D,1)
[例2.4] 给定系统1
25.03
2)(2323++++++=s s s s s s s G ，求系统的零极点增益模型和状态空间模型，并求其
单位脉冲响应及单位阶跃响应。 解：
num=[1 2 1 3];den=[1 0.5 2 1];
sys=tf(num,den) %系统的传递函数模型
Transfer function:
s^3 + 2 s^2 + s + 3
-----------------------------
s^3 + 0.5 s^2 + 2 s + 1
sys1=tf2zp(num,den) %系统的零极点增益模型 sys1 =
sys2=tf2ss(sys) %系统的状态空间模型模型；或用[a,b,c,d]=tf2ss(num,den)形式impulse(sys2) %系统的单位脉冲响应
step(sys2) %系统的单位阶跃响应
[例3.1] 对下面系统进行可控性、可观性分析。
解：
a=[-1 -2 2;0 -1 1;1 0 -1];b=[2 0 1]';c=[1 2 0]
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;1 0 0 0;0 1 0 0;0 0 1 0];
B=[1;0;0;0];C=[1 7 24 24];D=[0];
[z,p,k]=ss2zp(A,B,C,D,1);
Flagz=0;
n=length(A);
for i=1:n
if real(p(i))>0
Flagz=1;
end
end
disp('系统的零极点模型为');z,p,k
系统的零极点模型为
if Flagz==1
disp('系统不稳定');
else disp('系统是稳定的');
end
运行结果为:
系统是稳定的
step(A,B,C,D) ％系统的阶跃响应
资源与环境工程学院2008级硕士研究生《MatLab及其应用》试题注意，每题的格式均须包含3个部分
a. 程序（含程序名及完整程序）：
b. 运行过程：
c. 运行结果：
(1)求解线性规划问题：
7 4
4
3
5
7 4
2
1
3 2
1
3
2
1
3
2
1
3
2
1
≥
-
≤
-
+
≤
+
-
=
-
+
+
+
-
=
x,
x
x x
x
x
x
x
x
x
x .t.s
x
x
x
Z
min
问各x i分别取何值时，Z有何极小值。（10分）答：fprintf('线性规划问题求解 \n');
f = [-4;1;7];
A = [3,-1,1;1,1,-4;];
b = [4,-7]';
Aeq = [1,1,-1];
beq = [5]';
lb = [0,0,];
ub = [];
x = linprog(f,A,b,Aeq,beq,lb,ub);
x
z = f' * x;
fprintf('MIN z = %f \n' , z);
运行结果：线性规划问题求解
Optimization terminated successfully.
x =
2.2500
6.7500
4.0000
MIN z = 25.750000
(2)编写一个函数，使其能够产生如下的分段函数：??
?
??≤≤<-≤=x x x x x x f 65.0620.251.525.0)(，，，
，
并调用此函数，绘制曲线范围的，在2)()(2][0+?+=x f x f x 。（10分）
答：function y=f(x) if x<=2 y=0.5*x; else if x>6 y=0.5;
else y =1.5-0.25*x; end end
运行结果 x=2
f(x)=1
x = 0:0.05:2;
y = diag(A2(x)'*A2(x+2)); plot(x,y); xlabel('\bfx'); ylabel('\bfy');
(3) 将一个屏幕分4幅，选择合适的步长在右上幅与左下幅绘制出下列函数的图形。（10分）
①]2
2[)cos(π
π，，-∈x x （曲线图）； ②4)y 2,-4x (-242),(22
22≤≤≤≤+=；
y x y x f （曲面图）。
答：>> subplot(2,2,2)
>> ezplot('(cos(x))^(1/2)',[-pi/2 pi/2]) >> ylabel('y') >>subplot(2,2,3) >> x=-2:0.5:2; >> y=-4:1:4;
>>ezsurfc('x^2/2^2+y^2/4^2')
(4) A 是一个維度m×n 的矩阵. 写一段程序, 算出A 中有多少个零元素（10分） 答：>> A= input ('请输入一个矩阵') [m,n]= size(A); sig=0; for i=1:m for j=1:n if A(i,j)==0 sig = sig+1; end end end
请输入一个矩阵[0 1 2;1 0 2; 0 0 0] A =
0 1 2 1 0 2 0 0 0 >> sig sig = 5
(5) 向量]a ,,a ,[a A n 11 . 写一段程序, 找出A 中的最小元素（10分） 答：A= input ('请输入一个向量') [m,n]=sizeA min =A(1,n); for i=1:n if A(1,i) 请输入一个向量[1 2 3 -5 2 0] A = 1 2 3 -5 2 0 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.20 2.72 3.44 4.13 5.38] y=[2.39 3.83 6.07 6.39 6.84] plot(x,y); xlabel('pH'); ylabel('吸附量') 图1 x=[0.5 0.75 1.0 1.25 1.5 ] y=[6.05 6.19 5.33 4.69 4.02] plot(x,y); xlabel('秸秆用量g'); ylabel('吸附量') 图2 通过数据图，得到比较理想的实验条件pH和秸秆用量，接下来做动力学和等温线。 > x=[0.167 0.5 1 2 3 4 5 8 ] y=[0.033201 0.086059 0.169779 0.322061 0.480769 0.644122 0.809061 1.269841] plot(x,y); xlabel('时间t'); ylabel('时间/吸附量') 图3 x=[0.2363 0.15496 0.13619 0.12906 0.13373 0.13315] y=[0.25218 0.04707 0.02014 0.01267 0.00881 0.00706] plot(x,y); xlabel('1/吸附量'); ylabel('1/平衡浓度') 图4 x=[0.62654 0.80977 0.86585 0.8892 0.87377 0.87564] y=[0.59829 1.3273 1.69589 1.89737 2.05503 2.15149] plot(x,y); xlabel('Lg 吸附量'); ylabel('Lg 平衡浓度') 图5 d,总结：从图1和图2，分析看可以得到比较理想的对于本次实验的pH 和秸秆用量。后面实验是在前面的基础上得到的。图3是吸附动力学反应速率图，从图中可以看到线性拟合程度很好，符合二级反应速率方程。图4和图5是吸附等温线作图，看以看出图4的线性拟合较图5的好，说明符合Langmuir 吸附等温模型。 [例2.1] 已知SISO 系统的状态空间表达式为(2-3)式，求系统的传递函数。 A=[0 1 0;0 0 1;-4 -3 -2];B=[1;3;-6];C=[1 0 0];D=0; [num,den]=ss2tf(a,b,c,d,u) [num,den]=ss2tf(A,B,C,D,1) [例2.2] 从系统的传递函数(2-4)式求状态空间表达式。 num =[1 5 3]; den =[1 2 3 4]; [A,B,C,D]=tf2ss(num,den) [例2.3] 对上述结果进行验证编程。 %将[例2.2]上述结果赋值给A 、B 、C 、D 阵； A =[-2 -3 -4;1 0 0; 0 1 0]； B =[1;0;0]； C =[1 5 3]；D=0； [num,den]=ss2tf(A ，B ，C ，D,1) [例 2.4] 给定系统1 25.03 2)(2 323++++++=s s s s s s s G ，求系统的零极点增益模型和状态空间模型，并求其单位脉冲响应及单位阶跃响应。 解： num=[1 2 1 3];den=[1 0.5 2 1]; sys=tf(num,den) %系统的传递函数模型 Transfer function: s^3 + 2 s^2 + s + 3 ----------------------------- s^3 + 0.5 s^2 + 2 s + 1 sys1=tf2zp(num,den) %系统的零极点增益模型 sys1 = sys2=tf2ss(sys) %系统的状态空间模型模型；或用[a,b,c,d]=tf2ss(num,den)形式impulse(sys2) %系统的单位脉冲响应 step(sys2) %系统的单位阶跃响应 [例3.1] 对下面系统进行可控性、可观性分析。 解： a=[-1 -2 2;0 -1 1;1 0 -1];b=[2 0 1]';c=[1 2 0] 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;1 0 0 0;0 1 0 0;0 0 1 0]; B=[1;0;0;0];C=[1 7 24 24];D=[0]; [z,p,k]=ss2zp(A,B,C,D,1); Flagz=0; n=length(A); for i=1:n if real(p(i))>0 Flagz=1; end end disp('系统的零极点模型为');z,p,k 系统的零极点模型为 if Flagz==1 disp('系统不稳定'); else disp('系统是稳定的'); end 运行结果为: 系统是稳定的 step(A,B,C,D) ％系统的阶跃响应。 1、使用下列哪一个函数可以产生单位矩阵（ B ） A.zeros B.eye C.rand D.diag 2、下列哪一个函数是求模函数（ D ） A.rem B.sign C.fix D.mod 3、使用下列哪一个函数可以交换矩阵左右对称位置上的元素（ A ） A.fliplr B.flipdim C.flipud D.find 4、使用下列哪一个函数可以比较字符串，且比较时忽略字符的大小写（ D ） A.strncmp B.strcmp C.strncmpi D.strcmpi 5、要利用图形方式显示元胞数组，则应该使用下列哪一个函数（ B ） A.cellfun B.cellplot C.celldisp D.cell2mat 6、下列哪一个函数可以获取结构字段的数据（ B ） A.fieldnames B.getfield C.setfield D.rmfield 7、执行下列哪一条命令后，图形窗体的轴将显示坐标网格线（ A ） A.grid on B.hold on C.grid off D.hold off 8、进行格式化绘图时，使用哪一个函数可以添加图例（ B ） A.title B.legend https://www.wendangku.net/doc/9d7394580.html,bel D.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.surf B.plot C.mesh D.plot3 1、>> A=[1 2 3;4 5 6];B=[2 5;8 3] （2分） B = 2 5 8 3 2、假设向量A=[9 8 7 6 5 4 3 2 1 0] >> A([1:3:5]) （2分） ans = 9 6 >> A([1 3 5]) （2分） ans = 9 7 5 3、>> A=ones(2,2);A(:)=1:4; >> A*A’（2分） ans = 10 14 14 20 >> B=A.*A’ (2分) B = 1 6 6 16 4、使用三元组法，将下列满阵转变为稀疏矩阵 15 0 0 22 0 -15 0 11 3 0 0 0 S= 0 0 0 -6 0 0 91 0 0 0 0 0 0 0 0 28 0 0 解： >> ir=[1 4 2 2 1 3 5 1]; >> jc=[1 1 2 3 4 4 4 6]; >> data=[15 91 11 3 22 -6 28 -15]; >> s=sparse(ir,jc,data,5,6) s = (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(1:24,4,6); >> A(:,[2 3 4])=[ ] A = 1 17 21 2 18 22 3 19 23 4 20 24 6、使用函数struct创建一个结构。 此结构名为Student； 有三个字段，分别为name、age、grade； 有两条记录，分别为’Way’、23、3和’Deni’、21、1 解： >> Student=struct('name',{'Way','Deni'},'age',{23,21},'grade',{3,1}) 7、绘出下幅图 >> plot([1 4 2 8 5]) 三、写出使以下这段文字成为字符串的MATLAB指令。注意保持这段文字的格式。（10分） ‘Are you a teacher?’ ‘No,I’m a student.’ 解： >> a='''Are you a teacher?'''； >> b='''No,I''m a student.'''； >> c=strvcat(a,b) 或 >> char('''Are you a teacher?''','''No,I''m a student.''') 四、建立一个字符串向量sqTKghEad，统计字符串中大写字母的个数，然后再删除大写字母（10分） 解： >> a='sqTKghEad' a = sqTKghEad >> b=find(a>='A'&a<='Z') b = 3 4 7 >> c=length(b) c = 3 >> a(b)=[] a = sqghad （3分） 五、输入一个百分制成绩，要求输出成绩等级A、B、C、D、E。其中90分~100分为A，80分~89分为B，70分~79为C，60分~69分为D，60分以下为E。 要求：1）用if语句实现 2）对不合理的成绩应输出出错信息“输入的成绩不合理”（若成绩出现小数，则只能是“.5”）（16分） 解： clear clc x=input(’请输入一个百分制成绩：’); if rem(x*10,5)==0&&x>=0&&x<=100 if (x>=90&&x<=100) str=[’成绩等级为’,'A']; elseif (x>=80&x<90) str=[’成绩等级为’,'B']; elseif (x>=70&x<80) str=[’成绩等级为’,'C']; elseif (x>=60&x<70) str=[’成绩等级为’,'D']; elseif (x>=0&x<60) str=[’成绩等级为’,'E']; end else str='输入的成绩不合理'; end disp(str) 六、思考题：要绘制出如图所示的图形，请正确填写下列空格。 subplot( 6 , 4 , [2 3 6 7] );plot(1:10);grid on; subplot( 6 , 4 , [10 11 12 14 15 16] );plot(peaks);grid on; subplot( 6 , 4 , [5 9] );plot(membrane);grid on; subplot( 6 , 4 , [17 18 19 21 22 23] );surf(membrane);grid on; 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、求以下变量的值，并在MATLAB中验证。( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ；a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 8 2、下列运算是否合法，为什么？如合法， 结果是多少？ >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 2 5 8 11 >> result4=b*d result4 = 31 22 22 40 49 13 >> result5=[b;c']*d result5 = 31 22 22 40 49 13 -5 -8 7 >> result6=a.*b result6 = 2 8 -3 4 1 5 30 >> result7=a./b result7 = 0.5000 0.5000 -3.0000 4.0000 1.6667 1.2000>> result8=a.c Attempt to reference field of non-structure array. >> result9=a.\b result9 = 2.0000 2.0000 -0.3333 0.2500 0.6000 0.8333 >> result10=a.^2 result10 = 1 4 9 16 25 36 >> result11=2.^a result11 = 2 4 8 16 32 64 3、用MATLAB求解下面的的方程组。 (1) ? ? ? ? ? ? ? ? ? ? ? ? - = ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? - - - - 1 7 4 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 4 3 2 1 x x x x >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B (2) ? ? ? ? ? ? ? = - + + = - - = - + + = + + 5 6 5 3 3 3 3 2 8 2 1 w z y x w y x w z y x z y x >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5] >> x2=inv(A1)*B2 4、已知 ? ? ? ? ? ? ? ? ? ? ? ? - - - - = 13 2 3 1 5 11 2 2 2 3 15 9 2 1 2 7 A MATLAB的GUI 程序设计 Chapter 8: Design of MATLAB of GUI program GUI(Graphical User Interfaces):由各种图形对象组成的用户界面，在这种用户界面下，用户的命令和对程序的控制是通过“选择”各种图形对象来实现的。 基本图形对象分为控件对象和用户界面菜单对象，简称控件和菜单。 一. 控件对象及属性(Object and its attributes of controller)) 1. GUI控件对象类型(The mode of controller object) 控件对象是事件响应的图形界面对象。当某一事件发生时，应用程序会做出响应并执行某些预定的功能子程序（Callback）. 控件对象及其功能：（表7—1） 2. 控件对象的描述(Description of controller object) MATLAB中的控件大致可分为两种，一种为动作控件，鼠标点击这些控件时会产生相应的响应。一种为静态控件，是一种不产生响应的控件，如文本框等。 每种控件都有一些可以设置的参数，用于表现控件的外形、功能及效果，既属性。属性由两部分组成：属性名和属性值，它们必须是成对出现的。 （1）按钮(Push Buttons)：执行某种预定的功能或操作； （2）开关按钮(Toggle Button)：产生一个动作并指示一个二进制状态（开或关），当鼠点击它时按钮将下陷，并执行callback（回调函数）中指定的内容，再次点击，按钮复原，并再次执行callback 中的内容； （3）单选框(Radio Button)：单个的单选框用来在两种状态之间切换，多个单选框组成一个单选框组时，用户只能在一组状态中选择单一的状态，或称为单选项； （4）复选框(Check Boxes)：单个的复选框用来在两种状态之间切换，多个复选框组成一个复选框组时，可使用户在一组状态中作组合式的选择，或称为多选项； （5）文本编辑器(Editable Texts)：用来使用键盘输入字符串的值，可以对编辑框中的内容进行编辑、删除和替换等操作； （6）静态文本框(Static Texts)：仅仅用于显示单行的说明文字； （7）滚动条(Slider)：可输入指定范围的数量值； result5 = ( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ；a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 8 2、下列运算是否合法，为什么如合法， 结果是多少 >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 2 58 11 >> result4=b*d result4 = 31 22 22 40 49 13 31 22 22 40 49 13 -5 -8 7 >> result6=a.*b result6 = 2 8 -3 415 30 >> result7=a./b result7 = >> result8= Attempt to reference field of non-structure array. >> result9=a.\b result9 = >> result10=a92 result10 = 1 4 9 16 25 36 >> resultl 1=29a result11 = 2 4 8 16 32 64 >> result5=[b;c']*d 3、用MATLAB求解下面的的方程组。 1、求以下变量的值，并在MATLAB^验证。 1 2 x1 3 2 x2 11 5 x3 2 1 3 x4 >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5] >> x2=inv(A1)*B2 7 2 1 2 9 15 3 2 2 2 11 5 1 3 2 13 (1)求矩阵A的秩(rank) (2)求矩阵 A 的行列式(determinant) (3)求矩阵 A 的逆(inverse) (4)求矩阵 A 的特征值及特征向量 (eigenvalue and eigenvector) >> A3=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> r=rank(A3) >> b=inv(A3) >> a=det(A3) >> [V,D]=eig(A3) 10 n 10 查看y 的值) m1=0; for m=-10:10 m仁m1+2^m; end m1 m1 = 6、求分段函数的值。 用if 语句实现，算出下列表中x 对应的y 值。 x=input('enter x='); if x<0 y=x A2+x-6; elseif x>=0&&x<5 y=xA2-5*x+6; else y=xA2-x-1; end y 7、分别用if 和switch 语句实现，将百分 制成绩转换为成绩等级A、B、C、D、E。 其中90~1 00分为A，80~89 分为B，70~79 分为C，60~69 分为D，60 分以下为E。 对超出百分制范围的成绩，给出错误提示 信息。 if 结构程序： x=input('please enter score='); if x>=90&&x<=100 9 2 10 disp('A') 7 2 9 15 (1) 2 2 1 3 4 7 1 0 A 4、已知 2n 2 10 29 第一套 华东交通大学2015-2016学年第一学期Matlab期末考试 一、填空题（4*6=24分） 1、在MATLAB命令窗口中的“>>”标志为MATLAB的提示符，“│”标志为提示符。 2、符号表达式sin(2*a+t)+m 中独立的符号变量为__ _。 3、在通常情况下，左除x=a\b是的解，右除x=b/a是的解，一般情况下，。 4、为了使两个plot的图形在同一个坐标显示，可以使用_____命令进行图形保持；可以使用_ __命令为图形添加网格。 5、倘若要是对x进行赋值，从5到25，间隔是0.1,为；倘若要是对x进行赋值，从2到20，中间间隔100个点,为 6、A=[1,2;3,4];B=[1,0;0,1]; A*B= _ , A.*B=_ _； 二、选择题（4*6=24分） 7、如果x=1: 2 : 8,则x(1)和x(4)分别是( ) A．1，8 B．1, 7 C．2, 8 D．2, 7 8、运行如下程序后，输入９回车，命令窗口（command windows）显示的结果为( ) c=input('请输入一个字符','s'); if c>='A' & c<='Z' disp(setstr(abs(c)+abs('a')-abs('A'))); elseif c>='a'& c<='z' disp(setstr(abs(c)- abs('a')+abs('A'))); elseif c>='0'& c<='9' disp(abs(c)-abs('0')); else disp(c); end A. 9 B. 8 C. 7 D. 6 9、MATLAB表达式2*2^3^2的结果是() A．128 B．4096 C. 262144 D．256 10、在循环结构中跳出循环，执行循环后面代码的命令为( ) (A) return (B) break (C) continue (D) keyboard 11、在图形指定位置加标注命令是（） A. title(x,y,?y=sin(x)?); B. xlabel(x,y,?y=sin(x)?); C. text(x,y,?y=sin(x)?); D. legend(x,y,?y=sin(x)?); 12、下列哪个变量的定义是不合法的（）。 (A) abcd-3 (B) xyz_3 (C) abcdef (D) x3yz 三、程序题（52分） 13（10分）、在同一窗口绘制6个子图，已知t=0.1:0.1:2*pi，每个子图分别绘制（其中y1=sin(t),y2=cos(2t),y3=sin(t)+cos(t),y4=3t, y5=2ln(t),y6=et）； 14（10分）、某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price来表示)： price<200 没有折扣 200≤price<500 3%折扣 500≤price<1000 5%折扣 1000≤price<2500 8%折扣 2500≤price<5000 10%折扣 复习题 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= 1、求以下变量的值，并在MATLAB 中验证。 ( 1 ) a = 1 : 2 : 5 a = 1 3 5 ( 2 ) b = [ a' , a' , a' ；a ] b = 1 1 1 3 3 3 5 5 5 1 3 5 ( 3 ) c = a + b ( 2 , : ) c = 4 6 8 2、下列运算是否合法，为什么？如合法，结果 是多少？ >> result2=a*b Error using * Inner matrix dimensions must agree. >> result3=a+b result3 = 3 6 2 5 8 11 >> result4=b*d result4 = 31 22 22 40 49 13 >> result5=[b;c']*d result5 = 31 22 22 40 49 13 -5 -8 7 >> result6=a.*b result6 = 2 8 -3 4 1 5 30 >> result7=a./b result7 = 0.5000 0.5000 -3.0000 4.0000 1.6667 1.2000 >> result8=a.c Attempt to reference field of non-structure array. >> result9=a.\b result9 = 2.0000 2.0000 -0.3333 0.2500 0.6000 0.8333 >> result10=a.^2 result10 = 1 4 9 16 25 36 >> result11=2.^a result11 = 2 4 8 16 32 64 3、用MATLAB 求解下面的的方程组。 (1)????????????-=?????????????????????????----01741323151122231592127 4321x x x x >> A=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13] >> B=[4 7 -1 0] >> B=B' >> x=inv(A)*B (2)???????=-++=--=-++=++56533332821 w z y x w y x w z y x z y x >> A1=[1 1 1 0;1 2 1 -1;2 -1 0 -3;3 3 5 -6] >> B2=[1;8;3;5] >> x2=inv(A1)*B2 4、已知????????????----=1323151122231592127A (1)求矩阵A 的秩(rank) MATLAB期末考试试卷及其参考答案 一、填空题（每空1分，20分） 1、MATLAB常用操作界面包括、工作空间窗口、 、、内存数组编辑器、M文件编辑/调试器、帮助导航/浏览器、图形窗口等。 2、MATLAB中Inf或inf表示、NaN或nan表示、 nargout表示。 3、MATLAB中逗号主要用作； 用作输入量与输入量之间的分隔符；。 4、工作空间浏览器主要用于内存变量的、和。 5、MATLAB实现将全下标转换为单下标的指令为、据单下标换 算出全下标的指令为。 6、二维数组的标识有、、“逻辑1”标识。 7、在一个元胞数组A中寻访第2行第3列元胞元素用；寻访数 组第2行第3列元胞中所存的内容用。 8、4、MATLAB中clf用于、clc用于、 clear用于。 二、简答题（每题5分，共20分） 1、简述MATLAB历史指令窗的主要作用。 2、简述空数组的功用。 3、简述MATLAB函数的基本结构。 4、简述绘制二维图形的一般步骤。 三、阅读程序并回答问题（每题4分，共28分） 1、写出下列指令运行结果。 A=zeros(2,4); A(:)=1:8; s=[2 3 5]; A(s) Sa=[10 20 30]' A(s)=Sa 2、写出下列指令运行结果。 A=reshape(1:16,2,8) reshape(A,4,4) s=[1 3 6 8 9 11 14 16]; A(s)=0 3、写出下列指令运行结果。 A=[1,2;3,4]; B=[-1,-2;2;1]; S=3; A.*B A*B S.*A S*B 4、下面的函数主要完成什么功能？ function f=factor(n) if n<=1 f=1; else f=factor(n-1)*n; end 5、写出下列指令运行结果。 ch=‘ABc123d4e56Fg9’; subch=ch(1:5) revch=ch(end:-1:1) k=find(ch>=‘a’&ch<=‘z’); 实验二一维二维数组的创建和寻访 一、实验目的 1、掌握一维数组、二维数组创建和寻访的几种方法。 2、区别数组运算和矩阵运算的差别。 3、熟悉执行数组运算的常用数组操作函数。 4、掌握数组运算中的关系和逻辑操作及常用的关系、逻辑函数。 5、掌握“非数”、“空”数组在MA TLAB中的应用。 二、实验主要仪器与设备 装配有MA TLAB7.6软件的计算机 三、预习要求 做实验前必须认真复习第三章MATLAB的数值数组及向量化运算功能。 四、实验内容及实验步骤 1、一维数组的创建方法有哪几种？举例说明。 答：一维数组的创建方法有： ①递增/递减型一维数组的创建：冒号生成法：x=a：inc：b 线性（或对数）定点法：x=linspace(a,b,n)，x=logspace(a,b,n) ②逐个元素输入法：如x=[0.1，sin(pi/5)，-exp(-3)，-2*pi] ③运用MA TLAB函数生成法：例ones，rand等。 2、输入以下指令，并写出运行结果。本例演示：数组元素及子数组的各种标识和寻访格式；冒号的使用；end的作用。 A=zeros(2,6) %创建(2×6)的全零数组 A(:)=1:12 %赋值号左边：单下标寻访(2×6) 数组A的全部12个元素 %赋值号右边：拥有12个元素的一维数组 A(2,4) %双下标：A数组的第2行第4列元素 A(8) %单下标：数组A的第8个元素 A(: , [1,3]) %双下标：显示A的“第1列和第3列上全部行的元素” A([1, 2, 5, 6]') %单下标：把A数组第1,2,5,6个元素排成列向量 A(: , 4:end) %双下标：显示A的“从第4起到最后一列上全部行的元素” %在此end用于“列标识”，它表示“最后一列” A(2,1:2:5)=[-1, -3, -5] %把右边的3个数分别赋向A数组第2行的第1,3,5个元素位置 B=A([1, 2, 2, 2], [1, 3, 5]) %取A数组的1,3,5列的第1行元素作为B的第1行 %取A数组的1,3,5列的第2行分别作为B的第2,3,4行 L=A<3 %产生与A维数相同的“0，1”逻辑数组 A(L)=NaN %把逻辑1标识的位置上的元素赋为“非数” 运行结果： A = 0 0 0 0 0 0 0 0 0 0 0 0 第1章MATLAB基础 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 答：在MATLAB操作桌面上有五个窗口。在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口。在独立窗口的Desktop菜单中选择Dock...项就可以将独立的窗口重新放置到桌面上。 1.5 如何启动M文件编辑/调试器？ 答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选 定的内容拷贝到剪贴板中、将选定内容直接拷贝到M 文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path 菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB 运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB 中有几种获得帮助的途径？ 答：（1）帮助浏览器：选择view 菜单中的Help 菜单项或选择Help 菜单中的MATLAB Help 菜单项可以打开帮助浏览器。 （2）help 命令：在命令窗口键入“help ” 命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息。 （3）lookfor 命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数。 （4）模糊查询：输入命令的前几个字母，然后按Tab 键，就可以列出所有以这几个字母开始的命令和函数。 注意：lookfor 和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help 命令显示详细信息。 第2章 MATLAB 矩阵运算基础 2.1 在MATLAB 中如何建立矩阵?? ? ? ??194375，并将其赋予变量a ？ 答：在Command Window 窗口输入操作: 实验一：Matlab操作环境熟悉 一、实验目的 1．初步了解Matlab操作环境。 2．学习使用图形函数计算器命令funtool及其环境。 二、实验内容 熟悉Matlab操作环境，认识命令窗口、内存工作区窗口、历史命令窗口；学会使用format命令调整命令窗口的数据显示格式；学会使用变量和矩阵的输入，并进行简单的计算；学会使用who和whos命令查看内存变量信息；学会使用图形函数计算器funtool，并进行下列计算： 1．单函数运算操作。 求下列函数的符号导数 (1) y=sin(x); (2) y=(1+x)^3*(2-x); 求下列函数的符号积分 (1) y=cos(x); (2) y=1/(1+x^2); (3) y=1/sqrt(1-x^2); (4) y=(x-1)/(x+1)/(x+2); 求反函数 (1) y=(x-1)/(2*x+3); (2) y=exp(x); (3) y=log(x+sqrt(1+x^2)); 代数式的化简 (1) (x+1)*(x-1)*(x-2)/(x-3)/(x-4); (2) sin(x)^2+cos(x)^2; (3) x+sin(x)+2*x-3*cos(x)+4*x*sin(x); 2．函数与参数的运算操作。 从y=x^2通过参数的选择去观察下列函数的图形变化 (1) y1=(x+1)^2 (2) y2=(x+2)^2 (3) y3=2*x^2 (4) y4=x^2+2 (5) y5=x^4 (6) y6=x^2/2 3．两个函数之间的操作 求和 (1) sin(x)+cos(x) (2) 1+x+x^2+x^3+x^4+x^5 乘积 (1) exp(-x)*sin(x) matlab练习题和答案 控制系统仿真实验 Matlab部分实验结果 目录 实验一 MATLAB基本操 作 ............................................................................................ 1 实验二 Matlab编程 .................................................................................................... 5 实验三Matlab底层图形控制 (6) 实验四控制系统古典分析.............................................................................................12 实验五控制系统现代分析 . (15) 实验六 PID控制器的设 计 ...........................................................................................19 实验七系统状态空间设计.............................................................................................23 实验九直流双闭环调速系统仿真 . (25) 实验一 MATLAB基本操作 1 用MATLAB可以识别的格式输入下面两个矩阵 1233，，1443678，i，，,,2357,,,,2335542，i,,,, A,1357B,,,2675342， i,,3239,,,,189543，，,,1894，， 再求出它们的乘积矩阵C，并将C矩阵的右下角2×3子矩阵赋给D矩阵。赋值完成后，调 用相应的命令查看MATLAB工作空间的占用情况。 A=[1,2,3,3;2,3,5,7;1,3,5,7;3,2,3,9;1,8,9,4]; MATLAB 期末上机考试试题带答案版 姓名： 学号： 成绩： 1.请实现下图： 50 100 150 200 250 -1 -0.8-0.6-0.4-0.200.20.40.60.81x y y=sin(x) x=linspace(0,8*pi,250); y=sin(x); plot(x,y) area(y,-1) xlabel('x') ylabel('y') title('y=sin(x)') 2.请实现下图： x=linspace(0,2*pi,100); y1=sin(x); subplot(2,2,1) plot(x,y1,'k--') grid on xlabel('x') ylabel('y') title('sin(x)') legend('y=sin(x)') y2=cos(x); subplot(2,2,2) plot(x,y2,'r--') grid on xlabel('x') ylabel('y') title('cos(x)') legend('y=cos(x)') y3=tan(x); subplot(2,2,3) plot(x,y3,'k-') grid on xlabel('x') ylabel('y') title('tan(x)') legend('y=tan(x)') y4=cot(x); subplot(2,2,4) plot(x,y4) grid on xlabel('x') ylabel('y') title('cot(x)') legend('y=cot(x)') 3.解方程组： a=[3 2 1;1 -1 3;2 4 -4]； b=[7;6;-2] ； x=a\b 4.请实现下图： 2 4 6 8 10 12 14 -1-0.8-0.6-0.4-0.200.20.40.60.81x y '自己的姓名' sin(x)sin(2*x) x=linspace(0,4*pi,1000); y1=sin(x); 1、完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。(2) 建立一个字符串向 量，删除其中的大写字母。 （1） m=100:999; n=find(mod(m,21)==0); length(n) ans = 43 （2）>> ch='Maybe One Day' p=find(ch>='A'&ch<='Z') ch(p)=[] ch = Maybe One Day 2、自行产生一个5行5列的数组，分别得到最中间的三行三列矩阵、右下角2行2列矩阵， 奇数行矩阵、奇数列矩阵、奇数行奇数列矩阵。 >> t=rand(5)%生成矩阵 A=t(2:4,2:4)%中间三行散列矩阵 B=t(4:5,4:5)%右下角两行两列矩阵 C=t(1:2:end,:)%奇数行矩阵 D=t(:,1:2:end)%奇数列矩阵 E=t(1:2:end,1:2:end)%奇数行列矩阵 3、求方程组的根 syms x y z [X Y Z]=solve('x+4*y-3*z=2','2*x+5*y-z=11','x+6*y+z=12',x,y,z) 4、已知矩阵A=[1 2;3 4],运行指令B1=A.^(0.5), B2=A^(0.5), 可以观察到不同运算方法所得结果不同。（1）请分别写出根据B1, B2恢复原矩阵A的程序。（2）用指令检验所得的两个恢复矩阵是否相等(利用norm(…,’fro’)指令,误差矩阵F-范数，接近eps量级，认为实际相等)。 5、先运行clear,format long,rng('default'),A=rand(3,3)，然后根据A写出两个矩阵：一个对角 阵B，其相应元素由A的对角元素构成；另一个矩阵C，其对角元素全为0，而其余元素与对应的A阵元素相同（diag指令的使用）。 >> format long >> rand('twister',1) >> A=rand(3,3) A = 0.417022004702574 0.302332572631840 0.186260211377671 一、填空 1、编写M 函数文件时必须以关键字 开头。 2、命令clc 的含义是 ，命令clear 的含义是 。 3、若数组11[,,],[,,]n n a a αβββ== ,则两数组相乘的命令为： ,计算 数 组α的k 次幂的命令为： 。 4、Matlab 的最基本的运算是： ， ， 。 5、设矩阵111011212-?? ?= ? ??? A ,求矩阵A 的行列式的Matlab 命令为 ；求矩阵A 的行最简形式的命令为 。 6、设矩阵123456312?? ?= ? ??? A ,输入矩阵A 的Matlab 命令是 或 。 7、假定A 是一个10阶方阵，选取矩阵A 的第三行第二列的指令是 ，选取矩阵A 第四列和第七列的指令是 ，删除矩阵A 的第二行的指令是 。 8、Matlab 可以输入字母、汉字，但是M 文件中标点符号必须在 状态下输入。 9、若x=[2,3,4],y=[1,2;4,5；2,3];则plot(y)画出 条曲线。 10、若A ＝[2,3,4;4,5,6;6,7,8]，则A .^2= 11、命令subplot(m,n,k)的功能是把图形窗口分为 个子图，并把第 个子图作为当前图形窗口。 12、A ＝502103?? ??? ，则rank （A ）的结果是 。 二、选择题 1、在Matlab 操作中要列出当前内存变量的名称、大小、类型等信息应使用（ ）命令。 (A) List (B) Browse (C) Who (D) Whos 2、下列M 文件中命名错误的是：（ ）。 (A) as.m (B) yuli4.m (C) 4na.m (D) r45u.m 3、在matlab 中，若x=1:0.1:2,则正确输入表达式2sin ln x x e x π++的命令是：（）。 (A) x.^2+e.^x*sin π+lnx (B) x.^2+exp(x).*sin(pi)+ln(x) (C) x.^2+exp(x).*sin(pi)+log(x) (D) x.^2+exp(x).*sin(π)+lnx 4、若A ＝512143-?? ??? ，则A （1,2)和A(1,:)结果是：（ ）。 (A) -1和5 -1 2 (B) -1和5 1 (C) 1和5 -1 2 (D) 1和5 1 5、下面那个运算符为逻辑运算符（ ）。 (A) + (B) < (C) * (D) & 6、下面哪一种结构表示分支结构（ ） (A) for-end (B) while-end (C) if-else-end (D) switch-end 7、下列叙述正确的是：（ ） (A) M 文件中的％的含义是标明注释。 (B) 在matlab 所输入的命令后如果输入逗号，则不显示执行结果。 (C) matlab 的变量必须先定义其类型才可使用。 一个不错的Matlab的gui界面设计实例 %非常漂亮的日历， function CalendarTable; % calendar 日历 % Example: % CalendarTable; S=datestr(now); [y,m,d]=datevec(S); % d is day % m is month % y is year DD={'Sun','Mon','Tue','Wed','Thu','Fri','Sat'}; close all figure; for k=1:7; uicontrol(gcf,'style','text',... 'unit','normalized','position',[0.02+k*0.1,0.55,0.08,0.06],... 'BackgroundColor',0.6*[1,1,1],'ForegroundColor','b',... 'String',DD(k),'fontsize',16,'fontname','times new roman'); end h=1; ss='b'; qq=eomday(y,m); for k=1:qq; n=datenum(y,m,k); [da,w] = weekday(n); if k==d; ss='r'; end uicontrol(gcf,'style','push',... 'unit','normalized','position',[0.02+da*0.1,0.55-h*0.08,0.08,0.06],... 'BackgroundColor',0.6*[1,1,1],'ForegroundColor',ss,... 'String',num2str(k)); ss='b'; if da==7; h=h+1; 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个项目可供投资者选择。每个项目可以重复投资，根据专家经验，对每个项目投资总额不能太高，且有个上限。这些项目所需要的投资额已经知道，在一般情况下，投资一年后各项目所得利润也可估计出来（见表一）， 表一: 投资项目所需资金及预计一年后所得利润(单位：万元) WORD 格式 一、填空题（每空2分，总共30分） 1.已知A=[096;130];B=[143;150];写出下列各指令运行 的结果。 A&B 的运行结果ans=； A./B 的运行结果ans=。 方 阵的 命 令 为z e r o s (4 )1方阵的命 令为ones(3)。 3.A=rand(2,5);b=size(A);c=length(A)；则b 和c 的值分别为 和[2,5],[5]。 4.标点符号分号可以使命令行不显示运算结果，百分号用来表 示该行为注释行。 5.用if 判断语句判断80≥x>60,在MATLAB 中if 语言后的判断 应写为ifx<=80&x>60。 共10页 6.P,Q分别是个多项式的系数矢量，求P对应的多项式的积分（对应的常数项为K），使用的命令是polyint(P,K)；求P/Q的解，商和余数分别保存在k和r，使用的 命令是[k,r]=deconv(P,Q)； 7.为了使两个plot的图形在同一个坐标显示，可以使用holdon命令进行图形 保持；可以使用gridon命令为图形添加网格。 8.MATLAB的工作空间中有三个变量v1,v2,v3，写出把它们保存到文件my_data.mat 中的指令savemy_data；写出把my_data.mat文件中的变量读取到MATLAB 工作空间内的指令loadmy_data。 二、选择题（每空2分，总共20分） 1.下列哪个变量的定义是不合法的A (A)abcd-3(B)xyz_3(C)abcdef(D)x3yz 2.下列哪条指令是求矩阵的行列式的值C (A)inv(B)diag(C)det(D)eig 3.在循环结构中跳出循环，执行循环后面代码的命令为B (A)return(B)break(C)continue(D)keyboard 4.清空Matlab工作空间内所有变量的指令是C (A)clc(B)cls(C)clear(D)clfMATLAB数学实验练习题
MATLAB实验题答案
matlab中GUI设计
MATLAB实验题答案
Matlab 期末考试题库(共12套卷)
matlab数学实验复习题(有标准答案)
MATLAB实验题答案
MATLAB期末考试卷和参考答案..
MATLAB实验题目及答案
matlab基础练习题及答案讲解
Matlab实验指导书(含答案)汇总
matlab练习题和答案
MATLAB期末上机试题带答案
MATLAB实验上机易错题汇总
matlab期中测验
Matlab的gui界面设计实例练习
matlab数学实验练习题
MATLAB期末考试试卷及答案(20200412004949)