2013年春季学期
MATLAB 课程考查题
一、必答题:
1.matlab常见的数据类型有哪些?各有什么特点?
(1)numeric(数值):包括整数和浮点数,缺省状态下matlab将所有的数看
作是双精度的浮点数。
(2)char(字符):字符串必须用单引号括起来,是由多个字符组成的字符行
向量。
(3)logical (逻辑):仅有TRUE和FALSE两个值,一般用于关系运算或逻辑
运算。
(4)cell (单元、细胞、元胞):可以把不同类型的数据纳入到一个变量中,不
要求相同的数据类型,下标用大括号括起来。
(5)struct (结构):结构是包含已命名“数据容器”或字段的数组。结构中的
字段可以包含任何数据。
2.MATLAB中有几种帮助的途径?
(1)帮助浏览器:选择view菜单中的Help菜单项或选择Help菜单中的
MATLAB Help菜单项可以打开帮助浏览器;
(2)help命令:在命令窗口键入“help”命令可以列出帮助主题,键入“help
函数名”可以得到指定函数的在线帮助信息;
(3)lookfor命令:在命令窗口键入“lookfor 关键词”可以搜索出一系列与给
定关键词相关的命令和函数
(4)模糊查询:输入命令的前几个字母,然后按Tab键,就可以列出所有以
这几个字母开始的命令和函数。
3.Matlab常见的哪三种程序控制结构及包括的相应的语句?
(1)顺序结构:数据的输入-A=input(提示信息,选项);
数据的输出-disp(输出项)
程序的暂停-pause(延迟秒数)
(2)选择结构:if 条件语句组;else 语句组;end
Switch 表达式case 值1 语句组1;……otherwise 语句组n;
end
(3)循环结构:while 条件语句组;end
for(条件);end
do 语句组while 循环语句组;end
4.命令文件与函数文件的主要区别是什么?
答:(1)命令文件是一系列命令的组合,函数文件的第一行必须用function 说明;
(2)命令文件没有输入参数,也不用返回参数,函数文件可以接受输入参数,也可以返回参数;
(3)命令文件处理的变量为工作空间变量,函数文件处理的变量为函数内部的局部变量,也可以处理全局变量。
5.什么是全局变量和局部变量?
答:全局变量是没有定义存储类型的外部变量,其作用域是从定义点到程序结束,省略了存储类型符,系统将默认为是自动型
局部变量也称为内部变量。局部变量是在函数内作定义说明的。其作用域仅限于函数内,离开该函数后再使用这种变量是非法的。
如果一个函数内的变量没有特别声明,那么这个变量只在函数内部使用,即为局部变量,如果两个或多个函数公用一个变量(或者说在子程序中也要用到主程序中的变量,注意不是参数),那么可以用global来将它声明为全局变量。
全局变量的使用可以减少参数传递,合理利用全局变量可以提高程序执行效率。如果需要用到其他函数的变量,就要利用在主程序与子程序中分别声明全局变量的方式实现变量的传递,否则函数体内使用的都为局部变量。
6.数据插值和曲线拟合区别?
答:有这样一组数据(Xi,Yi),i=1,2,3,。。。。。,寻找Y与X之间的关系Y=F(X),就要根据这组数据求出其近似关系Y=G(X)。
插值与曲线拟合实际上是求这个近似关系的两种方法,不同的是,插值方法
求出的插值多项式要求所有的数据点(Xi,Yi)都在曲线上,而拟合求出的插值多项式只要反应数据的基本趋势就可以了,并不要求所有的数据点都在拟合曲线上。
但是,利用插值方法更多的依赖于插值点,比如选择插值点前后N个数值。。如果插值点很多,就不好利用插值法了,最好用曲线拟合,这样求出的曲线具有一定的函数关系,只要把要求的数据点代入,就可以得出结果。
7.数值运算与符号运算的区别?
答:数值运算:求出具体的数值,不含符号。(如解方程,求出未知数=1.565656 ,不是未知数=ab+c)数值计算是MATLAB最基本、最重要的功能,是MATLAB最具代表性的特点。MATLAB在数值计算过程中以数组和矩阵为基础。数组是MATLAB运算中的重要数据组织形式。
符号运算:带有符号变量、表达式的运算称为抽象计算即符号计算,MATLAB 借助于MAPLE符号数学特长,建立并发展了强大的符号运算功能。结果用符号表示。
8.GUI开发环境中提供了哪些方便的工具?各有什么用途?
答:开发工具及用途:
布局编辑器(Layout Edtor)-------在图形窗口中加入及安排对象。布局编辑器是可以启动用户界面的控制面板,上述工具都必须从布局编辑器中访问,用guide命令可以启动,或在启动平台窗口中选择GUIDE来启动布局编辑器。几何排列工具(Alignment Tool)-----调整各对象相互之间的几何关系和位置。
属性编辑器 (Property Inspector)-----查询并设置属性值。
对象浏览器Object Browser)-----用于获得当前MATLAB图形用户界面程序中所有的全部对象信息,对象的类型,同时显示控件的名称和标识,在控件上双击鼠标可以打开该控件的属性编辑器。
菜单编辑器(Menu Editor)-----建立窗口菜单条的菜单和任何构成布局的弹出菜单。
9.简述Simulink建立子系统的两种方法及其子系统封装的4个选项
答:建立子系统有两种方法:通过Subsystem模块建立子系统和通过已有的模块建立子系统。两者的区别是:前者先建立子系统,再为其添加功能模块;后者先选择模块,再建立子系统。
子系统的封装主要就是对这4页参数进行设置:Icon、Parameters、Initialization和Documentation。
10.分别采用函数fix、floor、ceil和round对A=[ -1.95,-1.25,1.45,1.85]取整,并
说明这些取整函数的区别。
答:A=[ -1.95,-1.25,1.45,1.85];
B=fix(A)
B =[-1 -1 1 1]
Fix函数为取截尾取整。
C=floor(A)
C =[ -2 -2 1 1]
Floor函数为取不超过x 的最大整数(高斯取整)。
D=ceil(A)
D = [ -1 -1 2 2]
Ceil函数为取大于x的最小整数。
E=round(A)
E = [-2 -1 1 2]
Round函数为取四舍五入取整后的整数。
11.求A=[1,2,3;4,5,6;7,8,9]矩阵的秩,迹,特征值和特征向量、范数和方阵的条
件数;并求A和a= magic(5)矩阵行列式,逆矩阵和伪逆矩阵,并说明逆矩阵和伪逆矩阵区别。
答: A的秩:rank(A);ans =2
A的迹:trace(A);ans =15
[V,D]=eig(A)
A的特征向量:
V =
-0.2320 -0.7858 0.4082
-0.5253 -0.0868 -0.8165
-0.8187 0.6123 0.4082
A的特征值为D的对角线元素:
D =
16.1168 0 0
0 -1.1168 0
A的范数:norm(A);ans =16.8481(2阶范数)
方阵的条件数:rcond(A);ans =1.5420e-018
A的矩阵行列式:det(A);ans =6.6613e-016
a的矩阵行列式:a=magic(5); det(a);ans = 5.0700e+006
A的逆矩阵:inv(A)
Warning: Matrix is close to singular or badly scaled.
Results may be inaccurate. RCOND = 1.541976e-018.
ans =
1.0e+016 *
-0.4504 0.9007 -0.4504
0.9007 -1.8014 0.9007
-0.4504 0.9007 -0.4504
A的伪逆矩阵:pinv(A)
ans =
-0.6389 -0.1667 0.3056
-0.0556 0.0000 0.0556
0.5278 0.1667 -0.1944
a的逆矩阵:inv(a)
ans =
-0.0049 0.0512 -0.0354 0.0012 0.0034
0.0431 -0.0373 -0.0046 0.0127 0.0015
-0.0303 0.0031 0.0031 0.0031 0.0364
0.0047 -0.0065 0.0108 0.0435 -0.0370
0.0028 0.0050 0.0415 -0.0450 0.0111
a的伪逆矩阵:pinv(a)
ans =
-0.0049 0.0512 -0.0354 0.0012 0.0034
0.0431 -0.0373 -0.0046 0.0127 0.0015
-0.0303 0.0031 0.0031 0.0031 0.0364
0.0047 -0.0065 0.0108 0.0435 -0.0370
0.0028 0.0050 0.0415 -0.0450 0.0111
同维数的单位阵,就称A为可逆矩阵(或者称A可逆),并称B是A的逆矩阵,简称逆阵。(此时的逆称为凯利逆)
矩阵A可逆的充分必要条件是|A|≠0。
伪逆矩阵是逆矩阵的广义形式。由于奇异矩阵或非方阵的矩阵不存在逆矩阵,但可以用函数pinv(A)求其伪逆矩阵。基本语法为X=pinv(A),X=pinv(A,tol),其中tol为误差,pinv为pseudo-inverse的缩写:max(size(A))*norm(A)*eps。函数返回一个与A的转置矩阵A' 同型的矩阵X,并且满足:AXA=A,XAX=X.此时,称矩阵X为矩阵A的伪逆,也称为广义逆矩阵。pinv(A)具有inv(A)的部分特性,但不与inv(A)完全等同。
如果A为非奇异方阵,pinv(A)=inv(A),但却会耗费大量的计算时间,相比较而言,inv(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=pascal(4)
A =
1 1 1 1
1 2 3 4
1 3 6 10
1 4 10 20
B=randn(4)
B =
-1.0689 0.3252 -0.1022 -0.8649
-0.8095 -0.7549 -0.2414 -0.0301
-2.9443 1.3703 0.3192 -0.1649
1.4384 -1.7115 0.3129 0.6277
A*B
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
A.*B
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
B/A
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
B./A
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
B\A
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
B.\A
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
a=6;A.^a
ans =
1 1 1 1
1 64 729 4096
1 729 46656 1000000
1 4096 1000000 64000000
a=6;A^a
ans =
1200089 4011200 9133680 17222240
4011200 13407257 30529056 57564976
9133680 30529056 69516457 131078976
17222240 57564976 131078976 247160297
上述运算中涉及到点乘和乘的概念,
(1)点乘表示两个矩阵对应位置元素相乘,所以这两个矩阵应该是尺寸等大的。例如:
>> a = [1 2 3;2 3 4];
>> b = [1 2 3;2 3 4];
>> a.*b
ans = [1 4 9 4 9 16]
>> c = [1 2 3;1 2 3;1 2 3];
>> a.*c
??? Error using ==> times
Matrix dimensions must agree.
这里a无法点乘c,因为这两个不是相同尺寸的矩阵。
(2)而乘*,这个和我们在数学中学到的矩阵相乘是一个意义。
例如:>>a*c
ans = [6 12 18 9 18 27]
“.^”和“^”区别实际上就是点乘与乘的区别。因为求幂实际上也是乘法;
(矩阵的点除(./)与点乘(.*)类似.)。
然后还有左除(\)和右除(/)的区别:
A*X=B的解为X=A\B(左除)
X*B=A的解为X=A/B(右除)
13.分别采用不循环语句和循环语句For和while计算级
数:S=1+6+62+63+···+633 ,并且说明循环语句For和while之间的区别。
答:不用循环语句:
a=0:1:33;
b=6;
m=b.^a;
s=sum(m(:))
>>s =5.7302e+025
For循环:
s=0;
for i=0:1:33
s=s+6^i;
end
s
>>s = 5.7302e+025
While循环:
s=0;i=0;
while(i<=33)
s=s+6^i;
i=i+1;
end
disp(s)
>>s=5.7302e+025
For循环,执行过程是依次将矩阵的各列元素赋给循环变量,然后执行循环
体语句,直至各列元素处理完毕。
而while 循环,其执行过程为:若条件成立,则执行循环体语句,执行后再判断条件是否成立,如果不成立则跳出循环。
14. 采用M 函数文件编写特殊矩阵函数Hilbter ,Hilbter 矩阵特点是每个元素
1
(,)1A i j i j =+-,并且调用matlab 中的函数Hilb ,进行对比。
答:Matlab 中的函数Hilb 如下:
function H = hilb(n,classname)
%HILB Hilbert matrix.
% HILB(N) is the N by N matrix with elements 1/(i+j-1),
% which is a famous example of a badly conditioned matrix. % See INVHILB for the exact inverse.
%
% HILB(N,CLASSNAME) produces a matrix of class CLASSNAME.
% CLASSNAME must be either 'single' or 'double' (the default). %
% This is also a good example of efficient MATLAB programming % style where conventional FOR or DO loops are replaced by % vectorized statements. This approach is faster, but uses % more storage.
%
% See also INVHILB.
% Copyright 1984-2005 The MathWorks, Inc.
% $Revision: 5.10.4.4 $ $Date: 2005/11/18 14:14:30 $
if nargin < 2
classname = 'double';
else % nargin == 2
if ~strcmpi(classname,'double') && ~strcmpi(classname,'single') error('MATLAB:hilb:notSupportedClass',...
'CLASSNAME must be ''double'' or ''single''.'); end
end
% I, J and E are matrices whose (i,j)-th element
% is i, j and 1 respectively.
J = 1:cast(n,classname);
J = J(ones(n,1),:);
I = J';
E = ones(n,classname);
H = E./(I+J-1);
调用hilb函数,若n=4时,
hilb(4)
ans =
1.0000 0.5000 0.3333 0.2500
0.5000 0.3333 0.2500 0.2000
0.3333 0.2500 0.2000 0.1667
0.2500 0.2000 0.1667 0.1429
编写一个四阶的hilbter矩阵:
function Hilbter4
i=[1,1,1,1;2,2,2,2;3,3,3,3;4,4,4,4];
j=[1,2,3,4;1,2,3,4;1,2,3,4;1,2,3,4];
a=i+j-1;
b=1./a
>>b =
1.0000 0.5000 0.3333 0.2500
0.5000 0.3333 0.2500 0.2000
0.3333 0.2500 0.2000 0.1667
0.2500 0.2000 0.1667 0.1429
和hilb函数一样。
15.分别求多项式a(x)=x4+8x3-10和多项式b(x)= 2x2-x+3的乘以和除以结果,
同时求这两个多项式乘积的导函数和两个多项式除法的导函数。
答:相乘:
P1=[1,8,0,0,-10];P2=[2,-1,3];
Q=conv(P1,P2)
>>Q =[2,15,-5,24,-20,10,-30]
则Q=2x6+15x5-5x4+24x3-20x2+10x-30;
导函数:p=polyder(Q)
>>p =
12 75 -20 72 -40 10
相除:
P1=[1,8,0,0,-10];P2=[2,-1,3];
[Q,r]=deconv(P1,P2)
>>Q =[ 0.5000 4.2500 1.3750]//P1除以P2的商式
r =[0 0 0 -11.3750 -14.1250]// P1除以P2的余式
则相除的多项式为:0.5x2-7.125x-12.75
相除的导函数:
x=[0.5,-7.125,-12.75]
x =
0.5000 -7.1250 -12.7500
p=polyder(x)
>>p =
1.0000 -7.1250
16.已知多项式P(x)=x4+8x3,分别用polyval和polyvalm函数计算这多项式在
x=[-1,1.2;2,-1.8]处的值,它们区别是什么?
答:polyval函数:
p=[1 8 0 0 0];x=[-1,1.2;2,-1.8];
polyval(p,x)
>>ans =
-7.0000 15.8976
80.0000 -36.1584
polyvalm函数:
p=[1 8 0 0 0];x=[-1,1.2;2,-1.8];
polyvalm(p,x)
>>ans =
-50.5840 50.6496
84.4160 -84.3504
两者的区别在于前者是代数多项式求值,而后者是矩阵多项式求值。
17.用2次和6次多项式拟合在[0,2错误!未找到引用源。]区间的函数y=cos(x)
*exp(x/20),并将原始数据和拟合结果绘制在一个图中,并且进行图形标注。答:程序:
x=0:pi/16:2*pi;
y=sin(x).*exp(-x/10);
a=polyfit(x,y,6);
b=polyfit(x,y,2);
[p1,len]=poly2str(a,'x')
[p2,len]=poly2str(b,'x')
p1=a(1)*x.^6+a(2)*x.^5+a(3)*x.^4+a(4)*x.^3+a(5)*x.^2+a(6)*x+a(7);
p2=b(1)*x.^2+b(2)*x+b(3);
plot(x,y,'-',x,p1,'*',x,p2,'+')
title('三次样条插值');
xlabel('Variable X');
ylabel('Variable Y');
求得:
>>p1 =
0.00037099 x^6 - 0.010803 x^5 + 0.10393 x^4 - 0.3593 x^3 + 0.11004
x^2
+ 0.91509 x + 0.0055543
len =73
p2 =
0.0026462 x^2 - 0.2225 x + 0.7347
len = 36
18.有一正弦衰减数据y=sin(x).*exp(-x/10),其中x=0:pi/5:4*pi,用三次样条法进
行插。
答:>> x0=0:pi/5:4*pi;
>> y0=sin(x0).*exp(-x0/10);
>> x=0:pi/20:4*pi;
>> y=spline(x0,y0,x);
>> plot(x0,y0,'or',x,y,'b')
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 =
43.0000
-11.5556
17.6667
-10.4444
求逆法:
a=[2,1,-3,1;1,-2,0,7;0,2,1,-1;1,4,-1,-2]; b=[11,-7,5,0]';
x=inv(a)*b
>>x =
43.0000
-11.5556
17.6667
-10.4444
伪逆法:
a=[2,1,-3,1;1,-2,0,7;0,2,1,-1;1,4,-1,-2]; b=[11,-7,5,0]';
x=pinv(a)*b
>>x =
43.0000
-11.5556
17.6667
-10.4444
Lu分解:
a=[2,1,-3,1;1,-2,0,7;0,2,1,-1;1,4,-1,-2]; b=[11,-7,5,0]';
[L,U]=lu(a)
X=U\(L\b)
>>L =
1.0000 0 0 0 0.5000 -0.7143 1.0000 0 0 0.5714 0.3846 1.0000 0.5000 1.0000 0 0 U =
2.0000 1.0000 -
3.0000 1.0000 0 3.5000 0.5000 -2.5000 0 0 1.8571
4.7143 0 0 0 -1.3846
X =
43.0000
-11.5556
17.6667
-10.4444
qr分解:
a=[2,1,-3,1;1,-2,0,7;0,2,1,-1;1,4,-1,-2]; b=[11,-7,5,0]';
[Q,R]=qr(a)
X=R\(Q\b)
>>Q =
0.8165 -0.0705 -0.3994 -0.4109 0.4082 -0.5643 0.6700 0.2568 0 0.4232 0.6124 -0.6678
0.4082 0.7053 0.1287 0.5650
R =
2.4495 1.6330 -2.8577 2.8577 0 4.7258 -0.0705 -5.8544 0 0 1.6818
3.4212 0 0 0 0.9246 X =
43.0000
-11.5556
17.6667
-10.4444
20.设
1
,
1t
y t
e
ππ
-
=-≤≤
+在同一图形窗口采用子图的形式绘制不同图形,条形
图、阶梯图、杆图和对数坐标图。答:t=-1*pi:pi/20:pi;
y=1./(1+exp(-1.*t));
subplot(2,2,1);bar(t,y,'g'); title('bar(t,y,''g'')'); subplot(2,2,2);stairs(t,y,'b'); title('stairs(t,y,''b'')'); subplot(2,2,3);stem(t,y,'k'); title('stem(t,y,''k'')'); subplot(2,2,4);loglog(t,y,'y');
title('loglog(t,y,''y'')');
一、 选答题(在下列题中选答4题):
1..符号运算(必选两道)
(1) 求微分方程222/()/(2)dy dt
t y t =+的通解 (2) 求代数方程1
122a x x +=
+-的解 (3) 求函数极限2lim
(1)x x x x →∞+- (4) 求不定积分sin 24x
xtd t ?
2. 求解方程的数值解(必选两道)
(1) 0-0x =0.5=x x sin(x)-10在附近的根
解:(1) 建立函数文件funx.m 。
function fx=funx(x)
fx=x-10.^x+2;
(2) 调用fzero 函数求根。
z=fzero('funx',0.5) z =
0.3758
(2) 222x -0.3c o s (y )=0
3y -0.6c o s (x )=0非线性方程组在(0.5,0.5) 附近的数值解
解:(1) 建立函数文件myfun.m。
function q=myfun(p)
x=p(1);
y=p(2);
q(1)=x-0.6*sin(x)-0.3*cos(y);
q(2)=y-0.6*cos(x)+0.3*sin(y);
(2) 在给定的初值x0=0.5,y0=0.5下,调用fsolve函数求方程的根。
x=fsolve('myfun',[0.5,0.5]',optimset('Display','off'))
x =
0.06354
0.3734
5. 有初始状态为0的二阶微分方程''0.2'0.40.2()
++=,其中u(t)是单位阶
x x x u t
跃函数,建立系统并仿真,并用scope模块观察结果。
解:自定义一个函数,从workspace获得数据.
函数定义的m文件:
function Ts=yuejie t=0:0.1:20;
y=heaviside(t);
Ts=[t', y'];
积分模块,设置Ts的属性
Ts的数据由运行Ts=yuejie 来获得
结果如图:
转换函数
x''+0.2x'+0.4x=0.2u(t) 两边进行变换s^2X(s)+0.2sX(s)+0.4X(s)=0.2U(s) 整理得G(s)=X(s)/U(s)=0.2/(s^2+0.2s+0.4)
如图包括设置转换函数属性
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('请输入商品价格:');
if price<200
y=price;
elseif (price>=200&price<500)
y=0.97*price;
elseif (price>=500&price<1000)
y=0.95*price;
elseif (price>=1000&price<2500)
y=0.92*price;
elseif (price>=2500&price<5000)
y=0.90*price; elseif (price>=5000) y=0.86*price; end
y
基于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期末作业 学院:机电工程学院 专业:机械制造设计及其自动化 学号: 班号: 姓名:
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期末作业 学院:机电工程学院 专业:机械制造设计及其自动化 学号: 班号: 姓名: 我去人也就有人!为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扼腕入站内信不存在向你偶同意调剖沙
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
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分)
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
哈工大机械设计电算实验(matlab2010)
机械设计电算实验 一:普通V带传动设计内容和任务 1、普通V带传动设计内容 给定原始数据:传递的功率P,小带轮转速n,传动比i及工作条件。 设计内容:带型号,基准长度Ld,根数Z,传动中心距a,带轮基准直径dd1、dd2,带轮轮缘宽度B,初拉力F0,和压轴力FQ。 2、CAD任务: (1)编制V带传动设计程序框图。 (2)编制V带传动设计原程序。 (3)按习题或作业中数据运行程序,要求对每一组数据各按三种V带型号计算,对每一种带型号选三种小带轮直径进行计算并输出所有结果。
二、变量标识符 kka 3 工作情况系数K A n r/min 4 小带轮转速n 1 5 工作机载荷特性载荷类型 6 原动机类型原动机类型 7 工作年限a年b班 8 带型号V带带型 dd1 mm 9 小带轮基准直径d d1 dd2 mm 10 大带论基准直径d d2 11 带速v v m/s vmax m/s 12 最大带速v max a0 mm 13 预估中心距a ld mm 14 基准长度L d 15 实际中心距 a aa mm r 度 16 小带轮包角α 1 17 V带根数Z z 18 基本额定功率P po kw pp0 kw 19 功率增量ΔP ki 20 传动比系数K i
三、程序框图
四、源程序与其说明 程序说明: 本程序用Matlab2010b软件编制,主要针对机械设计大作业上的题型设计。使用时只要打开m 文件,并点击运行,按照提示进行即可。首先输入原始数据,然后根据自己的需要选择带型,中心距即可得到设计结果,无需再查找资料,方便高效,计算过程如有错误会进行提示,并返回到输入处进行改正。而且该程序可以直接计算下一带轮直径或者计算下一带型,比较方便。源程序如下(先复制到记事本,再新建一个m文件,粘贴) clear all; disp('欢迎使用本程序,请输入V带传动设计的原始数据'); p=input('电动机工作功率(kw) P='); n=input('电动机满载转速(r/min) nm='); i=input('第一级传动比 i1='); a=input('请输入最短工作工作年限 a年b班 a='); b=input(' b='); disp('是否反复起动、正反转频繁或工作条件恶劣'); ka1=input('是请输入1,否请输入0。请输入:'); disp('原动机类型:'); disp('I类原动机包括普通笼型交流电动机,同步电动机,'); disp(' 直流电动机(并激),n>=600r/min的内燃机') disp('II类原动机包括交流电动机(双笼型、滑环式、单相、大转差率),'); disp(' 直流电动机(复激、串激),单缸发动机,n<=600r/min的内燃机') d1=input('请选择原动机的类型,输入1或2。请输入:'); disp('载荷类型:'); disp('1、载荷平稳;2、载荷变动小;3、载荷变动较大;4、载荷变动很大。'); e1=input('请选择载荷类型请输入:'); c=[3 710 38 80;2.2 940 28 60;3 960 38 80;4 960 38 80]; if p==3&&n==710 d=38; e=80;
2013年春季学期 MATLAB 课程考查题 一、必答题: 1.matlab常见的数据类型有哪些?各有什么特点? (1)numeric(数值):包括整数和浮点数,缺省状态下matlab将所有的数看 作是双精度的浮点数。 (2)char(字符):字符串必须用单引号括起来,是由多个字符组成的字符行 向量。 (3)logical (逻辑):仅有TRUE和FALSE两个值,一般用于关系运算或逻辑 运算。 (4)cell (单元、细胞、元胞):可以把不同类型的数据纳入到一个变量中,不 要求相同的数据类型,下标用大括号括起来。 (5)struct (结构):结构是包含已命名“数据容器”或字段的数组。结构中的 字段可以包含任何数据。 2.MATLAB中有几种帮助的途径? (1)帮助浏览器:选择view菜单中的Help菜单项或选择Help菜单中的 MATLAB Help菜单项可以打开帮助浏览器; (2)help命令:在命令窗口键入“help”命令可以列出帮助主题,键入“help 函数名”可以得到指定函数的在线帮助信息; (3)lookfor命令:在命令窗口键入“lookfor 关键词”可以搜索出一系列与给 定关键词相关的命令和函数 (4)模糊查询:输入命令的前几个字母,然后按Tab键,就可以列出所有以 这几个字母开始的命令和函数。 3.Matlab常见的哪三种程序控制结构及包括的相应的语句? (1)顺序结构:数据的输入-A=input(提示信息,选项); 数据的输出-disp(输出项) 程序的暂停-pause(延迟秒数)
(2)选择结构:if 条件语句组;else 语句组;end Switch 表达式case 值1 语句组1;……otherwise 语句组n; end (3)循环结构:while 条件语句组;end for(条件);end do 语句组while 循环语句组;end 4.命令文件与函数文件的主要区别是什么? 答:(1)命令文件是一系列命令的组合,函数文件的第一行必须用function 说明; (2)命令文件没有输入参数,也不用返回参数,函数文件可以接受输入参数,也可以返回参数; (3)命令文件处理的变量为工作空间变量,函数文件处理的变量为函数内部的局部变量,也可以处理全局变量。 5.什么是全局变量和局部变量? 答:全局变量是没有定义存储类型的外部变量,其作用域是从定义点到程序结束,省略了存储类型符,系统将默认为是自动型 局部变量也称为内部变量。局部变量是在函数内作定义说明的。其作用域仅限于函数内,离开该函数后再使用这种变量是非法的。 如果一个函数内的变量没有特别声明,那么这个变量只在函数内部使用,即为局部变量,如果两个或多个函数公用一个变量(或者说在子程序中也要用到主程序中的变量,注意不是参数),那么可以用global来将它声明为全局变量。 全局变量的使用可以减少参数传递,合理利用全局变量可以提高程序执行效率。如果需要用到其他函数的变量,就要利用在主程序与子程序中分别声明全局变量的方式实现变量的传递,否则函数体内使用的都为局部变量。 6.数据插值和曲线拟合区别? 答:有这样一组数据(Xi,Yi),i=1,2,3,。。。。。,寻找Y与X之间的关系Y=F(X),就要根据这组数据求出其近似关系Y=G(X)。 插值与曲线拟合实际上是求这个近似关系的两种方法,不同的是,插值方法
1. 表1 用三次样条方法插值计算0-90 度内整数点的sin 值和0-75 度内整数点的正切值,然后用5 次多项式拟合方法计算相同的函数。 a(度)0 15 30 45 60 75 90 Sin(a)0 0.2588 0.5000 0.7071 0.8660 0.9659 1.0000 tan(a)0 0.2679 0.5774 1.0000 1.7320 3.732 解:分别对应的程序如下: 正弦函数: x = pi*(0:90)/180; y = sin(x); xx = pi*(0:.25:90)/180; yy = spline(x,y,xx); plot(x,y,'o',xx,yy) 正切函数: x = pi*(0:75)/180; y = tan(x); xx = pi*(0:.25:75)/180; yy = spline(x,y,xx); plot(x,y,'o',xx,yy)
正弦拟合: figure x=pi*(0:15:90)/180; y=[0,0.2588,0.5,0.7071,0.866,0.9659,1.0]; xx=pi*(1:0.05:90)/180; p2=polyfit(x,y,5); yy=polyval(p2,xx); plot(x,y,'-ro',xx,yy); 正切拟合: figure x=pi*(0:15:75)/180; y=[0,0.2679,0.5774,1,1.732,3.732]; xx=pi*(1:0.05:75)/180; p2=polyfit(x,y,5); yy=polyval(p2,xx); plot(x,y,'-ro',xx,yy); legend('描点显示','五次拟合')
春季学期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)
a=[6,9,3;2,7,5]; b=[2,4;4,6;1,8]; a*b ans = 51 102 37 90 14. A=[4,9,2;7,6,4;3,5,7]; B=[37;26;28]; X=inv(A)*B X = -0.5118 4.0427 1.331815. syms x solve('x^3-7*x^2+2*x+4=0') 15. syms x solve('x^3-7*x^2+2*x+40=0') ans = -2 4 5 16. x=8; y=(x-1)*(x-2)*(x-3)*(x-4) y = 840 17. syms x y=3*x^3+13*x^2+6*x+8; A=factor(y); B=A/(x+4) B =
3*x^2 + x + 218. 18. syms x y=sqrt(4*x^4-12*x^3-14*x^2+5*x+9); df=diff(y) df = -(- 16*x^3 + 36*x^2 + 28*x - 5)/(2*(4*x^4 - 12*x^3 - 14*x^2 + 5*x + 9)^(1/2)) syms x f=4*x^4-12*x^3-14*x^2+5*x+9; nf=int(f) nf = (4*x^5)/5 - 3*x^4 - (14*x^3)/3 + (5*x^2)/2 + 9*x 19. 20. t=0.1:0.1^5*pi:2*pi; x=sin(3*t)*cos(t); y=xin(3*t)*sin(t); plot(x,y) 21. x=0:0.01:10; y1=exp(-0.1*x); y2=exp(-0.2*x); y3=exp(-0.5*x); plot(x,y1,'b',x,y2,'R',x,y3,'k') xlabel('x') ylabel('y') title('Sine and Cosine Curve') gtext('y=exp(-0.1*x)') gtext('y=exp(-0.2*x)') gtext('y=exp(-0.5*x)')
一.实验内容及目的 1.实验内容 一个单相合成电动势中含有明显的5次谐波和7次谐波,其表达式如下所示。采用短距方法,验证削弱谐波电动势的效果。 E=sin(Vs)+0.2sin(Vs)*sin(5Vs)+0.1sin(Vs)*sin(7Vs) 2实验目的 ①了解削弱谐波电动势的方法; ②探究采用短距绕组削弱谐波电动势的方法; ③探究采用增加槽数削弱谐波电动势的方法; ④探究采用斜槽削弱谐波电动势的方法; ⑤加深对绕组谐波电动势产生原因的理解。 二.实验要求及要点描述 1.削弱谐波电动势的方法 分析要点 ①采用屏幕图形方式直观显示; ②利用MATLAB 编程实现; ③用INSPACE 命令对画出合成波形; ④用计算公式绘制斜槽因数曲线。 三.基本知识及实验方法描述 基本知识 由谐波电动势公式v wv v v Nk f E φφ44.4=可见,通过减少v φ或wv k ,可以降低v E φ。常用的方法有: (1)采用短距绕组 适当选择线圈的节距,使得莫一次谐波的节距因数等于零或者接近于零,即可达到消除或削弱该次谐波的目的,如果要消除v 次谐波,只要选择 001 18090?=k y v τ,就能使该次谐波的节距因数0)90(sin 1 ==o pv y v k γ 从消除谐波的观点来看,式中的k 可选为任意整数;但从尽可能不削弱基波的 角度考虑,应当选用接近于整距的短节距,即选择12-=v k ,此时τ)11(1v y -=。也就是说,为消除第v 次谐波,应当选用比整距短τv 1的短距线圈。 由于三相绕组的线电压间不会出现三次谐波,所以选择三相绕组节距时,主要
应当考虑如何减小5次和7次谐波,故 1y 可选为5/6(=0.833)左右。 (2)采用分布绕组 对于分布绕组,每极每相槽数q 越多,抑制谐波的电动势的效果越好。但是q 增多,意味着总槽数增多,这将使电机的成本提高。q>6时,高次谐波分布因数的下降已不太显著。 (3)采用斜槽 采用斜槽后,由于同一根导体的各个小段在磁场中的位置互不相同,所以与直槽时相比,导体中的感应电动势将有所削弱。一般最常用的,是使得整个斜槽一个槽距。
地震工程大作业 哈工大李金平 弹性反应谱 原创性声明,主程序由本人独立编写完成,支持各种检验。所选地震动RSN2345CHICHI.AT2 peer索取号:2345 : 010203040506070 -0.25 -0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 0.25 时间 (s) 加 速 度 ( m / s 2 )
位移反应谱 局部放大 1234 5678910 周期T (s) 位移 (m m ) 1 1.02 1.04 1.06 1.08 1.1 1.12 9.2 9.259.39.359.49.459.5 9.559.69.659.7周期T (s) 位移 (m m )
求比值 1 2 3 4 567 8 9 10 0.940.950.960.970.980.991 1.011.021.031.04 周期T (s) 比值
绝对加速度反应谱 局部放大 0.2 0.4 0.6 0.8 1 1. 2 1.4 1.6 1.8 2 00.10.20.30.40.50.60.70.80.91周期T (s) S a (m /s 2) -0.25 -0.2 -0.15 -0.1 -0.05 00.050.1 0.15 0.2 0.25 0.150.2 0.250.3 0.35 0.4 0.45周期T (s) S a (m /s 2 )
求比值进行比较 速度反应谱 1 2 3 4 567 8 9 10 0.940.950.960.970.980.991 1.011.021.031.04 周期T (s) 比值 1 2 3 4 567 8 9 10 00.010.020.030.040.050.060.07 0.080.090.1周期T (s) 速度 (m / s )
2014年春季学期 MATLAB 课程考查题 姓名: 学号:11208 学院:机电工程学院 专业:机械设计制造及其自动化 一.必答题(80分) 1. 如何设置当前目录和搜索路径,在当前目录上的文件和在搜索路径上的文件 有何区别? 答:设置当前目录和搜索路径:在File菜单中选择SetPath选项,之后选择AddFolder增加目录。当前工作目录是指MATLAB运行文件时的目录,只有在当前工作目录或搜索路径下的文件、函数可以被运行或调用。 2. 创建符号变量和符号表达式有哪几种方法? 答:(1)符号变量: x = sym(‘x’) 创建x为符号变量,默认复数区域 x = sym(‘x’, ‘real’) 创建实数的符号变量x x = (‘x’, ‘positive’) 创建正数的符号变量x x = sym('x', 'clear')创建一个没有额外属性的纯形式上的符号变量x s=sym(‘ab’,’flag’) 创建flag数域(复数,实数,正数)符号变量名s,内容ab (2)符号表达式: ①直接法:>> x=sym('x'); >> a=sym('a'); >> b=sym('b'); >> f=sin(b*x)+exp(-a*x) ②整体定义法:f=sym(‘expression’) ③字符串符号表达式:f=‘expression’
3. GUIDE提供哪些常用的控件工具,各有什么功能?(5分) 答:按钮(Push Buttons) :通过鼠标单击按钮可以执行某种预定的功能或操作; 静态文本框(Static Texts):仅用于显示单行的说明文字. 文本编辑器(Editable Texts):用来使用键盘输入字符串的值,可以对编辑框中的内容进行编辑、删除和替换等操作; 单选按钮(Radio Button):单个的单选框用来在两种状态之间切换,多个单选框组成一个单选框组时,用户只能在一组状态中选择单一的状态,或称为单选项; 滚动条(Slider):可输入指定范围的数量值,通过移动滚动条来改变指定范围内的数值输入,滚动条的位置代表输入数值。 复选框(Check Boxes):单个的复选框用来在两种状态之间切换,多个复选框组成一个复选框组时,可使用户在一组状态中作组合式的选择,或称为多选项; 列表框(List Boxes):列表框显示列表项,并能够选择其中的一项或多项。 弹出式菜单(Popup Menus):让用户从一列菜单项中选择一项作为参数输入。 组合框(面板)panel:组合框是图形窗口中的一个封闭区域,它把相关联的控件组合在一起。 开关按钮(Toggle Button):产生一个一个二进制状态动作(开或关),当鼠点击它时按钮将下陷,并执行callback(回调函数)中指定的内容,再次点击,按钮复原,并再次执行callback 中的内容。 坐标轴axes:坐标轴可以设置关于外观和行为的参数。 按钮组button group:按钮组类似于组合框,但是它可以响应单选按钮以及开关按钮的高级属性。 4. 建立不同数据类型的变量:包括双精度的数值型(numeric)、符号型(sym)、结 构体(structure)、单元型(cell)、字符串型(char),函数句柄(function handle)、逻辑型(logical),并且采用whos列出各变量名,变量的数组维数,占用字节数,变量的类别。(4分)
哈尔滨理工大学 2007-2008学年 第二学期考试试题 A 卷 考试科目:MATLAB 语言及应用 考试时间:120分钟 试卷总分100分 考试班级:自动化06-5、6、7、8 一、单项选择题(每题2分,共20分) 1、下列关于help help 命令,说法正确的是( )。 A .打开帮助浏览器; B .显示help 命令的帮助注释区; C .显示help 命令帮助注释区的H1行; D .显示H1行含有help 字符的所有命令。 2、有一个2行2列的细胞数组c ,则c(2)是指( )。 A .第1行第2列的元素内容 B .第2行第1列的元素内容 C .第1行第2列的元素 D.第2行第1列的元素 3、计算三个多项式s1、s2和s3的乘积,则算式为( )。 A.conv(s1,s2,s3) B.s1*s2*s3 C.conv(conv(s1,s2),s3) D.conv(s1*s2*s3) 4、mesh 函数是( )。 A.绘制三维曲面图 B.绘制三维网线图并添加平行于z 轴的边框线 C.绘制三维表面图 D.绘制三维网线图 5、已知x 为一个向量,计算ln(x)的MATLAB 命令是( )。 A .ln(x) B .log(x) C .Ln(x) D .lg10(x) 6、在MA TLAB 中( )用于括住字符串。 A ., B .; C .’’ D .”” 7、已知x=linspace(1,5,9),则x 有( )个元素。 A .10 B .11 C.9 D .12 8、已知数组12 3456789a ?? ??=?? ???? ,则a(:,end)是指( )。 A .第三列元素 B .第一行元素 C .所有元素 D .第三行元素 9、下面switch 结构,正确的是( c )。 A .switch a B. switch a C.switch a D.switch a case a>1 case a=1 case 1 case =1 10、下列变量名中,合法的是( a )。