哈工大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语言功能时使用。

5.什么是全局变量和局部变量？

全局变量：作用域是在整个MATLAB工作空间,即全程有效，所有函数都可以对其进行存取和修改，定义全局变量是函数之间传递信息的手段。

局部变量：函数文件变量一般是局部变量，与其他函数和matlab工作空间是隔离的。即在一个函数文件中定义的变量不能被另一个函数文件引用。如果在若干函数中，把某个变量定义为全局变量，那么这些函数可以共用这个变量。

6.数据插值和曲线拟合区别？

曲线拟合研究的是如何寻找平滑曲线以最好的表现带噪声的测量数据，但不要求拟合曲线穿过这些测量数据点；插值是研究如何平滑的估算出基准数据之间的其他点的函数值，所以插值所得曲线一定要穿过基准数据。

7.数值运算与符号运算的区别？

数值运算：

1. 其运算的元素是数值；

2. 在运算过程中必须先对变量进行赋值，然后才能参加运算；

3. 其结果以数值形式出现。

符号运算：

1. 其运算的元素是符号；

2. 在运算过程中无须对变量进行赋值就可参加运算；但是必须先定义符号变量；

3. 其结果以符号形式出现。

8.GUI开发环境中提供了哪些方便的工具？各有什么用途？

布局编辑器(Layout Edtor)——在图形窗口中创建及布置图形对象；

几何排列工具(Alignment Tool)——调整各对象之间的几何关系和位置；

属性查看器(Property Inspector)——查询并设置对象的属性值;

对象浏览器(Object Browser)——用于获得当前Matlab窗口中图形对象句柄的分级排列；

菜单编辑器(Menu Editor)——建立和编辑主菜单和图形对象鼠标右键菜单Tab顺序编辑器（Tab Order Editor）——用于设置当用户按下键盘上的Tab

键时，对象被选中的先后顺序。

9.简述Simulink建立子系统的两种方法及其子系统封装的4个选项

1、使用Subsystems模块库中的Subsystem模块建立子系统，首先构成系统的整体模型，然后编辑空的子系统内的模块。（注意，对于多输入与多输出子系统而言，需要使用Sources模块库中的In1输入虚模块与Sinks模块库中的Out1输出虚模块来实现。）

2、首先选择能够完成一定功能的一组模块，然后选择Simulink模型创建编辑器中Edit菜单下的Create Subsystem，即可建立子系统并将这些模块封装（到此子系统中，Simulink自动生成子系统的输入与输出端口。）

3、Icon项Parameters Initialization标签页Documentation标签页

10.分别采用函数fix、floor、ceil和round对A=[ -1.95,-1.25,1.45,1.85]取整，并

说明这些取整函数的区别。

fix函数：结尾取整

Ceil函数：向上取整

Floor函数：向下取整

Round函数：四舍五入

11.求A=[1,2,3;4,5,6;7,8,9]矩阵的秩，迹，特征值和特征向量、范数和方阵的条

件数；并求A和a= magic(5)矩阵行列式，逆矩阵和伪逆矩阵，并说明逆矩阵

和伪逆矩阵区别。

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

rank=rank(A)

trace=trace(A)

[V,D]=eig(A)

norm=norm(A)

cond=cond(A)

结果：rank = 2

trace = 15

V =

-0.2320 -0.7858 0.4082 -0.5253 -0.0868 -0.8165 -0.8187 0.6123 0.4082 D =

16.1168 0 0 0 -1.1168 0 0 0 -0.0000 norm =16.8481

cond =3.8131e+16

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

a=magic(5)

detA=det(A)

invA=inv(A)

pinvA=pinv(A)

deta=det(a)

inva=inv(a)

pinva=pinv(a)

如果矩阵A不是一个方阵，或者A是一个非满秩的方阵时，矩阵A没有逆矩阵，但可以找到一个与A的转置矩阵A‘同型的矩阵B，使得：A·B·A=A B·A·B=B 此时称矩阵B为矩阵A的伪逆，也称为广义逆矩阵。

12.已知：a=6; A=pascal(4); B=randn(4);求A*B， A.*B, B/A，B./A, B\A,

B.\A ,A.^a，A^a,并且说明数组和矩阵在上述运算中的区别。

a=6; A=pascal(4); B=randn(4);

A*B

A.*B

B/A

B./A

B\A

B.\A

A.^a

A^a

ans =

-3.3843 -0.7710 0.2884 -0.4321

-5.7672 -3.9198 1.6239 1.0912

