MATLAB期末考试试题(卷)与答案

一、填空题（每空2分，总共30分）

1.已知A=[0 9 6;1 3 0];B=[1 4 3;1 5 0];写出下列各指令运行的结果。

A &

B 的运行结果ans= ；

A ./

B 的运行结果ans= 。

2. 产生4阶全0方阵的命令为 zeros(4) ；产生3阶全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语言后的判断应写为 if x<=80&x>60 。

6.P, Q分别是个多项式的系数矢量，求P对应的多项式的积分（对应的常数项为K），使用的命令是 polyint(P,K) ；求P/Q的解，商和余数分别保存在k和r，使用的命令是 [k,r]=deconv(P,Q) ；

7.为了使两个plot的图形在同一个坐标显示，可以使用 hold on 命令进行图形保持；可以使用 grid on 命令为图形添加网格。

8.MATLAB的工作空间中有三个变量v1, v2, v3，写出把它们保存到文件my_data.mat 中的指令 save my_data ；写出把my_data.mat文件中的变量读取到MATLAB 工作空间内的指令 load my_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) clf

5.用round函数四舍五入对数组[2.48

6.39 3.93 8.52]取整，结果为 C

(A) [2 6 3 8] (B) [2 6 4 8] (C) [2 6 4 9] (D) [3 7 4 9]

6.已知a=2:2:8, b=2:5，下面的运算表达式中，出错的为 C

(A) a'*b (B) a .*b (C) a*b (D) a-b

7.角度[]60

x，计算其正弦函数的运算为D

45

=

30

(A) SIN（deg2rad(x)） (B) SIN(x) (C) sin(x) (D) sin(deg2rad(x))

8.下面的程序执行后array的值为 ( A )

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.

9．i=2; a=2i; b=2*i; c=2*sqrt(-1); 程序执行后；a, b, c的值分别是多少？

(A)a=4, b=4, c=2.0000i (C)

(B)a=4, b=2.0000i, c=2.0000i

(C)a=2.0000i, b=4, c=2.0000i

(D) a=2.0000i, b=2.0000i, c=2.0000i

10. 求解方程x4-4x3+12x-9 = 0 的所有解

(A)1.0000, 3.0000, 1.7321, -1.7321

(B)1.0000, 3.0000, 1.7321i, -1.7321i

(C)1.0000i, 3.0000i, 1.7321, -1.7321

(D)-3.0000i, 3.0000i, 1.7321, -1.7321

三、写出程序的执行结果或写出给定要求的指令（总共35分）

1.写出执行以下代码后C，D，E的值 (6分)

A=[1,2,3;4:6;7:9];

C=[A;[10,11,12]],

D=C(1:3,[2 3])

E=C(2,[1 2])

2.写出执行以下代码后，MATLAB命令窗口上显示的x矩阵的值 (5分)

x=[0,1,0,2,0,3,0,4];

for k=1:8

if x(k)==0

x(k)=k;

else

x(k)=2*k+1;

end

end

disp(x);

3.创建符号函数并求解，要求写出步骤和运行结果（7分）

(1)创建符号函数f=ax2+bx+c

(2)求f=0的解

4. 求解以下线性方程组，要求写出程序代码和运行结果（5分）

2x1- 3x2+ x3+2x4=8

x1+3x2+ x4=6

x1- x2+ x3+8x4=1

7x1+ x2-2x3+2x4=5

5．绘制函数曲线，要求写出程序代码(12分)

(1)在区间[0:2π]均匀的取50个点，构成向量π

(2)在同一窗口绘制曲线y1=sin(2*t-0.3); y2=3cos(t+0.5)；要求y1曲线为红色点划线，标记点为圆圈；y2为蓝色虚线，标记点为星号

四、使用MATLAB语言进行编程（15分）

打印出所有的水仙花数。所谓“水仙花数”，是指一个三位数，其各位数字立方之和等于该数本身。

