您的位置：360文档中心› MTLB实验练习题计算机南邮MATLAB数学实验大作业答案

MTLB实验练习题计算机南邮MATLAB数学实验大作业答案

MTLB实验练习题计算机南邮MATLAB数学实验大作业答案
MTLB实验练习题计算机南邮MATLAB数学实验大作业答案

“M A T L A B ”练习题

要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。

1、求230x e x -=的所有根。（先画图后求解）（要求贴图）

>> solve('exp(x)-3*x^2',0)

ans =

-2*lambertw(-1/6*3^(1/2))

-2*lambertw(-1,-1/6*3^(1/2))

-2*lambertw(1/6*3^(1/2))

2、求下列方程的根。

1) 5510x x ++=

a=solve('x^5+5*x+1',0);a=vpa(a,6)

a =

1.10447+1.05983*i

-1.00450+1.06095*i

-.

-1.00450-1.06095*i

1.10447-1.05983*i

2）

1

sin0

2

x x-=至少三个根

>> fzero('x*sin(x)-1/2', 3)

ans =

2.9726

>> fzero('x*sin(x)-1/2',-3)

ans =

-2.9726

>> fzero('x*sin(x)-1/2',0)

ans =

-0.7408

3）2

sin cos0

x x x

-=所有根

>> fzero('sin(x)*cos(x)-x^2',0) ans =

>> fzero('sin(x)*cos(x)-x^2',0.6)

0.7022

3、求解下列各题：

1）3

0sin lim

x x x

x ->-

>> sym x;

>> limit((x-sin(x))/x^3)

ans =

1/6

2） (10)cos ,x y e x y =求

>> sym x;

>> diff(exp(x)*cos(x),10)

ans =

(-32)*exp(x)*sin(x)

3）2

1/2