-6.7793 -10.8330 4.2172 4.3328

-4.9821 -23.2219 8.3812 9.9203

ans =

-1.0689 0.3252 -0.1022 -0.8649

-0.8095 -1.5099 -0.7243 -0.1202

-2.9443 4.1109 1.9152 -1.6488

1.4384 -6.8461 3.1286 1

2.5541

ans =

-5.7707 9.4959 -6.2804 1.4863

0.3558 -3.1462 2.7419 -0.7610

-18.5572 32.8440 -23.1637 5.9327

16.6463 -34.1498 25.8257 -6.8838

ans =

-1.0689 0.3252 -0.1022 -0.8649

-0.8095 -0.3775 -0.0805 -0.0075

-2.9443 0.4568 0.0532 -0.0165

1.4384 -0.4279 0.0313 0.0314

ans =

-0.6348 -1.6245 -2.9494 -4.6346

-1.1688 -2.9031 -5.9575 -10.7762

1.7679 6.3637 16.2994 33.0876

-1.0201 -0.9924 -1.6780 -3.3918

ans =

-0.9356 3.0751 -9.7807 -1.1562

-1.2353 -2.6493 -12.4251 -133.1057

-0.3396 2.1893 18.7966 -60.6505

0.6952 -2.3371 31.9633 31.8620

ans =

1 1 1 1

1 64 729 4096

1 729 46656 1000000

1 4096 1000000 64000000

ans =

1200089 4011200 9133680 17222240

4011200 13407257 30529056 57564976

9133680 30529056 69516457 131078976

17222240 57564976 131078976 247160297

数组运算是指数组对应元素之间的运算,也称点运算.矩阵的乘法、乘方和除法有特殊的数学含义，并不是数组对应元素的运算，所以数组乘法、乘方和除法的运算符前特别加了一个点。

(1)对于乘法、乘方和除法等三种运算，矩阵运算与数组运算的运算符及含义都不同：矩阵运算按线性变换定义，使用通常符号；数组运算按对应元素运算定义，使用点运算符；

(2)数与矩阵加减、矩阵除法在数学是没有意义的，在MATLAB中为简便起见，定义了这两类运算

13.分别采用不循环语句和循环语句For和while计算级

数:S=1+6+62+63+···+633，并且说明循环语句For和while之间的区别。

While 语句

s=0;

i=0;

while i<34

s=s+6^i;

i=i+1;

end

S

For 语句：

s=0;

for i=(0:1:33)

s=s+6^i;

end

S

不循环语句：

i=(0:1:33);

sum(6.^i)

for 和while 的区别：while 循环语句与for 循环语句不同的是，前者是以条件的满足与否来判断循环是否结束的，而后者则是以执行次数是否达到指定值来判断的。

14. 采用M 函数文件编写特殊矩阵函数Hilbter ，Hilbter 矩阵特点是每个元素

1(,)1

A i j i j =+-，并且调用matlab 中的函数Hilb ，进行对比。 A=zeros(6,6)

A =

0 0 0 0 0 0

0 0 0 0 0 0

0 0 0 0 0 0

0 0 0 0 0 0

0 0 0 0 0 0

0 0 0 0 0 0

>> for i=1:1:6

for j=1:1:6

A(i,j)=1./(i+j-1);

end

end

>> A

A =

1.0000 0.5000 0.3333 0.2500 0.2000 0.1667

0.5000 0.3333 0.2500 0.2000 0.1667 0.1429

0.3333 0.2500 0.2000 0.1667 0.1429 0.1250

0.2500 0.2000 0.1667 0.1429 0.1250 0.1111

0.2000 0.1667 0.1429 0.1250 0.1111 0.1000

0.1667 0.1429 0.1250 0.1111 0.1000 0.0909

>> hilb(6)

ans =

1.0000 0.5000 0.3333 0.2500 0.2000 0.1667

0.5000 0.3333 0.2500 0.2000 0.1667 0.1429

0.3333 0.2500 0.2000 0.1667 0.1429 0.1250

0.2500 0.2000 0.1667 0.1429 0.1250 0.1111

0.2000 0.1667 0.1429 0.1250 0.1111 0.1000

0.1667 0.1429 0.1250 0.1111 0.1000 0.0909

两者相同。

15.分别求多项式a(x)=x4+8x3-10和多项式b(x)= 2x2-x+3的乘以和除以结果，同时

求这两个多项式乘积的导函数和两个多项式除法的导函数。

syms x