一、填空题

1、

ans =

0 1 1

1 1 0

ans =

0 2.2500 2.0000

1.0000 0.6000 NaN

2、 zeros(4) ones(3)

3、 [ 2 5 ] [ 5 ]

4、；%

5、 if x<= 80 & x>60

6、 polyint(P,K) [k,r]= deconv (P,Q)

7、 hold on grid on

8、 save my_data Load my_data

二、选择题

1、 A

2、 C

3、 B

4、 C

6、 C

7、 D

8、 A

9、 C

10、 A

三、写出程序的执行结果或写出给定要求的指令

1、

c =

0 + 2.0000i

D =

2 3

5 6

8 9

E =

4 5

2、

1 5 3 9 5 13 7 17

>> syms a b c x

>> s = a*x^2 + b*x + c;

>> solve(s)

ans =

[ 1/2/a*(-b+(b^2-4*a*c)^(1/2))] [ 1/2/a*(-b-(b^2-4*a*c)^(1/2))]

4、

>> syms a b c x

>> syms x1 x2 x3 x4

>> clear

>> syms x1 x2 x3 x4

>> f = 2*x1 + 3*x2 + x3 + 2*x4 -8; >> g = x1 + 3*x2 + x4 - 6;

>> h = x1 - x2 + x3 + 8*x4 - 1; >> i = 7*x1 + x2 - 2*x3 + 2*x4 -5; >> [x1,x2,x3,x4] = solve(f,g,h,i)

x1 =

13/17

x2 =

320/187

x3 =

212/187

x4 =

19/187

5、

t = linspace(0, 2*pi, 50);

plot(t,y1,'r-o' ,t,y2,'b--*');

四、使用MATLAB语言进行编程

for k = 100:999

a = fix(k/100);

b = rem(fix(k/10),10);

c = rem(k,10);

if a.^3 + b.^3 + c.^3 == k

fprintf( “%u,\t\t”, k ) ;

end

end

三相异步电动机Matlab仿真

中国石油大学胜利学院综合课程设计总结报告 题目：三相异步电机直接启动特性实验模型 学生姓名：潘伟鹏 系别：机械与电气工程系 专业年级： 2012级电气工程专业专升本2班 指导教师：王铭

2013年 6 月 27日

一、设计任务与要求 普通异步电动机直接起动电流达到额定电流的6--7倍，起动转矩能达到额定转矩的1.25倍以上。过高的温度、过快的加热速度、过大的温度梯度和电磁力，产生了极大的破坏力，缩短了定子线圈和转子铜条的使用寿命。但在电网条件和工艺条件允许的情况下，异步电动机也可以直接启动。本次课程设计通过MATLAB软件建模模拟三相异步电动机直接启动时的各个元器件上的电量变化。 参考： 电力系统matlab仿真类书籍 电机类教材 二、方案设计与论证 三相异步电动机直接起动就是利用开关或接触器将电动机的定子绕组直接接到具有额定电压的电网上。 由《电机学》知三相异步电动机的电磁转矩M与直流电动机的电磁转矩有相似的表达形式。它们都与电机结构（表现为转矩常数）和每级下磁通有关，只不过在三相异步电动机中不再是通过电枢的全部电流，而是点数电流的有功分量。三相异步电机电磁转矩的表达式为： （1-1）式中——转矩常数 ——每级下磁通 ——转子功率因数 式(1-1)表明，转子通入电流后，与气隙磁场相互作用产生电磁力，因此，反映了电机中电流、磁场和作用力之间符合左手定则的物理关系，故称为机械特性的物理表达式。该表达式在分析电磁转矩与磁通、电流之间的关系时非常方便。 从三相异步电动机的转子等值电路可知， （1-2） （1-3）将式(1-2)、(1-3)代入(1-1)得：

《Matlab语言》期末考试试题样例

