Matlab考试题库附加答案

Matab考试题库附加答案

一、辨析题（共2小题，每题 5 分，共计 10 分）

二、问答题（共3小题，每题 5 分，共计 15 分）

三、计算题（共7小题，每题 5 分，共计 35 分）

四、编程题（共4小题，每题 10 分，共计 40 分）

填空题（计算部分）

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

18.已知A=[1 3;4 6];C=[A,eye(size(A));A ,ones(size(A))] 则Ｃ＝

_____________.

[1 3 1 0；4 6 0 1；1 3 1 1；4 6 1 1]_ 19. A=[1,2,3;4,5,6]; A(4)=__________, A(3,2)=__________________

5_ 无答案

20.A=[1,2;3,4];B=[1,0;0,1]; A*B=__________, A.*B=____________.

[1 2；3 4]

[1 0；0 4]

21. A=[1 2 3;6 5 4]; reshape(A,3,2)=______________.

[1 5；6 3；2 4]

22. A=[1,2;3,1]; B=[1,0;0,1]; A~=B=___________.

[0 1；1 0]

23. A=[1,2,3;1,4,9;1,7,27]; triu(A,-1)=_____.提取矩阵上三角矩阵的函数

[1 2 3；1 4 9；0 7 27]

24. A=[4,15,-45,10,6;56,0,17,-45,0]; find(A>=10 & A<=20)=___ __________.

3 6 7

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

A &

B 的运行结果_______；

A ./

B 的运行结果_______。

ans=[0 1 1;1 1 0]

ans= [0 2.25 2;1 0.6 NaN]

26. 产生4阶全0方阵的命令为_______；

产生3阶全1方阵的命令为。

zeros(4)

ones(3)

27. A=rand(2,5);b=size(A);c=length(A)；则b和c的值分别为_______。

2 5和 5

28. 用if判断语句判断80≥ x >60,在MATLAB中if语言后的判断应写为_______。

if x<=80 & x>60

29.为了使两个pl ot的图形在同一个坐标显示，可以使用_______命令进行图形保持；可以使用_______命令为图形添加网格。

hold on

grid on

30. MATLAB的工作空间中有三个变量v1, v2, v3，写出把它们保存到文件my_data.mat中的指令_______；

save my_date

31. MATLAB提供了两种运行方式，是按_______和把_______文件copy到matlab里，回车。

F5

*.m

32．MATLAB的含义为_______。

矩阵实验室（Matrix Laboratory）33.MATLAB中，都包括哪些工具箱，举例不少于6个：_______

偏微分方程工具箱、最优化工具箱、插值运算工具箱、数理统计工具箱、信号处理工具箱、通信工具箱，滤波设计工具箱。

34、在MATLAB中，plot命令用于什么，。画曲线

35、size命令的功能，确定_______。阵列的长度

36、a=[1 2 2]，求a的转置_______。a=[1;2;2]

37、在线形图型格式的设置中，字符r表示什么颜色。红色

38. 在MATLAB中,矩阵的秩由函数_______来实现。rank

39．在MATLAB中,单位矩阵由函数_______来实现。eye

40. 在MATLAB中,二维饼图由函数_______来实现。pie

41. 在线形图型格式的设置中，字符b表示什么颜色。蓝色

42. MATLAB表达式3*2^3^2的结果是_______。192

43.条形图由函数_______来实现。bar

44. 倘若要是对x进行赋值，从5到25，间隔是0.1, x=5:0.1:25;

45. 倘若要是对x进行赋值，从2到20，中间间隔100个点,

x=linspace(2,20,100);

选择题（计算部分）

1. 下列变量中的最大数是()

A．eps B．realmin C．realmax D．-pi

2. 下列表示可以作为MATLAB的变量的是()

A．abcd- 2 B．xyz_2# C．@h D．X_1_a

3. 如果x=1: 2 : 8,则x(1)和x(4)分别是()

A．1，8 B．1, 7 C．2, 8 D．2, 7

4. 下列变量中比0大的最小数是()

A．i B．realmin C．realmax D．pi

5. 在MATLAB中下列数值的表示不正确的是（）

A．+99 B．1.3e-5 C．2-3*e^2 D．3-2*pi

c d b b b

6. 下列表示可以作为MATLAB的变量的是()

A．abcd-2 B．xyz_2a C．@h D．3hour

7. 如果x=1:-2:-8,则x(1)和x(5)分别是()