p1=[1 8 0 0 -10];

p2=[2 -1 3];

w = conv(p1,p2);

k = polyder(p1,p2);

[Q,r] = deconv(p1,p2);

[p,q] = polyder(p1,p2);

a=poly2str(p1,'x')

b=poly2str(p2,'x')

ab=poly2str(w,'x')

dab=poly2str(k,'x')

e=poly2str(Q,'x')

de=poly2str(p,'x')

g=poly2str(r,'x')

dg=poly2str(q,'x')

16.已知多项式P(x)=x4+8x3，分别用polyval和polyvalm函数计算这多项式在

x=[-1,1.2;2,-1.8]处的值，它们区别是什么？

p=[1,8,0,0,0];

x=[-1,1.2;2,-1.8];

polyval(p,x)

polyvalm(p,x)

x可以是一个矩阵或者一个向量,在这两种情况下，polyval计算在X中任意元素处的多项式p

的估值，计算在x=-1、1.2、2、-1.8的值，结果以矩阵的形式表示出来。而polyvalm 运算时把x=A矩阵直接带入多项式进行运算，进行矩阵运算，结果也是一个矩阵。

17.用2次和6次多项式拟合在[0,2]区间的函数y=cos(x) *exp(x/20)，并将原始

数据和拟合结果绘制在一个图中,并且进行图形标注。

x=0:0.1:2;

y=cos(x) .*exp(x/20);

p1=polyfit(x,y,2);

y1=polyval(p1,x);

p2=polyfit(x,y,6);

y2=polyval(p2,x);

hold on;

plot(x,y,'or--');

plot(x,y1,'b--');

plot(x,y2,'g*');

xlabel('x');

ylabel('y');

legend('原始数据','二次多项式拟合','六次多项式拟合')

18.有一正弦衰减数据y=sin(x).*exp(-x/10)，其中x=0:pi/5:4*pi，用三次样条法进

行插。

x=0:pi/5:4*pi;

y=sin(x).*exp(-x/10);

xi=0:0.1:4*pi;

yi_spline=interp1(x,y,xi,'spline');

plot(x,y,'ro',xi,yi_spline,'b-');

title('三次样条插值');

19.求解恰定方程组，采用左除法、求逆法、lu 和qr 分解法求数值解

2x+y-3z+w=11

x-2y+7w=-7

2y+z-w=5

x+4y-z-2w=0

A=[2,1,-3,1;1,-2,0,7;0,2,1,-1;1,4,-1,-2];

B=[11;-7;5;0];

x=A\B

x=pinv(A)*B

[L,U]=lu(A);

x=U\(L\B)

[Q,R]=qr(A);

x=R\(Q\B)

20. 设1,1t

y t e ππ-=-≤≤+在同一图形窗口采用子图的形式绘制不同图形，条形图、阶梯图、杆图和对数坐标图。

x=-pi:0.1:pi;

y=1./(1+exp(-x));

subplot(2,2,1);bar(x,y,'g');

title('bar(x,y,''g'')');axis([0,7,-2,2]);

subplot(2,2,2);stairs(x,y,'b');

title('stairs(x,y,''b'')');axis([0,7,-2,2]);

subplot(2,2,3);stem(x,y,'k');

title('stem(x,y,''k'')');axis([0,7,-2,2]);

subplot(2,2,4)

loglog(x,y)

二、 选答题（在下列题中选答4题）：

1..符号运算（必选两道）

(1) 求微分方程222/()/(2)dy dt t y t =+的通解

syms t y ;

dsolve('Dy = (t^2+y^2)/2/t^2')

ans =

-t*(1/(C4 + log(t)/2) - 1)

(2) 求代数方程1122

a x x +=+-的解 [x]=solve('1/(x+2)+a=1/(x-2)')

x =

(2*(a*(a + 1))^(1/2))/a

-(2*(a*(a + 1))^(1/2))/a

2. 求解方程的数值解（必选两道）

(1) 0-0x =0.5=x x sin(x)-10在附近的根

function yp =text2(x)

%UNTITLED5 Summary of this function goes here

% Detailed explanation goes here

yp = x-sin(x)-10^x;

end

z = fzero('text2',0.5)

z =

NaN

(2) 222x -0.3c o s (y )=0

3y -0.6c o s (x )=0非线性方程组在(0.5,0.5) 附近的数值解

function q=myfun1(p)

x=p(1);

y=p(2);

q(1)=2*x^2-0.3*cos(y);