一、填空题：（每空2分）大约20空 （注：写出完成要求所需的命令语句或题目给定语句执行后输出的 结果） 1. Matlab 中用来表示不确定值的变量是 NaN ；默认变量i 表示 虚数单位 。 2. Matlab 中符号 ； （分号） 可以使命令行不显示运行的结果，表示该行为注释行的符号为 % （百分号） 。 3. Matlab 将工作空间所有变量保存在默认目录下的命令是 Save 。(命令复习) 4. 骰子会掷出1到6的随机数字，用函数rand 和取整函数(Important)模拟骰子某次投掷出的数字的命令为 ceil(6*rand) ，模拟两 个骰子某次投掷的数字组合之和情况的命令为 。 5. 生成一个含有10个元素的行向量x ，其元素在[-1,1]之间均匀分布的的命令为 x=2*rand(1,10)-1 ，使用向量x 的后5个元素作为对角线 元素生成对角阵A 的命令为 A=diag(x(end-4:end)) ，将A 的主对 角线下方的元素（不含主对角线）全部填充为2014的命令为 B=ones(5); 2014*B-triu(B)+A 6. 为函数22 110at t x +=分布定义内联函数命令是 ，定义匿名函数的命令是 7. 在循环结构中使用 break 语句可以终止当前循环。 8. Matlab 用于绘制统计直方图的命令是 hist ；用于绘制y 轴为对数坐标的单对数坐标图的命令为 semilog 。