A．1，-8 B．-1,-7 C．-1,-8 D．1,-7

8. MATLAB表达式2*2^3^2的结果是()

A．128 B．4096 C. 262144 D．256

9. 下列哪个变量的定义是不合法的（）。

(A) abcd-3 (B) xyz_3 (C) abcdef (D) x3yz

10. 指出下列错误的指令（）。

(A) syms a b; (B) syms a, b; (C) syms(…a', …b'); (D) syms(…a','b','positive');

b d a A B

11.下列哪条指令是求矩阵的行列式的值（）。

(A) inv (B) diag (C) det (D) eig

inv是求逆矩阵

diag是矩阵对角元素的提取和创建对角阵

sqrt(x) ——求x的平方根

abs(x)——求x的绝对值

det(a)——求行列式的值

eig——计算矩阵A的特征值和特征向量的函数

clf; 用来清除图形的命令

12 .清空Matlab 工作空间内所有变量的指令是（）。

(A) clc (B) cls (C) clear (D) clf

13、下列变量名中————是合法的。

(A) char_1 ; (B) x*y ; (C) x\y ; (D) end

14.以下哪个说法是正确的（）。

(A)Matlab 进行数值计算的表达精度与其指令窗口中的数值显示精度相同。

(B)Matlab 指令窗口中显示的数值有效位数不能超过7 位。

(C)输入二维数值数组时，需要用到逗号和分号，它们可以在中文状态下输入。

(D)历史指令窗口所记录的内容与diary 指令所产生“日志”内容是不同的。

15. 产生四维单位矩阵的语句为().

A.ones(4)

B.eye(4)

C.zeros(4)

D.rand(4) >> eye(4)

ans =

1 0 0 0

0 1 0 0

0 0 1 0

0 0 0 1 >> ones(4)

ans =

1 1 1 1

1 1 1 1

1 1 1 1

1 1 1 1

>> zeros(4)

ans =

0 0 0 0

0 0 0 0

0 0 0 0

0 0 0 0 >> rand(4)

ans =

0.8147 0.6324 0.9575 0.9572

0.9058 0.0975 0.9649 0.4854

0.1270 0.2785 0.1576 0.8003

0.9134 0.5469 0.9706 0.1419

C C a

D b

16.求解线性方程组Ax=b,当det(A)≠0时，方程的解是().

A. A\b

B.b/A

C.b\A

D.A/b

17. 运行如下程序后，输入９回车，命令窗口（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

18. 运行下面的代码后，输入：１回车，２回车，１回车，则输出结果为（）a=input('a=?');

b=input('b=?');

c=input('c=?');

d=b*b-4*a*c;

x=[(-b+sqrt(d))/(2*a),(-b-sqrt(d))/(2*a)];

disp(['x1=',num2str(x(1)),',x2=',num2str(x(2))]);

Ａ. x1=num2str(1) x2=num2str(1)

B. x1=num2str(-1) x2=num2str(1)

C. x1= -1 x2=-1

D. x1=1 x2=-1

19. 清除工作空间（wordspace）的命令是（）

Ａ. clc; B. clear; C. clf; D.delete;

20. 运行如下程序后，输入９回车，命令窗口（command windows）显示的结果为( )

x=input('请输入x的值:');

if x==10

y=cos(x+1)+sqrt(x*x+1);

else

y=(3^(1/2))*sqrt(x+sqrt(x));

end

y

A. 9

B. 8

C. 7

D. 6

a a c

b d

21. 运行如下程序后, 命令窗口（command windows）显示的结果为( b)

s=0;

a=[12,13,14;15,16,17;18,19,20];

for k=a

s=s+k;

end

disp(s');

A. 144

B. 39 48 57

C.145

D. 45 48 51

22. 运行如下程序后, 命令窗口（command windows）显示的结果为(b )

k=0;

for n=100:200

if rem(n,21)~=0 %R=rem(X,Y),求余数函数,X,Y应该为正数

k=k+1;

continue

end

break;

end

A.105

B. 5

C.4

D.101

23. 图形窗口（Figure）显示网格的命令是（b）

A. axis on

B. grid on

C. box on

D. hold on

24. 已知函数文件如下，则factor(4)=(c )

function f=factor(n)

if n<=1

f=1;

else

f=factor(n-1)*n;

end

A. 4

B. 12

C. 24

D.48

25. 运行如下程序后, 命令窗口（command windows）显示的结果为( d)

A=[13,-56,78; 25,63,-735; 78,25,563; 1,0,-1];

y=max(max(A))

A. y=564

B.y=9

C.y=-735

D.y=563

26. 在图形指定位置加标注命令是（c）

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)?); %添加图例的标注,

