当前位置：搜档网 › 数值计算方法实验报告(西安工程大学附有程序)

数值计算方法实验报告(西安工程大学附有程序)

数值分析

上

机

实

验

报

告

学院：理学院

专业：数学与应用数学

学号：41108040110

姓名：吕炳林

实验报告1

例1

实验程序

>>I0 = exp(-1)*pond(‘x.^0.*exp(x.^2)’,0,1) I0 =

0.5381

>>vpa(I0,10)

实验结果：

ans =

0.5380795164

例1.9

实验程序：

function p=p(n,x)

a(1)=13;

fori=1:n+1

a(i+1)=2*a(i)+3;

end

S=a(n+1);

for j=1:n

S=x*S+a(n+1-j);

end

p=S;

End

在命令窗口进行计算：

实验结果

>>p(100,0.5)

ans =

600

>>p(150,13)

ans =

1.0995e+213

习题3

实验程序

Function y=y(n)

Y=28

fori=1:n

y=y-1/100*sqrt(783);

End

在命令窗口进行计算

实验结果

y(100)

ans=

0.0179

y(500)

ans=

-111.9107

实验报告2

习题3

实验原理：一次插值：2111

212201

()x x x x l x y y x x x x --=+-- 二次插值：

111121*********()()()()()()()()()()()()()k k k k k k k k k k k k k k k k k k k k k x x x x x x x x x x x x l x y y y x x x x x x x x x x x x +-+--+--+-++-+------=++------

实验运行环境：

本实验采用Matlab 编写。

实验程序：

线性插值：

x=0.4:0.1:0.8;

f=[-0.916291,-0.693147,-0.510826,-0.357765,-0.223144];

format long

interp1(x,f,0.54)

实验结果：

ans =

-0.620218600000000

二次插值：

x=0.54;

a=[0.4,0.5,0.6];

b=[-0.916291,-0.693147,-0.510826];

l=b(1)*(x-a(2))*(x-a(3))/((a(1)-a(2))*(a(1)-a(3)));

m=b(2)*(x-a(1))*(x-a(3))/((a(2)-a(1))*(a(2)-a(3)));

n=b(3)*(x-a(1))*(x-a(2))/((a(3)-a(1))*(a(3)-a(2)));

y=l+m+n

实验结果：

y =

-0.61531984000000

习题21

实验运行环境：

本实验采用Matlab 编写

实验程序1：（输入函数y ）

for x=-4.5:4.5

y=1/(x^2+1)

end

显示结果

y=0.04705882352941；

y=0.07547169811321；

y=0.13793103448276；

y=0.30769230769231；

y=0.80000000000000；

实验程序2（求()h I x 的值）

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

a=[x-0.5,x+0.5];

y=[1/(1+(x-0.5)^2),1/(1+(x+0.5)^2)];

I=y(1)*(x-a(2))/(a(1)-a(2))+y(2)*(x-a(1))/(a(2)-a(1))

显示结果

当分别输入 4.5, 3.5, 2.5, 1.5,0.5±±±±±时

I=0.0486 I=0.0794 I=0.1500I=0.3500 I =0.7500

习题24

[]

[]3333336.7593(0.30) 4.8810(0.25)10.0169(0.30)10.9662(0.25)0.25,0.302.7117(0.39) 1.9098(0.30) 6.1075(0.39) 6.9544(0.30)0.30,0.39() 2.8647(0.45) 2.2422(0.39)10.4186(0.45x x x x x x x x x x S x x x x ----+-+-∈----+-+-∈=----+-[][]

33)10.9662(0.39)0.39,0.451.6817(0.53) 1.3623(0.45)8.3958(0.53)9.1087(0.45)0.45,0.53x x x x x x x ????????+-??∈??----+-+-?∈??

(2)

[]

[]333336.2697(0.25)10(0.3)10.9697(0.25)0.25,0.303.4831(0.39) 1.5956(0.3) 6.1138(0.39) 6.9518(0.30)0.30,0.39() 2.3933(0.45) 2.8622(0.39)10.4186(0.45)11.1903(0.39)0.3x x x x x x x x x S x x x x x x --+-+-∈----+-+-∈∴=----+-+-∈[][]

39,0.452.1467(0.53)8.3987(0.53)9.1(0.45)0.45,0.53x x x x ????????????--+-+-?∈??实验报告3

习题22

实验程序

x=[19,25,31,38,44]