二、解答题：大约七道题目 1. 在[0,π2]区间上绘制如下图所示的sin(x)与sin(x 2)的曲线，并如图为图像添加标题以及为坐标轴添加标注。 2.绘制上面曲面的图形（提示：该曲面由上下两部分组成：） 3.创建命令式的M 文件计算圆形的面积S ，提示用户输入圆形的半径r ，并输出圆形的面积。（注意输入输出都要对用户添加一定的提示语） %Function Circle 4.请利用matlab 的符号工具箱做下面的运算，并写出相应的命令。 (1) 求极限x e x L x x -+=→10)1(lim ； (2)求积分：dx x x ?∞+02sin ； Syms x; f=((1+x)^(1/x)-exp(1)/x; limit(f,x,0) 2222, [2,2], [1,1]z x y x y =+∈-∈-

Matlab 期末考试题库(共12套卷)

第一套 华东交通大学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%折扣

西安石油大学MATLAB期末试卷

实用文档 标准文案 一、填空题。（30分，每小空1分） 1. MATLAB 变量命名，需要满足3条规则，分别是 、 、 2. 将变量值进行最小整数比的有理格式命令： 3. 标点符号_______可以使命令行不显示运算结果， 用来表示该行为注释行。 4. x 为0～π，步长为0.1π的向量，创建x 语句为 。 5. 输入矩阵A=[1 2 3 4;-2 -3 -4 -5;0 6 7 8]，使用全下标方式用 取出元素“-5”，使用单下标方式用_______取出元素“-5”。 6. 对多行进行注释的快捷键为 对多行去掉注释的快捷键为 7. 设x 是一维数组，x 的倒数第3个元素表示为_______； 8. 设y 为二维数组，要删除y 的第34行和48列，可使用命令_______;_______； 9. 在while 表达式, 语句体, end 循环语句中，表达式的值为____时表示循环条件为真，语句体将被执行，否则跳出该循环语句； 10. 要从键盘读入一个字符串并赋值给变量x ，且给出提示“What is your name?”，应使用命令_________； 11．已知A 、B 和C 均为n*n 矩阵，且存在于WorkSpace 中，要产生矩阵 () 1 ()C A C D A B B --= ，可用命令________ _，计算 12. 已知A=[123；450；789]；B=[103；150；012]；写出下列各指令运行的结果。A+B ；A.*B ；A==B ； 13. 已知A 是矩阵，求A 的对角矩阵B 的命令是_______，求A 的下三角矩阵C 的命令是_______。 14. 生成正态分布4?6随机矩阵A 的命令是 ，删除矩阵A 的第一列和第三列命令是 __________________ 15. 已知A=[1 5 3;4 2 6; 11 1 0]; A(7)=__________。

哈工大matlab期末考精彩试题精彩试题及问题详解(95分)分解

春季学期MATLAB期末作业 学院：机电工程学院 专业：机械制造设计及其自动化 学号： 班号： 姓名：

2013年春季学期 MATLAB 课程考查题 姓名: 学号： 学院：机电学院 专业：机械制造 一、必答题： 1.matlab常见的数据类型有哪些？各有什么特点? 常量：具体不变的数字 变量：会根据已知条件变化的数字 字符串：由单引号括起来的简单文本 复数：含有复数的数据 2.MATLAB中有几种帮助的途径？ （1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器； （2）help命令：在命令窗口键入“help”命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息； （3）lookfor命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数 （4）模糊查询：输入命令的前几个字母，然后按Tab键，就可以列出所有以这几个字母开始的命令和函数。 注意：lookfor和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help命令显示详细信息。 3.Matlab常见的哪三种程序控制结构及包括的相应的语句？ 1.顺序结构：数据输入A=input(提示信息,选项) 数据输出disp(X) 数据输出fprintf(fid,format,variables) 暂停pause 或 pause(n) 2.选择结构： If语句： if expression （条件）

statements1（语句组1） else statements2（语句组2） End Switch 语句： switch expression （表达式） case value1 （表达式1） statement1（语句组1） case value2 （表达式2） statement2（语句组2） ... ... case valuem （表达式m） statementm（语句组m） otherwise statement （语句组） end 3.循环结构： For循环： for variable=expression(循环变量) statement（循环体） end While循环： while expression （条件<循环判断语句>） statement（循环体） end 4.命令文件与函数文件的主要区别是什么？

MATLAB复习题.doc

MATLAB期末考试题目 ------电气二班 一．选择题。 1．下列哪个变量的定义是不合法的（A） A. abcd-3 B. xyz_3 C. abcdef D. x3yz 2．在同一个图形窗口中画出三行两列的图形并指定当前可编辑图形为第3个图 的正确的 MATLAB命令是： A.subplot(3,2,1)B subplot(3,3,) C.subplot(3,2,3 ) D.subplot(1,3,2) 3.x=[1234],y=[1111],x.*y的结果是：(B) A．1234 B.出错信息 C.4 3 2 1 D.1111 4．角度x=[30 45 60]，计算其正弦函数的运算为（D） A. sin(deg2rad(x)) B. sin(x) C. sin(x) D. sin(deg2rad(x)) 5．i=2;a=2i;b=2*i;c=2*sqrt(-1);执行程序后，a,b,c的值分别是多少？（C） A.a=4,b=2,c=2.0000i B.a=4,b=2.0000i,c=2.0000i C.a=2.0000i,b=4,c=2.0000i D.a=2.0000i,b=2.0000i,c=2.0000i 6. legend()命令是(A) A．在图中标注“图例” B.．在图中加上网格线C．在图中加上标题 D．打印图形 7.命令plot(x,y,′m:*′),其中字符串“m:o”中,第一个字符“m”(A) A. 表示曲线为m符号构成的线 B. 表示曲线颜色为紫色 C. 表示曲线x坐标的最大值为m变量的值 D. 表示曲线y坐标的最大值为m变量的值 8. 下列可作为matlab合法变量名的是(D) A、合计 B、123 C、@h D、xyz_2a 9. 使用语句t=0:7生成的是()个元素向量(A) A、8 B、7 C、6 D、5 10. 输入字符串或选项时，要用(D)括住 A、() B、【】 C、{} D、‘’ 11 . 清空 Matlab 工作空间内所有变量的指令是（A）。 A.clc B. cls C.clear D. clf

matlab 上机作业1,中国石油大学

第一次上机实习 1．编写一个程序ex1_1.m，对输入的方阵计算其伴随矩阵。 a)要求第一行帮助行注明这个函数的用途，能在matlab命令窗口可以用help或lookfor 查到该帮助； b)对输入的变量进行检测，如果输入的不是方阵或者矩阵元素不是数，则报错,退出程 序。 function A=ex1_1(X) % 建立函数A=ex1_1(X)，求矩阵X的伴随矩阵 [a,b]=size(X); if a~=b||~isnumeric(X) || ~ismatrix(X) error('输入的不是方阵或者矩阵元素不是数!'); end for i=1:a for j=1:b A(i,j)=(-1)^(i+j)*X(i,j)*det(X([1:i-1,i+1:end],[1:j-1,j+1:end])); end end A=A'; 2．已知表达式y=6*x^5+4*x^3+2*x^2-7x+10，x的范围是[0，100]，使用三阶拟合和五阶的方法得出多项式的表达式，并编程在图中绘制出原曲线、三阶拟合和五段拟合的曲线。 clc clear all x=linspace(0,100,2000); p1=[6 0 4 2 -7 0]; y=polyval(p1,x); p3=polyfit(x,y,3); p5=polyfit(x,y,5); y1=polyval(p3,x); y2=polyval(p5,x); plot(x,y,'r-o') hold on grid on plot(x,y1,'b-.',x,y2,'y--*'); grid on 三条曲线