27.下列哪个函数为插值函数（b）

A. P=polyfit(X,Y,3)

B. Y1=interp1(X,Y,X1,'method')

C. [Y,I]=sort(A,dim)

D. R=corrcoef(X)

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

(A)a=4, b=4, 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

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

(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

30、在循环结构中跳出循环，但继续下次循环的命令为 。c (A) return; (B) break ; (C) continue ; (D) keyboard

31. 用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]

32. 已知a=2:2:8, b=2:5，下面的运算表达式中，出错的为 c (A) a' *b (B) a .*b (C) a*b (D) a-b

辨析题：

1 MATLAB 有几种建立矩阵的方法？各有什么优点？。

2 MATLAB 提供编程语言，请辨析M 文件编程中的命令文件与函数文件两种形式。

3 辨析matlab 中的拟合与插值有什么区别

4. 判断下列一段程序的错误，并进行改正。

x=-2:0.2:2;

y=x.^2; plot(x,y)

str1 = num2str(min(x));% 数字转化为字符串 str2 = num2str(max(x));% 数字转化为字符串 out = [Value of f from 'str1' to 'str2']; xlabel(out);

out = ['Value of f from ' str1 ' to ' str2];

或out = ['Value of f from str1 to str2'];

5根据下面两个矩阵

??

??

?

?????=??????????=987654321,136782078451220124B A 执行下列的矩阵运算命令，并回答有关的问题

（1） A.*B 和 A*B 将分别给出什么结果，它们是否相同，为什么？ （2） 得出A.^B 、A/B 及A\B 的结果，并分别解释它们的物理意义。

问答题：

1 MATLAB系统包括哪些窗口，各个窗口功能是什么？

2 分析下面程序运行后s1,n1,s2,n2,m的值。

s1=0;s2=0;n1=0;n2=0；

x=[1,-4,-8,3,10,-9,7,-3,10, 8,-5,-2,2,0 ];

m=length(x);

for i=1:m

if x(i)<0

s1=s1+x(i);

n1=n1+1;

else

s2=s2+x(i);

n2=n2+1;

end

end

s1,n1,s2,n2,m

请回答s1,n1,s2,n2,m的结果。

3 已建立一个简单的SIMULINK模型，如下图，用信号发生器产生一个仿真图，请对下图

的符号和功能进行解释。

3 A = [1 2; 3 4]; B = [1 2; 1 4];

写出下面几条指令的运行结果：

C = A*B

D = A+B

E = A.*B

答：

C =

3 10

7 22

D =

2 4

4 8

E =

1 4

3 16

4.mysum=0；i=1；

while(i<=100)

mysum=mysum+i；

i=i+1；

end

mysum

答：

mysum =

5050

5.在MATLAB窗口中，用哪两种方法启动Simulink？答：

在MATLAB命令窗口中输入“simulink”再回车，或者单击工具栏中的SIMULINK 图标

6. 请修改下面的程序，让他们没有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

解：程序如下：a=[1 2 3;4 5 6;7 8 9];

b=(a<=8&a>=2);a=a.*b

兽血沸腾黑岩https://www.wendangku.net/doc/b011551351.html,/

7 编写一个函数，使其能够产生如下的分段函数：

并调用此函数，绘制x=[0,+2]范围内的f(x)*f(x+2)。

答：function y=f(x)

if x<=2

y=0.5*x;

else if x>6

y=0.5;

else y =1.5-0.25*x;

end

end

运行结果x=2

f(x)=1

x = 0:0.05:2;