y=[19.0,32.3,49.0,73.3 97.8]

[p,s]=polyfit(x,y,2)

实验结果 p =

0.0497 0.0193 0.6882

s =

R: [3x3 double]

df: 2

normr: 0.1145

习题23。

实验程序

t=[0,0.9,1.9,3.0,3.9,5.0];

s=[0,10,30,50,80,110];

[P,S]=polyfit(t,s,5)

显示结果

P =

-0.5432 6.4647 -26.5609 46.1436 -13.2601 -0.0000

S =

R: [6x6 double]

df: 0

normr: 1.2579e-012

实验结果：

习题 24.

实验程序：>> x=0:5:55;

-13.2601x 46.1436x -26.5609x 6.4647x -0.5432x 2

345++

y=[0,1.27,2.16,2.86,3.44,3.87,4.15,4.37,4.51,4.58,4.62,4.64]; >> [P,S]=polyfit(x,y,5)

显示结果

P =

0.0000 -0.0000 0.0002 -0.0084 0.2851 0.0082 S =

R: [6x6 double]

df: 6

normr: 0.0487

实验结果：实验报告4

习题4

实验程序

先求精确值：

x=sym('x');

f=exp(-x);

I=int(f,x,0,1)

I =

1-exp(-1)

利用simpson 公式得；

a=0;b=1;

S=(b-a)/6*[exp(-a)+4*exp(-(a+b)/2)+exp(-b)]

显示结果S = 0.6323

实验结果：其误差为|S-I|=|0.6323+exp(-1)-1|≈1.7944e -004

习题8

实验程序；

ifnargin<5

Eps=1E-6;

end

m=1;

h=(b-a);

err=1;

j=0;

T=zeros(4,4);

T(1,1)=h*(limit(f,a)+limit(f,b))/2;