0(17x e dx ?精确到位有效数字）

>> sym x;

>> vpa((int(exp(x^2),x,0,1/2)),17)

0.

4）4

2

254x dx x

+? >> sym x;

>> int(x^4/(25+x^2),x)

ans =

125*atan(x/5) - 25*x + x^3/3

5）求由参数方程arctan x y t

??=?=??dy dx 与二阶导数22d y dx 。

>> sym t;

>> x=log(sqrt(1+t^2));y=atan(t);

>> diff(y,t)/diff(x,t)

ans =

1/t

6）设函数y =f (x )由方程xy +e y = e 所确定，求y ′(x )。

>> syms x y;

f=x*y+exp(y)-exp(1); >> -diff(f,x)/diff(f,y) ans =

-y/(x + exp(y))

7）

sin2

x

e xdx

+∞-

?

>> syms x;

>> y=exp(-x)*sin(2*x);

>> int(y,0,inf)

ans =

2/5

8）08

x=展开（最高次幂为）

>> syms x

f=sqrt(1+x);

taylor(f,0,9)

ans =

- (429*x^8)/32768 + (33*x^7)/2048 - (21*x^6)/1024 + (7*x^5)/256 - (5*x^4)/128 +

x^3/16 - x^2/8 + x/2 + 1

9） 1sin

(3)(2)x

y e

y =求

>> syms x y;

>> y=exp(sin(1/x));

>> dy=subs(diff(y,3),x,2)

dy =

-0.5826

10）求变上限函数2

x x ?对变量x 的导数。

>> syms a t;

>> diff(int(sqrt(a+t),t,x,x^2))

Warning: Explicit integral could not be found.

ans =

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

4、求点(1,1,4)到直线L :

31

102

x y z --==

- 的距离 >> M0=[1,1,4];M1=[3,0,1];M0M1=M1-M0;

v=[-1,0,2];

d=norm(cross(M0M1,v))/norm(v)

d =

1.0954

5、已知22

()2(),

x f x μσ--=

分别在下列条件下画出()f x 的图形：（要求贴图）

(1)1,011σμ=时＝，－，，在同一坐标系里作图

>> syms x;

>> fplot('(1/sqrt(2*pi))*exp(-((x)^2)/2)',[-3,3],'r')

>> hold on

>> fplot('(1/sqrt(2*pi))*exp(-((x-1)^2)/2)',[-3,3],'y')

>> hold on

>> fplot('(1/sqrt(2*pi))*exp(-((x+1)^2)/2)',[-3,3],'g')

>> hold off

(2)0,124μσ＝时＝，，，在同一坐标系里作图。

>> syms x;

fplot('(1/sqrt(2*pi))*exp(-((x)^2)/2)',[-3,3],'r')

hold on

fplot('(1/(sqrt(2*pi)*2))*exp(-((x)^2)/(2*2^2))',[-3,3],'y') hold on

fplot('(1/(sqrt(2*pi)*4))*exp(-((x)^2)/(2*4^2))',[-3,3],'g') hold off

6、画下列函数的图形：（要求贴图）

（1）

sin

020

cos

02 4

x u t

t

y u t

u

t

z

?

?=

≤≤

?

=

?

≤≤?

?=

?

>> ezmesh('u*sin(t)','u*cos(t)','t/4',[0,20,0,2]) (2) sin()03,03

z xy x y

=≤≤≤≤

>> x=0:0.1:3;y=x;

[X Y]=meshgrid(x,y);

Z=sin(X*Y);

>> mesh(X,Y,Z)

（3）

sin(3cos)

02 cos(3cos)

02 sin

x t u

t

y t u

u

z u

π

π=+

?

≤≤

?

=+

?

≤≤

?=

?

ezmesh('sin(t)*(3+cos(u))','cos(t)*(3+cos(u))','sin(u)',[0,2*pi,0,2*pi])

7、 已知422134305,203153211A B -????

? ?=-=-- ? ? ? ?-????

，在MATLAB 命令窗口中建立A 、B 矩阵并对其进

行以下操作：

(1) 计算矩阵A 的行列式的值det()A

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

>> det(A)

ans =

-158

(2) 分别计算下列各式：1122,*,.*,,,,T A B A B A B AB A B A A ---

>> A=[4,-2,2;-3,0,5;1,5,3];B=[1,3,4;-2,0,-3;2,-1,1];

>> 2*A-B

ans =

7 -7 0

-4 0 13

0 11 5

>> A*B

ans =

12 10 24

7 -14 -7

-3 0 -8

>> A.*B

ans =

4 -6 8

6 0 -15

2 -5 3

>> A*inv(B)

ans =

-0.0000 -0.0000 2.0000 -2.7143 -8.0000 -8.1429 2.4286 3.0000 2.2857 >> inv(A)*B

ans =

0.4873 0.4114 1.0000

0.3671 -0.4304 0.0000

-0.1076 0.2468 0.0000

>> A*A

ans =

24 2 4

-7 31 9

-8 13 36

>> A'

ans =

4 -3 1

-2 0 5

2 5 3

>>

8、在MATLAB中分别利用矩阵的初等变换及函数rank、函数inv求下列矩阵的秩：

(1)

1632

3540,

11124

A

-

??

?

=-

?

?

--

??

求 rank(A)=?

>> A=[1,-6,3,2;3,-5,4,0;-1,-11,2,4];

>> rank(A)

ans =

3

(2) 35011

200,10201

20

2B ?? ?

?

= ? ???

求1B -。 >> B=[3,5,0,1;1,2,0,0;1,0,2,0;1,2,0,2]

>> inv(B)

ans =

2.0000 -4.0000 -0.0000 -1.0000

-1.0000 2.5000 0.0000 0.5000

-1.0000 2.0000 0.5000 0.5000

0 -0.5000 0 0.5000

9、在MATLAB 中判断下列向量组是否线性相关，并找出向量组1(1132),T α=

234(1113),(5289),(1317)T T T ααα=--=-=-中的一个最大线性无关组。

>> a1=[1 1 3 2]'

a2=[-1 1 -1 3]'

a3=[5 -2 8 9]'

a4=[-1 3 1 7]'

A= [a1, a2 ,a3 ,a4] ;[R jb]=rref(A) a1 =

1

1

3

2

a2 =

-1

1

-1

3

a3 =

5

8

9

a4 =

-1

3

1

7

R =

1.0000 0 0 1.0909 0 1.0000 0 1.7879 0 0 1.0000 -0.0606 0 0 0 0 jb =

1 2 3

>> A(:,jb)

1 -1 5

1 1 -2

3 -1 8

2 3 9

10、在MATLAB中判断下列方程组解的情况，若有多个解，写出通解。

（1）

1234

1234

1234

1234

420

20 3720 31260 x x x x

x x x x

x x x x

x x x x

-+-=?

?--+=

?

?

++-=?

?--+=?

一：

>> A=[1,-1,4,2;1,-1,-1,2;3,1,7,-2;1,-3,-12,6]; >> rank(A)

ans =

3

>> rref(A)

ans =

1 0 0 0

0 1 0 -2

0 0 1 0

0 0 0 0

二：

>> A=[1,-1,4,2;1,-1,-1,2;3,1,7,-2;1,-3,-12,6]; >> format rat

n=4;

RA=rank(A)

RA =

3

>> if(RA==n)

fprintf('%方程只有零解')

else

b=null(A,'r')

end

b =

0 2

1 >> syms k

X=k*b

X =

2*k

k

(2)

123

123

123

123

234

245 38213 496

x x x

x x x

x x x

x x x

++=?

?-+=-?

?

+-=?

?-+=-?

>> A=[2 3 1;1 -2 4;3 8 -2;4 -1 9]; b=[4 -5 13 -6]';

B=[A b];

>> n=3;

>> RA=rank(A)

RA =

2

>> RB=rank(B)

RB =

2

rref(B)

ans =

1 0

2 -1 0 1 -1 2 0 0 0 0 0 0 0 0 >> format rat

if RA==RB&RA==n %判断有唯一解

X=A\b

elseif RA==RB&RA

X=A\b %求特解

C=null(A,'r') %求AX=0的基础解系

else X='equition no solve' %判断无解

end

Warning: Rank deficient, rank = 2, tol = 8.9702e-015. X =

3/2

-1/2

C =

-2

1

1

11、求矩阵

211

020

413

A

-??

?

= ?

?

-??

的逆矩阵1

A-及特征值和特征向量。

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

>> a1=inv(A)

a1 =

-3/2 1/2 1/2 0 1/2 0 -2 1/2 1 >> [P,R]=eig(A)

P =

-985/1393 -528/2177 379/1257 0 0 379/419 -985/1393 -2112/2177 379/1257 R =

-1 0 0 0 2 0 0 0 2 A的三个特征值是：

Matlab上机实验答案

Matlab上机实验答案 实验一 MATLAB运算基础 1. 先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。 >> z1=2*sin(85*pi/180)/(1+exp(2)) z1 = >> x=[2 1+2i; 5]; >> z2=1/2*log(x+sqrt(1+x^2)) z2 = - + + -

>> a=::; >> z3=(exp.*a)-exp.*a))./2.*sin(a++log(+a)./2) (>> z33=(exp*a)-exp*a))/2.*sin(a++log(+a)/2)可以验证z3==z33，是否都为1） z3 = Columns 1 through 5 + + + + + Columns 6 through 10 + + + + + Columns 11 through 15 + + + + + Columns 16 through 20 + + + + +

Columns 21 through 25 + + + + + Columns 26 through 30 + + + + + Columns 31 through 35 + + + + + Columns 36 through 40 + + + + + Columns 41 through 45 + + + + + Columns 46 through 50