y = diag(A2(x)'*A2(x+2));

plot(x,y);

xlabel('\bfx');

ylabel('\bfy');

0.2

0.4

0.6

0.8

1 1.2

1.4

1.6

1.8

2

00.050.10.150.20.250.30.350.40.45

0.5x

y

8（1） 将一个屏幕分4幅，选择合适的步长在右上幅与左下幅绘制出下列函数的图形。 ①

,x

（曲线图）；②f(x,y)=(-2<=x<=2,-4<=x<=4)

（曲

面图）。

答：subplot(2,2,2);

ezplot('(cos(x))^(1/2)',[-pi/2 pi/2]); ylabel('y'); subplot(2,2,3); x=-2:0.5:2; y=-4:1:4;

ezsurfc('x^2/2^2+y^2/4^2')

-1

01

0.5

1

x

(cos(x))1/2

y

x x 2/22+y 2/42

y

(2) A 是一个維度m×n的矩阵. 写一段程序, 算出A中有多少个零元素。答：>> A= input ('请输入一个矩阵')

[m,n]= size(A);

sig=0;

for i=1:m

for j=1:n

if A(i,j)==0

sig = sig+1;

end

end

end

请输入一个矩阵[0 1 2;1 0 2; 0 0 0]

A =

0 1 2

1 0 2

0 0 0

>>sig

sig =

5

9 向量A=[]. 写一段程序, 找出A中的最小元素。

答：A= input ('请输入一个向量')

[m,n]=size(A)

min =A(1,n);

for i=1:n

if A(1,i)

min=A(1,i)

end

end

请输入一个向量[1 2 3 -5 2 0]

A =

1 2 3 -5 2 0

min =

-5

10 写出执行以下代码后C，D，E的值A=[1,2,3;4:6;7:9];

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

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

E=C(2,[1 2])

C =

1 2 3

4 5 6

7 8 9

10 11 12

D =

2 3

5 6

8 9

E =

4 5

11 写出执行以下代码后，MATLAB命令窗口上显示的x矩阵的值。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);

1 5 3 9 5 13 7 17

12 创建符号函数并求解，要求写出步骤和运行结果

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

(2)求f=0的解

syms a x b c;

f=a*x^2+b*x+c;

solve(f)

ans =

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

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

13 求解以下线性方程组，要求写出程序代码和运行结果。

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

x1+3x2+ x4=6

x1-x2+x3+8x4=1

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

解：程序代码：a=[2 -3 1 2;1 3 0 1;1 -1 1 8;7 1 -2 2];

b=[8 6 1 5]';

ra=rank(a);

rb=rank([a b]);

det(a);

xx=a\b

运行结果：

xx =

3.1079

1.3278

7.9502

-1.0913

14. (1) 写出关系运算符中的等号、不等号、小于号、大于号、小于等于号和

大于等于号；

(2) 写出逻辑操作中逻辑“与”、逻辑“或”及逻辑“非”的符号；

(3) 并用语句表达“如果a等于b且c等于0就让d 等于3，否则如果

a大于b且c=1让d等于0，其余情况让d等于3”；

解：（1）== ~= < > <= >=

（2）& | ~

（3）if a==b&c==0

d=3

else if a>b&c= =1

d=0

else

d=3

end

15. 设A=

和B=是两个10行10列的矩阵（数组），试说

明命令A*B, A/B, A\B, A.*(B.^A), A./B, A.\B和sin(A).*exp(B)的涵义。

解：A,B的矩阵乘法，矩阵A乘以B的逆，矩阵B乘以A的逆，B的每个元素的对应A次幂乘以A对应的每个元素

A的每个元素除以B对应的每个元素，B的每个元素除以A对应的每个元素，A的每个元素的正弦与B对应的每个元素的指数求积，

16. 设A为14*17数组，B为一个行数大于3的数组，请给出

（1）删除A的第4、8、12三列的命令；（2）删除B的倒数第3行的命令；

解：(1)A(:,[4 8 12])=[]

(2) B(end-2)=[]

编程简题

1 利用MATLAB数值运算，求解线性方程组(将程序保存为？？？.m文件)

2 求下列联立方程的解

3x+4y-7z-12w=4

5x-7y+4z+ 2w=-3

x +8z- 5w=9

-6x+5y-2z+10w=-8

求系数矩阵的秩；求出方程组的解。

解：（1）