放大图（原曲线与五阶拟合完全重合）： 3．读如图像5.tif（450*450），做fourier变换， a)将变换之后所得到的矩阵只保留前10行10列，后10行10列的数据，其他位置都 设置为0，然后再作反fourier变换，观察结果； 原图：

基于MATLAB的微机保护算法仿真

本科毕业设计(论文) 题目：基于MATLAB的微机保护算法仿真 学生姓名： 学号：05843420 专业班级：电气工程及其自动化1班 指导教师：

基于MATLAB的微机保护算法仿真 摘要 基于MATLAB软件，运用Simulink工具完成一种继电保护微机保护数据采集和处理系统。主要基于两点法、突变量算法、对称分量选相法等传统的微机保护算法，搭建MATLAB的仿真模型，根据采集的输入电气量的采样数据进行分析、运算和判断，以实现相应的继电保护功能。本文对MATLAB软件如何应用于微机保护做了详细说明，并运用MATLAB的动态仿真工具对电力系统中的故障以及微机保护的算法进行了仿真分析，说明了MATLAB在微机继电保护中应用的可行性。同时本文对工程中常用微机保护算法进行了原理层面的分析，并运用模型对各种算法的性能进行了仿真和研究，有很重要的现实意义。 关键词：MATLAB；微机保护；算法；采样数据

Microcomputer protection algorithm based on MATLAB simulation Abstract Based on MATLAB software, using Simulink tools perform one computer protection relay data acquisition and processing systems. Mainly based on two-point method, the amount of mutation algorithms, symmetric phase selector and other traditional computer protection algorithms to build MATLAB simulation model, based on the amount of collected samples of the input electrical data analysis, calculation and judgment, in order to achieve the corresponding relay protection. In this paper, MATLAB software how to apply a detailed description of computer protection, and the use of MATLAB dynamic simulation tool for power system failures and computer protection algorithms for the simulation analysis, illustrates the application of MATLAB in the feasibility of relay protection . This paper also commonly used in engineering computer protection algorithms theory level of analysis, and the use of models for the performance of various algorithms and simulation studies, there is a very important practical significance. Keywords：MATLAB；Microcomputer Protection；Algorithm; Sampling Data

MATLAB《数学软件及应用》公选课期末考试题目