+ + + + + Columns 51 through 55 + + + + + Columns 56 through 60 + + + + + Column 61 + (4) 2 2 4 2 01 112 2123 t t z t t t t t ?≤< ? =-≤< ? ?-+≤< ? ，其中t=0:: >> t=0::; >> z4=(t>=0&t<1).*(t.^2)+(t>=1&t<2).*(t.^2-1)+(t>=2&t<3).*(t.^ 2-2.*t+1) z4 =

数学软件MATLAB实验作业

数学软件与数学实验作业 一.《数学软件》练习题（任选12题，其中19－24题至少选2题）： 3.对下列各式进行因式分解. (1). syms x y >> factor(x^5-x^3) (2). syms x y >> factor(x^4-y^4) (3). syms x >> factor(16-x^4) (4). syms x >> factor(x^3-6*x^2+11*x-6) (5). syms x y >> factor((x+y)^2-10*(x+y)+25) (6). syms x y >> factor(x^2/4+x*y+y^2) (7). syms x y a b >> factor(3*a*x+4*b*y+4*a*y+3*b*x) (8). syms x >> factor(x^4+4*x^3-19*x^2-46*x+120) 5.解下列方程或方程组. (1).solve('(y-3)^2-(y+3)^3=9*y*(1-2*y)') (2). solve('3*x^2+5*(2*x+1)') (3). solve('a*b*x^2+(a^4+b^4)*x+a^3*b^3','x') (4). solve('x^2-(2*m+1)*x+m^2+m','x') (5). [x,y]=solve('4*x^2-9*y^2=15','2*x-3*y=15') 6.计算极限. (1). syms x f=(exp(x)-exp(-x))/sin(x); limit(f,x,0) (2) syms x >> f=(x/(x-1)-1/log(x)); >> limit(f,x,1) (3). syms x >> f=(1-cos(x))/x^2; >> limit(f,x,0)

MATLAB实验练习题(计算机)-南邮-MATLAB-数学实验大作业答案

“”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> ('(x)-3*x^2',0) = -2*(-1/6*3^(1/2)) -2*(-11/6*3^(1/2)) -2*(1/6*3^(1/2)) 3、求解下列各题： 1）30 sin lim x x x x ->- >> x;