q(2)=3*y^2-0.6*cos(x);

end

x=fsolve('myfun1',[0.5,0.5]',optimset('Display','off'))

x =

0.3691

0.4319

(3) 常微分方程数值解0f 0y'+(1.2+sin10t)y=0,t 0,t 5,(t )1y ===

(4) 采用牛顿-科特斯法数值求解定积分20sin 1cos x x I dxdy x

π

=+? 3. 采用GUI 设计工具设计用户界面，要求设计一个带有两个按钮和一个坐标轴的图形用户界面，当用鼠标点击两个按钮时，分别在坐标轴内画cylinder 和peaks 两个图形。 4. 2222sin()

x y z x y +=?+绘制，绘图区域[-8,8][-8,8]，分别绘制它的网络图、带等高

线的空间曲面、屏蔽边界面的空间曲面和表面图

[x,y]=meshgrid([-8:.5:8]);

z=sin(x.^2+y.^2)./(x.^2+y.^2);

mesh(x,y,z)

subplot(2,2,1), plot3(x,y,z)

title('plot3 (x,y,z)')

subplot(2,2,2), mesh(x,y,z)

title('mesh (x,y,z)')

subplot(2,2,3), surf(x,y,z)

title('surf (x,y,z)')

subplot(2,2,4), surf(x,y,z), shading interp

title('surf (x,y,z), shading interp')

x x x u t

++=，其中u(t)是单位阶跃函数，建立系统并仿真，并用scope模块观察结果。

6. 编程：某商场对顾客所购买的商品实行打折销售，标准如下(商品价格用price来表示)：

① price<200 没有折扣

② 200≤price<500 3%折扣

③ 500≤price<1000 5%折扣

④ 1000≤price<2500 8%折扣

⑤ 2500≤price<5000 10%折扣

⑥ 5000≤price 14%折扣

输入所售商品的价格，求其实际销售价格。

price=input('请输入商品价格');

a=price;

switch fix(a/100)

case {0,1}

rste=0;

case {2,3,4}

rate=0.03;

case num2cell(5:9)

rate=0.05;

case num2cell(10:24)

rate=0.08;

case num2cell(25:49)

rate=0.1;

otherwise

rate=0.14; end

a=a.*(1-rate)

MATLAB实验题答案

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 仿真 实验报告

基于MATLAB的电机学计算机辅助分析与仿真 实验报告 班级： 学号： 姓名： 完成时间：

一、实验内容 1.1单相变压器不同负载性质的相量图 通过MATLAB 画出单相变压器带感性，阻性，容性三种不同性质负载的变压器向量图 1.2感应电机的S T -曲线 通过MATLAB 画出三相感应电动机的转矩转差率曲线 二、实验要求 2.1单相变压器不同负载性质的相量图 根据给定的仿真实例画出负载相位角30,0,302-=j 三种情况下得向量图，观察电压大小与相位的关系，了解总结负载性质不同对向量图的影响 2.2感应电机的S T -曲线 根据给定的实例，画出3.1~3.1-=s 的S T -曲线，了解感应电机临界转差率的大小和稳定工作区间的大小，给出定性分析 三、实验方法 3.1单相变压器不同负载性质的相量图 1.单相变压器不同负载性质的相量图 (1)先画出负载电压'2U 的相量； (2)根据负载的性质和阻抗角画出二次电流（规算值）的相量 (3)在2U 上加上一个与电流方向相同的压降，其大小为二次电流规算值'2I 与二次漏电阻规算值'2R 之积；再加上一个超前电流方向?90的压降，其大小为二次电流'2I 规算值与二次漏电抗规算值'2χ之积； (4)根据上一步结果连线，得出'2E ； (5)超前'2E 方向?90画出m Φ； (6)根据励磁电阻与电抗的大小得出励磁阻抗角，并超前m Φ一个励磁阻抗角的大小得出m I 的方向； (7)根据平行四边形法则，做出'2I -与m I 的和，即为1I ； (8)根据'21E E =得出1E ，并得出1E -。

(9)在1E -上加上一个与电流方向相同的压降，其大小为一次电流1I 与一次漏电阻1R 之积；再加上一个超前电流方向?90的压降，其大小为一次电流1I 与一次漏电抗1χ之积； (10) 根据上一步结果连线，得出1U ； 3.2感应电机的S T -曲线 实验采用matlab 对转矩转差率曲线进行仿真。 由转矩转差率关系公式知， 2212 2122 1)()(x c x s r c r s r U m T s s +++?Ω= 只有s 为自变量，其他参数均为已知。 编程时，先取s 在0.01-1.3正区间的S T -，进行绘图；再取相应负区间对S T -绘图；最后加入（0,0） 四、实验源程序(1分) 4.1单相变压器不同负载性质的相量图 见附录 4.2感应电机的T-S 曲线 %T-S 曲线绘制 %定义常量 R2 = 0.04; R1 = 0.06; M1 = 3; U1 = 380; W = 2*pi*1485/60; X1 = 0.27; X2 = 0.56; C = 1+X1/16.4; %画出s=0.01~1.3的T-S 曲线 s = 0.01:0.01:1.3; T=ones(1,length(s));

《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实验题答案

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

(完整word版)哈工大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）

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期末考试题试题及答案(95分)分解,推荐文档

建议收藏下载本文，以便随时学习！ 春季学期MATLAB期末作业 学院：机电工程学院 专业：机械制造设计及其自动化 学号： 班号： 姓名： 我去人也就有人！为UR扼腕入站内信不存在向你偶同意调剖沙

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）建议收藏下载本文，以便随时学习！我去人也就有人！为UR扼腕入站内信不存在向你偶同意调剖沙

MATLAB实验题目及答案

实验二一维二维数组的创建和寻访 一、实验目的 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

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实验指导书(含答案)汇总

实验一：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编程

H a r b i n I n s t i t u t e o f T e c h n o l o g y 材料力学上机作业 课程名称：材料力学 设计题目：应力状态分析 院系：机电学院 班级： 分析者： 学号： 指导教师：张桂莲 设计时间：2013年6月18日 哈尔滨工业大学

材料力学上机课设计说明书 一， 设计题目 题目7 应力状态分析 输入： 1. 平面应力状态输入： x y xy σστ（,,）；某截面方位角α 2. 空间应力状态输入： ,x y z xy yz zx σσστττ（,,,,） 输出： 1. 输出主应力123σσσ（,,） 2. 最大切应力（13 max 132 σσττ-== ） 3. 如为平面应力状态则需要输出方位角α斜截面上的应力αα στ、及主方向角*σα 4. 画出应力圆示意图 二， 程序计算设计过程 1. 平面应力状态分析 对于任意平面应力状态，有 max min σσ = 2x y σσ+± 主应力为： 1max 23min ,0,σσσσσ=== 并且由 2tan 2xy x y στασσ=- 可求得主应力方向角 1 3 σσ αα、。 对于任意一个方位角α，有： = cos 2sin 222 sin 2cos 22 x y x y xy x y xy αασσσσσατα σστατα +++ +-=- +

从而，输入任意角α，即可求得该截面的应力状态 ααστ、 并且 ααστ、都是关于α的函数，上式即为应力圆的参数方程，参数为α。 将α从0到pi 取一系列的值，则可以求出一系列的ααστ、，在坐标系中找 到对应点，连接即可作出应力圆。 2. 三向应力状态分析 解特征方程 32 1230I I I σσσ-+-= 即可求出主应力123σσσ、、 其中： 123|||||||| x y z x yx y zy z xz xy y yz z zx x x yx zx xy y zy xz yz z I I I σσσστστσττστστσστττστττσ=-+??????=++ ? ? ? ???????? ?= ? ??? 再由 13 max 132 σσττ-== 可求得最大切应力。 求解三向应力圆： 三个圆121323C 、C 、C 的圆心分别为： 231312122313,0,0,0222C C C σσσσσσ+++?????? ? ? ??????? 、、 半径非别为： 23 13 12 122313r = ,r = ,r = 2 2 2 σσσσσσ--- 由此可以求出三个应力圆的方程，从而作出三向应力圆。 三， 程序代码 reg=input('选择应力状态方式(1或2)：');%1表示平面应力状态，2表示空间应力状态 if reg==1 %选择平面应力状态分析 %输入已知量，应力单位为MPa ，转角单位为rad

MATLAB实验上机易错题汇总

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

哈工大_Matlab__2013年春季学期《MATLAB语言及应用》试题答案

2013年春季学期 《MATLAB语言及应用》课程试卷 姓名： 学号： 学院： 专业： 必答题 1.常用的matlab界面由哪些窗口组成，各有什么主要作用？(4分) (1)菜单和工具栏功能：【File】菜单主要用于对文件的处理。【Edit】菜单主 要用于复制、粘贴等操作，与一般Windows程序的类似，在此不作详细介绍。【Debug】菜单用于调试程序。【Desktop】菜单用于设置主窗口中需要打开的窗口。【Window】菜单列出当前所有打开的窗口。【Help】菜单用于选择打开不同的帮助系统。 (2)命令窗口功能：用于输入命令并显示除图形以外的所有执行结果 (3)历史命令窗口功能：主要用于记录所有执行过的命令 (4)当前工作目录窗口功能：对文件和目录进行操作 (5)工作空间窗口功能：查看、载入和保存变量 2.如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么 区别？(2分) 方法一：在MATLAB命令窗口中输入editpath或pathtool命令或通过【File】/|【SetPath】菜单，进入“设置搜索路径”对话框，通过该对话框编辑搜索路径。 方法二：在命令窗口执行“path(path,…D:\Study ?)”，然后通过“设置搜索路径”对话查看“D:\Study”是否在搜索路径中。 方法三：在命令窗口执行“addpath D:\Study- end”，将新的目录加到整个搜索路径的末尾。如果将end改为begin，可以将新的目录加到整个搜索路径的开始。 区别：当前文件目录是正在运行的文件的目录，显示文件及文件夹的详细信息，且只有将文件设置为当前目录才能直接调用。搜索路径中的文件可以来自多个 不同目录，在调用时不用将其都设置为当前目录，为同时调用多个文件提供 方便。 3.有几种建立矩阵的方法？各有什么优点？(4分) 1.在命令窗口中直接输入优点：适合输入不规则和较小的矩阵 2.通过语句和函数生成矩阵优点：适合输入规则且较大的矩阵 3.通过M文件来建立矩阵，从外部数据文件中导入矩阵 优点：方便创建和导入大型矩阵 4.说明break语句、continue语句和return语句的用法。(3分)

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

MATLAB实验题答案

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)???????=-++=--=-++=++5 65333328 21w 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

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期末上机试题带答案

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选修课试题

