《Matlab 软件实训》
一、 训练目的
1、 熟悉Matlab 工具软件的开发环境
2、 掌握Matlab 的一些常用命令
3、 掌握Matlab/Simulink 建模仿真的方法
二、 训练内容
1、 利用命令进行各种相关数学计算；
2、 编写M 文件实现各种相关数学计算；
3、 利用Simulink 建模并仿真；
三、 训练任务
1、 解方程：Ax=B ，其中，
A=?????????
???------5407395012818053，B=????
?
?
??????-6120 本题解决思路：
这是一个简单的数组运算，核心点在于数组的输入方法和一些基本的数组运算基础，本题中需要注意的是左除与右除的区别。
命令如下
>> A=[-3 5 0 8;1 -8 2 -1;0 -5 9 3;-7 0 -4 5] A =
-3 5 0 8 1 -8 2 -1 0 -5 9 3 -7 0 -4 5
>> B=[0;2;-1;6] B =
0 2 -1 6
>> X=A\B
X =
-0.6386
-0.4210
-0.3529
0.0237
2、数值运算
已知某班的5名学生的三门成绩如下：
序号 1 2 3 4 5
高数78 89 64 73 68
外语83 77 80 78 70
电路82 91 78 82 68
试写出有关命令，先分别找出三门课的最高分及其学生序号，然后找出三门课总分的最高分及其学生序号。
本题解决思路：
详细的分析题目，可以发现本题就是一个多组求最大值并输出的问题。有三个核心点，第一点是求最大值问题。第二点是如何将多组最大值问题合并在一个程序中解出，第三点是输入和输出的方法。
最大值我采用了MA TLAB中的max命令，其实也可以采用自己编写M文件的方法，但是程序越精简运行越快，我选择了自带命令。第二个问题我采用了一个for和if配合的循环判断的语句。第三个问题我采用disp输出，并用num2str将字符串转换输出。
程序如下
a=input('高数成绩');
b=input('英语成绩');
c=input('电工成绩');
for i=1:3
if i==1
[y1,t1]=max(a,[],2);
else if i==2
[y2,t2]=max(b,[],2);
else
[y3,t3]=max(c,[],2);
end
end
end
s=a+b+c;
[y4,t4]=max(s,[],2);
disp(['高数最高分是',num2str(t1),'号',num2str(y1),'分'])
disp(['英语最高分是',num2str(t2),'号',num2str(y2),'分'])
disp(['电工最高分是',num2str(t3),'号',num2str(y3),'分'])
disp(['总分最高分是',num2str(t4),'号',num2str(y4),'分'])
高数成绩[78 89 64 73 68]
英语成绩[83 77 80 78 70]
电工成绩[82 91 78 82 68] 高数最高分是2号89分 英语最高分是1号83分 电工最高分是2号91分 总分最高分是2号257分 >>
3、 建立M 文件，计算积分y=
?
20
sin π
xdx ；
M 文件如下
function [ y ] = J( a,b ) syms x
y=int(sin(x),x,a,b) end
在命令行中输入J(0,pi/2) 运算得出结果 >> J(0,pi/2) y = 1
ans = 1
4、 建立M 文件，在[0 2π]范围内绘制二维曲线图y=)5cos(*)sin(x x
本题解决思路：
在M 文件中要有输入变量和输出量，分析本题得知，本题的x 并不是输入文件，他是一个中间变量，真正的输入变量其实是x 的输入范围，定义变量为(a,b),确定这个以后直接调用MATLAB 中的内置命令求解。
M 文件如下
function [ y ] = H( a,b ) x=linspace(a,b); y=sin(x).*cos(5*x); plot(x,y); end
编辑好后输入x 的范围并调运M 文件“H ”输入H （0,2*pi ） >> H(0,2*pi)
ans =
Columns 1 through 9
0 0.0603 0.1019 0.1098 0.0746 -0.0050 -0.1216 -0.2603 -0.4005
Columns 10 through 18
-0.5187 -0.5926 -0.6040 -0.5424 -0.4069 -0.2068 0.0388 0.3033 0.5558
Columns 19 through 27
0.7652 0.9042 0.9530 0.9022 0.7544 0.5240 0.2355 -0.0792 -0.3851
Columns 28 through 36
-0.6482 -0.8398 -0.9402 -0.9407 -0.8450 -0.6678 -0.4330 -0.1705 0.0882
Columns 37 through 45
0.3139 0.4836 0.5827 0.6070 0.5625 0.4638 0.3316 0.1893 0.0594
Columns 46 through 54
-0.0401 -0.0978 -0.1108 -0.0845 -0.0313 0.0313 0.0845 0.1108 0.0978
Columns 55 through 63
0.0401 -0.0594 -0.1893 -0.3316 -0.4638 -0.5625 -0.6070 -0.5827 -0.4836
Columns 64 through 72
-0.3139 -0.0882 0.1705 0.4330 0.6678 0.8450 0.9407 0.9402 0.8398
Columns 73 through 81
0.6482 0.3851 0.0792 -0.2355 -0.5240 -0.7544 -0.9022 -0.9530 -0.9042
Columns 82 through 90
-0.7652 -0.5558 -0.3033 -0.0388 0.2068 0.4069 0.5424 0.6040 0.5926
Columns 91 through 99
0.5187 0.4005 0.2603 0.1216 0.0050 -0.0746 -0.1098 -0.1019 -0.0603
Column 100
-0.0000
5、利用Simulink建立仿真模型，并给出建模步骤。
本题思路:
确定输入函数的类型，确定中间运算过程，确定最终结果显示器
结果截图如下：
四、个人总结
MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，MATLAB的最突出的特点就是简洁。
MATLAB相对于其他的一些编程软件有许多的优点：
一、语言简洁紧凑，使用方便灵活，库函数极其丰富。
二、运算符丰富。
三、MATLAB既具有结构化的控制语句（如for循环、while循环），又有面向对
象编程的特性。
四、语法限制不严格，程序设计自由度大。
五、程序的可移植性很好，基本上不做修改就可以在各种型号的计算机和操作
系统上运行。
六、MATLAB的图形功能强大。
用madlab创建矩阵时，方法有两种：第一、可以直接依次输入矩阵各行各列的元素，但矩阵元素必须用[ ]括住，矩阵元素必须用逗号或空格分隔，在[ ]内矩阵的行与行之间必须用分号分隔。第二、用MATLAB函数创建矩阵。MATLAB可以进行矩阵的加减、乘除的元素，求可逆矩阵、转置矩阵，求矩阵的特征值，求线性方程组等等。 MATLAB 的功能是非常强大的，MATLAB不仅有强大的运算功能，它还有强大的绘图功能，我对它的了解也仅仅就是一点点，或许说还没有入门。比如说它含有丰富的内建函数，例如数学函数中的三角函数、复函数、多项式函数、数据分析函数的求平均值、最大最小值、排序等，以及逻辑/选择函数如if－else等，还有用来模拟随机发生事件的随机函数。
这些我都不了解。
虽学习MATLAB的时间虽然很短，但却让我了解到了它的强大和它的功能。我想就算时间足够，老师也不能把所有的都讲解给我们，因为一个软件的功能需要我们自己不断的去摸索，老师也不可能知道所有。老师只是个指路的明灯，最终的学习还是要靠自己。
而且在摸索过程中，我们能够发现和体会学习的快乐！
五、学习到的技能
基本的命令:
1、!dir 可以查看当前工作目录的文件。 !dir& 可以在dos状态下查看。
2、who 可以查看当前工作空间变量名， whos 可以查看变量名细节。
3、功能键：
功能键快捷键说明
方向上键 Ctrl+P 返回前一行输入
方向下键 Ctrl+N 返回下一行输入
方向左键 Ctrl+B 光标向后移一个字符
方向右键 Ctrl+F 光标向前移一个字符
Ctrl+方向右键 Ctrl+R 光标向右移一个字符
Ctrl+方向左键 Ctrl+L 光标向左移一个字符
home Ctrl+A 光标移到行首
End Ctrl+E 光标移到行尾
Esc Ctrl+U 清除一行
Del Ctrl+D 清除光标所在的字符
Backspace Ctrl+H 删除光标前一个字符
Ctrl+K 删除到行尾
Ctrl+C 中断正在执行的命令
4、clc可以命令窗口显示的内容，但并不清除工作空间。
二、函数及运算
1、运算符：
＋：加，－：减， *：乘， /：除， \：左除 ^：幂，‘：复数的共轭转置，（）：制定运算顺序。
2、常用函数表：
sin( ) 正弦（变量为弧度）
Cot( ) 余切（变量为弧度）
sind( ) 正弦（变量为度数）
Cotd( ) 余切（变量为度数）
asin( ) 反正弦（返回弧度）
acot( ) 反余切（返回弧度）
Asind( ) 反正弦（返回度数）
acotd( ) 反余切（返回度数）
cos( ) 余弦（变量为弧度）
exp( ) 指数
cosd( ) 余弦（变量为度数）
log( ) 对数
acos( ) 余正弦（返回弧度）
log10( ) 以10为底对数
acosd( ) 余正弦（返回度数）
sqrt( ) 开方
tan( ) 正切（变量为弧度）
realsqrt( ) 返回非负根
tand( ) 正切（变量为度数）
abs( ) 取绝对值
atan( ) 反正切（返回弧度）
angle( ) 返回复数的相位角
atand( ) 反正切（返回度数）
mod(x,y) 返回x/y的余数
sum( ) 向量元素求和
3、其余函数可以用help elfun和help specfun命令获得。
4、常用常数的值：
pi 3.1415926…….
realmin 最小浮点数，2^-1022
i 虚数单位
realmax 最大浮点数，（2－eps）2^1022
j 虚数单位
Inf 无限值
eps 浮点相对经度＝2^-52
NaN 空值
三、数组和矩阵：
1、构造数组的方法：增量发和linspace(first,last,num)first和last为起始和终止数，num为需要的数组元素个数。
2、构造矩阵的方法：可以直接用[ ]来输入数组，也可以用以下提供的函数来生成矩阵。ones( ) 创建一个所有元素都为1的矩阵，其中可以制定维数，1，2….个变量zeros() 创建一个所有元素都为0的矩阵
eye() 创建对角元素为1，其他元素为0的矩阵
diag() 根据向量创建对角矩阵，即以向量的元素为对角元素
magic() 创建魔方矩阵
rand() 创建随机矩阵，服从均匀分布
randn() 创建随机矩阵，服从正态分布
randperm() 创建随机行向量
horcat C=[A,B]，水平聚合矩阵，还可以用cat(1,A,B)
vercat C=[A;B]，垂直聚合矩阵, 还可以用cat(2,A,B)
repmat(M,v,h) 将矩阵M在垂直方向上聚合v次，在水平方向上聚合h次blkdiag（A，B）以A，和B为块创建块对角矩阵
length 返回矩阵最长维的的长度
ndims 返回维数
numel 返回矩阵元素个数
size 返回每一维的长度，[rows,cols]=size(A)
reshape 重塑矩阵，reshape(A,2,6),将A变为2×6的矩阵，按列排列。rot90 旋转矩阵90度，逆时针方向
fliplr 沿垂轴翻转矩阵
flipud 沿水平轴翻转矩阵
transpose 沿主对角线翻转矩阵
ctranspose 转置矩阵，也可用A’或A.’，这仅当矩阵为复数矩阵时才有区别
inv 矩阵的逆
det 矩阵的行列式值
trace 矩阵对角元素的和
norm 矩阵或矢量的范数，norm（a，1），norm（a，Inf）…….
normest 估计矩阵的最大范数矢量
chol 矩阵的cholesky分解
cholinc 不完全cholesky分解
lu LU分解
luinc 不完全LU分解
qr 正交分解
kron（A，B） A为m×n，B为p×q，则生成mp×nq的矩阵，A的每一个元素都会乘上B，并占据p×q大小的空间
rank 求出矩阵的刺
pinv 求伪逆矩阵
A^p 对A进行操作
A.^P 对A中的每一个元素进行操作
四、数值计算
1、线性方程组求解
（1）AX=B的解可以用X＝A\B求。XA=B的解可以用X= A/B求。如果A是m×n的矩阵，当m＝n时可以找到唯一解，m
（2）AX=b, A＝L×U，[L,U]=lu(A), X=U\(L\b),即用LU分解求解。
（3）QR（正交）分解是将一矩阵表示为一正交矩阵和一上三角矩阵之积，A＝Q×R[Q,R]=chol(A), X=Q\(U\b)
（4）cholesky分解类似。
2、特征值
D＝eig（A）返回A的所有特征值组成的矩阵。[V,D]=eig(A),还返回特征向量矩阵。
3、A＝U×S×UT，[U,S]=schur(A).其中S的对角线元素为A的特征值。
4、多项式Matlab里面的多项式是以向量来表示的，其具体操作函数如下：
conv 多项式的乘法
deconv 多项式的除法，【a，b】＝deconv（s），返回商和余数
poly 求多项式的系数（由已知根求多项式的系数）
polyeig 求多项式的特征值
Polyfit（x，y，n）多项式的曲线拟合，x，y为被拟合的向量，n为拟合多项式阶数。
polyder 求多项式的一阶导数，polyder（a，b）返回ab的导数
[a,b]＝polyder（a，b）返回a/b的导数。
polyint 多项式的积分
polyval 求多项式的值
polyvalm 以矩阵为变量求多项式的值
residue 部分分式展开式
roots 求多项式的根（返回所有根组成的向量）
注：用ploy（A）求出矩阵的特征多项式，然后再求其根，即为矩阵的特征值。
5、插值常用的插值函数如下：
griddata 数据网格化合曲面拟合
Griddata3 三维数据网格化合超曲面拟合
interp1 一维插值
（yi=interp1(x,y,xi,’method’)Method=nearest/linear/spline/pchip/cubic Interp2 二维插值zi=interp1(x,y,z,xi,yi’method’),bilinear
Interp3 三维插值
interpft 用快速傅立叶变换进行一维插值，help fft。
mkpp 使用分段多项式
spline 三次样条插值
pchip 分段hermit插值
6、函数最值的求解
fminbnd（‘f’，x1，x2，optiset（,））求f在 x1和x2之间的最小值。Optiset 选项可以有‘Display’+‘iter’/’off’/’final’,分别表示显示计算过程/不显示/只显示最后结果。fminsearch求多元函数的最小值。fzero（‘f’，x1）求一元函数的零点。X1为起始点。同样可以用上面的选项。
五、图像绘制：
1、基本绘图函数
plot 绘制二维线性图形和两个坐标轴
plot3 绘制三维线性图形和两个坐标轴
fplot 在制定区间绘制某函数的图像。fplot（‘f’，区域，线型，颜色）loglog 绘制对数图形及两个坐标轴（两个坐标都为对数坐标）semilogx 绘制半对数坐标图形
semilogy 绘制半对数坐标图形
2、可以用subplot（3，3，1）表示将绘图区域分为三行三列，目前使用第一区域。
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、标点符号; %—用来表示该行为注释行。 可以使命令行不显示运算结果， 2、x 为0 ～4pi ，步长为 的向量，使用命令 x=0:*pi:4*pi 创建。 3、输入矩阵A=[1323?57569 ] ，使用全下标方式用A(2,2)取出元素“-5 ”，使用单下标方式用A(5)取出元素“-5 ”。 4、符号表达式sin(2*a+t)+m 中独立的符号变量为 t 。 5、M 脚本文件和M 函数文件的主要区别是M 脚本文件没有函数定义和M 函数文件有函数定义_______。 6. 设x 是一维数组，x 的倒数第3个元素表示为x(_end-2_) 设y 为二维数组，要删除y 的第34行和48列，可使用命令y(34,:)=[] ； y(:,48)=[]； 7. 将变量x 以Ascii 文本格式存储到文件，应使用命令 save _x ； 8. 在while 表达式, 语句体, End 循环语句中，表达式的值 非零 时表示循环条件为真，语句体将被执行，否则跳出该循环语句； 9.要从键盘读入一个字符串并赋值给变量x ，且给出提示“Who is she ”，应使用命令x=input(‘Who is she’,’s’) ； 10．设A=(a ij )m?n 和B=(b ij )m?n 和C=(c ij )m?n 均为m*n 矩阵，且存在于WorkSpace 中，要产生矩阵D=(a ij ?c ij b ij c ij ?)m?n ，可用命令D=(A-C)/B.^C ， 计算|(A ′B )?1| 可用命令det(inv(A’*B) 11. 在MATLAB 命令窗口中的“>>”标志为MATLAB 的 命令行 提示符，“│”标志为 输入 提示符。 12.已知A=[1 2 3；4 5 0；7 8 9]；B=[1 0 3；1 5 0；0 1 2]；写出下列各指令运行的结果。 A+B ； A.*B ； A==B ； ans= [2,2,6;5,10,0;7,9,11] ans= [1,0,9;4,25,0;0,8,18] ans= [1,0,1;0,1,1;0,0,0] 13.已知A 是矩阵，求A 的对角矩阵函数是 diag ； 求A 的下三角矩阵函数是 tril 。 的程序文件和Simulink 模型文件的扩展名分别是.m 、.mdl 最基本的绘图函数为plot() 16. A=[1,2,3;4,5,6]; A(:,[1,3])=[];A= [2；5] 17. fix= -1 , round= -2 . 18.已知A=[1 3;4 6];C=[A,eye(sizeA);A ,ones(sizeA)] 则Ｃ＝ [1 3 1 0；4 6 0 1；1 3 1 1；4 6 1 1]
一、填空题 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的值分别为和。
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___。
第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 函数名”可以得到
第一套 华东交通大学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、编制一个解数论问题的函数文件：取任意整数，若是偶数，则用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、标点符号; _______可以使命令行不显示运算结果， %——用来表示该行为注释行。 2、x 为0 ～4pi ，步长为0.1pi 的向量，使用命令_______创建。 x=0:0.1*pi:4*pi 3、输入矩阵A= 错误！未找到引用源。 ，使用全下标方式用A(2,2) 取出元素“-5 ”，使用单下标方式用_______取出元素“-5 ”。 A(5) 4、符号表达式sin(2*a+t)+m 中独立的符号变量为_______。 t 5、M 脚本文件和M 函数文件的主要区别是M 脚本文件没有函数定义和M 函数文件有函数定义_______。 6. 设x 是一维数组，x 的倒数第3个元素表示为_______； 设y 为二维数组，要删除y 的第34行和48列，可使用命令_______; _______； x(_end-2_) y(34,:)=[] y(:,48)=[] 7. 将变量x 以Ascii 文本格式存储到文件fname.txt ，应使用命令_________ _； save _x 8. 在while 表达式, 语句体, End 循环语句中，表达式的值__ __时表示循环条件为真，语句体将被执行，否则跳出该循环语句； 非零 9.要从键盘读入一个字符串并赋值给变量x ，且给出提示“Who is she?”，应使用命令_________； x=input(‘Who is she?’,’s’)_ 10．设A=错误！未找到引用源。和B=错误！未找到引用源。 和C=错误！未找到引用源。均为m*n 矩阵，且存在于WorkSpace 中，要产生矩阵 D=错误！未找到引用源。 ，可用命令________ _， 计算错误！未找到引用源。 可用命令________; D=(A-C)/B.^C det(inv(A’*B ) 11. 在MATLAB 命令窗口中的“>>”标志为MATLAB 的_______ 提示符，“│”标志为_______提示符。 命令行 输入 12.已知A=[1 2 3；4 5 0；7 8 9]；B=[1 0 3；1 5 0；0 1 2]；写出下列各指令运行的结果。 A+B ； A.*B ； A==B ； ans= [2,2,6;5,10,0;7,9,11] ans= [1,0,9;4,25,0;0,8,18] ans= [1,0,1;0,1,1;0,0,0] 13.已知A 是矩阵，求A 的对角矩阵函数是_______， 求A 的下三角矩阵函数是_______ 。 diag tril 14.MATLAB 的程序文件和Simulink 模型文件的扩展名分别是_______、 。 .m .mdl 15.MATLAB 最基本的绘图函数为_______ 。 plot() 16. A=[1,2,3;4,5,6]; A(:,[1,3])=[];A=__________________ [2；5] 17. fix(-1.5)=___ ________, round(-1.5)=__ _______________. -1 -2
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 。
一、选择题共10小题 二、填空题共10空 三、看程序填空 四、编程题 3题 填空题（计算部分） 1、标点符号; _______可以使命令行不显示运算结果， %——用来表示该行为注释行。 2、x为0～4pi，步长为0.1pi的向量，使用命令_______创建。 x=0:0.1*pi:4*pi 3、输入矩阵A=，使用全下标方式用A(2,2) 取出元素“-5”，使用单下标方式用_______取出元素“-5”。 A(5) 4、符号表达式sin(2*a+t)+m中独立的符号变量为_______。 t 5、M脚本文件和M函数文件的主要区别是M脚本文件没有函数定义和M函数文件有函数定义_______。 6. 设x是一维数组，x的倒数第3个元素表示为_______； 设y为二维数组，要删除y的第34行和48列，可使用命令_______; _______； x(_end-2_) y(34,:)=[] y(:,48)=[] 7. 将变量x以Ascii文本格式存储到文件fname.txt，应使用命令_________ _； save _x 8. 在while 表达式, 语句体, End 循环语句中，表达式的值__ __时表示循环条件为真，语句体将被执行，否则跳出该循环语句； 非零 9.要从键盘读入一个字符串并赋值给变量x，且给出提示“Who is she?”，应使用命令_________； x=input(‘Who is she?’,’s’)_ 10．设A=和B=和C=均为m*n矩阵，且存在于WorkSpace中，要产生矩阵D=，可用命令________ _， 计算可用命令________; D=(A-C)/B.^C
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)
第1章 MATLAB 基础 MATLAB 操作桌面有几个窗口如何使某个窗口脱离桌面成为独立窗口又如何将脱离出去的窗口重新放置到桌面上 答：在MATLAB 操作桌面上有五个窗口。在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口成为独立窗口的Undock 按钮，点击Undock 按钮就可以使该窗口脱离桌面成为独立窗口。在独立窗口的Desktop 菜单中选择Dock...项就可以将独立的窗口重新放置到桌面上。 如何启动M 文件编辑/调试器 答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M 文件编辑/调试器将被启动。在命令窗口中键入edit 命令时也可以启动M 文件编辑/调试器。 存储在工作空间中的数组能编辑吗如何操作 答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 命令历史窗口除了可以观察前面键入的命令外，还有什么用途 答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M 文件中。 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别 答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path 菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB 运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 在MATLAB 中有几种获得帮助的途径 答：（1）帮助浏览器：选择view 菜单中的Help 菜单项或选择Help 菜单中的MATLAB Help 菜单项可以打开帮助浏览器。 （2）help 命令：在命令窗口键入“help ” 命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息。 （3）lookfor 命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数。 （4）模糊查询：输入命令的前几个字母，然后按Tab 键，就可以列出所有以这几个字母开始的命令和函数。 注意：lookfor 和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help 命令显示详细信息。 第2章 MATLAB 矩阵运算基础 在MATLAB 中如何建立矩阵??? ???194375，并将其赋予变量a 答：在Command Window 窗口输入操作: >> a=[5 7 3;4 9 1] 有几种建立矩阵的方法各有什么优点 答：（1）直接输入法，如a=[1 2 3 4]，优点是输入方法方便简捷； （2）通过M 文件建立矩阵，该方法适用于建立尺寸较大的矩阵，并且易于修改； （3）由函数建立，如y=sin(x)，可以由MATLAB 的内部函数建立一些特殊矩阵； （4）通过数据文件建立，该方法可以调用由其他软件产生数据。 在进行算术运算时，数组运算和矩阵运算各有什么要求 答：进行数组运算的两个数组必须有相同的尺寸。进行矩阵运算的两个矩阵必须满足矩阵运算规则，如矩阵a 与b 相乘（a*b ）时必须满足a 的列数等于b 的行数。 数组运算和矩阵运算的运算符有什么区别 答：在加、减运算时数组运算与矩阵运算的运算符相同，乘、除和乘方运算时，在矩阵运算的运算符前加一个点即为数组运算，如a*b 为矩阵乘，a.*b 为数组乘。 计算矩阵??????????897473535与??? ?? ?????638976242之和。 答： 求??? ???+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x 的共轭转置。 答： 计算??????=572396a 与?? ? ???=864142b 的数组乘积。 答：
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.请实现下图： 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);
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 考试试题(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]写到文件里，写完后文件看起来是这样的 1 2 3 4 5 6 7 8 9 5.试从Yahoo网站上获得微软公司股票的2008年9月的每日收盘价。 6.编写M文件，从Yahoo网站批量读取至在2008年9月份的每日收盘价（提示：使用字符串函数）。 7. 将金牛股份（000937）2005年12月14日至2006年1月10日的交易记录保存到Excel中，编写程序将数据读入MATLAB中，进一步将数据读入Access数据库文件。 8.已知资产每日回报率为，标准差为，资产现在价值为亿，求5％水平下资产的10天在险价值(Var)。 =[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、求下列联立方程的解