Matlab考题题整理带答案
MATLAB考试试题(1)
产生一个1x10的随机矩阵，大小位于（-55），并且按照从大到小的顺序排列好！（注：要程序和运行结果的截屏）
答案：
a=10*rand(1,10)-5;
b=sort(a,'descend')
1.请产生一个100*5的矩阵，矩阵的每一行都是[12345]
2.已知变量：A=’ilovematlab；’B=’matlab’请,找出：
（A）B在A中的位置。
（B）把B放在A后面，形成C=‘ilovematlabmatlab’
3.请修改下面的程序，让他们没有for循环语句！
A=[123;456;789];
[rc]=size(A);
fori=1:1:r
forj=1:1:c
if(A(i,j)>8|A(i,j)<2)
A(i,j)=0;
end
end
end
4.请把变量A=[123;456;789]写到文件里(output.xls)，写完后文件看起来是这样的
123456789
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=[12345],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=[34-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)/(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-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.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）在同一图形窗口绘制曲线与，要求曲线为黑色点画线，曲线为红色虚线圆圈；并在图中恰当位置标注两条曲线的图例；给图形加上标题“y1andy2。”
五、（本题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）创建符号函数symsx
（2）求该符号函数对的微分；
（3）对趋向于求该符号函数的极限；
（4）求该符号函数在区间上对的定积分；
（5）求符号方程的解。
四、（本题20分）编写MATALAB程序，完成下列任务（将程序保存为test04.m文件）：
（1）在区间上均匀地取20个点构成向量；
（2）分别计算函数与在向量处的函数值；
（3）在同一图形窗口绘制曲线与，要求曲线为黑色点画线，曲线为红色虚线圆圈；并在图中恰当位置标注两条曲线的图例；给图形加上标题“y1andy2。”
五、（本题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=[010;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阵；
A=[-2-3-4;100;010]；B=[1;0;0]；C=[153]；D=0；
[num,den]=ss2tf(A，B，C，D,1)
[例2.4]给定系统
32
s2ss3
G(s)，求系统的零极点增益模型和状态空间模型，并求其
32
s0.5s2s1
单位脉冲响应及单位阶跃响应。
解：
num=[1213];den=[10.521];
sys=tf(num,den)%系统的传递函数模型Transferfunction:
s^3+2s^2+s+3
-----------------------------
s^3+0.5s^2+2s+1
sys1=tf2zp(num,den)%系统的零极点增益模型sys1=
sys2=tf2ss(sys)%系统的状态空间模型模型；或用[a,b,c,d]=tf2ss(num,den)形式impulse(sys2)%系统的单位脉冲响应
step(sys2)%系统的单位阶跃响应
[例3.1]对下面系统进行可控性、可观性分析。
解：
a=[-1-22;0-11;10-1];b=[201]';c=[120]
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=[1;0;0;0];C=[172424];D=[0];
[z,p,k]=ss2zp(A,B,C,D,1);
Flagz=0;
n=length(A);
fori=1:n
ifreal(p(i))>0
Flagz=1;
end
end
disp('系统的零极点模型为');z,p,k
系统的零极点模型为
ifFlagz==1
disp('系统不稳定');
elsedisp('系统是稳定的');
end
运行结果为:
系统是稳定的
step(A,B,C,D)％系统的阶跃响应
资源与环境工程学院2008级硕士研究生《MatLab及其应用》试题注意，每题的格式均须包含3个部分
a.程序（含程序名及完整程序）：
b.运行过程：
c.运行结果：
(1)求解线性规划问题：
minZ 4x
1 x
2
7x
3
s.t. x
1 x
2
x
3
5
3x
1 x
2
x
3
4
x 1 x
2
4x
3
7
x, 1 x
2
问各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('MINz=%f\n',z);
运行结果：线性规划问题求解
Optimizationterminatedsuccessfully.
x=
2.2500
6.7500
4.0
MINz=25.750000
0.5x，x2 (2)编写一个函数，使其能够产生如下的分段函数：f(x)1.50.25x2x6
，，
0.56x
，并调用此函数，绘制在x[0，2]范围的f(x)f(x2)曲线。（10分）
答：functiony=f(x)
ifx<=2
y=0.5*x;
elseifx>6
y=0.5;
elsey=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 分）
22
xy
cos(x，x，（曲线图）；②(-2x2,-4y4)
①)[]f(x,y)（曲面图）。
；
22
22
24
答：>>subplot(2,2,2)
>>ezplot('(cos(x))^(1/2)',[-pi/2pi/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;
fori=1:m
forj=1:n
ifA(i,j)==0
sig=sig+1;
endendend
请输入一个矩阵[012;102;000]
A=
012
102
000>>sig
sig=
5
(5)向量A[a,a,,a]
11n.写一段程序,找出A中的最小元素（10分）
答：A=input('请输入一个向量')
[m,n]=sizeA min=A(1,n); fori=1:n
ifA(1,i) WORD格式 min=A(1,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.396.84] plot(x,y); xlabel('pH'); ylabel('吸附量') 图 1 x=[0.50.751.01.251.5]y=[6.05 6.195.334.694.02] plot(x,y); xlabel('秸秆用量g'); ylabel('吸附量') 图2 通过数据图，得到比较理想的实验条件pH和秸秆用量，接下来做动力学和等温线。 >x=[0.1670.5123458] y=[0.0332010.0860590.1697790.3220610.4807690.6441220.809061 10.] plot(x,y); xlabel('时间t'); ylabel('时间/吸附量') 图 3 x=[0.23630.154960.136190.129060.133730.13315]y=[0.25218 0.047070.020140.012670.008810.00706] plot(x,y); xlabel('1/吸附量'); ylabel('1/平衡浓度') 图 4 x=[0.626540.809770.865850.88920.873770.87564] y=[0.598291.32731.695891.897372.055032.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=[010;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阵； A=[-2-3-4;100;010]；B=[1;0;0]；C=[153]；D=0； [num,den]=ss2tf(A，B，C，D,1) [例2.4]给定系统 32 s2ss3 G(s)，求系统的零极点增益模型和状态空间模型， 32 s0.5s2s1 并求其单位脉冲响应及单位阶跃响应。 解： num=[1213];den=[10.521]; sys=tf(num,den)%系统的传递函数模型 Transferfunction: s^3+2s^2+s+3 ----------------------------- s^3+0.5s^2+2s+1 sys1=tf2zp(num,den)%系统的零极点增益模型sys1= sys2=tf2ss(sys)%系统的状态空间模型模型；或用[a,b,c,d]=tf2ss(num,den)形式 impulse(sys2)%系统的单位脉冲响应 step(sys2)%系统的单位阶跃响应 [例3.1]对下面系统进行可控性、可观性分析。 解： a=[-1-22;0-11;10-1];b=[201]';c=[120] 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=[1;0;0;0];C=[172424];D=[0]; [z,p,k]=ss2zp(A,B,C,D,1); Flagz=0; n=length(A); fori=1:n ifreal(p(i))>0 Flagz=1; end end disp('系统的零极点模型为');z,p,k 系统的零极点模型为 ifFlagz==1 disp('系统不稳定'); elsedisp('系统是稳定的'); end 运行结果为: 系统是稳定的 step(A,B,C,D)％系统的阶跃响应。 2．一单位反馈控制系统，若要求： （1）跟踪单位斜坡输入时系统的稳态误差为2。 （2）设该系统为三阶，其中一对复数闭环极点为-1±j1 求满足上述要求的开环传递函数。 利用MATALAB 完成以下操作： （1）求控制系统的单位阶跃响应、单位斜坡响应、单位加速度响应； （2）绘制控制系统的根轨迹，并判断系统的稳定性； （3）绘制控制系统的奈奎斯特图和Bode 图。 解：由于控制系统为单位反馈系统 则 H （s ）=1 由条件一可知，输入信号R(s)= 21s 时，稳态误差为： 2) (10)()(1)(==→+=s sG s s H s G s sR Ess 由条件二可知，系统为三阶且有。两个极点s1=-1+j1 ；s2=-1-j1 则可得，开环传递函数为 ) 222(1)(++=s s s s G （1）求控制系统的单位阶跃响应： G=tf([1],[1, 2, 2, 0]); Step(G) title('单位阶跃响应'); grid (2)控制系统的单位斜坡响应：z0=0; p0=[0, -1+j,-1-j]; k0=1; [num0,den0]=zp2tf(z0,p0,k0); g0=tf(num0,den0) g=feedback(g0,1); num=g.num{1}; den=[g.den{1},0]; sys1=tf(num,den) step(sys1) (3)控制系统的单位加速度响应：z0=0; p0=[0, -1+j,-1-j]; k0=1; [num0,den0]=zp2tf(z0,p0,k0); g0=tf(num0,den0) g=feedback(g0,1); num=g.num{1}; den=[g.den{1} 0 0]; sys2=tf(num,den) step(sys2) 一、填空题 1、求可逆矩阵A的逆矩阵的指令是 inv(A) 2、A=[0 2 3 0],B=[4 5 0 0],则X=A&B的显示结果为 0 1 0 0 3、求矩阵A的转置矩阵的指令是： A’ 4、x=[10 -20 30]; a=sign(x);显示结果为 1 -1 1 5、Notebook是Matlab的文字处理工具，启动Notebook有两种方法，即一 种是在 Word 中启动，另一种是在 Matlab 中启动。 6、Matlab7.3提供了非常强大的帮助系统，包括帮助命令、帮助 窗口、HTML格式帮助、pdf格式帮助和帮助演示等。 7、clc命令和clear命令的区别是 clc清除命令窗口所有显示内容； Clear清除所有显示的变量 . 8、命令窗口中的type命令的功能为显示文件的内容 . 9、matlab7.3提供了三种程序调试的工具，包括直接检测调试、专用调试命 令和使用 M文件编辑/调试器。 10、matlab除了最常用的主函数和子函数外，还提供了嵌套函数、私有函数、 重载函数和匿名函数。 11、在matlab7.3中，创建句柄图形对象的命令格式为____________ h_obj=funname(‘PropertyName’,PropertyValue,……)_ 12、计算5*e3*|cos175°|语句为 5*exp(3)*abs(cos(175*pi/180)). 13、创建随机矩阵的指令用rand,创建稀疏矩阵用sparse 14、为图形添加网格使用grid on，使用legend命令添加图例。 1、MATLAB常用操作界面包括命令窗口、工作空间窗口（浏览器）、 命令历史窗口、当前目录窗口、内存数组编辑器、M文件编辑/调试器、 帮助导航/浏览器、图形窗口等。 2、MATLABInf或inf表示无穷大、NaN或nan表示不是一个数、nargout表 示函数输出宗量数目。 3、MATLAB中逗号主要用作要显示计算结果的指令与其后指令的分隔； 用作输入量与输入量之间的分隔符；用作数组元素分隔符号。 4、工作空间浏览器主要用于内存变量的查阅、保存和编辑。 5、MATLAB实现将全下标转换为单下标的指令为Sub2ind、据单下标换算出全 下标的指令为Ind2sub。 6、二维数组的标识有“全下标”标识、“单下标”标识、“逻辑1”标识。 7、在一个元胞数组A中寻访第2行第3列元胞元素用A(2,3)；寻访数组第 2行第3列元胞中的内容用A{2,3}。 8、MATLAB中clf用于清除图形窗、clc用于清除指令窗中显示内容、clear 用于清除MATLAB工作空间中保存的变量。 1.已知A=[0 9 6;1 3 0];B=[1 4 3;1 5 0];写出下列各指令运行的结果。 A & B 的运行结果ans= ； A ./ B 的运行结果ans= 。 2. 产生4阶全0方阵的命令为；产生3阶全1方阵的命令为。 3. A=rand(2,5);b=size(A);c=length(A)；则b和c的值分别为和。 第1章MATLAB概论 1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？ 答：MATLAB语言有一下特点： 1、起点高。 2、人机界面适合科技人员。 3、强大而简易的作图功能。 4、智能化程度高。 5、功能丰富，可扩展性强。 1.2 MATLAB系统由那些部分组成？ 答：MATLAB系统主要由开发环境、MATLAB数学函数库、MMATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？ 答：在安装MATLAB时，安装内容由选择组件窗口中各复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即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 函数名”可以得到 第2章符号运算 习题2及解答 1 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度” 对象，还是“符号”符号对象 3/7+; sym(3/7+; sym('3/7+'); vpa(sym(3/7+) 〖目的〗 不能从显示形式判断数据类型，而必须依靠class指令。 〖解答〗 c1=3/7+ c2=sym(3/7+ c3=sym('3/7+') c4=vpa(sym(3/7+) Cs1=class(c1) Cs2=class(c2) Cs3=class(c3) Cs4=class(c4) c1 = c2 = 37/70 c3 = c4 = Cs1 = double Cs2 = sym Cs3 = sym Cs4 = sym 2 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认 为是自由符号变量. sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)') 〖目的〗 理解自由符号变量的确认规则。 〖解答〗 symvar(sym('sin(w*t)'),1) ans = w symvar(sym('a*exp(-X)'),1) ans = a symvar(sym('z*exp(j*th)'),1) ans = z 5求符号矩阵???? ??????=3332 31 232221 131211 a a a a a a a a a A 的行列式值和逆，所得结果应采用“子表达式置换”简洁化。 〖目的〗 理解subexpr 指令。 〖解答〗 A=sym('[a11 a12 a13;a21 a22 a23;a31 a32 a33]') DA=det(A) IA=inv(A); [IAs,d]=subexpr(IA,d) A = [ a11, a12, a13] [ a21, a22, a23] [ a31, a32, a33] DA = a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31 IAs = [ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 - a13*a22)] [ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d*(a11*a23 - a13*a21)] [ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 - a12*a21)] d = 1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31) 8（1）通过符号计算求t t y sin )(=的导数 dt dy 。（2）然后根据此结果，求- =0t dt dy 和2 π = t dt dy 。 〖目的〗 diff, limit 指令的应用。 如何理解运行结果。 〖解答〗 syms t ones表示1矩阵 zeros表示0矩阵 ones(4)表示4x4的1矩阵 zeros(4)表示4x4的0矩阵 zeros(4,5)表示4x5的矩阵 eye(10,10)表示10x10的单位矩阵 rand(4,5)表示4x5的伴随矩阵 det(a)表示计算a的行列式 inv(a)表示计算a的逆矩阵 Jordan(a)表示求a矩阵的约当标准块rank(a)表示求矩阵a的秩 [v,d]=eig(a)对角矩阵 b=a’表示求a矩阵的转置矩阵 sqrt表示求平方根 exp表示自然指数函数 log自然对数函数 abs绝对值 第一章 一、5（1） b=[97 67 34 10;-78 75 65 5;32 5 -23 -59]; >> c=[97 67;-78 75;32 5;0 -12]; >> d=[65 5;-23 -59;54 7]; >> e=b*c e = 5271 11574 -11336 664 1978 3112 （2）a=50:1:100 二、1 、x=-74; y=-27; z=(sin(x.^2+y.^2))/(sqrt(tan(abs(x+y)) )+pi) z = -0.0901 2、a=-3.0:0.1:3.0; >> b=exp(-0.3*a).*sin(a+0.3) 3、x=[2 4;-0.45 5]; y=log(x+sqrt(1+x.^2))/2 y = 0.7218 1.0474 -0.2180 1.1562 4、a*b表示a矩阵和b矩阵相乘 a.*b表示a矩阵和b矩阵单个元素相乘A(m,n)表示取a矩阵第m行，第n列 A(m,:)表示取a矩阵第m行的全部元素A(:,n)表示取a矩阵的第n列全部元素 A./B表示a矩阵除以b矩阵的对应元素， B.\A等价于A./B A.^B表示两个矩阵对应元素进行乘方运算 A.^2表示a中的每个元素的平方 A^2表示A*A 例：x=[1,2,3]; y=[4,5,6]; z=x.^y z= 1 3 2 729 指数可以是标量（如y=2）.底数也可以是标量（如x=2） 5、a=1+2i; >> b=3+4i; >> c=exp((pi*i)/6) c = 0.8660 + 0.5000i d=c+a*b/(a+b) d = 1.6353 + 1.8462i 第二章 二、4、（1） y=0;k=0; >> while y<3 k=k+1; y=y+1/(2*k-1); end >> display([k-1,y-1/(2*k-1)]) ans = %1、编制一个解数论问题的函数文件：取任意整数，若是偶数，则用2除，否则乘3加1，重复此过程，直到整数变为1。 function f=NO_1(X); Y(1)=X;k=1; while (X~=1) k=k+1; if (mod(X,2)==0) X=X/2; else X=3*X+1; end Y(k)=X; end plot(Y,'b.') end % 2、编制程序产生一个数组，满足：a1=1，a2=1，从第三个元素开始，每个元素等于前两个元素的和，直到数组的前后两个元素的比值比小于1e-4，并且以红色点线的形式画出这个数组。 clear; A(1)=1; A(2)=1; i=3; Z=1; while (abs(Z)>=1e-4) A(i)=A(i-1)+A(i-2); Z=A(i-1)/A(i); i=i+1; end plot (A,'r.') % 3、编写一个函数，能够产生分段函数。function y=test_3_1(X) if (X<=2) y=*X; elseif (X>6) y=; else y=调用分段函数，绘制曲线。 clear; i=1; for j=0::2; x(i)=j;y(i)=test_3_1(j)*test_3_1(j+2); i=i+1; end plot(x,y) % 4、在2pi周期内画正弦函数曲线，并加注坐标轴标识和标题，然后在3pi/4，pi，5pi/4处分别加入带箭头的说明性文本，最后加注图例。 clear; t = 0:pi/50:2*pi; n = length(t); y = sin(t); plot(t,y,'-bo','linewidth',1) xlabel('X');ylabel('Y'); title('正弦函数曲线'); text,, ' \leftarrow 3pi/4','FontSize',18); text,, ' \leftarrow 3pi/4','FontSize',18); text,,' \leftarrow 3pi/4','FontSize',18); hleg1 = legend('sin(x)'); % 5、A为任意一个n*m矩阵，写程序来计算A 中有多少个零元素，并输出个数。 A=input('输入一个矩阵 A = ') n=length(find(A==0)) % 6、A为任意一个向量，写程序找出A中的最小元素，并且输出这个最小元素。 A=input('输入一个向量 A = '); x=length(A); i=1; y=A(i); while (i 第1章 MATLAB概论 1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？ MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MATLAB系统由那些部分组成？ MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？ 在安装MATLAB时，安装内容由选择组件窗口中个复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M文件编辑/调试器？ 在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB中有几种获得帮助的途径？ 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、简述MA TLAB历史指令窗的主要作用。 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) 2. 用MATLAB 语句输入矩阵A 和B 3．假设已知矩阵A ，试给出相应的MATLAB 命令，将其全部偶数行提取出来， 赋给B 矩阵，用magic(8)A =命令生成A 矩阵，用上述命令检验一下结果是不是正确。 4．用数值方法可以求出∑=++++++==63 63622284212i i S Λ，试不采用循环的 形式求出和式的数值解。由于数值方法是采用double 形式进行计算的，难以保证有效位数字，所以结果不一定精确。试采用运算的方法求该和式的精确值。 5．选择合适的步距绘制出下面的图形。 （1）)/1sin(t ，其中)1,1(-∈t ； （2）)tan(sin )sin(tan t t -，其中),(ππ-∈t 6. 试绘制出二元函数2 2 2 2 )1(1)1(1),(y x y x y x f z +++ +-= =的三维图和三 视图 7. 试求出如下极限。 （1）x x x x 1)93(lim +∞ →； （2）1 1lim 0-+→→xy xy y x ； （3）2 2)()cos(1lim 2 2 220 0y x y x e y x y x +→→++- 8. 已知参数方程? ??-==t t t y t x sin cos cos ln ，试求出x y d d 和3 /2 2d d π=t x y 9. 假设?-=xy t t e y x f 0 d ),(2 ，试求2222 22y f y x f x f y x ??+???-?? 10. 试求出下面的极限。 （1）??????-++-+-+-∞→1)2(1 161141121lim 2222n n Λ； （2）)131211( lim 2 222π πππn n n n n n n ++++++++∞ →Λ 11. 试求出以下的曲线积分。 （1）?+l s y x d )(22，l 为曲线)sin (cos t t t a x +=，)cos (sin t t t a y -=， )20(π≤≤t 。 MATLAB试卷 一、选择题（每空2分，总共20分） 1.下列哪个变量的定义是不合法的（） (A) abcd-3 (B) xyz_3 (C ) abcdef (D) x3yz 2.下列哪条指令是求矩阵的行列式的值（） (A)inv (B) diag (C ) det (D) eig 3.在循环结构中跳出循环，执行循环后面代码的命令为（） (A) return (B) break (C) continue (D) keyboard 4.已知a=2:2:8, b=2:5，下面的运算表达式中，出错的 为（） (A) a'*b (B) a .*b (C) a*b (D) a-b 5.用round函数四舍五入对数组[2.48 6.39 3.93 8.52]取整，结果为（） (A) [2 6 3 8] (B) [2 6 4 8] (C) [2 6 4 9] (D) [3 7 4 9] 6. 下面的程序执行后array的值为( ) for k=1:10 if k>6 break; else array(k) = k; end end (A) array = [1, 2, 3, 4, 5, 6] (B) array = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10] (C) array =6 (D) array =10. 7.下列关于脚本文件和函数文件的描述中不正确的是( ) A. 去掉函数文件第一行的定义行可转变成脚本文件; B. 函数文件可以在命令窗口直接运行; C. 脚本文件可以调用函数文件; D. 函数文件中的第一行必须以function开始; 8.对应MATLAB提供的绘制平面曲线、离散序列、三维曲线、三维网格曲线、等高线图的函数为( ) A. plot, stem, mesh, plot3, surf B. plot, stem, plot3, mesh, contour C. plot, hist, mesh, plot3, contour D. plot, hist, plot3, meshc, contour 9.设X=[0:0.5*pi:pi], Y=cos(X)且Z=sin(X)，那么Y.*Z, Y*Z’和cat(1,Y,Z)的结果分别为( ) A. [0 0 0] 0 [1 0 -1; 0 1 0] B. [0 0 0] 0 [1 0 -1 0 1 0] C. [0 0 0] [0 1 0; 0 0 0; 0 -1 0] [1 0 -1; 0 1 0] D. [0 1 0; 0 0 0; 0 -1 0] [0 0 0] [1 0 -1 0 1 0] 10.使用下列哪一条指令可以将图形窗体分割成二行三列，并且将第二行第二列的绘图区域设置为当前的绘图区域（） A.subplot(3,2,2) B.subplot(2,3,2) C.subplot(3 ,2,5) D.subplot(2,3,5) MATLAB 基本运算 1．在MATLAB 中如何建立矩阵?? ?? ??194375，并将其赋予变量a ？>>a=[573;491] 2．在进行算术运算时，数组运算和矩阵运算各有什么要求？ 进行数组运算的两个数组必须有相同的尺寸。进行矩阵运算的两个矩阵必须满足矩阵运算规则，如矩阵a 与b 相乘（a*b ）时必须满足a 的列数等于b 的行数。 3．数组运算和矩阵运算的运算符有什么区别？ 在加、减运算时数组运算与矩阵运算的运算符相同，乘、除和乘方运算时，在矩阵运算的运算符前加一个点即为数组运算，如a*b 为矩阵乘，a.*b 为数组乘。 4．计算矩阵??????????897473535与??????????638976242之和。>>a=[535;374;798]; >>b=[242;679;836]; >>a+b ans = 7 779 1413151214 5．计算??????=572396a 与?? ????=864142b 的数组乘积。>>a=[693;275]; >>b=[241;468]; >>a.*b ans = 12 36384240 6．“左除”与“右除”有什么区别？ 在通常情况下，左除x=a\b 是a*x=b 的解，右除x=b/a 是x*a=b 的解，一般情况下，a\b ≠b/a 。 7．对于B AX =，如果??????????=753467294A ，???? ??????=282637B ，求解X 。>>A=[492;764;357]; >>B=[372628]’; >>X=A\B X = -0.5118 4.0427 1.3318 8．已知：???? ??????=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。 >>a=[123;456;789]; >>a.^2 ans = 1 4916 253649 6481 >>a^2 ans = 30 364266 81961021261509．[]7.0802.05--=a ，在进行逻辑运算时，a 相当于什么样的逻辑量。 相当于a=[11011]。 10．在sin(x )运算中，x 是角度还是弧度？ 在sin(x)运算中，x 是弧度，MATLAB 规定所有的三角函数运算都是按弧度进行运算。 Matlab 基础练习题 常量、变量、表达式 1、 MATLAB 中，下面哪些变量名是合法的？（ ） （A ）_num （B ）num_ （C ）num- （D ）-num 2、 在MA TLAB 中，要给出一个复数z 的模，应该使用（ ）函数。 （A ）mod(z) （B ）abs(z) （C ）double(z) （D ）angle(z) 3、 下面属于MATLAB 的预定义特殊变量的是？（ ） （A ）eps （B ）none （C ）zero （D ）exp 4、 判断：在MA TLAB 的内存工作区中，存放一个英文字符 'a' 需要占用1个字节，存放 一个中文字符‘啊’需要占用2个字节。（ 错，都是2个字节 ） 5、 判断：MA TLAB 中，i 和j （ 对 ） 6、 判断：MA TLAB 中，pi 代表圆周率，它等于3.14。（ 错，后面还有很多位小数 ） 7、 在MA TLAB 中，若想计算的5 1)3.0sin(21+= πy 值，那么应该在MA TLAB 的指令窗中 输入的MA TLAB 指令是__y1=2*sin(0.3*pi)/(1+sqrt(5))_。 8、 在MA TLAB 中，a = 1，b = i ，则a 占_8__个字节，b 占_16_个字节，c 占________字 节。 9、 在MA TLAB 中，inf 的含义是__无穷大__，nan 的含义是__非数（结果不定）___。 数组 1、 在MA TLAB 中，X 是一个一维数值数组，现在要把数组X 中的所有元素按原来次序 的逆序排列输出，应该使用下面的（ ）指令。 （A ）X[end:1] （B ）X[end:-1:1] （C ）X (end:-1:1) （D ）X(end:1) 2、 在MA TLAB 中，A 是一个字二维数组，要获取A 的行数和列数，应该使用的MATLAB 的命令是（ ）。 （A ）class(A) （B ）sizeof(A) （C ）size(A) （D ）isa(A) 3、 在MATLAB 中，用指令x=1:9生成数组x 。现在要把x 数组的第二和第七个元素都 赋值为0，应该在指令窗中输入（ ） （A ）x([2 7])=(0 0) （B ）x([2,7])=[0,0] （C ）x[(2,7)]=[0 0] （D ）x[(2 7)]=(0 0) 4、 在MA TLAB 中，依次执行以下指令：clear;A=ones(3,4); A(:)=[-6:5];这时， 若在指令窗中输入指令b=A(:,2)'，那么，MATLAB 输出的结果应该是（ ） （A ）b = -3 -2 -1 （B ）b = -2 -1 0 1 （C ）b = -5 -1 3 （D ）b = -5 -2 1 4 5、 在MA TLAB 中，A = 1:9，现在执行如下指令L1 = ~(A>5)，则MATLAB 的执行结果应 该是L1 =___ 1 1 1 1 1 0 0 0 0___。 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.如何理解“矩阵是MATLAB最基本的数据对象”？ 答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。 因此，矩阵是MATLAB最基本、最重要的数据对象。 2.设A和B是两个同维同大小的矩阵，问： (1)A*B和A.*B的值是否相等？ 答：不相等。 (2)A./B和B.\A的值是否相等？ 答：相等。 (3)A/B和B\A的值是否相等？ 答：不相等。 (4)A/B和B\A所代表的数学含义是什么？ 答：A/B等效于B的逆右乘A矩阵，即A*inv(B)，而B\A等效于B矩阵的逆左乘A矩阵，即inv(B)*A。 3.写出完成下列操作的命令。 (1)将矩阵A第2~5行中第1, 3, 5列元素赋给矩阵B。 答：B=A(2:5,1:2:5); 或B=A(2:5,[1 3 5]) (2)删除矩阵A的第7号元素。 答：A(7)=[] (3)将矩阵A的每个元素值加30。 答：A=A+30; (4)求矩阵A的大小和维数。 答：size(A); ndims(A); (5)将向量t的0元素用机器零来代替。 答：t(find(t==0))=eps; (6)将含有12个元素的向量x转换成34 矩阵。 答：reshape(x,3,4); (7)求一个字符串的ASCII码。 答：abs(‘123’); 或double(‘123’); (8)求一个ASCII码所对应的字符。 答：char(49); 4. 下列命令执行后，L1、L2、L3、L4的值分别是多少？ A=1:9;B=10-A;... L1=A==B; L2=A<=5; L3=A>3&A<7; L4=find(A>3&A<7); 答：L1的值为[0, 0, 0, 0, 1, 0, 0, 0, 0] L2的值为[1, 1, 1, 1, 1, 0, 0, 0, 0] L3的值为[0, 0, 0, 1, 1, 1, 0, 0, 0] L4的值为[4, 5, 6] 5. 已知 23100.7780414565532503269.5454 3.14A -????-??=????-?? 完成下列操作： (1) 取出A 的前3行构成矩阵B ，前两列构成矩阵C ，右下角32?子矩阵构 成矩阵D ，B 与C 的乘积构成矩阵E 。 答：B=A(1:3,:); C=A(:,1:2); D=A(2:4,3:4); E=B*C; (2) 分别求E 第4章数值运算 习题 4 及解答 1 根据题给的模拟实际测量数据的一组t和)(t y试用数值差分diff 或数值梯度gradient指令计算)(t y'曲线绘制 y'，然后把)(t y和)(t 在同一张图上，观察数值求导的后果。（模拟数据从prob_data401.mat获得） 〖目的〗 ●强调：要非常慎用数值导数计算。 ●练习mat数据文件中数据的获取。 ●实验数据求导的后果 ●把两条曲线绘制在同一图上的一种方法。 〖解答〗 （1）从数据文件获得数据的指令 假如prob_data401.mat文件在当前目录或搜索路径上 clear load prob_data401.mat （2）用diff求导的指令 dt=t(2)-t(1); yc=diff(y)/dt; %注意yc的长度将比y短1 plot(t,y,'b',t(2:end),yc,'r') （3）用gradent求导的指令（图形与上相似） dt=t(2)-t(1); yc=gradient(y)/dt; plot(t,y,'b',t,yc,'r') grid on 〖说明〗 ● 不到万不得已，不要进行数值求导。 ● 假若一定要计算数值导数，自变量增量dt 要取得比原有数据相对误差高1、2个量级 以上。 ● 求导会使数据中原有的噪声放大。 2 采用数值计算方法，画出dt t t x y x ? =0 sin )(在]10 ,0[区间曲线，并计算)5.4(y 。 〖提示〗 ● 指定区间内的积分函数可用cumtrapz 指令给出。 ● )5.4(y 在计算要求不太高的地方可用find 指令算得。 〖目的〗 ● 指定区间内的积分函数的数值计算法和cumtrapz 指令。 ● find 指令的应用。 〖解答〗 dt=1e-4; t=0:dt:10; t=t+(t==0)*eps; f=sin(t)./t; s=cumtrapz(f)*dt; plot(t,s,'LineWidth',3) ii=find(t==4.5); s45=s(ii) s45 = MATLAB期末上机试题带答案 MATLAB 期末上机考试试题带答案版 姓名： 学号： 成绩： 1.请实现下图： x 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.请实现下图： matlab练习题和答案 控制系统仿真实验 Matlab 部分实验结果 目录 实验一MATLAB基本操作..................................................................... .......................1 实验二Matlab 编程..................................................................... ...............................5 实验三Matlab 底层图形控制..................................................................... .....................6 实验四控制系统古典分析...................................................................... .......................12 实验五控制系统现代分析..................................................................... . (15) 实验六PID 控制器的设计..................................................................... ......................19 实验七系统状态空间设计...................................................................... .......................23 实验九直流双闭环调速系统仿真..................................................................... . (25) 实验一MATLAB基本操作 学习中心/函授站_ 姓名学号 西安电子科技大学网络与继续教育学院 2015学年上学期 《MATLAB与系统仿真》期末考试试题 （综合大作业） 考试说明： 1、大作业于2015年4月3日公布，2015年5月9日前在线提交； 2、考试必须独立完成，如发现抄袭、雷同、拷贝均按零分计。 3、程序设计题（三（8，10））要求写出完整的程序代码，并在matlab软件环境调试并运行通过，连同运行结果一并附上。 一、填空题（1’ ×25=25’） 1、Matlab的全称为MATrix LABoratory 。 2、在Matlab编辑器中运行程序的快捷键是：F5 。 3、Matlab的工作界面主要由以下五个部分组成，它们分别是：Command window(指令窗口、Current directory（当前路径浏览器）、Workspace(工作台窗口) 、Command history （指令历史窗口）和Launch pad（发布平台）。 4、在Matlab中inf表示：无穷大；clc表示：清除当前Command window区域内容；more表示：显示命令窗口输出一个整个萤幕的时间；who表示：返回当前工作空间中所有变量名；whos表示：返回当前工作空间中所有变量，以及它们的名字、尺寸、所占字节数、属性等信息。 5、在Matlab命令窗口中运行命令Simulink 可以打开Simulink模块库浏览器窗口。 6、求矩阵行列式的函数：det ；求矩阵特征值和特征向量的函数eig 。 7、Matlab预定义变量ans表示：当前未定义变量和变量名的答案；eps表示： 浮点相对精度（其值是从1.0到下一个最大浮点数之间的差值）；nargin表示：函 9、建立结构数组或转换结构数组的函数为：structMATLAB试题解答
matlab试卷合集附答案
matlab练习题及答案
matlab课后习题解答第二章
matlab课后答案解析完整版
matlab考试题及答案
matlab课后习题及答案详解
MATLAB期末考试试卷及其参考答案
matlab习题及答案
MATLAB试卷
matlab习题及答案2
matlab基础练习题(带答案)
MATLAB期末上机试题带答案
matlab课后习题与答案
matlab课后习题答案
MATLAB期末上机试题带答案
matlab练习题和答案
MATLAB试卷+答案
- matlab课后习题答案1到6章
- MATLAB课后习题集附标准答案
- matlab课后习题及答案详解
- matlab课后习题及答案详解
- matlab课后习题答案第四章
- matlab课后习题及答案详解
- MATLAB习题及参考答案
- matlab习题及答案
- matlab课后习题与答案
- matlab课后习题答案
- MATLAB数据分析课后习题答案
- matlab课后习题及答案
- matlab课后习题解答第二章
- MATLAB基础教程-薛山第二版-课后习题答案
- MATLAB程序设计与应用课后习题答案
- Matlab习题与答案2
- matlab程序设计第三章课后习题答案
- matlab课后习题答案-(附图)
- matlab基础练习题(带答案)
- 《MATLAB》课后习题答案