>> a=[3 4 -7 -12

5 -7 4 2 ;

三相异步电动机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试卷合集附答案

一、填空题 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, >> figure('Color',[1,0,0],'WindowButtonDownFcn','text(0.5,0.5,''Left Button Pressed'')');axis off 2, (1)默认属性 >> x=linspace(0,600,100000);y=(log(x+sqrt(1+x.*x)))/2;line(x,y) 句柄操作 >> x=linspace(0,600,100000); y=(log(x+sqrt(1+x.*x)))/2; line(x,y,'linewidth',3,'linestyle','-.','color','r'); text(300,3.5,'文字标注') (2)默认属性 >> t=linspace(0,20,500);x=t.*t;y=5.*t.*t;line(x,y) 句柄操作 >> t=linspace(0,20,500); x=t.*t;y=5.*t.*t;line(x,y); line(x,y,'linewidth',3,'linestyle','-.','color','r'); text(250,1600,'文字标注') 3, (1) >> x=linspace(0,50,1000); [x,y]=meshgrid(x); z=x.*x+y.*y-5.*sin(x.*y); axes('view',[-37.5,30]); surf(x,y,z); light('position',[10,20,4000]); shading interp (2) >> x=linspace(0,50,1000); [x,y]=meshgrid(x); z=y.*y.*y; axes('view',[-37.5,30]); surf(x,y,z); light('position',[30,20,1300000]); shading interp 4, >> x=-2*pi:0.01:2*pi; y1=sin(x); y2=cos(x); axes('position',[0.1,0.6,0.2,0.2]); plot(x,y1); axes('position',[0.6,0.6,0.2,0.2]); plot(x,y2); axes('position',[0.1,0.1,0.2,0.2]); fplot('tan(x)',[-1.5,1.5]); axes('position',[0.6,0.1,0.2,0.2]);

中南大学材料学院科学计算与MATLAB考试题库

练习题 1.求函数在指定点的数值导数 x=sym('x'); >> y=[x x.^2 x.^3;1 2*x 3*x.^2;0 2 6*x]; >> x=1; >> eval(diff(y)) ans = 1 2 3 0 2 6 0 0 6 >> x=2; >> eval(diff(y)) ans = 1 4 12 0 2 12 0 0 6 >> x=3; >> eval(diff(y)) ans = 1 6 27 0 2 18 0 0 6 2.求下列函数导数 （1） x=sym('x'); >> y=x^10+10^x+(log(10))/log(x); >> diff(y) ans = 10*x^9+10^x*log(10)-2592480341699211/1125899906842624/log(x)^2/x （2） x=sym('x');

>> y=log(1+x); >> x=1; >> eval(diff(y,2)) %在x=1的条件下对y表达式求两次导数后导函数的值 ans = -0.2500 3.用数值方法求下列积分 首先先讲一下trapz的用法，如下题 t=0:0.001:1; >> y=t; >> trapz(t,y) ans = 0.5000 （1） >> x=1:0.01:5; >> y=(x.^2).*sqrt(2*x.^2+3); >> trapz(x,y) ans = 232.8066 (2) x=pi/4:0.01:pi/3; >> y=x./(sin(x).^2); >> trapz(x,y) ans = 0.3810 第三题拟合曲线题 x=[0:0.1:1]; >> y=[-0.447 1.978 3.28 6.16 7.08 7.34 7.66 9.56 9.48 9.30 11.2]; >> a=polyfit(x,y,2); >> x=[0.05:0.2:1.05]; >> y=a(3)+a(2)*x+a(1)*x.^2 %注意x要在y前先赋值，不然y不会运行为最新的x对呀的y值 y =

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

中南大学matlab题目

1 求函数在指定点的导数值 () 23 2 123,1,2,3 026 x x x f x x x x x == >> syms x >> a=[x x^2 x^3;1 2*x 3*x^2;0 2 6*x]; >> f=det(a); >> diff(f,1) ans = 6*x^2 >> diff(f,2) ans = 12*x >> diff(f,3) ans = 12 2 符号法求下列函数的导数或积分 1）y=x10+10x+log x 10,求y’ f=('x^10+10^x+log(10)/log(x)') f = x^10+10^x+log(10)/log(x) >> diff(f) ans = 10*x^9+10^x*log(10)-log(10)/log(x)^2/x 2）y=ln(1+x), 求y’’∣ x=1 f=('log(1+x)/log(e)') f = log(1+x)/log(e) >> diff(f,1,2) ans = -1/(1+x)^2/log(e) 3) y=e x/cosx,求y’ f=('exp(x)/cos(x)') f = exp(x)/cos(x)

