南邮2013MATLAB数学实验答案

第一次练习

教学要求：熟练掌握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((216*x-sin(216*x))/x^3) ans =

366935404/3

limit((216*x-sin(216*x))/x^3,inf) ans = 0 1.2 cos

1000

x

mx

y e =，求''y syms x

diff(exp(x)*cos(216*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/(216^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(216*x),10) ans =

-356485076957717053044344387763*cos(216*x)*exp(x)-3952323024277642494822005884*sin(216*x)*exp(x)

1.6 0x =的泰勒展式(最高次幂为4).

syms x

taylor((216/1000+x)^(1/3),5,x) ans =

-(9765625*451^(1/2)*500^(1/2)*x^4)/82743933602 +(15625*451^(1/2)*500^(1/2)*x^3)/91733851

-(125*451^(1/2)*500^(1/2)*x^2)/406802 + (451^(1/2)*500^(1/2)*x)/216 +(451^(1/2)*500^(1/2))/500 1.7 Fibonacci 数列{}n x 的定义是121,1x x ==12,(3,4,

)n n n x x x n --=+=用循环语句编

程给出该数列的前20项（要求将结果用向量的形式给出）。 x=[1,1]; for n=3:20

x(n)=x(n-1)+x(n-2); end x x=

Columns 1 through 10

1 1

2

3 5 8 13 21 3

4 5

5 Columns 11 through 20

89 144 233 377 610 987 1597 2584 4181 6765

1.8 对矩阵211020411000A m ?? ?- ?= ?

?- ?

?

?，求该矩阵的逆矩阵，特征值，特征向量，行列式，

计算6

A ，并求矩阵,P D （D 是对角矩阵），使得1

A PDP

-=。

A=[-2,1,1;0,2,0;-4,1,216/1000];inv(A)

ans =

0.4107 0.0223 -0.4554 0 0.5000 0 1.8215 -0.4554 -0.9107 eig(A) ans =

-0.5490 + 1.3764i -0.5490 - 1.3764i 2.0000 det(A) ans =

4.3920

[P,D]=eig(A)

P = %特征向量

0.3245 - 0.3078i 0.3245 + 0.3078i 0.2425 0 0 0.9701 0.8944 0.8944 0.0000 D =

-0.5490 + 1.3764i 0 0 0 -0.5490 - 1.3764i 0 0 0 2.0000 P*D^6*inv(P) %A^6的值 ans =

15.3661 12.1585 + 0.0000i -5.8531 0 64.0000 0 23.4124 -5.8531 + 0.0000i -1.6196

1.9 作出如下函数的图形（注：先用M 文件定义函数，再用fplot 进行函数作图）：

1202

()12(1)1

2

x x f x x x ?

≤≤??=??-<≤??

m 文件:

function y=fenduan(x)

if x<=1/2 y=2*x

else x<=1 y=2-2*x

end end

执行函数：fplot('fenduan',[0,1]); grid on title('1.9') 得下图：

0.0.0.0.0.0.0.0.0.第1.9题图

1.10 在同一坐标系下作出下面两条空间曲线（要求两条曲线用不同的颜色表示）

（1）cos sin x t y t z t =??=??=? （2）2cos 2sin x t y t z t =??

=??=?

t=-15:0.01:15; x1=cos(t); y1=sin(t); z1=t;

plot3(x1,y1,z1); hold on

x2=cos(2*t); y2=sin(2*t); z2=t;

plot3(x2,y2,z2,'m'); grid on

title('1.10') 得下图：

第1.10题图

1.11 已知422134305,203153211A B m -???? ? ?

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

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

对其进行以下操作：

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

(2) 分别计算下列各式：1122,*,.*,,,,T A B A B A B AB A B A A --- 解：A=[4,-2,2;-3,0,5;1,5*216,3]; B=[1,3,4;-2,0,3;2,-1,1];

det(A) ans =

-117288 2*A-B ans =

7 -7 0

-4 0 7 0 2161 5 A*B ans =

12 10 12

7 -14 -7 -9013 0 13537 A.*B ans =

4 -6 8

6 0 15 2 -4510 3 A*inv(B) ans =

1.0e+003 *

-0.0000 0 0.0020 0.0000 0.0016 0.0001 1.0311 -0.9016 -1.4167 inv(A)*B ans =

0.3463 0.5767 0.5383 0.0005 -0.0006 -0.0005 -0.1922 0.3460 0.9230 A*A ans =

24 9012 4

-7 22556 9 -13523 13528 22561 A' ans =

4 -3 1

-2 0 4510 2 5 3

1.12

已知22

()2()x f x μσ--=分别在下列条件下画出)(x f 的图形：

(1)/600m σ=,μ分别为0,1,1-(在同一坐标系上作图); (2)0μ=,σ分别为1,2,4,/100m (在同一坐标系上作图). (1)x=-4:0.1:4;

h=inline('1/sqrt(2*pi)/s*exp(-(x-mu).^2/(2*s^2))'); y1=h(0,216/400,x);y2=h(-1,216/400,x);y3=h(1,216/400,x); plot(x,y1,'b',x,y2,'m',x,y3,'y')

grid on

title('1.11(1)')

第1.12题

(2) z1=h(0,1,x);z2=h(0,3,x);z3=h(0,5,x); z4=h(0,216/150,x);

plot(x,z1,x,z2,'y',x,z3,'m',x,z4, 'g') grid on

title('1.11(2)')

z1=h(0,1,x);z2=h(0,3,x);z3=h(0,5,x); z4=h(0,216/150,x);

第1.12题

1.13 作出2

4

z mx y =+的函数图形。

x=-10:0.2:10;y=x;

[X Y]=meshgrid(x,y);Z=216*X.^2+Y.^4; mesh(X,Y,Z); title('1.13')

第1.13题

1.14对于方程5

0.10200

m

x x -

-=，先画出左边的函数在合适的区间上的图形，借助于软件中的方程求根的命令求出所有的实根,找出函数的单调区间,结合高等数学的知识说明函数为什么在这些区间上是单调的,以及该方程确实只有你求出的这些实根。最后写出你做此题的体会。 解：作图程序：（注：x 范围的选择是经过试探而得到的）

x=-1.7:0.02:1.7;y=x.^5-216/200*x-0.1;

plot(x,y);grid on; title('1.14')

第1.14题

由图形观察，在x=-1.5,x=0,x=1.5附近各有一个实根 solve('x^5-216/200*x-0.1') ans =

-1.4516870267499636199995749888894 -0.022172950190557703188753952167919 1.4627751059480654637229232196174

1.4573364935933870280941533926624*i + 0.0055424354962279297327028641499658 0.0055424354962279297327028641499658 - 1.4573364935933870280941533926624*i 三个实根的近似值分别为：-1.4517,-0.0222,1.4628

由图形可以看出，函数在区间(,1)-∞-单调上升，在区间(1,1)-单调下降，在区间(1,)∞单调上升。 syms x

diff('x^5-216/200*x-0.1',x)

结果为5*x^4-27/25

solve('5*x^4-216/200') ans =

1/5*3^(1/2)*15^(1/4) -1/5*3^(1/2)*15^(1/4) 1/5*i*3^(1/2)*15^(1/4) -1/5*i*3^(1/2)*15^(1/4) vpa(ans) ans =

.68173161988049962110253262670818 -.68173161988049962110253262670818 .68173161988049962110253262670818*i -.68173161988049962110253262670818*i 得到两个实根：-0.6817与0.6817

可以验证导函数在)9745

.0,(--∞内为正，函数单调上升 导函数在)9745.0,9745

.0(-内为负，函数单调下降 导函数在),9745

.0(∞内为正，函数单调上升 根据函数的单调性，最多有3个实根。

1.15 求23m 0x e x -=的所有根。（先画图后求解）（要求贴图）

作图命令：（注：x 范围的选择是经过试探而得到的） x=-5:0.001:15;y=exp(x)-3/1000*216*x.^2; plot(x,y);grid on; title('1.15(1)') 得到下图

进一步细化

6

第1.15题图

x=-0.05:0.0001:0.05;y=exp(x)-(3*216*x.^2)/1000; plot(x,y);grid on; title('1.15(2)')

x=10:0.001:15;y=exp(x)-(3*216*x.^2)/1000; plot(x,y);grid on; title('1.15(3)')

第1.15题

图

6

第1.15题图

可看出解在-0.02，0.02，13附近，进一步求得 fzero('exp(x)-3*216*x^2',0.02) ans =

0.0194

fzero('exp(x)-3*216*x^2',-0.02) ans =

-0.0190

fzero('exp(x)-3*216*x^2',13) ans =

13.0391

第二次练习

教学要求：要求学生掌握迭代、混沌的判断方法，以及利用迭代思想解决实际问题。

2.1 设11

()/23n n

n m x x x x +?

=+???=?，数列{}n x 是否收敛？若收敛，其值为多少？精确到8位有效

数字。

解：程序代码如下（m=216）：

f=inline('(x+216/x)/2'); x0=3; for i=1:20; x0=f(x0);

end

1 37.500000

2 21.630000

3 15.808065

4 14.735988

5 14.696990

6 14.696938

7 14.696938

8 14.696938

9 14.696938

10 14.696938

11 14.696938

12 14.696938

13 14.696938

14 14.696938

15 14.696938

16 14.696938

17 14.696938

18 14.696938

19 14.696938

20 14.696938由运行结果可以看出，，数列{}n x收敛，其值为14.696938。

2.2求出分式线性函数

2

12

1

(),()

x x m

f x f x

x m x m

-+

==

++

的不动点，再编程判断它们的迭代序

列是否收敛。

解：取m=1000.

（1）程序如下：

f=inline('(x-1)/(x+1000)');

x0=2;

for i=1:20;

x0=f(x0);

end

运行结果：

1,0.000998004 11,-0.001001

2,-0.000999001 12,-0.001001

3,-0.001001 13,-0.001001

4,-0.001001 14,-0.001001

5,-0.001001 15,-0.001001

6,-0.001001 16,-0.001001

7,-0.001001 17,-0.001001

8,-0.001001 18,-0.001001

9,-0.001001 19,-0.001001

10,-0.001001 20,-0.001001

由运行结果可以看出，，分式线性函数收敛，其值为-0.001001。易见函数的不动点为-0.001001（吸引点）。

（2）程序如下：

f=inline('(x+1000000)/(x+1000)');

x0=2;

for i=1:20;

x0=f(x0);

fprintf('%g,%g\n',i,x0);

end

运行结果：

1,998.006 11,618.332

2,500.999 12,618.302

3,666.557 13,618.314

4,600.439 14,618.309

5,625.204 15,618.311

6,615.692 16,618.31

7,619.311 17,618.311

8,617.929 18,618.31

9,618.456 19,618.31

10,618.255 20,618.31

由运行结果可以看出，，分式线性函数收敛，其值为618.31。易见函数的不动点为618.31（吸引点）。

2.3 下面函数的迭代是否会产生混沌？（56页练习7（1））

1202

()12(1)12

x x f x x x ?

≤≤??=??-<≤??

解：程序如下：

f=inline('1-2*abs(x-1/2)'); x=[]; y=[]; x(1)=rand();

y(1)=0;x(2)=x(1);y(2)=f(x(1)); for i=1:100; x(1+2*i)=y(2*i); x(2+2*i)=x(1+2*i); y(2+2*i)=f(x(2+2*i)); end

plot(x,y,'r'); hold on; syms x;

ezplot(x,[0,1/2]); ezplot(f(x),[0,1]); axis([0,1/2,0,1]); hold off 运行结果：

x

1 -

2 abs(x - 1/2)

2.4 函数()(1)(01)f x x x x α=-≤≤称为Logistic 映射，试从“蜘蛛网”图观察它取初值

为00.5x =产生的迭代序列的收敛性，将观察记录填人下表，若出现循环，请指出它的周期．（56页练习8）

解：当α=3.3时，程序代码如下： f=inline('3.3*x*(1-x)'); x=[]; y=[]; x(1)=0.5;

y(1)=0;x(2)=x(1);y(2)=f(x(1)); for i=1:1000; x(1+2*i)=y(2*i); x(2+2*i)=x(1+2*i); y(1+2*i)=x(1+2*i); y(2+2*i)=f(x(2+2*i)); end

plot (x,y,'r'); hold on;

syms x;

ezplot(x,[0,1]);

ezplot(f(x),[0,1]);

axis([0,1,0,1]);

hold off运行结果：

x

当α=3.5时，上述程序稍加修改，得：

x

当α=3.56时，得：

x

当α=3.568时，得：

x

当α=3.6时，得：

x

当α=3.84时，得：

x

2.5 对于Martin 迭代，取参数,,a b c 为其它的值会得到什么图形？参考下表（取自63页练

习13）

解：取m=10000；迭代次数N=20000；

在M-文件里面输入代码：

function Martin(a,b,c,N)

f=@(x,y)(y-sign(x)*sqrt(abs(b*x-c)));

g=@(x)(a-x);

m=[0;0];

for n=1:N

m(:,n+1)=[f(m(1,n),m(2,n)),g(m(1,n))];

end

plot(m(1,:),m(2,:),'kx');

axis equal

在命令窗口中执行Martin（10000，10000，10000，20000），得：

x 104

执行Martin（-10000，-10000，10000，20000），得：

x 104

执行Martin（-10000，10，-10000，20000），得：

执行Martin（10，10，0.5，20000），得：

执行Martin（10，10000，-10000，20000），得：

执行Martin（100，1000，-10，20000），得：

执行Martin（-1000，17，4，20000），得：

2.6 能否找到分式函数

2ax b

cx dx e

+++（其中,,,,a b c d e 是整数）,使它产生的迭代序列（迭代

(。如果迭代收敛,那么迭代的初值与收敛的速度有什么关系.写出你做此题的体会. 提示：教材54页练习4的一些分析。

若分式线性函数()ax b

f x cx d

+=+为()f x 的不动点，因此

=

化简得：(2)(0c b d a -+-=。

若,,,a b c d 为整数，易见2,b c d a ==。

取满足这种条件的不同的,,,a b c d 以及迭代初值进行编。 解：

2ax b cx dx e +++迭代收敛到指定的数3902，则3

902

为2ax b cx dx e

+++的不动点，所以

3

902=

e

d c b a +++3

2

3

3902)902(902，解得a=e,d=0,b=216c

取m=216；根据上述提示，取：a=e=1,b=216,c=1,d=0,程序如下：

f=inline('(x+216)/(x^2+1)');

x0=1;

for i=1:100; x0=f(x0);

fprintf('%g %g\n',i,x0); end

结果如下

1 451.5

2 0.00663958

3 901.967

4 0.002217413

5 901.998

6 0.0022173

7 901.99

8 8 0.0022173

9 901.998 10 0.0022173

12 0.0022173

13 901.998

14 0.0022173

15 901.998

16 0.0022173

17 901.998

……

95 901.998

96 0.0022173

97 901.998

98 0.0022173

99 901.998

100 0.002217

初值为-1，结果为

1 450.5

2 0.00666416

3 901.967

4 0.00221742

5 901.998

6 0.0022173

7 901.998

8 0.0022173

9 901.998

10 0.0022173

11 901.998

12 0.0022173

13 901.998

14 0.0022173

15 901.998

16 0.0022173

17 901.998

18 0.0022173

19 901.998

……

95 901.998

96 0.0022173

97 901.998

98 0.0022173

99 901.998

100 0.0022173

初值为1000，结果为

1 0.001216

2 901.999

4 901.998

5 0.0022173

6 901.998

7 0.0022173

8 901.998

9 0.0022173 10 901.998 11 0.0022173 12 901.998 13 0.0022173 14 901.998 15 0.0022173 16 901.998 17 0.0022173 18 901.998 19 0.0022173 20 901.998 ……

93 0.0022173 94 901.998 95 0.0022173 96 901.998 97 0.0022173 98 901.998 99 0.0022173 100 901.998

第三次练习

教学要求：理解线性映射的思想，会用线性映射和特征值的思想方法解决诸如天气等实际问题。 3.1 对???

?

??=3124A ，T

T x x )2,1(),()0(2)0(1=，求出}{n x 的通项. >> syms n

>> A=sym('[4,2;1,3]');x=[1;2];[P,D]=eig(A) %没有sym 下面的矩阵就会显示为小数

P = [ -1, 2] [ 1, 1] D = [ 2, 0] [ 0, 5]

>> An=P*D^n*inv(P) An =

[ 2^n/3 + (2*5^n)/3, (2*5^n)/3 - (2*2^n)/3]

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实验练习题(计算机)-南邮-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实验练习题(计算机) 南邮 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

高等数学实验1 函数与极限 - 参考答案

高等数学实验1 函数与极限 参考答案 一．用MA TLAB 计算： 1． 433sin log 210.235 π +- sin(3*pi/5)+log(21)/log(3)-0.23^4+452^(1/3)-sqrt(43) ans = 4.8365 2．2ln 645 1.2374cos 48 -?+π 4*cos(4*pi/7)+3*2.1^8/(sqrt (645))-log(2) ans = 43.0950 二. 用MA TLAB 计算： 设向量(1,2,3,4,5)x =，求 1．sin 2y x x =+ clear >> x=[1,2,3,4,5]; >> y=sin(x)+2*x y = 2.8415 4.9093 6.1411 7.2432 9.0411 2．2 3sin z x x x =- z=3*x.*sin(x)-x.^2 z =1.5244 1.4558 -7.7299 -25.0816 -39.3839 3． ()2 cos 2ln(21) x x u e x = -+ u=((cos(2.*x)).^2+(sin(x)+1).^(1/2))./(exp(x)-log(2.*x+1)) u = 0.9448 0.3130 0.1097 0.0098 0.0062

三．用MA TLAB 绘图： x 1 2.1 3 3.9 5.3 6.1 6.9 8 9.1 y 1.01 3.98 8.99 16.01 25.41 37.01 48.89 63.89 81.21 clear >> x=[1,2.1,3,3.9,5.3,6.1,6.9,8,9.1]; >> y=[1.01,3.98,8.99,16.01,25.41,37.01,48.89,63.89,81.21]; >> plot(x,y) >> hold on >> plot(x,y,'s') 2.作出函数2y x =与3 y x = x ∈[-3，3]的图象； clear >> hold off >> fplot('x^2',[-3,3]) >> fplot('x^2',[-3,3],'r') hold on >> fplot('x^3',[-3,3],'g') 3．在同一坐标系作出下列函数的图形，并用不同颜色表示。 （1）sin y x x =+ （2）cos y x x =+ clear >> hold off >>fplot('x+sin(x)',[-5,5],'s') >> hold on >> fplot('x+cos(x)',[-5,5],'r') >> hold off >>ezplot('(y-(x+sin(x)))*x*y',[-5,5]) >> hold on >> fplot('x+cos(x)',[-5,5],'r') 4．作下列函数图形：

河南城建学院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数学实验答案(全)

第一次练习 教学要求：熟练掌握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)

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)

数学实验1-3章习题答案

2 一元微积分实验 2.1 曲线绘图 练习题2.1 会出下列常见曲线的图形（其中a=1,b=2,c=3）. 1. 立方抛物线3x y = syms x y; >> ezplot('y=x^(1/3)',[-5,5]) >> title('y=x^(1/3)') -5 -4 -3 -2 -1 01 2 3 4 5 x y y=x (1/3) 2. 高斯曲线2 x e y -= syms x y; >> ezplot('y=exp(-x^2)',[-5,5]) >> title('y=exp(-x^2)')

-5 -4 -3 -2 -1 01 2 3 4 5 x y y=exp(-x 2) 3. 笛卡尔曲线 2 2 213,13t at y t at x +=+= )3(33axy y x =+ syms x y; >> ezplot('x^3+y^3=3*x*y',[-2,2]) >> title('x^3+y^3-3*x*y=0') -2 -1.5-1-0.5 00.51 1.5 2 x y x 3+y 3-3*x*y=0 4. 蔓叶线 ).(1,132 2 322x a x y t at y t at x -=+=+= syms x y; >> ezplot('y^2*(1-x)=x^3',[-10,10]) >> title('y^2=x^3/(1-x)')

-10 -8 -6 -4 -2 02 4 6 8 10 x y y 2=x 3/(1-x) 5. 摆线 ).cos 1(),sin (t b y t t a x -=-= syms t; >> x=t-sin(t); >> y=2-2*cos(t); >> ezplot(x,y) 1 2 345 6 00.511.522.533.54x y x = t-sin(t), y = 2-2 cos(t) 6. 星形线 )(sin ,cos 3 23 23 23 3 a y x t a y t a x =+== syms t; >> x=cos(t)^3; >> y=sin(t)^3; >> ezplot(x,y)

南邮数学实验问题详解

第一次练习题 1、求032=-x e x 的所有根。 >>x=-5:0.01:5;y=exp(x)-3*x.^2;plot(x,y);grid on >> fsolve('exp(x)-3*x.^2',-1) Equation solved. fsolve completed because the vector of function values is near zero as measured by the default value of the function tolerance, and the problem appears regular as measured by the gradient. ans = -0.4590

>> fsolve('exp(x)-3*x.^2',1) Equation solved. fsolve completed because the vector of function values is near zero as measured by the default value of the function tolerance, and the problem appears regular as measured by the gradient. ans = 0.9100 >> fsolve('exp(x)-3*x.^2',4) Equation solved. fsolve completed because the vector of function values is near zero as measured by the default value of the function tolerance, and the problem appears regular as measured by the gradient.

MAAB数学实验第二版答案胡良剑

数学实验答案 Chapter 1 Page20,ex1 (5) 等于[exp(1),exp(2);exp(3),exp(4)] (7) 3=1*3, 8=2*4 (8) a为各列最小值，b为最小值所在的行号 (10) 1>=4,false, 2>=3,false, 3>=2, ture, 4>=1,ture (11) 答案表明：编址第2元素满足不等式(30>=20)和编址第4元素满足不等式(40>=10) (12) 答案表明：编址第2行第1列元素满足不等式(30>=20)和编址第2行第2列元素满足不等式(40>=10) Page20, ex2 (1)a, b, c的值尽管都是1，但数据类型分别为数值，字符，逻辑，注意a与c相等，但他们不等于b (2)double(fun)输出的分别是字符a,b,s,(,x,)的ASCII码 Page20,ex3 >> r=2;p=0.5;n=12;

>> T=log(r)/n/log(1+0.01*p) Page20,ex4 >> x=-2:0.05:2;f=x.^4-2.^x; >> [fmin,min_index]=min(f) 最小值最小值点编址 >> x(min_index) ans = 0.6500 最小值点 >> [f1,x1_index]=min(abs(f)) 求近似根--绝对值最小的点f1 = 0.0328 x1_index = 24 >> x(x1_index) ans = -0.8500

>> x(x1_index)=[];f=x.^4-2.^x; 删去绝对值最小的点以求函数绝对值次小的点>> [f2,x2_index]=min(abs(f)) 求另一近似根--函数绝对值次小的点 f2 = 0.0630 x2_index = 65 >> x(x2_index) ans = 1.2500 Page20,ex5 >> z=magic(10) z = 92 99 1 8 15 67 74 51 58 40 98 80 7 14 16 73 55 57 64 41 4 81 88 20 22 54 56 63 70 47 85 87 19 21 3 60 62 69 71 28

离散数学实验指导书及其答案

实验一命题逻辑公式化简 【实验目的】加深对五个基本联结词（否定、合取、析取、条件、双条件）的理解、掌握利用基本等价公式化简公式的方法。 【实验内容】用化简命题逻辑公式的方法设计一个表决开关电路。 实验用例：用化简命题逻辑公式的方法设计一个5人表决开关电路，要求3人以上（含3人）同意则表决通过（表决开关亮）。 【实验原理和方法】 （1）写出5人表决开关电路真值表，从真值表得出5人表决开关电路的主合取公式（或主析取公式），将公式化简成尽可能含五个基本联结词最少的等价公式。 （2）上面公式中的每一个联结词是一个开关元件，将它们定义成C语言中的函数。 （3）输入5人表决值（0或1），调用上面定义的函数，将5人表决开关电路真值表的等价公式写成一个函数表达式。 （4）输出函数表达式的结果，如果是1，则表明表决通过，否则表决不通过。 参考代码： #include int vote(int a,int b,int c,int d,int e) { //五人中任取三人的不同的取法有10种。 i f( a&&b&&c || a&&b&&d || a&&b&&e || a&&c&&d || a&&c&&e || a&&d&&e || b&&c&&d || b&&c&&e || b&&d&&e || c&&d&&e) return 1; e lse return 0; } void main() { i nt a,b,c,d,e; p rintf("请输入第五个人的表决值（0或1，空格分开）："); s canf("%d%d%d%d%d",&a,&b,&c,&d,&e); i f(vote(a,b,c,d,e)) printf("很好，表决通过！\n"); e lse printf("遗憾，表决没有通过！\n"); } //注：联结词不定义成函数，否则太繁