《数学软件及应用》期末考试题目 1.计算表达式 5cos(/5) pi +x 。 2. 在matlab 的命令窗口里完成如下计算，其中t 的值分别取-2,-1,0,1,2，表达式如下： 4/3) 2t y e π-=+。 3.计算23sin cos 0.5(1)x y x x ??=+??+?? 在x=1,2,3处的函数值 4.用语句a=fix(rand(4,4)*20)产生一个矩阵。解释下列语句的功能： (:,3)b a =， ([1,2],[3,4])c a =， a(:,5)=[ ] 5. 已知12536 4a ?? =??-??，8743 62b -??=???? 写出以下指令的结果并加以解释： a>=b a.*b 6.用创建数组的方法生成0~100以内3的倍数，将此数组赋给变量d 。（要求使用冒号法和linspace 两种方法） 7.创建一个数组，其元素分布在区间[0, 2]上。(使用冒号法和linspace) 8.自定义一个矩阵，取出它的第二行；将它第一行第二列的元素改为3. 9.求方程432234180x x x x ----=的所有根。 10.求方程323210x x x --+=的根。 11.对于B AX =，如果4 92764357A ????=? ????? ，372628B ????=??????，求解X 12.解线性方程组?? ? ????=+-+-=-+=++-=--+410256958424754 12743w z y x w z x w z y x w z y x 13.解非线性方程组221231 x xy y x y ?--=?-=? 14.求解微分方程2 /5)1(12+=+-x x y dx dy ； 15.求解微分方程2 )0(,1)0(),sin(32322='==++y y x y dx dy dx y d

中国石油大学(北京)自动控制原理胡寿松 课件PPT

实验四 控制系统的根轨迹分析 一. 实验目的: 1. 学习利用MATLAB 语言绘制控制系统根轨迹的方法。 2. 学习利用根轨迹分析系统的稳定性及动态特性。 二. 实验内容: 1. 应用MATLAB 语句画出控制系统的根轨迹。 2. 求出系统稳定时，增益K 的范围。 3. 实验前利用图解法画出系统的根轨迹，算出系统稳定的增益范围，与实测值相比较。 4. 应用SIMULINK 仿真工具，建立闭环系统的实验方块图进行仿真。观察不同增益下系统的阶跃响应，观察闭环极点全部为实数时响应曲线的形状；有共轭复数时响应曲线的形状。（实验方法参考实验二） 5. 分析系统开环零点和极点对系统稳定性的影响。 三. 实验原理: 根轨迹分析法是由系统的开环传递函数的零极点分布情况画出系统闭环根轨迹，从而确定增益K 的稳定范围等参数。假定某闭环系统的开环传递函数为 )164)(1() 1()()(2++-+=s s s s s K s H s G 利用MATLAB 的下列语句即可画出该系统的根轨迹。 b=[1 1]； ％确定开环传递函数的分子系数向量 a1=[l 0]； ％确定开环传递函数的分母第一项的系数 a2=[l -1]； ％确定开环传递函数的分母第二项的系数 a3=[l 4 16]； ％确定开环传递函数的分母第三项的系数 a=conv(al ，a2)； ％开环传递函数分母第一项和第二项乘积的系数 a=conv(a ，a3)； ％分母第一项、第二项和第三项乘积的系数 rlocus(b,a) ％绘制根轨迹，如图（4-l ）所示。 p=1.5i ； ％ p 为离根轨迹较近的虚轴上的一个点。 [k ，poles]=rlocfind(b ，a ，p) ％求出根轨迹上离p 点很近的一个根及所对应

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、简述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’);

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、建立结构数组或转换结构数组的函数为：struct

ADI(交替方向隐格式)求解二维抛物方程(含matlab程序)

ADI 法求解二维抛物方程 学校：中国石油大学（华东） 学院：理学院 姓名：张道德 时间：2013.4.27 1、ADI 法介绍 作为模型，考虑二维热传导方程的边值问题： （3.6.1）,0,,0(,,0)(,) (0,,)(,,)(,0,)(,,)0t xx yy u u u x y l t u x y x y u y t u l y t u x t u x l t ?=+<<>?? =??====? 取空间步长1 h M =，时间步长0t >，作两族平行于坐标轴的网线： ,,,0,1,, j k x x jh y y kh j k M ===== 将区域0,x y l ≤≤分割成2M 个小矩形。第 一个ADI 算法（交替方向隐格式）是Peaceman 和Rachford （1955）提出的。 方法： 由第n 层到第n+1层计算分为两步： （1） 第一步： 1 2,12 n j k xx yy u + 从n->n+ ，求u 对向后差分,u 向前差分，构造出差分格 式为： 1 (3.6.1)11 11 2 22 2 ,,1,,1, ,1,,1 2 2 1 222,,2 -22=2 1() n n n n n n n n j k j k j k j k j k j k j k j k n n x j k y j k h h h τδδ+ + + + +-+-+-+-+= +u u u u u u u u （+ ） u u （2） 第二步：1 2,12 n j k xx yy u + 从n+ ->n+1，求u 对向前差分,u 向后差分，构造出差分格 式为： 2 (3.6.1)1 1 1 1 111 222,,1,,1, ,1,,1 2 2 1 221 2,,2 -22=2 1() n n n n n n n n j k j k j k j k j k j k j k j k n n x j k y j k h h h τδδ+ + + +++++-+-++-+-+= +u u u u u u u u （+ ） u u 其中1211,1,,1,0,1,2,,()2 2 n j k M n n n τ+ =-=+ =+ 上表表示在t=t 取值 。 假定第n 层的,n j k u 已求得，则由1(3.6.1)求出1 2,n j k u + ，这只需按行 (1,,1)j M =- 解一些具有三对角系数矩阵的方程组；再由2 (3.6.1)求出

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程序)