>> diff(f) ans = exp(x)/cos(x)+exp(x)/cos(x)^2*sin(x) 4） function f=fun0(t) f=t*sin(t) int('fun0','0','pi') ans = 1/2*pi^2 5) 已知函数z=sin(xy), 计算 syms x y >> z=('sin(x*y)') z = sin(x*y) >> diff(diff(z,y,2),x) ans = -cos(x*y)*y*x^2-2*sin(x*y)*x 3 用数值方法求定积分 1） function f=fun(x) f=x.^2.*sqrt(2.*x.^2+3) quad('fun',1,5) ans = 232.8057 2） function f=fun(x) f=x./sin(x).^2 quad('fun',pi/4,pi/3) ans = 0.3835 4 已知数据[x,y]如下表，试求2次拟合多项式f(x)，然后求x=0.05,0.25,0.45,0.65,0.85,1.05各点的函数近似值，并绘出拟合曲线及求得

MATLAB试卷

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)

中南大学matlab课后习题(1)

第二章 1·求下列表达式的值。 (1）w=sqrｔ(２)*（1+0.34245*10＾-6） ｗ= 1.4１42 （2)a＝3.5;b=５；c=-9．8; x=（2＊pi＊a＋(ｃ+b）/(pi＋a＊ｂ*c)-exp（２))/(ｔａn（ｂ+c)+a)；x x ＝ ０.98２9 （3)a=３．32;b=-7.9； y＝２*pi*a＾２*[(1-ｐi/４)*b-(0.８33３-ｐi／4)*a］; y ｙ= -128.4２7１ (4）t=[2,１-3i；５，－０.65］； z=１／2*exp（2＊t)*ｌog（t+ｓｑrｔ(１＋t^２)); z z= 1.０e+00４＊ 0.0057 － 0.000７i 0．０04９ - 0.0０2７i 1.９8８4 －０.3６9６i 1.7706 - １.0５39i 2,已知a,ｂ,求下列表达式的值。 a=[-1,5,－４；０,7，８;3,6１,7];b=[８,3,-1;２,5，3;-3,2，0]; (1）a+6*b ans= 47２3－10 12 3726 -15 73 7 a^２-b+eye（3) ans= －１8 -２１7 17 22 53３ 109 ２1867 52６ (２)a＊b anｓ＝ 14 14 1６ -10 5１ 21 125 328 １80 a.*b anｓ = -8 １5 ４ 0 35 24 -9 1２2 0 b*a ans ＝ -１１０-15 7 22８５3 3 -1 28 (3）a/b anｓ =

1.２234 -0．925５ 2．9787 －0.９468 2.3511 -0．９574 ４.6170 3.8723 13．8９３6 ｂ\a ａns = -0．51０6 -8.617０ -1.1277 ０.７340 17.5７45 1.8085 －0.8８30 -21.21２8 0.4043 （４)[a,b] anｓ= －１ 5 -4 ８ 3 －１ 0 7 8 2 5 3 ３６1 7 -3 2 0 ［a([1,３],:);b^2] aｎs ＝ -１５ -4 3 61 7 73 3７ 1 17 37 １3 -20 1 ９ 3.已知ａ,完成下列操作。 a=[２3,１0,－0.778,0;4１,-45，65,５;32,5,０,32;6,-９.５4,5４,3.1４］; （1）输出a在[10，25]范围内的全部元素。 k=ｆｉｎd(ａ>1０＆a<2５） ａ(k） k = 1 ans＝ 23 (２)取出a前３行构成矩阵b，前两列构成矩阵c,右下角3＊2子矩阵构成矩阵ｄ,ｂ与c的乘积构成矩阵e。 ｂ＝a(1:3,:) b = ２3.0000 １０.0000 -0.77８0 0 41.0000 -45.0000 65．0０0０５.00００ ３2．０000 5.0000 0 32.0000 c=a（:,1:2) c = 23.00００ 10.０000 41.0０00 -4５.0000 32.00０0 ５.00００ ６.0000 -9.54０0 d=a(2:4,3：4） d = 65．000０ 5．0０0０ 0 32.00００ ５４.00０0 3．１４00 e=b*c e = 1.0e+003 * ０.9141 -0.2２39

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 期末上机考试试题带答案版 姓名： 学号： 成绩： 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 下面的MATLAB 语句中不正确的有： a) 2a ＝pi; b) record_1=3+4i c) a=2.0, d) c=1+6j 2. 已知水的黏度随温度的变化公式如下，其中a=0.03368，b=0.000221，计算温度t 为20，30，40度时的粘度分别是： 02 1at bt μμ=++0μ为0℃水的黏度，值为3 1.78510-?；a 、b 为常数，分别为0.03368、0.000221。 3. 请补充语句以画出如图所示的图形： [x,y]=meshgrid(-2:0.1:2, -2:0.1:2); Z=x.*exp(-x.^2-y.^2); ; a) Plot3(x,y,Z) b) plot3(x,y,Z) c) mesh(x,y,Z) d) plot3(x,y,z) 2 a) 0.4900 1.2501 0.8560 b) 0.8560 1.2501 0.4900 c) -0.6341 3.8189 -3.7749 d) 3.8189 -3.7749 2.8533 解释说明： >> x=0.5:0.5:3.0; >> y=[1.75,2.45,3.81,4.80,8.00,8.60]; >> a=polyfit(x,y,2)

