(完整word版)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

30

45

(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

5、 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

3、

>> 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课后答案-第九章

实验指导 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题目

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试卷合集附答案

一、填空题 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)

第二章 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试卷 一、选择题（每空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与科学计算样题 (加主观题答案)

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

中国石油大学(北京)自动控制原理胡寿松 课件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 ?=-+=-<

中南大学MATLAB课程设计

一、《MATLAB程序设计实践》Matlab基础 班级： 学号： 姓名： 表示多晶体材料织构的三维取向分布函数（f＝f（φ1，φ，φ2））是一个非常复杂的函数，难以精确的用解析函数表达，通常采用离散空间函数值来表示取向分布函数，Data.txt是三维取向分布函数的一个实例。由于数据量非常大，不便于分析，需要借助图形来分析。请你编写一个matlab程序画出如下的几种图形来分析其取向分布特征：（1）用Slice函数给出其整体分布特征； （2）用pcolor或contour函数分别给出(φ2＝0, 5, 10, 15, 20, 25, 30,

35 …90)切面上f分布情况(需要用到subplot函数)； (3) 用plot函数给出沿α取向线(φ1=0~90，φ＝45，φ2＝0)的f 分布情况。

流程图

解： （1）将文件Data.txt内的数据按照要求读取到矩阵f(phi1,phi,phi2)中，代码如下： fid=fopen('data.txt'); %读取数据文件Data.txt for i=1:18 tline=fgetl(fid); end phi1=1;phi=1;phi2=1;line=0; f=zeros(19,19,19); while ~feof(fid) tline=fgetl(fid); data=str2num(tline); line=line+1; if mod(line,20)==1 phi2=(data/5)+1; phi=1; else for phi1=1:19 f(phi1,phi,phi2)=data(phi1); end phi=phi+1; end end fclose(fid); 将以上代码保存为readtext.m文件并在MATLAB中运行，运行结果如下图所示：

MatLab实习心得体会(体会心得)

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

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