二阶常微分方程的中心差分求解 学校:中国石油大学（华东）理学院 姓名：张道德 一、 实验目的 1、 构造二阶常微分边值问题: 22,(),(), d u Lu qu f a x b dx u a u b αβ?=-+=<

11122 222222333222122112 100121012010012 00N N N u f q h h u f q h h h u f q h h h q u f h h ---???? ??+-???? ??? ???? ???????-+-? ?????? ???????????=-+? ?????? ???????????-???? ????????-+????? ?? ????? 可以看出系数矩阵为三对角矩阵，而对于系数矩阵为三对角矩阵的方程组可以用“追赶法”求解，则可以得出二阶常微分方程问题的数值解。 四、 举例求解 我们选取的二阶常微分方程边值问题为： 2 22242,01 (0)1,(1), x d u Lu x u e x dx u u e ?=-+=-<

MATLAB期末上机试题带答案

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实习心得体会 我从中国石油大学（北京）来到大庆油田勘探开发研究院实习，刚刚来到这里时我内心充满了激动和喜悦，对这个有着辉煌成绩的科研单位充满了期待。到了这里没多久我就参加培训中心组织的入院教育，对研究院的辉煌过去有了更深刻的了解，那时我就下定决心在研究院好好实习，多学习一些现场经验，为以后的工作打下良好的基础。 一个月的入院培训结束后，经过自己申请和中心安排我被分到了地震资料处理一室，实习期间地震一室的副主任陈志德是我的现场导师，陈主任有着丰富的理论知识和现场经验，能够成为他的学生我感到非常的高兴。最初他把我分到了高分辨率处理解释项目组，工作内容主要是跟着项目组成员用landmark做一些解释性的工作，期间我跟着张晶师傅学习到了如何加载地震数据和测井数据；如何应用测井数据做了合成地震记录，另外我还做了一些地层解释方面的工作。landmark实习大概进行了三个多月，从中我学习到了地震解释的一些基本流程，虽然以后从事的是地震资料处理方面的工作，但有解释方面的实习经历对于提高自己的能力也是很有好处的。 十二月份研究院举行了硕士研究生开题汇报，现场导师陈志德针对大庆地区扶扬油层中浅层高频噪声干扰严重的现状给我定出了小波阈值去噪的论文题目，并且小波变换也是我学校导师的研究方向。他根据我的研究方向给我的论文制定了如下三个目标： 1. 小波分解压制单频噪音