a = 0.4900 1.2501 0.8560 >> x1=[0.5:0.25:3.0]; >> y1=a(1)*x1.^2+a(2)*x1+a(3) >> plot(x,y,'*') >> hold on >> plot(x1,y1,'--r') 5. 求方程在x=0.5附近的根. 21x x += a) 0.6180 b) -1.1719e-25 c) －1 d) -1.6180 6. 用Newton-Cotes 方法计算如下积分 1 5x ? （a ）133.6625 (b) 23.8600 (c) 87.9027 (d) -1.6180 7. y=ln(1+x)，求x=1时y" a) -0.25 b) 0.5 c) -0.6137 d) -1.6137 8. 某公司用3台轧机来生产规格相同的铝合金薄板。取样测量薄板的 厚度，精确至‰厘米。得结果如下： 轧机1：0.236 0.238 0.248 0.245 0.243 轧机2：0.257 0.253 0.255 0.254 0.261 轧机3：0.258 0.264 0.259 0.267 0.262 计算方差分析结果，并判定各台轧机所生产的薄板的厚度有无显著的差异？ a) p ＝1.3431e-005，没有显著差异。 b) p ＝0.9688，没有显著差异。 c) p ＝0.4956，有显著差异。 d) p ＝0.9688，有显著差异。 22x y x y e x y e --?-=??-+=??

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

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.请实现下图：

中国石油大学(北京)自动控制原理胡寿松 课件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、作出函数y=x4-4x3+3x+5 （x [0,6]）的图形，用小红点标出其在[0,6]之间的最小值点，并在最小值点附近标出该最小值点的坐标值； 程序： function f=myfun(x) f=x.^4-4*x.^3+3*x+5; x=linspace(0,6,100); y=x.^4-4*x.^3+3*x+5; x1=fminbnd(@myfun,0,6) y1=myfun(x1) 结果： x1 = 2.9115 y1 = -13.1300 plot(x,y,x1,y1,'r*') text(x1,y1,'2.9115,-13.1300'); 0123456 2、某公司有一批以每桶2元购进的彩漆，为了获得较高的利润，希望以较高的价格卖出，但价格越高，售出量就越少，二者之间的关系由表一给出。于是打算增加广告投入来促销。而广告费与销售量的关系可由销售增长因子来描述。例如，投入3万元的广告费，销售因子为1.85，意味着做广告后的销售量将是未做广告销售量的1.85倍。根据经验，广告费与销售因子的关系如表2，现请你作出决策：投入多少广告费和售价为多少时所获得的利润最大?

表1 表2 彩漆的销售量 摘要 在经济学中，某种产品的销售量与产品自身的价格存在着负相关关系，即产品价格上升会导致产品的销售量减少，产品价格下降会导致产品的销售量增加。与此同时，广告宣传对产品的销售量也是影响深远的。对一个企业而言，广告费既不是越少越好，也不是多多益善。广告活动的规模和广告费用的大小，应与企业的生产和流通规模相适应，在发展中求节约。 为研究产品销售量与售价和广告费用的关系，我们收集了某售价与预期销售量和广告费与销售增长因子的一些数据（见附录一），并建立了预期销售量1y 与售价1x 的线性模型： 11^ 1333.54222 .50x y -= 销售增长因子2y 与广告费2x 的二次函数模型： 22^ 0409.00188 .1x y +=—2 2 0004.0x 利润p 与售价1x 和销售增长因子2x 的模型： 2 2 2211)0004.00409.00188 .1)(1333.54222 .50)(2(x x x x x p --+--= 关键字：预期销售量 广告费 销售增长因子 线性回归

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试卷+答案

学习中心/函授站_ 姓名学号 西安电子科技大学网络与继续教育学院 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、 构造二阶常微分边值问题: 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 ?=-+=-<