while ((err>Eps) & (j

j=j+1;

h=h/2;

s=0;

for p=1:m

x0=a+h*(2*p-1);

s=s+limit(f,x0);

end

0082.02851

.00084.00002.023++-x x x

T(j+1,1)=T(j,1)/2+h*s;

m=2*m;

for k=1:j

T(j+1,k+1)=T(j+1,k)+(T(j+1,k)-T(j,k))/(4^k-1); end

err=abs(T(j,j)-T(j+1,k+1));

end

I=T(j+1,j+1);

ifnargout==1

T=[];

end

然后在输入框输入：

>>syms x;

>> f=sym('exp(-x)*2/sqrt(3.1415926)')

f =

exp(-x)*2/sqrt(3.1415926)

>> [I,T]=romberg(f,0,1,4,1E-8)

I =

0.7133

T =

0.7717 0 0 0 0

0.7281 0.7135 0 0 0

0.7170 0.7133 0.7133 0 0

0.7142 0.7133 0.7133 0.7133 0

0.7135 0.7133 0.7133 0.7133 0.7133

实验结果I=0.7133

习题11

Romberg方法

利用先前的Romberg.m文件

实验程序:

>>syms y;

f=sym('1/y');

[I,T]=romberg(f,1,3,3,1E-5)

I =

1.0986

T =

1.3333 0 0 0 0

1.1667 1.1111 0 0 0

1.1167 1.1000 1.0993 0 0

1.1032 1.0987 1.0986 1.0986 0

1.0998 1.0986 1.0986 1.0986 1.0986

实验结果积分值为1.0986

(2)三点及五点Gauss公式

三点公式:

查表有 , :

先作变换,将积分区间变换到上,令y=t+2,则有:

于是就可以利用求积公式来解题.

%三点公式

f=inline('1/(t+2)')

f =

Inline function:

f(t) = 1/(t+2)

>> x=[0.7745966692,- 0.7745966692,0]

x =

0.7746 -0.7746 0

>> A=[0.555555556, 0.555555556,0.888888889]

A =

0.5556 0.5556 0.8889

>> I=0;

>>fori=1:length(x)

I=I+feval(f,x(i))*A(i);%三点Simpson公式

end

>> I

I =

1.0980

%五点公式

>> f=inline('1/(t+2)')

f =

Inline function:

f(t) = 1/(t+2)

>> x=[0.9061793,-0.9061793,0.5384693,- 0.5384693,0]

x =

0.9062 -0.9062 0.5385 -0.5385 0

>> A=[0.2369263,0.2369263, 0.4786287,0.4786287,0.5688889]

A =

0.2369 0.2369 0.4786 0.4786 0.5689

>> I=0;fori=1:length(x)

I=I+feval(f,x(i))*A(i);%三点Simpson公式

end

>> I

I =

1.0986

a=1;

b=3;

f=inline('1/y');

%m=2,n=1Gauss-Legendre复化求积公式

n=1;

m=4;

h=(b-a)/(m);

x=[0.5773502692 -0.5773502692];

A=[1 1];

k=1;

I(k)=0;

for j=1:n+1

temp=0;

fori=1:m

tx=a+(2*(i-1)+1)*h/2+h*x(j)/2;

temp=temp+feval(f,tx);

end

I(k)=I(k)+A(j)*temp;

end

I(k)=I(k)*h/2;

>> I

I =

1.0985

实验报告5 习题2实验程序：

(1)先用Euler方法：

>> F='x+y';

a=0;

b=1;

h=0.1;

n=(b-a)/h;

X=a:h:b;

Y=zeros(1,n+1);

Y(1)=0;

fori=2:n+1

x=X(i-1);

y=Y(i-1);

Y(i)=Y(i-1)+eval(F)*h; end

>>disp([X',Y']);

实验结果

ans =

0 1.0000

0.1000 1.1103

0.2000 1.2428

0.3000 1.3997

0.4000 1.5836

0.5000 1.7974

0.6000 2.0442

0.7000 2.3275

0.8000 2.6511

0.9000 3.0192

1.0000 3.4366

(2)用改进的Euler方法：实验程序

>> F='x+y';

>> Y1=zeros(1,n+1);

Y1(1)=0;

fori=2:n+1

x=X(i-1);

y=Y1(i-1);

ty=Y1(i-1)+eval(F)*h;

Y1(i)=Y1(i-1)+h/2*eval(F); x=X(i);

y=ty;

Y1(i)=Y1(i)+h/2*eval(F); end

>>disp([X',Y1']);

实验结果ans =

0 1.0000

0.1000 1.1103

0.2000 1.2428

0.3000 1.3997

0.4000 1.5836

0.5000 1.7974

0.6000 2.0442

0.7000 2.3275

0.8000 2.6511

0.9000 3.0192

1.0000 3.4366 习题6

（1）实验程序

>> F='x+y';

a=0;

b=1;

h=0.2;

n=(b-a)/h;

X=a:h:b;

Y=zeros(1,n+1);

Y(1)=1;

fori=1:n

x=X(i);

y=Y(i);

K1=h*eval(F);

x=x+h/2;

y=y+K1/2;

K2=h*eval(F);

x=x;

y=Y(i)+K2/2;

K3=h*eval(F);

x=X(i)+h;

y=Y(i)+K3;

K4=h*eval(F);

Y(i+1)=Y(i)+(K1+2*K2+2*K3+K4)/6; end

实验结果；

>>disp([X',Y']);

0 1.0000

0.2000 1.2428

0.4000 1.5836

0.6000 2.0442

0.8000 2.6510

1.0000 3.4365

(2)

实验程序

>> F='3*y/(1+x)';

a=0;

b=1;

h=0.2;

n=(b-a)/h;

X=a:h:b;

Y=zeros(1,n+1);

Y(1)=1;

fori=1:n

x=X(i);

y=Y(i);

K1=h*eval(F);

x=x+h/2;

y=y+K1/2;

K2=h*eval(F);

x=x;

y=Y(i)+K2/2;

K3=h*eval(F);

x=X(i)+h;

y=Y(i)+K3;

K4=h*eval(F);

Y(i+1)=Y(i)+(K1+2*K2+2*K3+K4)/6; end

实验结果

>>disp([X',Y']);

0 1.0000

0.2000 1.7275

0.4000 2.7430

0.6000 4.0942

0.8000 5.8292

1.0000 7.9960

实验六习题2 （1）实验程序

fa=’1-a*sin(a)’;

fc=’1-c*sin(c)’;

b=pi/2;

a=0;

R=1;

k=0；

while (R>5e-6) ;

c=(a+b)/2;

ifeval(fa)*eval(fc)>0;

a=c;

else

b=c;

end

R=b-a;

k=k+1;

end

x=c

实验结果

ans =1.41021

习题4实验程序

fa=’exp(a)+10*a-2’;

fc=’exp(c)+10*c-2’;

b=1;

a=0;

R=1;

k=0；

while (R>5e-6) ;

c=(a+b)/2;

ifeval(fa)*eval(fc)>0;

a=c;

else

b=c;

end

R=b-a;

k=k+1;

end

x=c

实验结果

x =0.090515136718750

（2）实验程序

f=input('请输入需要求解函数>>','s')

df=diff(f);

miu=2;

x0=input('input initial value x0>>');

k=0; max=100; R=eval(subs(f,'x0','x'))

while (abs(R)>1e-8)

x1=x0-miu*eval(subs(f,'x0','x'))/eval(subs(df,'x0','x'));

R=x1-x0;

x0=x1;

k=k+1;

if (eval(subs(f,'x0','x'))<1e-10);

break

end

if k>max;

ss=input('maybe result is error,choose a new x0,y/n?>>','s');

ifstrcmp(ss,'y')

x0=input('input initial value x0>>');

k=0;

else

break

end

x=x0

实验结果x = 0.090526468052644

习题7实验程序

x=sym('x');

f=sym('x^3-3*x-1');

df=diff(f,x);

FX=x-f/df

FX =x + (3*x - x^3 + 1)/(3*x^2 - 3)

Fx=inline(FX);

x0=2;

fori=1:10

disp(x0);

x0=feval(Fx,x0);

end

format short;

显示结果

1.888888888888889

1.879451566951567

1.879385244836671

1.879385241571817

实验结果

x0 =

1.879385241571817

习题13 实验程序

f=input('请输入需要求解函数>>','s')

df=diff(f);

miu=2;

x0=input('input initial value x0>>');

k=0; max=100; R=eval(subs(f,'x0','x'))

while (abs(R)>1e-8)

x1=x0-miu*eval(subs(f,'x0','x'))/eval(subs(df,'x0','x'));

R=x1-x0;

x0=x1;

k=k+1;

if (eval(subs(f,'x0','x'))<1e-10);

break

end

if k>max;

ss=input('maybe result is error,choose a new x0,y/n?>>','s');

ifstrcmp(ss,'y')

x0=input('input initial value x0>>');

k=0;

else

break

end

x=x0

实验结果

1072380529

实验七课题利用lu 分解求出方程组CijXi=bi(其中i=j=100,C 是一百阶方阵，b 是从1到100的列向量)解线形方程组，根据题目要求C 是一百阶方阵且满足下面约束条件

A=ones(100);

C=diag([diag(A,1)],-1)+diag([diag(A,1)],1)+diag([diag(A,2)],2)+diag([diag(A,2)],-2)+4*eye(100)

b=(1:100)’;

1234510.33043478,10.70242553,10.7237414,

10.72380529,10.72380529

x x x x x =====

实验分析利用matlab中的lu函数对方阵C进行L,U分解

[L,U]=lu(C)

根据Ly=b可以求出y=inv(L)*b

实验结果

y=(1.0000, 1.7500 , 2.4000 , 3.0000 , 3.6771 , 4.3475 ,

5.0087 ,5.6727 ,

6.3401 ,

7.0060 , 7.6713

8.3372,

9.0031 ,

9.6689 ,10.3347, 11.0005 , 11.6663 ,12.3321 , 12.9980 13.6638 , 14.3296 , 14.9954 , 15.6612 ,16.3270, 16.9928 ,

17.6586 , 18.3245 18.9903 , 19.6561 , 20.3219 , 20.9877 , 21.6535 , 22.3193 , 22.9851 , 23.6510 24.3168,

24.9826 , 25.6484 , 26.3142, 26.9800 , 27.6458 , 28.3117 ,

28.9775 29.6433 , 30.3091 , 30.9749 , 31.6407 , 32.3065,

32.9723 ,

33.6382 , 34.3040 34.9698 , 35.6356 36.3014 36.9672 37.6330 3 8.2988 38.9647 , 39.6305 40.2963,

40.9621 , 41.6279 , 42.2937 , 42.9595 ,

43.6254 , 44.2912 , 44.9570 45.6228 , 46.2886 ,46.9544 , 47.6202 , 48.2860 , 48.9519 , 49.6177 , 50.2835 50.9493 , 51.6151 ,

52.2809 , 52.9467 ,

53.6125 , 54.2784, 54.9442 , 55.6100 56.2758, 56.9416 , 57.6074 , 58.2732 , 58.9390 , 59.6049 , 60.2707 ,

60.9365 61.6023 , 62.2681 , 62.9339 , 63.5997 ,

64.2656, 64.9314 , 65.5972, 66.2630 66.9288)’

在根据Ｕx=y可以求出x=inv(U)*y则方程的解为

实验结果

(0.0898 , 0.2578 , 0.3832 , 0.4960 , 0.6236 , 0.7514 , 0.8751 ,

0.9996 ， 1.1251 , 1.2501 , 1.3750 , 1.5000 ,

1.6250 , 1.7500 , 1.8750 ,

2.0000 2.1250 , 2.2500 ,

2.3750 , 2.5000 , 2.6250 , 2.7500 , 2.8750 ,

3.0000 3.1250 ,

3.2500 , 3.3750 , 3.5000 , 3.6250 , 3.7500 , 3.8750 ,

4.0000

4.1250 , 4.2500 , 4.3750 , 4.5000 ,

4.6250 , 4.7500 , 4.8750 ,

5.0000 5.1250 , 5.2500 , 5.3750 ,

5.5000 , 5.6250 , 5.7500 ,

5.8750 ,

6.0000 6.1250 , 6.2500 , 6.3750 ,

6.5000 , 6.6250 , 6.7500 , 6.8750 ,

7.0000 7.1250 ,

7.2500 , 7.3750 , 7.5000 ,7.6250 , 7.7500 , 7.8750 , 8.0000 8.1 250 ,.2500 , 8.3750 , 8.5000 , 8.6250 , 8.7500 , 8.8750 , 9.000 0 , 9.1250 9.2500 , 9.3750 , 9.5000 , 9.6250 ,

9.7500 , 9.8750 , 10.0000 , 10.1250 10.2501 , 10.3749 ,

10.4997 ,10.6254 , 10.7507 , 10.8730 , 10.9991 , 11.1330 11.2471 ,

11.3490 , 11.5307 , 11.6932 , 11.5880 , 11.7596, 12.6645 ,

12.0160 9.9785 , 14.5334 , 18.8720)'

实验八

用Cholesky分解（题目和上实验的相同）

C=diag([diag(A,1)],-1)+diag([diag(A,1)],1)+diag([diag(A,2)],2)+diag([diag( A,2)],-2)+4*eye(100)

b=(1:100)’

西安交通大学接口技术实验报告

西安交通大学微型计算机接口技术实验报告班级:物联网姓名: 学号：

实验一基本I/O扩展实验一、实验目的 1、了解 TTL 芯片扩展简单 I/O 口的方法，掌握数据输入输出程序编制的方法； 2、对利用单片机进行 I/O 操作有一个初步体会。二、实验内容 74LS244 是一种三态输出的8 总线缓冲驱动器，无锁存功能，当G 为低电平时，Ai 信号传送到Yi，当为高电平时，Yi 处于禁止高阻状态。 74LS273 是一种8D 触发器，当CLR 为高电平且CLK 端电平正跳变时，D0——D7 端数据被锁存到8D 触发器中。实验原理图：三、实验说明利用74LS244 作为输入口，读取开关状态，并将此状态通过74LS273 再驱动发光二极管显示出来，连续运行程序，发光二极管显示开关状态。

四、实验流程图五、实验连线 1、244的cs连接到CPU地址A15，Y7—Y0连接开关K1-K8; 2、273的CS连接到CPU地址A14，Q7-Q0连接到发光二极管L1-L8; 3、该模块的WR,RD连接CPU的WR,RD,数据线AD7-AD0，地址线A7-A0分别与CPU的数据线AD7-AD0，地址线A7-A0相连接。

六、程序源代码（略）七、实验结果通过开关K01 到K08 可以对应依次控制LED 灯的L1 到L8 ，即当将开关Ki 上拨时，对应的Li 被点亮，Ki 下拨时，对应的Li熄灭。此外，如果将开关拨到AAH 时，将会产生LED 灯左移花样显示；如果开关拨到55H 时，将会产生LED 灯右移花样显示。七、实验心得通过本次实验，我了解了TTL 芯片扩展简单I/O 口的方法，同时也对数据输入输出程序编制的方法有一定的了解与掌握，对利用单片机进行I/O 操作有一个初步体会，实验使我对自己在课堂上学的理论知识更加理解，同时也锻炼了我的动手操作能力。

西安邮电大学课程设计报告

一. 实验目的及实验环境实验目的：目前，对于中小超市来说，需要采取一套行之有效的超市管理系统来管理超市的日常进货、销售、员工管理等信息。本次设计的超市管理系统就是想让用户以最少的手工业务处理，来完成超市中的日常工作，增加超市业务处理的效率。一个典型的超市管理系统，除了应该完成超市中必要的销售工作，如收银台，还应该具有后台管理的功能。实验环境：系统的体系结构：C/S结构前台：Java (Eclipse 3.5) 后台数据库：SQL Server 2008 二. 实验内容超市管理系统该系统应具备前台收银操作，还有后台对各种数据的管理。该系统应为多用户使用的，应以用户的权限来进入系统管理模块。友好的界面及方便简单的操作。数据存放在数据库中，系统对数据库可进行增、删、查、改的操作。三．方案设计系统开发的总体任务是实现Java与SQL Server 2008的组合应用，实现用JDBC 方式调用SQL数据库，且有一定的透明性。超市管理系统需要完成的功能主要有： 1、用户分类登录系统及权限管理。 2、超市收银台的销售管理方式。 3、销售商品信息的录入、查询、修改、删除。 4、供应商信息的录入，查询，修改、删除。 5、超市员工的职位基本信息的设定。 6、员工的基本信息的录入，查询，修改和删除。 7、仓库信息的录入与查看。 8、数据间的相互联系：销售商品的数量和库存量的基本联系；商品信息和供应商的相关数据联系；员工与职位间的联系等。 9、数据存放在数据库中，任何数据操作都是对数据库的数据的操作。

系统结构图（功能模块图）各个模块的功能描述： (1)、登录界面：本系统简单的设置有三种职位，用户使用密码登录，成功后给其相应的操作界面。 (2)、仓库管理：是仓库管理员和老板可使用的操作，主要是查看仓库的入库信息和添加入库信息。 (3)、管理商品信息：只有老板可用，主要是查看现有商品的信息、修改、删除及添加新的商品信息。 (4)、管理供应商信息：只有老板可用，主要是查看现有的供应商的信息、修改、删除及添加新供应商的信息。 (5)、管理员工信息：只有老板可用，主要是查看现有的员工信息、修改、删除及增加新员工。 (6)、查看销售纪录。四．测试数据及运行结果 1、Java与数据库SQL 2008的链接： String driverName="com.microsoft.sqlserver.jdbc.SQLServerDriver"; String dbURL="jdbc:sqlserver://localhost:1433;DatabaseName=SuperMarket"; 2、用户登陆(登陆界面，老板界面，仓库管理员界面，销售员界面)：

西安交通大学19年5月补考软件工程概论作业考核试题含主观题答案

西安交通大学19年5月补考《软件工程概论》作业考核试题（含主观题）一、单选题(共10 道试题,共30 分) 1.可行性分析中,系统流程图用于描述()。 A.当前运行系统 B.当前逻辑模型 C.目标系统 D.新系统正确答案:A 2.软件生存周期包括可行性分析和项目开发计划、需求分析、概要设计、详细设计、编码、()、维护等活动。 A.应用 B.测试 C.检测 D.以上答案都不正确正确答案:B 3.在软件开发和维护中所产生的一系列严重的问题通常称为软件危机,这些问题中相对次要的因素是()。 A.文档质量 B.开发效率 C.软件功能 D.软件性能正确答案:D 4.在软件的开发与维护过程中,用来存储.更新.恢复和管理一个软件的多版本,它是()工具。 A.文档分析 B.项目管理 C.成本估算 D.版本控制正确答案:D 5.()是计算机程序及其说明程序的各种文档。 A.软件 B.文档 C.数据 D.程序正确答案:A 6.软件生存周期中时间最长的阶段是()。 A.需求分析阶段 B.概要设计阶段 C.测试阶段

D.维护阶段正确答案:D 7.软件开发方法中的原型化方法是一种()型的设计过程。 A.自外向内 B.自顶向下 C.自内向外 D.自底向上正确答案:A 8.软件开发环境支持一定的()或按照一定的软件开发模型组织而成。 A.软件生存周期 B.软件过程 C.软件开发方法 D.软件开发模型正确答案:C 9.加工是对数据流图中不能再分解的基本加工的精确说明,下述()是加工的最核心。 A.加工顺序 B.加工逻辑 C.执行频率 D.激发条件正确答案:B 10.软件生存周期模型有多种,下列选项中,()不是软件生存周期模型。 A.螺旋模型 B.增量模型 C.功能模型 D.瀑布模型正确答案:C 二、判断题(共5 道试题,共10 分) 1.瀑布模型本质上是一种线性顺序模型。 A.对 B.错正确答案:A 2.需求分析要回答“软件必须做什么?” A.对 B.错正确答案:A 3.原型化方法适用于需求不确定性高的系统。 A.对

某厂排气管车间生产计划的优化分析实验报告

————————————————————————————————作者: ————————————————————————————————日期:

西安理工大学实验报告成绩第１页(共页) 课程: 实验日期: 年月日专业班号：组别：交报告日期：年月日姓名：学号:报告退发: （订正、重做）同组者: 教师审批签字: 实验名称：某厂排气管车间生产计划的优化分析实验内容: 1、制定利润最大的生产计划（1）建立线性规划数学模型；车间基本生产工人工时计算如下: S＝27×(22×8 – 2 – 4 ) =4590(小时/月）按产量1%的工废实际可用时间:4590×（1– 1% )=４544（小时/月) 每种设备月实际可利用工时(减去工废后的可用时间）如下表所示: 设备名称月可利用工时(小时／月）设备名称月可利用工时（小时/月) 1.平面铣床139２×( 1– 1% )＝1378 2.卧铣床1４00×（1 – 1% )=１386 2０82×( 1 – 1％）=２061 4．单面铣床７0０×( １–１% )=693 3．组合钻床 5.攻丝床2088×（1– 1% )＝２067 ６．精铣床1396×( 1 – 1%)=1３8２ 1376×（1– 1%）=１３62 8．摇臂钻床1384×( 1 –１% ）=137０7．扩孔钻床９.去毛刺 700×(１– 1% )=693 10.清洗机700×( 1 – 1%)=693 机每种排气管利润（减去料废后的利润）如下表所示：

品项目 1 2 3 4 5 6 7 8 毛坯价格 ① 98１7１０4１02 料废成本(① ×2％) １．96 2.0８ 1.88 2.２4 2.1２ 1.94２.08 2.０4辅料消耗 2 2 2 2 2 2 2 2 动能消耗１10 10 工具等消耗 1 １管理费用１.４５５１.099 1.21 1.44 1.188 1．2２65 1.3０8 1.56 税收15 16 14.8 １7 16.5 14.5 15.6 15．５售价150 １６ 0.1 1４9 １72 166 145．6 15７.8 1５５.8 利润(元/根）11．5８５11.921 13．１１1 ３.32 1 ３.192 １0．９ 335 13.812 １1.7 根据给定的条件及要求,设x1，x2,…ｘ8分别为第一种、第二种……第八种排气管的月产量，月产量单位为1000根/月。将目标函数取为每月的总产品利润,得出关于月计划生产的线性规划模型为: ｍaxＺ= 11５8５ｘ1＋1192１ｘ2+１3110 x3+1３32０x４＋13192 x5+10９33.5x ６+13812 x7＋１１7００x8 ４1.9ｘ1+４3.3ｘ2+４1.７x3+4８．1ｘ4+46．5ｘ５＋40.9 x6+47．8 x7+45.７x8 ≤ ４544 4.0 x1＋4．5 x2+４．8 x3+５.8 x４+５.2ｘ5+４.０ｘ６+４.６ｘ7+５.６x8 ≤1378 3．９x1＋4.5 x２+4．３x3+5．0ｘ4+4.9ｘ5+4.4 x６+5.1 x７+４.8 x８≤１38６ 5．9x1＋5．8ｘ２＋5.7 x3+６．3 x4+６.5 x５+6.0 x６＋6.6ｘ７+6．４ｘ8 ≤2０61 3.５ｘ１+3.0ｘ２+3.7ｘ3＋ 4.０ｘ４+３.８x5+3.0 x6+4.1 x7＋3.4ｘ8≤693 ５.８x1+6.2 x2+５．７ｘ3+6.４x４+6.3 x5+6．0ｘ６+6.５x7＋６.２x8 ≤20６7 5．5 x1＋５.7x２＋4.7ｘ3+6．0 x4+5.９x5+5．２ｘ6＋6.2 x7＋５.6ｘ８≤138２ 3.9ｘ1+3.8x2+ 4.０x3+４.1 x４+3.7 ｘ5+3.5 x６+4．1 x7+3.６x8 ≤13６2

西安交大数字图像处理第二次实验报告

数字图像处理第二次作业

摘要本次报告主要记录第二次作业中的各项任务完成情况。本次作业以Matlab 2013为平台，结合matlab函数编程实现对lena.bmp,elain1.bmp图像文件的相关处理：1.分别得到了lena.bmp 512*512图像灰度级逐级递减8-1显示，2.计算得到lena.bmp图像的均值和方差，3.通过近邻、双线性和双三次插值法将lena.bmp zoom到2048*2048，4. 把lena和elain 图像分别进行水平shear（参数可设置为1.5，或者自行选择）和旋转30度，并采用用近邻、双线性和双三次插值法zoom到2048*2048。以上任务完成后均得到了预期的结果。 1.把lena 512*512图像灰度级逐级递减8-1显示（1）实验原理：给定的lena.bmp是一幅8位灰阶的图像，即有256个灰度色。则K位灰阶图像中某像素的灰度值k(x,y)(以阶色为基准)与原图同像素的灰度值v(x,y)(以256阶色为基准)的对应关系为：式中floor函数为向下取整操作。取一确定k值，对原图进行上式运算即得降阶后的k位灰阶图像矩阵。（2）实验方法首先通过imread()函数读入lena.bmp得到图像的灰度矩阵I，上式对I矩阵进行灰度降阶运算，最后利用imshow()函数输出显示图像。对应源程序为img1.m。（3）处理结果 8灰度级

7灰度级 6灰度级 5灰度级

4灰度级 3灰度级 2灰度级

1灰度级（4）结果讨论：由上图可以看出，在灰度级下降到5之前，肉眼几乎感觉不出降阶后图像发生的变化。但从灰度级4开始，肉眼明显能感觉到图像有稍许的不连续，在灰度缓变区常会出现一些几乎看不出来的非常细的山脊状结构。随着灰度阶数的继续下降，图像开始出现大片的伪轮廓，灰度级数越低，越不能将图像的细节刻画出来，最终的极端情况是退化为只有黑白两色的二值化图像。由此可以得出，图像采样的灰度阶数越高，灰度围越大，细节越丰富，肉眼看去更接近实际情况。 2.计算lena图像的均值方差（1）实验原理对分辨率为M*N的灰度图像，其均值和方差分别为：（2）实验方法首先通过imread（）函数读入图像文件到灰度矩阵I中，然后利用 mean2函数和std2函数计算灰度矩阵（即图像）的均值和标准差，再由标准差平方得到方差。对应源程序：img1.m （3）处理结果均值me =99.0512，标准差st =52.8776，方差sf =2.7960e+03。（4）结果分析图像的均值可反应图像整体的明暗程度，而方差可以反应图像整体的对比度情况，方差越大，图像的对比度越大，可以显示的细节就越多。 3.把lena图像用近邻、双线性和双三次插值法zoom到2048*2048；（1）实验原理图像插值就是利用已知邻近像素点的灰度值来产生未知像素点的灰度值，以便由原始图

西南交通大学限修课数学实验题目及答案四

实验课题四曲面图与统计图第一大题：编程作下列曲面绘图：用平面曲线r=2+cos(t)+sin(t)，t∈(0,π)绘制旋转曲面 t=0:0.02*pi:pi; r=2+cos(t)+sin(t); cylinder(r,30) title('旋转曲面'); shading interp 用直角坐标绘制双曲抛物面曲面网线图，z2=xy (-3

axis off 用直角坐标绘制修饰过的光滑曲面曲面：z 4=sin(x )－cos(y ) x 与y 的取值在(-π,π) [x,y]=meshgrid(-pi:0.02*pi:pi); z4=sin(x)-cos(y); surf(x,y,z4); title('picture 4'); shading interp axis off 用连续函数绘图方法绘制曲面)2 s in (6522x y x z ++=,x ∈[-2pi,2pi], y ∈[-2pi,2pi],并作图形修饰。 ezsurf(@(x,y)(x^2+y^2+6*sin(2*x)),[-2*pi 2*pi -2*pi 2*pi]) title('picture 5'); shading interp axis off 第二大题：按要求作下列问题的统计图： x21是1—10的10维自然数构成的向量，y21是随机产生的10维整数向量，画出条形图。(提示bar(x,y)) x21=1:10; y21=randn(10,1); bar(x21,y21) 随机生成50维向量y22，画出分5组的数据直方图。（提示hist(y,n)）