中南大学matlab习题复习课
matlab期末复习题及答案
matlab期末复习题及答案1. MATLAB基础操作在MATLAB中，如何创建一个名为"myMatrix"的3x3矩阵，其元素分别为1, 2, 3, 4, 5, 6, 7, 8, 9？答案：在MATLAB中，可以通过直接输入矩阵元素来创建矩阵。
例如，要创建一个名为"myMatrix"的3x3矩阵，可以使用以下命令：```matlabmyMatrix = [1 2 3; 4 5 6; 7 8 9];```这样，MATLAB就会创建一个3x3的矩阵，其元素按照行顺序排列。
2. 矩阵运算给定两个矩阵A和B，其中A = [1 2; 3 4]，B = [5 6; 7 8]，计算矩阵A和B的和。
答案：在MATLAB中，可以使用加号(+)来计算两个矩阵的和。
对于给定的矩阵A和B，可以使用以下命令来计算它们的和：```matlabA = [1 2; 3 4];B = [5 6; 7 8];C = A + B;```执行上述命令后，矩阵C的结果将是：```C = [6 8;10 12];```3. 条件语句编写一个MATLAB脚本，判断一个给定的数n是否为质数，并输出相应的信息。
答案：在MATLAB中，可以使用if-else语句来判断一个数是否为质数。
以下是一个简单的脚本示例：```matlabn = input('请输入一个数：');if mod(n, 2) == 0 && n > 2disp('该数不是质数');elseif n == 1disp('1不是质数');elseisPrime = true;for i = 2:sqrt(n)if mod(n, i) == 0isPrime = false;break;endendif isPrimedisp('该数是质数');elsedisp('该数不是质数');endend```该脚本首先接收用户输入的数n，然后通过一个for循环检查n是否有除了1和它自身以外的因数，从而判断n是否为质数。
中南大学材料院matlab操作题集答案
的491操作题集操作题1 《Matlab与科学计算》P28 例2.1要求计算水在温度为0℃，20℃，40℃，60℃，80℃时的黏度，已知水的黏度随/(1+at+bt2)温度的变化公式为μ=μ为0℃水的黏度，值为1.758*10-3。
其中μ要求分别以长短格式显示计算结果。
操作题2 《Matlab与科学计算》P31 例2.2数值数组【1：5】和字符串转换示例。
操作题3《Matlab与科学计算》P44 例2.9比较用左除和右除法分别求解恰定方程的解。
操作题4 《Matlab与科学计算》 P48 例2.14计算矩阵magic(3)的指数，并比较不同函数的结果。
指数运算expm(a)操作题5 《Matlab与科学计算》 P50 例2.18计算下面矩阵[-149 -50 -154；537 180 546；-27 -9 -25]的特征值条件数。
特征值函数condeig(A)操作题6《Matlab与科学计算》P62 例2.29对矩阵a = 1 1 1 11 2 3 41 3 6 101 4 10 20 进行三角抽取使之分别变成0 0 0 0 1 1 1 01 0 0 0 及 123 41 3 0 0 1 3 6 101 4 10 0 1 4 10 20操作题7《Matlab与科学计算》P62 例2.30用3种方法建立向量v 为[1 2 6 20]多项式的伴随矩阵。
操作题8《Matlab与科学计算》P64 例2.31已知a=[2 1 -3 -1；3 1 0 7；-1 2 4 -2；1 0 -1 5]；求该矩阵的3次方及数组的3次方。
操作题9《Matlab与科学计算》P66 例2.32 数组逻辑运算演示。
已知a=[1:3；4：6；7：9]，b=[0 1 0；1 0 1；0 0 1]， x=5, y=ones （3）*5,求x<=a，ab=a&b，～b操作题10《Matlab与科学计算》P71 例2.41对多项式p =[2 -5 6 -1 9]求根，计算其在x=3的值，进行微分计算并写出其表达式。
matlab课后习题及答案
第一章 5题已知a=4.96，b=8.11，计算)ln(b a eba +-的值。
解：clear clc a=4.96; b=8.11;exp(a-b)/log(a+b) ans =0.0167 6题已知三角形的三边a=9.6,b=13.7, c=19.4,求三角形的面积。
提示：利用海伦公式area =))()((c s b s a s s ---计算，其中S=(A+B+C)/2. 解：clear clc a=9.6; b=13.7; c=19.4; s=(a+b+c)/2area=sqrt(s*(s-a)*(s-b)*(s-c)) s =21.3500 第二章 8题已知S=1+2+2^2+2^3+……+2^63,求S 的值 解：clear clc S=0;for i=0:1:63 S=S+2^i; end S S =1.8447e+019 9题分别用for 和while 循环结构编写程序，计算∑=-1001n 1n 2）（的值。
解：clear clc s=0;for n=1:100 s=s+(2*n-1); end s s =10000 clear clc n=1; s=0;while n<=100 s=s+(2*n-1); n=n+1; end s s =10000 第三章 2题在同一坐标下绘制函数x ,,2x-,2x xsin(x)在()∏∈,0x 的曲线。
解：clear clcx=0:0.2:pi; y1=x; y2=x.^2; y3=-(x.^2); y4=x.*sin(x);plot(x,y1,'-' ,x,y2,'-' ,x,y3,'-' ,x,y4,'-')0.511.522.53-10-8-6-4-202468109题用不同的线型和颜色在同一坐标内绘制曲线y1=2ex5.0 、y2=sin(2∏x )的图形。
matlab复习题标准答案
matlab复习题标准答案# MATLAB复习题标准答案一、选择题1. MATLAB中的矩阵默认是按什么顺序存储的？- A. 行主序- B. 列主序- C. 随机序- D. 行主序2. 下列哪个函数用于计算矩阵的行列式？- A. `det`- B. `diag`- C. `trace`- A. `det`3. MATLAB中如何创建一个5x5的单位矩阵？- A. `eye(5,5)`- B. `zeros(5,5)`- C. `ones(5,5)`- A. `eye(5,5)`4. 下列哪个命令用于生成0到1之间的随机数？- A. `rand()`- B. `randi()`- C. `randn()`- A. `rand()`5. MATLAB中如何计算向量的范数？- A. `norm(vector)`- B. `abs(vector)`- C. `sum(vector)`- A. `norm(vector)`二、填空题1. MATLAB中的`%`符号用于______。
- 答案：注释2. 要创建一个从10到1的递减向量，可以使用______。
- 答案：`1:-1:10`3. MATLAB中，`zeros(3,4)`会生成一个______。
- 答案：3行4列的全零矩阵4. 函数`max(A)`会返回矩阵A中的______。
- 答案：最大值5. `plot(x,y)`函数在MATLAB中用于______。
- 答案：绘制x和y的数据点的图形三、简答题1. 请简述MATLAB中矩阵的基本操作有哪些？- 答案：矩阵的基本操作包括矩阵的创建、转置、求逆、矩阵乘法、矩阵加法、矩阵减法等。
2. MATLAB中如何实现循环结构？- 答案：MATLAB中实现循环结构主要有两种方式，一种是使用`for`循环，另一种是使用`while`循环。
3. 请解释MATLAB中的函数文件和脚本文件的区别。
- 答案：函数文件通常包含返回值和参数，可以被其他脚本或函数调用；而脚本文件主要用于执行一系列命令，不包含返回值。
matlab课后习题与答案
习题二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 。
中南大学matlab课后习题（10）
中南大学matlab课后习题（10）Unit 1实验内容1.答：用help命令可以查询到自己的工作目录。
输入help命令：help <函数名>2.答：MATLAB的主要优点：通过例1-1至例1-4的验证，MATLAB的优点是MATLAB以矩阵作为数据操作的基本单位，使得矩阵运算变得非常简捷，方便，高效。
还提供了丰富的数值计算函数。
MATLAB绘图十分方便，只需输入绘图命令，MATLAB便可自动绘出图形。
3.答：INV(X) is the inverse of the square matrix X。
A warning message is printed if X is badly scaled or nearly singular. PLOT(X,Y) plots vector Y versus vector X. If X or Y is a matrix, then the vector is plotted versus the rows or columns of the matrix, whichever line up. If X is a scalar and Y is a vector, length(Y) disconnected points are plotted. PLOT(Y) plots the columns of Y versus their index. If Y is complex, PLOT(Y) is equivalent to PLOT(real(Y),imag(Y)).In all other uses of PLOT, the imaginary part is ignored. For vectors, MAX(X) is the largest element in X. For matrices,MAX(X) is a row vector containing the maximum element from each column. For N-D arrays, MAX(X) operates along the first non-singleton dimension. [Y,I] = MAX(X) returns the indices of the maximum values in vector I. If the values along the first non-singleton dimension contain more than one maximal element, the index of the first one is returned. ROUND(X) rounds the elements of X to the nearest integers. MAX(X,Y) returns an array the same size as X and Y with the largest elements taken from X or Y. Either one can be a scalar。
中南大学matlab课后答案-第十一章
第十一章实验指导1，(1)(2)2，M文件的函数为：function x=f(t)x=(8*1/(pi^2))*(cos(2*pi*t)+(1/9)*cos(3*2*pi*t)+(1/25)*cos(5*2*pi*t))3,4,5,5，限幅：function [sys,x0,str,ts]=xianSfun(t,x,u,flag,max,min) switch flag,case 0sizes= simsizes;sizes.NumContStates = 0;sizes.NumDiscStates = 0;sizes.NumOutputs = 1;sizes.NumInputs = 1;sizes.DirFeedthrough = 1;sizes.NumSampleTimes = 1;sys= simsizes(sizes);str =[];x0 =[];ts =[-1 0];case 3if u>maxsys=max;elseif u<minsys=min;elsesys=u;endcase{1,2,4,9}sys=[];otherwiseerror(['Unhandled flag=',num2str(flag)]); end绝对值：function [sys,x0,str,ts]=abs(t,x,u,flag)switch flagcase 0sizes= simsizes;sizes.NumContStates = 0;sizes.NumDiscStates = 0;sizes.NumOutputs = 1;sizes.NumInputs = 1;sizes.DirFeedthrough = 1;sizes.NumSampleTimes = 1;sys= simsizes(sizes);str = [];x0 =[];ts =[0 0];case 3if u>0sys=u;else u<0sys=-u;endcase {1,2,4,9}sys=[];otherwiseerror(['Unhandled flag=',num2str(flag)]);end思考练习：1，两个主要功能Simu仿真和Link连接即该软件；天价模块——模块编辑——模块的连接——模块的参数和属性设置——设置仿真参数——启动系统仿真与仿真结果分析2，适用sinks模块库中的Scope模型或其他显示模块来观察有关信号的变化曲线，也可以将结果存放在MATLAB中的工作空间，工艺后处理和使用。
MATLAB练习题和复习资料1
控制系统仿真实验部分实验结果目录实验一基本操作............................. 错误!未指定书签。
实验二编程............................... 错误!未指定书签。
实验三底层图形控制......................... 错误!未指定书签。
实验四控制系统古典分析..................... 错误!未指定书签。
实验五控制系统现代分析...................... 错误!未指定书签。
实验六控制器的设计........................ 错误!未指定书签。
实验七系统状态空间设计..................... 错误!未指定书签。
实验九直流双闭环调速系统仿真............... 错误!未指定书签。
实验一基本操作1 用可以识别的格式输入下面两个矩阵再求出它们的乘积矩阵C，并将C矩阵的右下角2×3子矩阵赋给D矩阵。
赋值完成后，调用相应的命令查看工作空间的占用情况。
[1,2,3,3;2,3,5,7;1,3,5,7;3,2,3,9;1,8,9,4];[1+4i,4,3,6,7,8;2,3,3,5,5,4+2i;2,6+7i,5,3,4,2;1,8,9,5,4,3];*B;(4:5,4:6);A 5x4 160B 4x6 384C 5x6 480D 2x3 962 选择合适的步距绘制出下面的图形，其中[-1:0.1:1];(1);3 对下面给出的各个矩阵求取矩阵的行列式、秩、特征多项式、范数、特征根、特征向量和逆矩阵。
[7.5,3.5,0,0;8,33,4.1,0;0,9,1031.5;0,0,3.7,19.3];[5,7,6,5;7,10,8,7;6,8,10,9;5,7,9,10];[1:4;5:8;9:12;13:16];[332,4;55,1,8;11,8,57;5131];(A)(B)(C)(D);(A);(B);(C);(D);(A);(B);(C);(D);(A);(B);(C);(D);[](A,'');[](B,'');[](C,'');[](D,'');(A);(B);(C);(D);4 求解下面的线性代数方程，并验证得出的解真正满足原方程。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
2020/3/3
11
2、利用非线性方程组的符号法求解如下方程组：
x1 0.7sin x1 0.2cos x2 0
x2 0.7 cos x1 0.2sin x2 0
在（0.5，0.5）附近的解。
（a）
a) x1=0.5265 x2=0.5079
b) x1=0.5079 x2=0.5265
ans =
x0=[0.5 0.5] >> fsolve('fc',x0)
0.5265 0.5079
2020/3/3
10
九、方程符号求解
1 求微分方程 y(0)=0的解
x
d2 y dx2
3
dy dx
x2
，y(1)=0，
（c）
a)x^4/5 - x^3/4 b)x^3/3 - x^2/3 c) x^4/3 - x^3/3 d) x^4/5- x^3/4
>> [x,y]=dsolve('Dx=y,Dy=x','t')
2020/3/3
9
非线性方程的符号注解
语法： fsolve(‘fun’, x0) （1）fun为所要求解的函数名，常以M文件形式给出； （2）X0为初始向量或矩阵
Fc.m function y=fc(x) y(1)=x(1)-0.7*sin(x(1))0.2*cos(x(2)); y(2)=x(2)0.7*cos(x(1))+0.2*sin(x(2)); y=[y(1) y(2)]
2020/3/3
8
【例
】 求微分方程
x
d2 y dx2
3
dy dx
x2
，y(1)=0，y(0)=0的解。
>> y=dsolve('x*D2y3*Dy=x^2','x')
%求微分方程的通解
y=
C3*x^4 - x^3/3 + C2
>> y=dsolve('x*D2y3*Dy=x^2','y(1)=0,y(0)=0','x') %求微分方程的特解
-x+10y-2x=7
-2x+10z=6 的解
2020/3/3
5
2020/3/3
6
九、方程符号求解
2020/3/3
7
符号常微分方程
MATLAB提供了dsolve命令，可以用于对符号常微分方程进行求解。
语法：
dsolve('eq', 'con', 'v') dsolve('eq1,eq2…', 'con1,con2…', 'v1,v2…')
（1）当y是因变量时，微分方程'eq'的表述规定为：
y的一阶导数
dy dx
或dy
dt
表示为Dy；
%求解微分方程 %求解微分方程组
y的n阶导数 dn y 或dn y 表示为Dny。
dxn dtn
（2）微分初始条件'con'应写成'y(a)=b，Dy(c)=d'的格式；当初始条件少于
微分方程数时，在所得解中将出现任意常数符C1、C2……，解中任意常数符的数目 等于所缺少的初始条件数。
c) x1=0.6265 x2=0.7079
d) x1=0.7079 x2=0.6265
2020/3/3
12
1 对于圆周率值，若取值为3.1415，则其
有数字为
（b）
a、3位
b、4位
c、5位 d、2位
2020/3/3
13
2020/3/3
14
2下面的MATLAB语句中正确的有：（cd）
a) a=sin(x)
2020/3/3
20
3 试用一个二次多项式拟合下列数据 x 0.5 1.0 1.5 2.0 y 1.75 2.45 3.81 4.80
2.5 8.00
3.0 8.60
（d）
a) 0.4600 1.8601 0.9730 b) 0.8400 2.9401 0.85608 c) 0.4900 1.8601 0.9730 d) 0.4900 1.2501 0.8560
(a) return; (b) break (c) continue (d) keyboad
2020/3/3
17
二、多项式及其运算
2020/3/3
18
二、多项式及其运算
1 求多项式 x4 12x3 25x 116 在x=3,8处的值
a)-52 -1730
（c）
b) -54 -1732
c) [4 19 3;21 75 54; 22 72 72]
d) [6 25 63;0 36 18; 1 15 0]
2020/3/3
23
（b） 2、求数组b=[1 2 3;4 5 6;7 8 9]的指数运算
b)Record-1=3+4i
c) a=sin(pi)
d) C=1+6j
2020/3/3
15
3 符合MATLAB变量命名规则的有a) var!205 b)VAR15（b）
c)my-var2015 d)2015_VAR
2020/3/3
16
4在循环结构中跳出循环，但继续下次循环
的命令为
（c）
c) -52 -1732
d) -54 -1730
2020/3/3
19
2 计算多项式x3+11x2+55x+125 在x=[1
2;3 4]的矩阵计算值
（a）
a)[192 287;416 585]
b)[190 285;406 582]
c)[194 292;422 580]
d)[192 288;420 580]
习题课
Any 函数
判断元素是否为非零元素any(v),如果v是 非零元素返回true(即1)否则返回flase(即 0)
如果A是矩阵，则把矩阵的列当做向量来处 理，函数返回每个列向量的逻辑值
2020/3/3
2
2020/3/3
3
Magic函数
2020/3/3
4
线性方程组求解
试求数值方程组 10x-y=9
2020/3/3
21
三、矩阵及数组代数运算
2020/3/3
22
三、矩阵及数组代数运算
1、已知A=[0 9 6;1 3 0;3 5 7];B=[1 4 3;1 5 0;2 5
9], 则A.*B
（b）
a) [21 75 54; 4 19 3; 22 72 72]
b) [0 36 18; 1 15 0; 6 25 63]
y=
x^4/3 - x^3/3
【例 】 求微分方程组
dx y, dy x dt dt
的解。
>> [x,y]=dsolve('Dx=y,Dy=x') x= C8*cos(t) + C7*sin(t) y= C7*cos(t) - C8*sin(t)
程序分析：默认的自由变量是t，C1、C2为任意常数，程序也可指定自由变量， 结果相同：