应用matlab设计一个含有随机噪音的一维正演模型,采用小波分解的方法压制噪音，最大程度的保留有效信号。 2. 小波分解压制随机噪音 应用matlab 设计一个含有随机噪音的二维正演模型，并利用小波变换阈值方法压制噪音，然后对软、硬阈值去噪方法进行比较，最后得出结论。 3. 对实际地震记录的噪音压制 在对理论模型进行噪音压制的基础上，对实际的二维数据进行了噪音压制，对比分析噪音压制前后的结果。 经过这一年半的实习我基本上达到了上述目标，完成了论文的预期效果。论文中编制的程序加以优化后就可以在实际生产上进行应用。 下面我简单谈一下我论文的具体实现过程： 我首先做了两个正演模型：一个为楔状砂体正演模型，另外一个为砂泥岩透镜体正演模型。模型的生成均是以褶积理论为基础的，用标准雷克子波与给定反射系数进行褶积，最后得到合成记录。合成记录生成后我又在其上加上了均值为零的随机分布的高频白噪声，如何去掉这些噪声就是本次论文的重点。我采用的是小波阈值去噪的方法，该方法的理论基础为：属于besov空间的信号在小波域内其能量主要集中在有限的几个系数中，而噪声的能量却分布于整个小波小波域内，因此经过小波变换后，信号的小波变换系数大于噪声的小波变换系数，于是可以找到一个合适的数λ作为阈值。当小波变换系数w小于该阈值时，认为这是的w主要时由噪声引起的；当w大于该阈值时，认为这是的w主要是由信号引起的，从而实现信噪分离。在论文研究过程中对于小波基的选择、尺度参数的选择

MATLAB期末考试试卷及答案(20200412004949)

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)clf

三相异步电动机Matlab仿真

中国石油大学胜利学院 综合课程设计总结报告 题目：三相异步电机直接启动特性实验模型 学生姓名：潘伟鹏 系别：机械与电气工程系 专业年级：2012 级电气工程专业专升本2 班指导教师：王铭 2013 年 6 月27 日

一、设计任务与要求 普通异步电动机直接起动电流达到额定电流的 6--7 倍，起动转矩能达到额定转矩的 1.25 倍以上。过高的温度、 过快的加热速度、 过大的温度梯度和电磁力， 产生了极大的破坏力， 缩短了定子线圈和转子铜条的使用寿命。 但在电网条件和工艺条件允许的情况下，异步电动机也可以直接启动。本次课程设计通过 MATLAB 软件建模 模拟三相异步电动机直接启动时的各个元器件上的电量变化。 电力系统 matlab 仿真类书籍 电机类教材 二、方案设计与论证 三相异步电动机直接起动就是利用开关或接触器将电动机的定子绕组直接接到具有额定电压的电网上。 由《电机学》知三相异步电动机的电磁转矩 M 与直流电动机的电磁转矩 有相似的表达形式。它 们都与电机结构(表现为转矩常数)和每级下磁通有关，只不过在三相异步电动机中不再是通过电枢的全部 电流，而是点数电流的有功分量。三相异步电机电磁转矩的表达式为： 1-1) 式中 ——转矩常数 ——转子功率因数 式(1-1)表明，转子通入电流后，与气隙磁场相互作用产生电磁力，因此，反映了电机中电流、磁场和 作用力 之间符合左手定则的物理关系，故称为机械特性的物理表达式。该表达式在分析电磁转矩与磁通、电流之间 的关系时非常方便。 从三相异步电动机的转子等值电路可知， (1-3) 将式(1-2)、 (1-3)代入 (1-1) 得： 每级下磁通