>> (((x))^3) = 1/6 2） (10)cos ,x y e x y =求 >> x; >> ((x)*(x),10) = (-32)*(x)*(x) 3）2 1/2 0(17x e dx ?精确到位有效数字） >> x; >> ((((x^2),0,1/2)),17) =

0.54498710418362222 4）4 2 254x dx x +? >> x; >> (x^4/(25^2)) = 125*(5) - 25*x + x^3/3 5）求由参数方程arctan x y t ??=? =??dy dx 与二阶导 数22 d y dx 。 >> t; >> ((1^2))(t); >> ()() = 1

6）设函数(x)由方程e所确定，求y′(x)。>> x y; *(y)(1); >> ()() = (x + (y)) 7） sin2 x e xdx +∞- ? >> x; >> ()*(2*x); >> (y,0) = 2/5

8） 08x =展开（最高次幂为） >> x (1); taylor(f,0,9) = - (429*x^8)/32768 + (33*x^7)/2048 - (21*x^6)/1024 + (7*x^5)/256 - (5*x^4)/128 + x^3/16 - x^2/8 + 2 + 1 9） 1sin (3)(2)x y e y =求 >> x y; >> ((1)); >> ((y,3),2) =

数学实验答案-1

1.（1） [1 2 3 4;0 2 -1 1;1 -1 2 5;]+(1/2).*([2 1 4 10;0 -1 2 0;0 2 3 -2]) 2. A=[3 0 1;-1 2 1;3 4 2],B=[1 0 2;-1 1 1;2 1 1] X=(B+2*A)/2 3. A=[-4 -2 0 2 4;-3 -1 1 3 5] abs(A)>3 % 4. A=[-2 3 2 4;1 -2 3 2;3 2 3 4;0 4 -2 5] det(A),eig(A),rank(A),inv(A) 求计算机高手用matlab解决。 >> A=[-2,3,2,4;1,-2,3,2;3,2,3,4;0,4,-2,5] 求|A| >> abs(A) ans = （ 2 3 2 4 1 2 3 2 3 2 3 4 0 4 2 5 求r（A） >> rank(A) ans =

4 求A-1 《 >> A-1 ans = -3 2 1 3 0 -3 2 1 2 1 2 3 -1 3 -3 4 求特征值、特征向量 >> [V,D]=eig(A) %返回矩阵A的特征值矩阵D 与特征向量矩阵V , V = - + + - - + - + - + - + D = { + 0 0 0 0 - 0 0 0 0 + 0 0 0 0 - 将A的第2行与第3列联成一行赋给b >> b=[A(2,:),A(:,3)'] b = 《 1 - 2 3 2 2 3 3 -2

1． a=round(unifrnd(1,100)) i=7; while i>=0 i=i-1; b=input('请输入一个介于0到100的数字：'); if b==a ￥ disp('You won!'); break; else if b>a disp('High'); else if b

MATLAB上机实验(答案)

MATLAB工具软件实验(1) （1）生成一个4×4的随机矩阵，求该矩阵的特征值和特征向量。程序： A=rand(4) [L,D]=eig(A) 结果： A = 0.9501 0.8913 0.8214 0.9218 0.2311 0.7621 0.4447 0.7382 0.6068 0.4565 0.6154 0.1763 0.4860 0.0185 0.7919 0.4057 L = -0.7412 -0.2729 - 0.1338i -0.2729 + 0.1338i -0.5413 -0.3955 -0.2609 - 0.4421i -0.2609 + 0.4421i 0.5416 -0.4062 -0.0833 + 0.4672i -0.0833 - 0.4672i 0.4276 -0.3595 0.6472 0.6472 -0.4804 D = 2.3230 0 0 0 0 0.0914 + 0.4586i 0 0 0 0 0.0914 - 0.4586i 0 0 0 0 0.2275 （2）给出一系列的a值，采用函数 22 22 1 25 x y a a += - 画一组椭圆。 程序： a=0.5:0.5:4.5; % a的绝对值不能大于5 t=[0:pi/50:2*pi]'; % 用参数t表示椭圆方程 X=cos(t)*a; Y=sin(t)*sqrt(25-a.^2); plot(X,Y) 结果： （3）X=[9,2,-3,-6,7,-2,1,7,4,-6,8,4,0,-2], （a）写出计算其负元素个数的程序。程序： X=[9,2,-3,-6,7,-2,1,7,4,-6,8,4,0,-2]; L=X<0; A=sum(L) 结果： A =

南邮课内实验-运筹学-线性规划-第一次0407

南邮课内实验-运筹学-线性规划-第一次0407

课内实验报告 课程名：运筹学 任课教师：邢光军 专业： 学号： 姓名：

/学年第学期 南京邮电大学管理学院

实验背景：某商场是个中型的百货商场，它对售货人员的需求经过统计分析如表1所示。 休息两天，并要求休息的两天是连续的，问应该如何安排售货人员的作息，既满足了工作需要，又使配备的售货人员人数最少？ 实验结果：一：问题分析和建立模型： 解：设xi表示星期i开始上班的售货人员数，建立如下求解模型：目标函数：Min f（x）=x1+x2+x3+x4+x5+x6+x7 约束条件：s.t. X3+x4+x5+x6+x7≥28 X1+x4+x5+x6+x7≥15 X1+x2+x5+x6+x7≥24 X1+X2+x3+x6+x7≥25 X1+X2+X3+x4+x7≥19 X1+X2+X3+X4+x5≥31 X2+X3+X4+X5+X6≥28 二：计算过程： 下面利用Spreadsheet来求解该问题： 在Excel2003版本中，单击“工具”栏中“加载宏”命令，在弹出的的“加载宏”对话框选择“规划求解”，在“工具”下拉菜单中会增加“规划求解”命

令，这样就可以使用了。 1、将求解模型及数据输入至Spreadsheet工作表中。 在工作表中的B1~H1单元格分别输入x1,x2,x3,x4,x5,x6,x7,B2~H2单元格分别表示决策变量的取值。B3~H10单元格数据为技术系数矩阵，I3~I10单元格值为目标函数及约束1~7不等式符号左边部分，如I3=SUMPRODUCT（B3：H3，B2：H2），即I3=1*x1+1*x2+1*x3+1*x4+1*x5+1*x6+1*x7,其余I4~I10含义雷同。K4~K10单元格数据为约束1~7不等式符号右端系数。（如图①） 图① 2、单击“工具”菜单中的“规划求解”命令，弹出“规划求解参数”对话框。在“规划求解参数”对话框中设置目标单元格为I3，选中“最小值”前的单选按钮，设置可变单元格为B2：H2。单击“规划求解参数”对话框中的“添加”按钮，打开“添加约束”对话框，单击单元格引用位置文本框，然后选定工作表的I4单元格，则在文本框中显示“$I$4”,选择“>=”的约束条件，在约束值文本框中输入K4单元格，则在文本框中显示“$K$4”。单击“添加”按钮，把所有的约束条件都添加到“规划求解参数”对话框的“约束”列表框中。其余6条约束不等式的输入方法雷同。按照同样的方法继续输入决策变量的非负约束、整数约束。（如图②） 图② 3、在“规划求解参数”对话框中单击“求解”按钮，弹出“规划求解结果”对话框，选中“保存规划求解结果”前的单选按钮，单击“确定”按钮，工作表中就显示规划求解的结果。（如图③）

山东建筑大学数学实验期末作业matlab

数学实验 期 末 作 业 学号： 班级： 姓名：

1. 求函数x x y 2sin 3=的5阶导数。 2. 使用sparse 命令描述? ? ???? ? ? ??30001 020******* 01020 10003。 3. 求解边值问题 1)0(,0)0(,34,43==+-=+=g f g f dx dg g f dx df 。 4. 建立函数1 2sin )(3-=x x f x 的M-文件，并计算)2(f 和)10(f 。 5. 计算二重积分dy dx x y ??211 0][。 6. 已知数列满足2,11 01=+= +a ka a k k ，求5a ，并要求最后结果分别以小数点后两位和有理数这两种数据显示格式输出。

7. 大约在1500年前，《孙子算经》中就记载了这个有趣的问题：“今有鸡兔同笼，上有三十五头，下有九十四足，问鸡兔各几何？”请根据你的思路编程求解。 8. 绘制以下方程所表示的图形。 （1）x x y -=23 2 （2）y z cos =绕z 轴的旋转曲面 （3）)40(,) 2sin(sin )]2cos(4[cos )]2cos(4[π<

10.根据中华人民共和国个人所得税法规定：公民的个人工资、薪金应依法缴纳个人所得税。所得税计算办法为：在每个人的月收入中超过2000元以上的部分应该纳税，这部分收入称为应纳税所得额。应纳税所得额实行分段累计税率，按下列税率表计算： 个人所得税税率表： 等级全月应纳税所得额税率（%） 1 不超过500元的部分 5 2 超过500元，不到2000元的部分10 3 超过2000元，不到5000元的部分15 4 超过5000元，不到20000元的部分20 5 超过20000元，不到40000元的部分25 6 超过40000元，不到60000元的部分30 7 超过60000元，不到80000元的部分35 8 超过80000元，不到100000元的部分40 9 超过100000元的部分45 若某人的工资是x元，试建立税款y与收入x之间的M-文件，并要求程序运行时可以告知操作者“please input the number of your wage”。

MATLAB实验练习题(计算机) 南邮 MATLAB 数学实验大作业答案

“MATLAB”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> solve('exp(x)-3*x^2',0) ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 2、求下列方程的根。 1) 5510x x ++= a=solve('x^5+5*x+1',0);a=vpa(a,6)

1.10447+1.05983*i -1.00450+1.06095*i -.199936 -1.00450-1.06095*i 1.10447-1.05983*i 2） 1 sin0 2 x x-=至少三个根 >> fzero('x*sin(x)-1/2', 3) ans = 2.9726 >> fzero('x*sin(x)-1/2',-3) ans = -2.9726 >> fzero('x*sin(x)-1/2',0) ans = -0.7408

3）2sin cos 0x x x -= 所有根 >> fzero('sin(x)*cos(x)-x^2',0) ans = >> fzero('sin(x)*cos(x)-x^2',0.6) ans = 0.7022 3、求解下列各题： 1）30sin lim x x x x ->- >> sym x; >> limit((x-sin(x))/x^3) ans = 1/6 2） (10)cos ,x y e x y =求 >> sym x; >> diff(exp(x)*cos(x),10) ans =

数学实验4答案

第四次练习题 1、 编程找出 5,1000+=≤b c c 的所有勾股数，并问：能否利用通项表示 },,{c b a ? >> for b=1: 995 a=sqrt((b+5)^2-b^2); if(a==floor(a)) fprintf('a=%i,b=%i,c=%i\n',a,b,b+5) end end a=15,b=20,c=25 a=25,b=60,c=65 a=35,b=120,c=125 a=45,b=200,c=205 a=55,b=300,c=305 a=65,b=420,c=425 a=75,b=560,c=565 a=85,b=720,c=725 a=95,b=900,c=905 >> for c=6:1000 a=sqrt(c^2-(c-5)^2); if(a==floor(a)) fprintf('a=%i,b=%i,c=%i\n',a,c-5,c) end end a=15,b=20,c=25 a=25,b=60,c=65 a=35,b=120,c=125 a=45,b=200,c=205 a=55,b=300,c=305 a=65,b=420,c=425 a=75,b=560,c=565 a=85,b=720,c=725 a=95,b=900,c=905 {a,b,c}={100*n^2-100*n+25,10*n^2-10*n,10*n^2-10*n+5} 2、编程找出不定方程 )35000(122<-=-y Dy x 的所有正整数解。（学号为单号的取D=2, 学号为双号的取D=5） D=2（学号为单号） >> for y=1:34999 x=sqrt(2*y^2-1); if(x==floor(x)) fprintf('x=%i,y=%i\n',x,y) end

河南城建学院MATLAB上机实验答案

一熟悉Matlab工作环境 1、熟悉Matlab的5个基本窗口 思考题： （1）变量如何声明，变量名须遵守什么规则、是否区分大小写。 答：变量一般不需事先对变量的数据类型进行声明，系统会依据变量被赋值的类型自动进行类型识别，也就是说变量可以直接赋值而不用提前声明。变量名要遵守以下几条规则：?变量名必须以字母开头，只能由字母、数字或下划线组成。 ?变量名区分大小写。 ?变量名不能超过63个字符。 ?关键字不能作为变量名。 ?最好不要用特殊常量作为变量名。 （2）试说明分号、逗号、冒号的用法。 分号：分隔不想显示计算结果的各语句；矩阵行与行的分隔符。 逗号：分隔欲显示计算结果的各语句；变量分隔符；矩阵一行中各元素间的分隔符。 冒号：用于生成一维数值数组；表示一维数组的全部元素或多维数组某一维的全部元素。 （3）linspace（）称为“线性等分”函数，说明它的用法。 LINSPACE Linearly spaced vector. 线性等分函数 LINSPACE(X1, X2) generates a row vector of 100 linearly equally spaced points between X1 and X2. 以X1为首元素，X2为末元素平均生成100个元素的行向量。 LINSPACE(X1, X2, N) generates N points between X1 and X2. For N < 2, LINSPACE returns X2. 以X1为首元素，X2为末元素平均生成n个元素的行向量。如果n＜2，返回X2。 Class support for inputs X1,X2: float: double, single 数据类型：单精度、双精度浮点型。 （4）说明函数ones（）、zeros（）、eye（）的用法。 ones（）生成全1矩阵。 zeros（）生成全0矩阵。 eye（）生成单位矩阵。 2、Matlab的数值显示格式

matlab数学实验报告5

数学实验报告 制作成员班级学号 2011年6月12日

培养容器温度变化率模型 一、实验目的 利用matlab软件估测培养容器温度变化率 二、实验问题 现在大棚技术越来越好，能够将温度控制在一定温度范围内。为利用这种优势，实验室现在需要培植某种适于在8.16℃到10.74℃下能够快速长大的甜菜品种。为达到实验所需温度，又尽可能地节约成本，研究所决定使用如下方式控制培养容器的温度：1，每天加热一次或两次，每次约两小时； 2，当温度降至8.16℃时，加热装置开始工作；当温度达到10.74℃时，加热装置停止工作。 已知实验的时间是冬天，实验室为了其它实验的需要已经将实验室的温度大致稳定在0℃。下表记录的是该培养容器某一天的温度 时间(h)温度（℃）时间（h）温度（℃）09.68 1.849.31 0.929.45 2.959.13 3.878.981 4.989.65 4.988.811 5.909.41 5.908.691 6.839.18 7.008.5217.938.92 7.938.3919.048.66 8.978.2219.968.43 9.89加热装置工作20.848.22 10.93加热装置工作22.02加热装置工作10.9510.8222.96加热装置工作12.0310.5023.8810.59 12.9510.2124.9910.35 13.889.9425.9110.18 三、建立数学模型 1，分析：由物理学中的傅利叶传热定律知温度变化率只取决于温度

差，与温度本身无关。因为培养容器最低温度和最高温度分别是：8.16℃和10.74℃；即最低温度差和最高温度差分别是：8.16℃和10.74℃。而且，16.8/74.10≈1.1467，约为1，故可以忽略温度对温度变化率的影响2， 将温度变化率看成是时间的连续函数，为计算简单，不妨将温度变化率定义成单位时间温度变化的多少，即温度对时间连续变化的绝对值（温度是下降的），得到结果后再乘以一系数即可。 四、问题求解和程序设计流程1）温度变化率的估计方法 根据上表的数据，利用matlab 做出温度-时间散点图如下： 下面计算温度变化率与时间的关系。由图选择将数据分三段，然后对每一段数据做如下处理：设某段数据为{(0x ,0y ),(1x ,1y ),(2x , 2y ),…,(n x ,n y )},相邻数据中点的平均温度变化率采取公式： 温度变化率=（左端点的温度-右端点的温度）/区间长度算得即：v( 2 1i i x x ++)=(1+-i i y y )/(i i x x - +1). 每段首尾点的温度变化率采用下面的公式计算：v(0x )=(30y -41y +2y )/(2x -0x )v(n x )=(3n y -41+n y +2+n y )/(n x -2-n x )

南邮MATLAB数学实验答案(全)

第一次练习 教学要求：熟练掌握Matlab 软件的基本命令和操作，会作二维、三维几何图形，能够用Matlab 软件解决微积分、线性代数与解析几何中的计算问题。 补充命令 vpa(x,n) 显示x 的n 位有效数字，教材102页 fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形 在下面的题目中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin lim x mx mx x →-与3 sin lim x mx mx x →∞- syms x limit((902*x-sin(902*x))/x^3) ans = 366935404/3 limit((902*x-sin(902*x))/x^3,inf) ans = 0 1.2 cos 1000 x mx y e =，求''y syms x diff(exp(x)*cos(902*x/1000),2) ans = (46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算 22 11 00 x y e dxdy +?? dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1) ans = 2.1394 1.4 计算4 2 2 4x dx m x +? syms x int(x^4/(902^2+4*x^2)) ans = (91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求 syms x diff(exp(x)*cos(902*x),10) ans = -356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x) 1.6 0x =的泰勒展式(最高次幂为4).

MATLAB上机实验练习题答案

数学建模 MATLAB上机实验练习题 1、给出一个系数矩阵A[234;541;132]，U=[123]，求出线性方程组的一个精确解。 2、给出两组数据x=[00.30.81.11.62.3]’y=[0.820.720.630.600.550.50]’，我们可以简单的 认为这组数据在一条衰减的指数函数曲线上，y=C1+C2e-t通过曲线拟合求出这条衰减曲线的表达式，并且在图形窗口画出这条曲线，已知的点用*表示。 3、解线性方程 4、通过测量得到一组数据： 5、已知一组测量值 6、从某一个过程中通过测量得到： 分别采用多项式和指数函数进行曲线拟合。 7、将一个窗口分成四个子窗口，分别用四种方法做出多峰函数的表面图（原始数据法，临近 插值法，双线性插值法，二重三次方插值法） 8、在同一窗口使用函数作图的方法绘出正弦、余弦、双曲正弦、双曲余弦。分别使用不同的 颜色，线形和标识符。 9、下面的矩阵X表示三种产品五年内的销售额，用函数pie显示每种产品在五年内的销售额

占总销售额的比例，并分离第三种产品的切片。 X=19.322.151.6 34.270.382.4 61.482.990.8 50.554.959.1 29.436.347.0 10、对应时间矢量t，测得一组矢量y t00.30.8 1.1 1.6 2.3 y0.50.82 1.14 1.25 1.35 1.40 采用一个带有线性参数的指数函数进行拟合，y=a0+a1e-t+a2te-t，利用回归方法求出拟合函数，并画出拟合曲线，已知点用圆点表示。 11、请创建如图所示的结构数组（9分） 姓名编号指标 江明顺071023身高：176，体重：82 于越忠060134身高：168，体重：74 邓拓050839身高：182，体重：77 12、创建如图所示的元胞数组。（9分） 13、某钢材厂从1990年到2010年的产量如下表所示，请利用三次样条插值的方法计算1999年该钢材厂的产量，并画出曲线，已知数据用‘*’表示。要求写出达到题目要求的MATLAB 操作过程，不要求计算结果。 年份19901992199419961998200020022004200620082010 产量（万吨）75.99591.972105.711123.203131.669150.697179.323203.212226.505249.633256.344 14、在一次化学动力学实验中，在某温度下乙醇溶液中，两种化合物反应的产物浓度与反应时间关系的原始数据如下，请对这组数据进行三次多项式拟合，并画出拟合曲线，已知数据如下。 time=[2.55.07.510.013.017.020.030.040.050.060.070.0] res=[0.290.560.771.051.361.522.002.272.813.053.253.56]

数学实验(MATLAB版韩明版)5.1,5.3,5.5,5.6部分答案

练习 B的分布规律和分布函数的图形，通过观1、仿照本节的例子，分别画出二项分布()7.0,20 察图形，进一步理解二项分布的性质。 解：分布规律编程作图：>> x=0:1:20;y=binopdf(x,20,; >> plot(x,y,'*') 图像： y x 分布函数编程作图：>> x=0::20; >>y=binocdf(x,20, >> plot(x,y) 图像： 《

1 x 观察图像可知二项分布规律图像像一条抛物线，其分布函数图像呈阶梯状。 2、仿照本节的例子，分别画出正态分布()25,2N的概率密度函数和分布函数的图形，通过观察图形，进一步理解正态分布的性质。 解：概率密度函数编程作图：>> x=-10::10; >> y=normpdf(x,2,5); >> plot(x,y) 图像：

00.010.020.030.040.050.060.070.08x y 分布函数编程作图：>> x=-10::10; >> y=normcdf(x,2,5); ~ >> plot(x,y) 图像：

01x y 观察图像可知正态分布概率密度函数图像像抛物线，起分布函数图像呈递增趋势。 3、设()1,0~N X ，通过分布函数的调用计算{}11<<-X P ,{}22<<-X P , {}33<<-X P . 解：编程求解： >> x1=normcdf(1)-normcdf(-1),x2=normcdf(2)-normcdf(-2),x3=normcdf(3)-normcdf(-3) x1 = x2 = ） x3 = 即：{}6827.011=<<-X P ，{}9545.022=<<-X P ，{}9973.033=<<-X P . 4、设()7.0,20~B X ,通过分布函数的调用计算{}10=X P 与{}10> x1=binopdf(10,20,,x2=binocdf(10,20,-binopdf(10,20, x1 = x2 =

Matlab上机实验答案 (1)

Matlab上机实验答案 实验一MATLAB运算基础 1. 先求下列表达式的值,然后显示MATLAB工作空间的使用情况并保存全部变量。 >> z1=2*sin(85*pi/200)/(1+exp(2)) z1 = 0.2375 >> x=[2 1+2i;-0.45 5]; >> z2=1/2*log(x+sqrt(1+x^2)) z2 = 0.7120 - 0.0253i 0.8968 + 0.3658i 0.2209 + 0.9343i 1.2041 - 0.0044i 2.9,,2.9, 3.0

>> a=-3.0:0.1:3.0; >> z3=(exp(0.3.*a)-exp(-0.3.*a))./2.*sin(a+0.3)+log((0.3+a)./2) (>> z33=(exp(0.3*a)-exp(-0.3*a))/2.*sin(a+0.3)+log((0.3+a)/2)可以验证z3==z33,是否都为1） z3 = Columns 1 through 5 0.7388 + 3.2020i 0.7696 + 3.2020i 0.7871 + 3.2020i 0.7920 + 3.2020i 0.7822 + 3.2020i Columns 6 through 10 0.7602 + 3.2020i 0.7254 + 3.2020i 0.6784 + 3.2020i 0.6206 + 3.2020i 0.5496 + 3.2020i Columns 11 through 20 0.4688 + 3.2020i 0.3780 + 3.2020i 0.2775 + 3.2020i 0.2080 + 3.2020i 0.0497 + 3.2020i

南邮课程设计实验报告

课程设计I报告 题目：课程设计 班级：44 姓名：范海霞 指导教师：黄双颖 职称： 成绩： 通达学院 2015 年 1 月 4 日

一：SPSS的安装和使用 在PC机上安装SPSS软件，打开软件： 基本统计分析功能包括描述统计和行列计算，还包括在基本分析中最受欢迎的常见统计功能，如汇总、计数、交叉分析、分类比较、描述性统计、因子分析、回归分析及聚类分析等等。具体如下： 1.数据访问、数据准备、数据管理与输出管理； 2.描述统计和探索分析：频数、描述、集中趋势和离散趋势分析、分布分析与查看、正态性检验与正态转换、均值的置信区间估计； 3.交叉表：计数；行、列和总计百分比；独立性检验；定类变量和定序变量的相关性测度； 4.二元统计：均值比较、T检验、单因素方差分析； 5.相关分析：双变量相关分析、偏相关分析、距离分析； 6.线性回归分析：自动线性建模、线性回归、Ordinal回归—PLUM、曲线估计； 7.非参数检验：单一样本检验、双重相关样本检验、K重相关样本检验、双重独立样本检验、K重独立样本检验； 8.多重响应分析：交叉表、频数表； 9.预测数值结果和区分群体：K-means聚类分析、分级聚类分析、两步聚类分析、快速聚类分析、因子分析、主成分分析、最近邻元素分析； 10. 判别分析； 11.尺度分析； 12. 报告：各种报告、记录摘要、图表功能（分类图表、条型图、线型图、面积图、高低图、箱线图、散点图、质量控制图、诊断和探测图等）； 13.数据管理、数据转换与文件管理； 二．数据文件的处理 SPSS数据文件是一种结构性数据文件，由数据的结构和数据的内容两部分构成，也可以说由变量和观测两部分构成。定义一个变量至少要定义它的两个属性，即变量名和变量类型其他属性可以暂时采用系统默认值，待以后分析过程中如果有需要再对其进行设置。在spss数据编辑窗口中单击“变量视窗”标签，进入变量视窗界面，即可对变量的各个属性进行设置。 1．创建一个数据文件数据 （1）选择菜单【文件】→【新建】→【数据】新建一个数据文件，进入数据编辑窗口。窗口顶部标题为“PASW Statistics数据编辑器”。 （2）单击左下角【变量视窗】标签进入变量视图界面，根据试验的设计定义每个变量类型。

Matlab编程与应用习题和一些参考答案

Matlab 上机实验一、二 3.求下列联立方程的解???????=+-+-=-+=++-=--+4 1025695842475412743w z y x w z x w z y x w z y x >> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10]; >> b=[4;4;9;4]; >> c=a\b 4．设???? ??????------=81272956313841A ，??????????-----=793183262345B ，求C1=A*B’;C2=A’*B;C3=A.*B,并求上述所有方阵的逆阵。 >> A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8]; >> B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7]; >> C1=A*B' >> C2=A'*B >> C3=A.*B >> inv(C1) >> inv(C2) >> inv(C3) 5．设 ????? ?++=)1(sin 35.0cos 2x x x y ，把x=0~2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线。 >> x=linspace(0,2*pi,101); >> y=cos(x)*(0.5+(1+x.^2)\3*sin(x)); >> plot(x,y,'r') 6．产生8×6阶的正态分布随机数矩阵R1, 求其各列的平均值和均方差。并求该矩阵全体数的平均值和均方差。 （mean var ） a=randn(8,6) mean(a) var(a) k=mean(a) k1=mean(k) i=ones(8,6) i1=i*k1 i2=a-i1 i3=i2.*i2 g=mean(i3) g2=mean(g) 或者 u=reshape(a,1,48); p1=mean(u)

matlab实验报告

Matlab实验报告 ——定积分的近似计算 学生姓名： 学号： 专业：数学与应用数学专业

数学实验报告 实验序号:1001114030 日期：2012年10月20日 班级应一姓名陈璐学号1001114030 实验名称：定积分的近似运算 问题背景描述： 利用牛顿—莱布尼茨公式虽然可以精确地计算定积分的值，但它仅适合于被积分函数的原函数能用初等函数表达出来的情形。如果这点办不到或不容易办到， 这就有必要考虑近似计算的方法。在定积分的很多应用问题中，被积函数甚至没 有解析表达式，可能只是一条实验记录曲线，或者是一组离散的采样值，这时只 能应用近似方法去计算相应的定积分。 实验目的： 本实验将主要研究定积分的三种近似计算算法：矩形法、梯形法、抛物线发。对于定积分的近似数值计算，Matlab有专门函数可用。 实验原理与数学模型： 1.sum（a）：求数组a的和。 2.format long：长格式，即屏幕显示15位有效数字。 3.double（）：若输入的是字符则转化为相应的ASCII码；若输入的是整型数之则转化为 相应的实型数值。 4.quad（）：抛物线法求数值积分。格式：quad（fun，a，b）。此处的fun是函数，并且

为数值形式，所以使用*、/、^等运算时要在其前加上小数点。 5.trapz():梯形法求数值积分。格式：trapz（x，y）。其中x为带有步长的积分区间;y为数 值形式的运算。 6.fprintf（文件地址，格式，写入的变量）：把数据写入指定文件。 7.syms 变量1变量2……：定义变量为符号。 8.sym（'表达式'）：将表达式定义为符号。 9.int（f，v，a，b）：求f关于v积分，积分区间由a到b。 10.subs（f，'x'，a）：将a的值赋给符号表达式f中的x，并计算出值。若简单地使用subs （f），则将f的所有符号变量用可能的数值代入，并计算出值。 实验所用软件及版本：Matlab 7.0.1

MTLB实验练习题计算机南邮MATLAB数学实验大作业答案

“M A T L A B ”练习题 要求：抄题、写出操作命令、运行结果，并根据要求，贴上运行图。 1、求230x e x -=的所有根。（先画图后求解）（要求贴图） >> solve('exp(x)-3*x^2',0) ans = -2*lambertw(-1/6*3^(1/2)) -2*lambertw(-1,-1/6*3^(1/2)) -2*lambertw(1/6*3^(1/2)) 2、求下列方程的根。 1) 5510x x ++= a=solve('x^5+5*x+1',0);a=vpa(a,6) a = 1.10447+1.05983*i -1.00450+1.06095*i -. -1.00450-1.06095*i

1.10447-1.05983*i 2） 1 sin0 2 x x-=至少三个根 >> fzero('x*sin(x)-1/2', 3) ans = 2.9726 >> fzero('x*sin(x)-1/2',-3) ans = -2.9726 >> fzero('x*sin(x)-1/2',0) ans = -0.7408 3）2 sin cos0 x x x -=所有根 >> fzero('sin(x)*cos(x)-x^2',0) ans = >> fzero('sin(x)*cos(x)-x^2',0.6)

0.7022 3、求解下列各题： 1）3 0sin lim x x x x ->- >> sym x; >> limit((x-sin(x))/x^3) ans = 1/6 2） (10)cos ,x y e x y =求 >> sym x; >> diff(exp(x)*cos(x),10) ans = (-32)*exp(x)*sin(x) 3）2 1/2 0(17x e dx ?精确到位有效数字） >> sym x; >> vpa((int(exp(x^2),x,0,1/2)),17)

MATLAB数学实验报告

Matlab 数学实验报告

一、实验目的 通过以下四组实验,熟悉MATLAB的编程技巧,学会运用MATLAB的一些主要功能、命令，通过建立数学模型解决理论或实际问题。了解诸如分岔、混沌等概念、学会建立Malthu模型和Logistic 模型、懂得最小二乘法、线性规划等基本思想。 二、实验内容 2.1实验题目一 2.1.1实验问题 Feigenbaum曾对超越函数y=λsin(πｘ)（λ为非负实数）进行了分岔与混沌的研究，试进行迭代格式x k+1=λsin(πx k),做出相应的Feigenbaum图 2.1.2程序设计 clear;clf; axis([0,4,0,4]); hold on for r=0:0.3:3.9 x=[0.1]; for i=2:150 x(i)=r*sin(3.14*x(i-1)); end pause(0.5) for i=101:150

plot(r,x(i),'k.'); end text(r-0.1,max(x(101:150))+0.05,['\it{r}=',num2str(r)]) end 加密迭代后 clear;clf; axis([0,4,0,4]); hold on for r=0:0.005:3.9 x=[0.1];

for i=2:150 x(i)=r*sin(3.14*x(i-1)); end pause(0.1) for i=101:150 plot(r,x(i),'k.'); end end 运行后得到Feigenbaum图

2.2实验题目二 2.2.1实验问题 某农夫有一个半径10米的圆形牛栏，长满了草。他要将一头牛拴在牛栏边界的桩栏上，但只让牛吃到一半草，问拴牛鼻子的绳子应为多长？ 2.2.2问题分析 如图所示，E为圆ABD的圆心，AB为拴牛的绳子，圆ABD为草场，区域ABCD为牛能到达的区域。问题要求区域ABCD等于圆ABC 的一半，可以设BC等于x,只要求出∠a和∠b就能求出所求面积。先计算扇形ABCD的面积，2a÷π×πx2=2aπ2，再求AB的面积，用扇形ABE的面积减去三角形ABE的面积即可。

相关主题
相关文档
最新文档