2015年春季学期 《MATLAB语言及应用》课程试卷 姓名：李清蔚 学号：11408010304 学院：机电学院 专业：机械设计制造及其自动化 (1)请概括Matlab中冒号(:)，逗号(,)，分号(;)的功能？(6分) ：用于生成一位数组，表示一位数组的全部元素或多维数组的某一维的全部元素 ，用于要显示计算结果的命令之间的分隔符，用于输入变量之间的分隔符，用于数组行元素之间的分隔符。 ;用于不显示计算结果命令行的结尾，用于不显示计算结果命令之间的分隔符，用于数组元素行之间的分隔符。 (2)a=[1 -6 0 0.45 0 100]，b=[0 1.5 4 0 0.5 -34]，请问在进行逻辑运算时，a、b分别相当于什么样的逻辑量？a与b的逻辑与、逻辑或都分别是什么？(4分) 元素方式的逻辑变量a&b=[0 1 0 0 0 1], a|b=[1,1,1.1,1] (3)有几种建立矩阵的方法？各有什么优点？(4分) 直接输入矩阵，eye()生成单位矩阵，通过M文件建立矩阵，通过外部数据文件导入矩阵 (4)生成一个7阶的魔方矩阵A，再采用相应的MATLAB命令，将其全部奇数行提取出来，赋给矩阵B，给出执行语句以及计算结果。(4分) A=magic(7) table=[1 3 5 7]; [n,len]=size(table); B=A(table(1),:); for i=2:len B=[B;A(table(i),:)]; End

B= 30 39 48 1 10 19 28 46 6 8 17 26 35 37 13 15 24 33 42 44 4 22 31 40 49 2 11 20 (5)生成5阶魔方矩阵A，再分别计算A中元素的平方、A矩阵的平方，给出执行语句以及计算结果。(5分) A=magic(5) B=A.^2 C=A^2 A = 17 24 1 8 15 23 5 7 14 16 4 6 13 20 22 10 12 19 21 3 11 18 25 2 9 B = 289 576 1 64 225 529 25 49 196 256 16 36 169 400 484 100 144 361 441 9 121 324 625 4 81 C = 1090 900 725 690 820 850 1075 815 720 765 700 840 1145 840 700 765 720 815 1075 850 820 690 725 900 1090 (6)生成矩阵 8412 6115 689 a ?? - ?? =-- ?? ?? - ??，计算a的行列式、逆矩阵，给出执行语句以 及计算结果。(5分) b = 784 c = 0.0753 0.0765 0.1429 0.0306 0.1837 0.1429 -0.0230 0.1122 0.1429