河南城建学院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的数值显示格式

思考题：

（1）3次执行exist（’pi’）的结果一样吗？如果不一样，试解释为什么？

>> pi

ans =

3.1416 >> sin(pi); >> exist('pi') ans =

5 >> pi=0;

>> exist('pi')

ans =

1

>> pi

pi =

>> clear

>> exist('pi')

ans =

5

>> pi

ans =

3.1416

答：3次执行的结果不一样。exist（）函数是返回变量搜索顺序的一个函数。在第一次

执行时返回5代表变量pi是由Matlab构建的变量。在第二次执行时已经通过赋值语句定义

了变量pi，返回1代表pi是工作空间变量。第三次执行前清除了工作空间，此时pi为系统

默认常量，和第一次执行时性质一样，所以又返回5。

（2）圆周率pi是系统默认常量，为什么会被改变为0。

pi=0 为赋值语句，此时pi不再是系统默认常量，而是定义的变量了。

二 MATLAB语言基础

1、向量的生成和运算

练习：使用logspace（）创建1~4π的有10个元素的行向量。

>> A=logspace(0,1.0992,10)

A =

1.0000 1.3247 1.7550

2.3249

3.0799

4.0801

5.4051 7.1603 9.4856 12.5661

2、矩阵的创建、引用和运算

（1）矩阵的创建和引用

练习：创建以下矩阵：A为3×4的全1矩阵、B为3×3的0矩阵、C为3×3的单位矩阵、D为3×3的魔方阵、E由C和D纵向拼接而成、F抽取E的2~5行元素生成、G由F经变

形为3×4的矩阵而得、以G为子矩阵用复制函数生成6×8的大矩阵H。

>> A=ones(3,4),B=zeros(3,3),C=eye(3,3),D=magic(3)

A =

1 1 1 1 1 1 1 1 1 1 1 1

B =

0 0 0 0 0 0

0 0 0

C =

1 0 0 0 1 0

0 0 1

D =

8 1 6

3 5 7

4 9 2

>> E=[C;D], F=E(2:5,:), G=reshape(F,3,4)

E =

1 0 0

0 1 0

0 0 1

8 1 6

3 5 7

4 9 2

F =

0 1 0

0 0 1

8 1 6

3 5 7

G =

0 3 1 1

0 1 5 6

8 0 0 7

>> H=repmat(G,2)

H =

0 3 1 1 0 3 1 1 0 1 5 6 0 1 5 6 8 0 0 7 8 0 0 7 0 3 1 1 0 3 1 1 0 1 5 6 0 1 5 6

8 0 0 7 8 0 0 7 2）矩阵运算

练习：1）用矩阵除法求下列方程组的解x=

>> A=[6 3 4;-2 5 7;8 -1 -3],B=[3;-4;-7]

A =

6 3 4 -2 5

7

8 -1 -3 B =

3 -

4 -7

>> x=A\B

x =

1.0200

-14.0000

9.7200

2）求矩阵的秩；

>> r=rank(A)

r =

3

3）求矩阵的特征值与特征向量>> [X,Lamda]=eig(A)

X =

0.8013 -0.1094 -0.1606 0.3638 -0.6564 0.8669 0.4749 0.7464 -0.4719 Lamda =

9.7326 0 0 0 -3.2928 0 0 0 1.5602

4）矩阵的乘幂（平方）与开方>> A^2

ans =

62 29 33

34 12 6

26 22 34

>> A1=sqrtm(A)

A1 =

2.2447 + 0.2706i 0.6974 - 0.1400i 0.9422 - 0.3494i -0.5815 + 1.6244i 2.1005 - 0.8405i 1.7620 - 2.0970i 1.9719 - 1.8471i -0.3017 + 0.9557i 0.0236 + 2.3845i 5）矩阵的指数与对数（以e为底）

>> Ae=expm(A)

Ae =

1.0e+004 *

1.0653 0.5415 0.6323

0.4830 0.2465 0.2876

0.6316 0.3206 0.3745

>> Ael=logm(A)

Ael =

1.7129 + 0.4686i 0.5305 - 0.2425i 0.5429 - 0.6049i 1.1938 +

2.8123i 0.3658 - 1.4552i -0.5514 -

3.6305i -0.0748 - 3.1978i 0.7419 + 1.6546i 1.8333 +

4.1282i 6）矩阵的提取（取右上三角）与翻转（逆时针转90度）

>> a=triu(A)

a =

6 3 4 0 5

7 0 0 -3 >> a1=rot90(A)

a1 =

4 7 -3 3

5 -1

6 -2 8

3、多维数组的创建及运算

练习：创建三维数组A ，第一页为，第二页为，第三页为。然后用reshape函数重排为数组B，B为3行、2列、2页。

>> a=[1 3;4 2],b=[1 2;2 1],c=[3 5;7 1]

>> A=cat(3,a,b,c)

A(:,:,1) = 1 3 4 2 A(:,:,2) =

1 2

2 1

A(:,:,3) =

3 5

7 1

>> B=reshape(A,3,2,2)

B(:,:,1) = 1 2 4 1 3 2

B(:,:,2) =

2 7

1 5

3 1

三 Matlab数值运算

1、多项式运算

练习：求的商及余多项式。

>> p1=conv([1 0 1],conv([1 3],[1 1]))

p1 =

1 4 4 4 3

>> [q r]=deconv(p1,[1 0 2 1])

q =

1 4

r =

0 0 2 -5 -1

2、多形式插值和拟合

有一组实验数据如附表1-1所示。请分别用拟合（二阶至三阶）和插值（线性和三次样条）的方法来估测X=9.5时Y的值

X 1 2 3 4 5 6 7 8 9 10

Y 16 32 70 142 260 436 682 1010 1432 1960 >> x=1:10;y=[16 32 70 142 260 436 682 1010 1432 1960];

>> p1=polyfit(x,y,1) p1 =

204.8000 -522.4000 >> y1=polyval(p1,9.5)

y1 =

1.4232e+003

>>

p2=polyfit(x,y,2),y2=polyval(p2,9.5)

p2 =

32.0000 -147.2000 181.6000

y2 =

1.6712e+003

>>

p3=polyfit(x,y,3),y3=polyval(p3,9.5) p3 =

2.0000 -1.0000 5.0000 10.0000 y3 =

1.6820e+003

>> y4=interp1(x,y,9.5) y4 =

1696

>> y5=spline(x,y,9.5) y5 =

1682

3、习题

（1）用函数roots 求方程的根

>> roots([1 -1 -1])

ans =

-0.6180

1.6180

（2），在n个节点（n不要太大，如取5~11）上用分段线性和三次

样条插值方法，计算m个插值点（m可取50~100）的函数值。通过数值和图形输出，将两种

插值结果与精度进行比较。适当增加n，再作比较。

>> x=linspace(0,2*pi,8),y=sin(x)

x =

0 0.8976 1.7952 2.6928 3.5904 4.4880 5.3856 6.2832

y =

0 0.7818 0.9749 0.4339 -0.4339 -0.9749 -0.7818 -0.0000

>>

xi=linspace(0,2*pi,100);y0=sin(xi);y1=interp1(x,y,xi);y2=interp1(x,y,xi,'spline'

);

>> plot(xi,y0,'*',xi,y1,'-.',xi,y2)

>> e1=y1-y0;e2=y2-y0; >> plot(xi,e1)

1

2

3

4

5

6

7

-0.1

-0.08-0.06-0.04-0.0200.020.040.060.080.1

>> plot(xi,e2)

01234567

-0.015

-0.01

-0.005

0.005

0.01

0.015

（3）大气压强p 随高度x 变化的理论公式为

，为验证这一公式，

测得某地大气压强随高度变化的一组数据如表所示。试用插值法和拟合法进行计算并绘图，看那种方法较为合理，且总误差最小。

高度/m 0 300 600 1000 1500 2000 压强/Pa

0.9689

0.9322

0.8969

0.8519

0.7989

0.7491

>> x=[0 300 600 1000 1500 2000]; p=[0.9689 0.9322 0.8969 0.8519 0.7989 0.7491]; >> xi=linspace(0,2000);p0=1.0332*exp(-(xi+500)/7756); >> p1=interp1(x,p,xi,'spline'); >> plot(xi,p0,'*',xi,p1) >> e1=p1-p0; >> e=sum(e1.^2) e =

1.8652e-005 拟合法：

>> x=[0 300 600 1000 1500 2000]; p=[0.9689 0.9322 0.8969 0.8519 0.7989 0.7491]; >> P=log10(p)

P =

-0.0137 -0.0305 -0.0473 -0.0696 -0.0975 -0.1255 >> p1=polyfit(x,P,1)

p1 =

-0.0001 -0.0137

>> b=p1(1)/0.4343,a=10.^p1(2)

b =

-1.2863e-004

a =

0.9689

>> xi=linspace(0,2000);p0=1.0332*exp(-(xi+500)/7756);

>> p2=polyval(p1,xi);P2=10.^p2;

>> e2=P2-p0;e=sum(e2.^2)

e =

1.8116e-005

四 Matlab数值运算

1、数值微积分

练习：瑞士地图如图所示，为了算出其国土面积，首先对地图作如下测量：以由西向东方向为X轴，由南到北方向为Y轴，选择方便的原点，并将从最西边界点到最东边界点在X轴上的区间适当划分为若干段，在每个分点的Y方向测出南边界点和北边界点的Y坐标Y1和Y2，根据地图比例尺知道18mm相当于40km，试由测量数据计算瑞士国土近似面积，与其精确值41228km2比较。

>>

x=[7,10.5,13,17.5,34,40.5,44.5,48,56,61,68.5,76.5,80.5,91,96,101,104,106.5,111.5 ,118,123.5,136.5,142,146,150,157,158];

>>

y1=[44,45,47,50,50,38,30,30,34,36,34,41,45,46,43,37,33,28,32,65,55,54,52,50,66,6 6,68];

>>

y2=[44,59,70,72,93,100,110,110,110,117,118,116,118,118,121,124,121,121,121,116,1 22,83,81,82,86,85,68];

>> X=x./18*40;Y1=y1./18*40;Y2=y2./18*40;

>> t1=trapz(X,Y1),t2=trapz(X,Y2), t=t2-t1

t1 =

3.3819e+004

t2 =

7.6328e+004

t =

4.2510e+004

>> expt=t-41228

expt =

1.2819e+003

2、习题

（4）利用梯形法和辛普森法求定积分的值，并对结果进行比较。如果积分区间改为-5~5结果有何不同？梯形积分中改变自变量x的维数，结果有何不同？

>> x=linspace(-3,3);y=exp(-x.^2/2);

>> t=(1/2*pi)*trapz(x,y)

t =

3.9267

>>

q=(1/2*pi)*quad('exp(-x.^2/2)',-3,3) q =

3.9268

>> x=linspace(-5,5);y=exp(-x.^2/2);

>> t=(1/2*pi)*trapz(x,y)

t =

3.9374

>>

q=(1/2*pi)*quad('exp(-x.^2/2)',-5,5) q =

3.9374

>> x=linspace(-3,3,150);y=exp(-x.^2/2);

>> t=(1/2*pi)*trapz(x,y)

t =

3.9268

（5）分别用矩形法、梯形法、辛普森法和牛顿-科茨4种方法近似计算定积分，取n=4，保留4位有效数字。

矩形法：

>> x=linspace(0,1);y=x./(x.^2+4); >> t=cumsum(y)*1/99;T=t(100) T =

0.1126

梯形法：

>> x=linspace(0,1);y=x./(x.^2+4); >> t=trapz(x,y) t =

0.1116

辛普森法：

>> q=quad('x./(x.^2+4)',0,1) q =

0.1116 牛顿-科茨法：

>> q=quadl('x./(x.^2+4)',0,1) q =

0.1116

五 Matlab符号运算

1、符号矩阵创建

练习：分别用sym和syms 创建符号表达式：，。

>> f1=sym('cos(x)+(-(sin(x)^2))^(1/2)')

f1 =cos(x)+(-(sin(x)^2))^(1/2)

>> syms y e t

>> f2=y/exp(-2*t)

f2 =y/exp(-2*t)

2、习题

（2）试创建以下2个矩阵：

6、符号表达式的变量替换

练习：（1）已知,按照自变量x和自变量a，对表达式f分别进行降幂排列。

>> f=sym('(a*x^2+b*x+c-3)^3-a*(c*x^2+4*b*x-1)')

f =(a*x^2+b*x+c-3)^3-a*(c*x^2+4*b*x-1)

>> f1=collect(f),f2=collect(f,'a')

f1 =

a^3*x^6+3*b*a^2*x^5+((c-3)*a^2+2*b^2*a+a*(2*(c-3)*a+b^2))*x^4+(4*(c-3)*b*a+b *(2*(c-3)*a+b^2))*x^3+((c-3)*(2*(c-3)*a+b^2)+2*b^2*(c-3)+a*(c-3)^2-a*c)*x^2+(3*( c-3)^2*b-4*b*a)*x+(c-3)^3+a

f2 =

a^3*x^6+3*(b*x+c-3)*x^4*a^2+(3*(b*x+c-3)^2*x^2-c*x^2-4*b*x+1)*a+(b*x+c-3)^3

8、符号方程的求解

练习：（1）求

>> f=sym('(x^2-1)/(x^2-3*x+2)');

>> limit(f,'x',2)

ans =NaN

（2）求函数f（x）=cos2x-sin2x的积分；求函数的导数。

>> f=sym('cos(2*x)-sin(2*x)');

>> int(f)

ans =1/2*sin(2*x)+1/2*cos(2*x)

>> g=sym('(exp(x)+x*sin(x))^(1/2)');

>> diff(g)

ans =1/2/(exp(x)+x*sin(x))^(1/2)*(exp(x)+sin(x)+x*cos(x))

（3）计算定积分

>> f=sym('sin(x)+2');

>> int(f,'x',0,pi/6)

ans =-1/2*3^(1/2)+1/3*pi+1

（4）求下列线性方程组的解

>> f1=sym('x+y+z=10');

>> f2=sym('3*x+2*y+z=14');

>> f3=sym('2*x+3*y-z=1');

>> g=solve(f1,f2,f3,'x','y','z')

g =

x: [1x1 sym]

y: [1x1 sym]

z: [1x1 sym]

>> g.x

ans =1

>> g.y

ans =2

>> g.z

ans =7

（5）求解当y（0）=2，z（0）=7时，微分方程组的解

>> [g_y,g_z]=dsolve('Dy-z=sin(x)','Dz+y=1+x','y(0)=2','z(0)=7','x')

g_y =cos(x)+6*sin(x)+1/2*sin(x)*x+1+x

g_z =-3/2*sin(x)+6*cos(x)+1+1/2*cos(x)*x

六 Matlab程序设计

1、程序流程控制结构

练习：（1）请把exp2.函数文件用while循环改写。

function s=exp3(x)

n=1;s=0;

while n<=x

s=s+n;

n=n+1;

end

s

（2）用公式求pi的近似值，直到最后一项的绝对值小于10-6为止，试编写其M脚本文件。

k=0;jspi=1;i=3;

while (1/i)>=10e-6

k=k+1;

if rem(k,2)==0

jspi=jspi+1/i;

else

jspi=jspi-1/i;

end

i=i+2;

end

p=4*jspi ,k

2、子函数和参数传递

练习：编写求矩形面积函数rect，当没有输入参数时，显示提示信息；当只输入一个参数时，则以该参数作为正方形的边长计算其面积；当有两个参数时，则以这两个参数为长和宽计算其面积。

function s=mianji(a,b)

switch nargin

case 0

error('没有输入参数')

case 1

s=a*a;

case 2

s=a*b;

end

3、习题

（3）编写一个函数project1.m，其功能是判断某一年是否为闰年。

function ryear(year)

s=0;

if rem(year,4)==0

s=s+1;

end

if rem(year,100)==0

s=s-1;

end

if rem(year,400)==0

s=s+1;

end

if s==1

fprintf('%4d 是闰年.\n',year)

else

fprintf('%4d 不是闰年.\n',year)

end

（4）编制一个函数，使得该函数能对输入的两个数值进行比较并返回其中的最小值。function c = bijiao(a,b)

if nargin==2

if a < b

c=a;

else

c=b;

end

else

error('输入参数不正确')

end

（6）观察以下循环语句，计算每个循环的循环次数和循环结束之后var的值。

var=1;

while mod(var,10)~=0

var=var+1

end

循环次数10，var=10。

var=2;

>> while var<=100

var=var^2;

end

循环次数4，var=256。

>> var=3;

>> while var>100

var=var^2;

end

循环次数0次，var=3。

七 Matlab数据可视化

1、二维图形绘制

练习：写出图A2的绘制方法。

y1=sin(x);y2=cos(x);

plot(x,y1,'r -',x,y2,'m --')

x=linspace(0,4*pi);

y1=sin(x);y2=cos(x);

plot(x,y1,'r --',x,y2,'m -')

ylabel('·ù?è'),xlabel('ê±??')

legend('sinx','cosx')

gtext('\leftarrowsinx')

gtext('\leftarrowcosx')

axis([0 16 -1 1])

xdate=0:0.5:16;ydate=0

line(xdate,ydate,'Color','k','Marker','.')

2、三维曲线和三维曲面绘制

练习：（1）绘制以上空间螺旋线的俯视图、左侧视图和前视图。z=0:0.1:6*pi;

x=cos(z);y=sin(z);

subplot(2,2,1);plot3(cos(z),sin(z),z);

title('èy???ú??')

subplot(2,2,2);plot(cos(z),sin(z));

title('??êóí?')

subplot(2,2,3);plot(cos(z),z);

title('×óêóí?')

subplot(2,2,4);plot(sin(z),z);

title('?°êóí?')

（2）设，求定义域x=[-2，2]，y=[-2，2]内的z值（网格取0.1）。请把z的值用网面图形象的表示出来，如图A3所示。

x=-2:0.1:2;y=x;

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

Z=X.^2.*exp(-(X.^2+Y.^2));

surf(X,Y,Z)

河南城建学院期末考试考查试题 工程热力学试卷答案

1．可逆过程一定是准静态过程； 对，因为一个可逆过程，首先应该是准平衡过程。 2．工质经一不可逆循环后其熵一定增加； 错，工质经一不可逆循环后其熵是不变的。 3．第二类永动机违反了热力学第一和第二定律； 错，第二类永动机并不违背热力学第一定律。 4．节流过程是一个等焓过程； 错，节流过程中，流体在孔口附近产生严重的扰动及涡流，处于不稳定状态，因而不能确定各截面的焓值，所以，尽管节流前后焓不变，但不能把节流过程理解为等焓过程。 5．制冷系数是大于1的数； 错，制冷系数是指制冷循环中制冷量与消耗功的比值，可能大于1，也可能小于1。 6．压缩机的余隙容积对压缩机压缩每公斤工质的理论功没有影响； 对，有余隙容积后，如生产增压比相同、质量相同的同种压缩气体，理论上所消耗的功与无余隙容积时相同。 7．熵减少的过程是可以发生的； 对，只对于孤立系统，熵必定增加或不变，但是对于其他非孤立系统，不一定熵增。 8．孤立热力系熵减少的过程是无法实现的； 对，根据熵增原理，孤立系统熵必定增加或不变。 9．热量不可能从低温热源传向高温热源； 错，根据热力学第二定律，热不可能自发的从低温物体转向高温物体。 10．当蒸汽的温度高于饱和温度时，称该蒸汽为过热蒸汽； 错，当蒸汽温度高于对应压力下的饱和温度时，该蒸汽为过热蒸汽。二．名二、名词解释：（15分，每题3分） 1. 混合气体的分压力； 答：在与混合气体温度相同的情况下，每一种组成气体都独自占据体积V时，组成气体的压力成为分压力。 2. “焓”； 答：在有关热工计算中时常有U出现，为了简化公式和简化计算，把它定义为焓，用符号H表示。pV+ 3. 第二类永动机； 答：从环境大气或海水吸收热量不断转化为机械功，这种单一热源下做功的动力机械称为第二类永动机。 4. 孤立热力系； 答：任何一个热力系（闭口系、开口系、绝热系、非绝热系），总可以将它连同与其相互作用的一切物体组成一个复合系统，该复合系统不再与外界有任何形式的能量交换和质量交换，我们称该复合系统为孤立系统。 5. 理想气体； 答：理想气体是一种实际上不存在的假想气体，其分子是些弹性的、不具体积的质点，分子间相互没有作用力。 三．（15分） 空气由初态压力为0.1MPa，温度20℃，经2级压缩机压缩后，压力提高到2MPa。若空气进入各级气缸的温度相同，且各级压缩过程的多变指数均为1.2，求最佳的中间压力为多少？并求生产1kg质量的压缩空气所消耗的理论功？求各级气缸的排气温度为多少？ 解：最佳中间压力MPa p p p447 .0 2 1.0 2 1 2 = ? = =； 理论比功：() KJ p p p p T R n n w n n n n g c 21 . 286 2 1.0 447 .0 2 20 15 . 273 10 06 . 287 1 2.1 2.1 2 1 2.1 1 2.1 3 1 2 3 1 1 2 1 = ? ? ? ? ? ? ? ? - ? ? ? ? ? ? + ? ? ? - = ? ? ? ? ? ? ? ? ? ? - ?? ? ? ? ? + ?? ? ? ? ? - = - - - -

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上机习题2

实验二 基本矩阵操作 实验目的： ① 掌握matlab 变量和数据操作； ② 掌握matlab 矩阵的创立、拆分及特殊矩阵； ③ 掌握matlab 运算，掌握matlab 在矩阵分析中的应用，掌握稀疏矩阵的存储方式和创建方法； ④ 掌握矩阵的基本运算和结构体数组的创建。 ⑤ 了解字符串处理函数。 实验要求:给出程序和实验结果。 实验内容： 1、利用列向量()1,2,3,,10T 建立一个范得蒙矩阵A ，并利用位于矩阵A 的奇数行偶数列的元素建立一个新的矩阵B ，须保持这些元素的相对位置不变。 2、矩阵的基本运算与点运算的区别。 3、给出矩阵的两种存储方式的联系和区别，这两种存储方式在实际应用中主要应用于具有什么特点的矩阵？ 4*．将字符串'very good'转换为等值的整数。 5．按水平和竖直方向分别合并下述两个矩阵： 6． 分别删除第5题两个结果的第2行。 7． 分别将第5题两个结果的第2行最后3列的数值改为[11 12 13]。 8． 分别查看第5题两个结果的各方向长度。 9*． 分别判断pi 是否为字符串和浮点数。 10．分别将第5题两个结果均转换为2?9的矩阵。 11．计算第5题矩阵A 的转置。 12．分别计算第5题矩阵A 和B 的A +B 、A .* B 和A \B 。 13．判断第5题矩阵A 和B 中哪些元素值不小于4。 14*．分别用函数strcat()和矩阵合并符合并如下字符串：' The picture is '和' very good '。 15*．创建字符串数组，其中元素分别为‘Picture ’和'Pitch '。 16*．在第14题结果中查找字符串'e'。 17*．在第15题结果中匹配字符串'Pi'。

河南城建学院编译原理试卷及答案

河南城建学院2010学年第一学期期末考试 《编译原理》试题（A卷） 一、填空题：（每空1分，共10分） 1、符号表项的组织常采用线性法、二分法和（散列法）。 2、整个编译过程可以划分成五个阶段：（词法分析）、语法分析阶段、（语义分析及中间代码生成）、（代码优化）和目标代码生成阶段。 3、对于文法G，仅含终结符号的句型称为（句子）。 4、逆波兰式ab+c+d*e-所表达式为（(a+b+c)*d-e）。 5、语言翻译常用的两种形式是（编译）和（解释）。 6、词法分析器输出的是单词符号，语法分析器输出的是（语法单元）。 二、选择题：（每空2分，共10分） 1、3型文法是（D）是语法分析使用的文法。 A．短语文法 B.上下文有关文法 C.上下文无关文法 D.正规文法 2、语法分析是依据语言的（A）规则进行的，中间代码产生是依据语言的（）规则进行的。 A．语法，推导B．语义，产生式C．语法，语义D．推导，产生式 3、错误“变量类型声明不一致”将在（ C ）阶段发现。 A．词法分析 B.语法分析 C.语义分析 D.目标代码生成 4、下列（D）不是数据空间的使用方法和管理方法 A．静态存储分配 B.栈式动态存储分配 C.堆式动态存储分配 D.段页式存储分配 三、计算题：（每题6分，共24分） 1、对给定正规表达式b*(d∣ad) (b∣ab)+构造其NFA M。

解：先用R+=RR*改造题设的正规表达式为b*(d ∣ad) (b ∣ab) (b ∣ab)*，然后构造其 NFA M ，如图 2、试给出下列语句的四元式序列： if (a<0∧b ＞5) X[1,1]＝=1; else X[3,2]＝0; 其中，X 是10×20的数组(每维下界为1)且按行存放；一个数组元素占用两个字节，机器按字节编址。 解： 100 (j<, a, 0, 102) 101 (j,_, _,110) 102 (j>,b,5,104) 103 (j, _,_,110) 104 (*,1,40,T1) 105 (*,1,2,T2) 106 (+,T1,T2,T3) 107 (?,X,42,T4) 108 ([ ]=,1, _,T4[T3]) 109 (j, _,_,115) 110 (*,3,40,T1) 111 (*,2,2,T2) 112 (+,T1,T2,T3) 113 (?,X,42,T4) 114 ([ ]=,0, _,T4[T3]) 115 b*(d | ad)(b | ab)(b | ab)* X Y (b | ab)* d | ad b* b | ab X 1 2 3 Y ab ad ε ε b ε ε d b | ab b X 41 2 3 5 Y ε ε b X 41d 26 ε ε b b 35Y 7 8 a d b a a 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 =

Matlab上机实验

Matlab 上机实验 一、 实验目的 1、 掌握绘制MATLAB 二维、三维和特殊图形的常用函数； 2、 熟悉并掌握图像输入、输出及其常用处理的函数。 二、 实验内容 1 绘制函数的网格图和等高线图。42 2cos cos y x ye x z +-= 其中x 的21个值均匀分布在[-5,5]范围，y 的31个值均匀分布在[0,10]，要求将产生的网格图和等高线图画在同一个图形窗口上。 2 绘制三维曲面图，使用纯铜色调色图阵进行着色，并进行插值着色 处理。?????===s z t s y t s x sin sin cos cos cos 230,20ππ≤≤≤≤t s 3 已知 ???????>++≤+=0),1ln(210,22x x x x e x y π 在-5<=x<=5区间绘制函数曲线。 4 已知y1=x2，y2=cos(2x)，y3＝y1*y2，其中x 为取值-2π~2π的等差数列(每次增加0.02π)，完成下列操作： a) 在同一坐标系下用不同的颜色和线型绘制三条曲线，给三条曲线添加图例； b) 以子图形式，分别用条形图、阶梯图、杆图绘制三条曲线，并分别给三个图形添加标题“y1=x^2”，“y2=cos(2x)”和“y3＝

y1*y2”。 5 在xy 平面内选择区域[][],,-?-8888 ，绘制函数 z =的三 种三维曲面图。 6 在[0,4pi]画sin(x),cos(x)(在同一个图象中); 其中cos(x)图象用红色小圆圈画.并在函数图上标注 “y=sin(x)”, “y=cos(x)” ,x 轴,y 轴,标题为“正弦余弦函数图象”. 7 分别用线框图和曲面图表现函数z=cos(x)sin(y)/y ，其中x 的取值为 [-1.5pi,1.5pi]，y=x ，要求：要有标题、坐标轴标签 8 有一组测量数据满足-at e =y ，t 的变化范围为0~10，用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线，并加入标题和图列框（用代码形式生成） 9 2 2y x xe z --=，当x 和y 的取值范围均为-2到2时，用建立子窗口 的方法在同一个图形窗口中绘制出三维线图、网线图、表面图和带渲染效果的表面图 10 x= [66 49 71 56 38]，绘制饼图，并将第五个切块分离出来。 11 用sphere 函数产生球表面坐标，绘制不通明网线图、透明网线图、表面图和带剪孔的表面图。 12 以自己的个人画像或照片(JPG)为对象，读入该图像并了解图像的信息，同时利用所学函数对其进行灰度、二值、旋转及缩放等处理，并以PNG 形式输出。

matlab上机实验指导书

MATLAB应用基础实验指导书

第一章 MATLAB及其工作环境介绍 (1) 1.1 MATLAB简介 (1) 1.2 MATLAB的工作环境介绍 (1) 1.3 MATLAB的基本管理命令 (4) 第二章 MATLAB的数值计算功能 (5) 2.1 变量与赋值语句 (5) 2.2 MATLAB矩阵 (5) 2.3 MATLAB表达式 (10) 2.4 MATLAB常用数学函数 (11) 2.5 矩阵的基本运算 (12) 2.6 数组运算 (16) 2.7 多项式及其运算 (17) 第三章 MATLAB程序设计入门 (19) 3.1 M文件 (19) 3.2 数据的输入输出 (21) 3.3 全局变量和局部变量 (23) 3.4 程序流程控制 (23) 第四章 MATLAB的符号运算功能 (28) 4.1 建立符号对象 (28) 4.2 符号算术运算 (29) 4.3 符号微积分运算 (32) 4.4 符号函数的可视化 (34) 第五章 MATLAB的可视化功能 (37) 5.1 二维图形 (37) 5.2绘制三维图形 (42) 5.3 特殊坐标图形 (44) 5.4 图形句柄 (45)

第一章 MATLAB及其工作环境介绍 1.1 MATLAB简介 MATLAB是matrix和laboratory前三个字母的缩写，意思是实验室矩阵。MATLAB 语言是一种广泛应用于工程计算及数值分析领域的新型高级语言，自1984年由美国MathWorks公司推向市场以来，经过十多年的发展与完善，MATLAB已发展成为由MATLAB语言、MATLAB工作环境、MATLAB图象处理系统、MATLAB数学函数库和MATLAB 应用程序接口五大部分组成的集数值计算、图形处理、程序开发为一体的功能强大的体系。MATLAB由“主包”和三十多个扩展功能和应用学科性的工具箱组成。 MATLAB具有以下基本功能： ●数值计算功能 ●符号计算功能 ●图形处理及可视化功能 ●可视化建模及动态仿真功能 MATLAB语言是以矩阵计算为基础的程序设计语言，语法规则简单易学。其指令格式与数学表达式非常相近，用MATLAB编写程序犹如在便笺上列写公式和求解，因而被称为“便笺式”的编程语言。另外，MATLAB还具有功能丰富和完备的数学函数库及工具箱，大量繁杂的数学运算和分析可通过调用MATLAB函数直接求解，大大提高效率，其程序编译和执行速度远远超过了传统的C和FORTRAN语言，因而用MATLAB 编写程序，往往可以达到事半功倍的效果。在图形处理方面，MATLAB可以给数据以二维、三维乃至四维的直观表现，并在图形色彩、视角、品性等方面具有较强的渲染和控制能力，使技术人员对大量原始数据的分析变得轻松和得心应手。 MATLAB的上述特点，使它深受工程技术人员及科技专家的欢迎，并成为应用学科计算机辅助分析、设计、仿真、教学等领域不可缺少的基础软件。目前MATLAB已成为国际上公认的最优秀的科技应用软件。 1.2 MATLAB的工作环境介绍 一、MATLAB的工作环境

河南城建学院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的数值显示格式

河南城建学院《中外城建史》试题B答案

平顶山工学院2005—2006学年第一学期期末考试 《中外城建史》试题B答案 供城市规划专业0513041/2 班使用 一、填空题（25分） 1、商城 2、汉朝 3、瓦子 4、秦 5、里弄 6、帝国主义共占城市 7、新巴比伦王国 8、帝国时期9、要塞型、城堡型、商业交通型10、田园城市、带形城市、广亩城市11、中心阁12、“功能城市”13、建筑十书14、盖迪斯 15、人的需要和以人为出发点的价值衡量是一切建设工作成功的关键 16、人的活动要求流动的、连续的综合空间 17、适度的集中、均衡的分散 二、选择题（10分） CAACC AACBA 三、简答题（40分） 1、1）唐长安城主要受儒家思想影响，儒家思想的“居中不偏”，“不正不威”的思想影响到城市布局。 长安城布局采用中轴线对称布局原则，宫城、皇城居中在城市最北面，其它设施则按照封建礼制顺序布局。道路是规整的方格网状系统，有集中的市，市和坊都有封闭的坊墙。 2）宋东京城的布局主要立足点在促进商业的交易，所以城市总体的布局没有严格按照儒家思想的学说。除了皇城、宫城内部布局采取轴线布局，外城的布局比较自由，大街道一般都与商业结合，形成商业街，有集中的市，也有分散的商业街道，市和坊无坊墙。 2、答：城市的产生是在手工业从农业分离后，剩余产品的出现，产生了私有制，才出现了以交换产品的场所-市。所以早期的城市是从农业居民点的基础上发展而来的。城市与农业的关系很紧密，因此大部分城市都靠近河流。秦汉以后国内有统一的市场，商品流转数量大，商业的运输主要是靠自然水运条件，所以在通航河道的重要渡口或交汇处往往形成较大的商业都会，因此这个时期城市的建设靠近河流主要是方便商业的运输。 3、城市的布局体现了明确的侵略意图，有的偏重于军事，有的偏重于商业贸易，有的偏重于交通运输；城市的布局体现了民族矛盾，中国居住区和外国居住区在建筑质量上，人口密度，建筑密度，绿地标准，公用设施水平上有显著差别；城市的布局有明显的商业性质，密集的道路网，众多的临界街坊，反映了资本主义城市以营利为目的的特征。 4、雅典卫城建立在城内一个陡峭的山顶上，用乱石在四周砌成挡土墙形成的大平台，山势险要，只有西面有个上下的通道。卫市的主要建筑是供奉城邦保护神雅典娜的帕提农神庙。卫城建筑的布局不是刻板简单的轴线关系，而是按照人们祭祀是行走的路线设计的，利用地形，把最好的角度展现给人们，使没段路程展现给人们面前的都是不同的建筑优美景象。 5、圣马可广场最突出的艺术成就是建筑群体的空间组合达到统一和谐的情调。大广场与靠海的小广场通过钟塔过渡，靠海的大广场和靠海的小广场用用台阶和一对狮子雕像进行过渡。同时广场四周建筑底层采用外廊式做法，使广场空间与建筑室内空间有个过渡。在艺术处理方面，高耸的直塔作为竖直线条，与周围这些建筑物众多的水平线条形成对比。总之，圣马可广场运用了不同空间的互迭手法，视觉上采用相似、对比的手法，使得建筑群体的组合和环境达到和谐统一的艺术高峰。 6、1）、进步性：提出了通过限制大城市规模来控制大城市恶性膨胀问题；通过在大城市中布置大面积的绿化带来改变城市环境质量；在城市中心设置公共活动场所，满足人们室外活动的需要；工业区，铁路干线放在城市边缘地带，避免干扰城市。 2）、霍华德没有认清资本主义社会的以营利为目的的本质，资产阶级不会把大量的资金投入毫无利润可赚的绿化设施和公共活动场地，更不会对城市人口进行疏散

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上机实验答案 (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

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)

MATLAB上机习题一

MATLAB上机习题一 请按以下步骤完成上机实验： 1）在FTP上下载“MATLAB上机习题一.doc”文件，所有习题列在该文件内； 2）在MATLAB中完成所有习题，并将屏幕截图粘贴到相应习题后面； 3）如果习题是问答题，请将答案写在题目后； 4）如果有的习题要求提供脚本文件，请将脚本文件内容拷贝到相应习题后； 5）将文件保存并重命名为“自己的学号-姓名”，例如“20110771-张三.doc”； 6）上传该文件到FTP的相关目录。 1. 运行MATLAB软件，拖放、关闭界面上的子窗口，并恢复到原始试图。 2. 采用鼠标及命令两种方式将桌面添加到MATLAB搜索路径列表的起始及最后位置。 3. 采用鼠标及命令方式将当前工作目录设置为桌面指向的文件夹。 4. 通过使用帮助确定内置变量ispc 的含义。 6. 观察MATLAB中关键字、字符串、注释的字体显示颜色。 5. 创建double类型的变量，并进行计算 1）a=87，b=190，计算a+b、a-b、a*b；

2）创建uint8类型的变量，数值与（1）中相同，进行相同的计算，观察计算结果与预想的是否一致。 6. 计算如下表达式： 1）()sin 60 2）3e 3）3cos 4π?? ??? 4）2 7562323336 +?-?? 7. 设u=2，v=3，计算： 1）4log uv v 2）() 2 2u e v v u +- 3 8. 计算如下表达式： 1）()()3542i i -+ 2）()sin 28i - 9. 判断下面语句的运算结果，并与MATLAB 计算结果做比较：

1）420< 2）420≤ 3）420== 4）4~20= 5）''''b B < 10. 设a=39，b=58，c=3，d=7，判断下面表达式的值与MATLAB 计算结果做比较： 1）a b > 2）a b < 3）&&a b b c >> 4）a d == 5）|a b c > 6）~~d 11. 判断下列变量名哪些是合法的MATLAB 变量： 1）fred 7）fred! 2）book_1 8）book-1 3）2ndplace 9）Second_Place 4）#1 10）No_1 5）vel_5 11）vel.5 6）tan 12）while 提示，可以使用isvarname 对上述变量名进行检验。 12. 编写脚本程序，命名为barbell.m ，完成如下计算： 如图1所示的杠铃，

河南城建学院《工程经济学》期末考试试卷

河南城建学院2011—2012学年第一学期期末考试（试） 《工程经济学》试题（A卷） 本套试卷共 3 页 答题卡 （请把相关试题的答案填写在答题卡相应位置） 一、单项选择题 二、多项选择题 一、单项选择题(共20题，每题2分，共40分。每题的备选项中，只有一个最符合题意，请把正确答案填写在答题卡相应位置。) 1.现金流入与现金流出之差称为（）。 A．现金流量 B．净现金流量 C．净现金 D．纯现金 2.资金时间价值是( )。 A．利息 B．资金随时间推移在经营活动中产生的增值（或减值） C．通货膨胀 D．利率 3.基准收益率应（）单位资金成本和单位投资的机会成本，这样才能使资金得到最有效的利用。 A．高于 B．低于 C．等于 D．不低于 4.根据对项目不同方案的敏感性分析，投资者应选择（）的方案实施。 A.项目盈亏平衡点高，抗风险能力适中 B.项目盈亏平衡点低，承受风险能力弱 C.项目敏感程度大，抗风险能力强 D.项目敏感程度小，抗风险能力强 5.年利率8％，按季度复利计息，则半年期实际利率为（）。 A．4.00% B．4.04% C．4.07% D．4.12% 6.盈亏平衡分析是基于( )。 A．收入等于成本 B．产量等于收入 C．收入等于可变成本 D．产量等于成本 7. 盈亏平衡分析适应( )。 A．财务评价 B．国民经济评价 C．一般经济评价 D．风险评价

8.单因素敏感性分析是指( )。 A．一个因素变化，其他因素不变 B．所有因素不变 C．评价指标变化 D．一个因素变化，其他因素同时变化 9.建设投资是由( )构成。 A．工程费用、工程建设其他费用、预备费 B．工程费用 C．工程建设其他费用、预备费 D．固定资产 10.可行性研究是建设程序中( )。 A .建设前期 B. 勘探设计期 C. 建设实施期 D. 建成投产期 11.清偿能力分析指标主要有：利息备付率、（）、偿债备付率。 A．资产负债率 B．利息偿还期 C．投资回收期 D．偿债额 12.将技术方案经济效果评价分为静态分析和动态分析的依据是( )。 A.评价方法是否考虑主观因素 B.评价指标是否能够量化 C.评价方法是否考虑时间因素 D.经济效果评价是否考虑融资的影响 13.若净现值为负数，表明该投资项目( )。 A．它的投资报酬率小于零，不可行 B．为亏损项目，不可行 C．它的投资报酬率不一定小于零，因此也有可能是可行方案 D．它的投资报酬率没有达到预定的贴现率，不可行 14.财务评价与国民经济评价共同点( )。 A．评价角度 B．费用与效益的划分 C．价格体系 D．评价的基础 15.敏感系数高，表示项目效益对该不确定因素的影响（）。 A．敏感程度低 B．敏感程度高 C．不会有影响 D．敏感程度与之无关 16.某施工企业希望从银行借款500万元，借款期限2年，期满一次还本。经咨询有甲、乙、丙、丁四家银行愿意提供贷款，年利率均为7%。其中，甲要求按月计算并支付利息，乙要求按季度计算并支付利息，丙要求按半年计算并支付利息，丁要求按年计算并支付利息。则对该企业来说，借款实际利率最低的银行是( ) A.甲 B.乙 C.丙 D.丁 17.反映技术方案资本金盈利水平的经济效果评价指标是 ( ) A.内部收益率 B.总投资收益率 C.资本积累率 D.资金本净利润率 18.反映企业在一定会计期间经营成果的报表是( ) A.利润表 B.资产负债表 C.现金流量表 D.所有者权益变动表 19.单因素敏感分析过程包括：1确定敏感因素；2确定分析指标；3选择需要分析的不确定性因素；4分析每个不确定因素的波动程度及其对分析指标可能带来的增减变化情况。正确的排列顺序是（）。 A.3241 B.1234 C.2431 D.2341 20.某化工建设项目设计年生产能力5万吨，预计年固定总成本为800万元，产品销售价格1500元/吨，产品销售税金及附加为销售收入的10%，产品变动成本1150元/吨，则该项目用生产能力利用率表示的盈亏平衡点是（）。 A．100% B．40% C．80% D．55% 二、多项选择题（共10题，每题3分，共30分。每题的备选项中，有2个或2个以上符合答案，请把正确答案填写在答题卡相应位置。） 1.折旧方法有（） A.平均年限法 B.双倍余额递减法 C.工作量法 D.年数总和法 E.以上都是 2.绘制现金流量图需要把握的现金流量的要素有（）

Matlab上机实验二

1. 对Runge 函数2()1/(125)R x x =+在区间[-1,1]作下列插值逼近，并和()R x 的图 像进行比较，并对结果进行分析。 （1） 用等距节点1,0.1,020i x ih h i =-+=≤≤，绘出它的20次Newton 插值 多项式的图像。 （2） 用节点21cos()(0,1,2...,20)42 i i x i π+==，绘出它的20次Lagrange 插值多项式的图像。

（3） 用等距节点1,0.1,020i x ih h i =-+=≤≤，绘出它的分段线性插值函数 的图像。 （4） 用等距节点1,0.1,020i x ih h i =-+=≤≤，绘出它的三次自然样条插值 函数的图像。

结果分析： 等距节点的Newton插值在边界附近呈现不收敛的状态，也就是所谓的Runge 现象。这是因为被插函数2 ()1/(125) R x x =+及其任意阶导数并不是一致有界的。 由Lagrange插值的余项 (n1) 1 () ()()()() (1)! n n n f R x f x L x x n ξ ω + + =-= + ，由于高阶导数值 非常大，在等距节点的高次插值中，余项的值很大，亦即误差特别大，表现为Runge现象。在实际中，除非被插函数的光滑性十分好，否则等距节点的高次插值一般不使用。 切比雪夫多项式零点插值没有出现Runge现象，这是因为对于Lagrange插 值的余项 (n1) 1 () ()()()() (1)! n n n f R x f x L x x n ξ ω + + =-= + ，切比雪夫多项式零点插值可以 将 1() n x ω +项控制在 1 2n 以下（ 1 [1,1] max() n x x ω + ∈- = 1 2n ），由此可得切比雪夫多项式零点 插值的误差范围 (n1) [1,1] () max() 2(1)! n n x f R x n ξ + ∞ ∈- ≤ + ，比等距节点的高次插值小很多，因此 对于被插函数的光滑性要求要比等距节点的插值也低很多。 分段低次线性插值是常用的一种插值方式。线性插值的优点在于计算量小，但是从图中可以看出，其光滑性不够。三次样条插值要求其一阶导数是连续的，因此其准确性好，光滑性也较优。

MATLAB上机题汇总(基本题,抽3题)合肥学院

1设置matlab 的工作环境，将工作目录设置为d:\work ，添加搜索目录d:\example 2在matlab 的命令窗口里完成如下计算，其中t 的值分别取-1,0,1，表达式如下： 4/3)2 t y e π-= 3自行产生一个5行5列的数组，得到最中间的三行三行矩阵。 4用magic 产生一个5*5的矩阵，将这个矩阵的第二行与第三行互换位置 5求方程组的根 x 1+4x 2-3x 3=2 2x 1+5x 2-x 3=11 x 1+6x 2+x 3=12 6已知：一个多项式的系数向量是p=[1 -6-72 -27]，求这个多项式的根。 7已经两个多项式的系数分别是：[1 2 3 4]和[1 4 9 16]，请求这两个多项式的乘积，及商和余数。 8给定一个多项式的根是[-5 -3+4i -3-4i],求原来的多项式 9 A=[2 3 4； 1 5 7； 6 2 5] 用什么函数，保证第一列排序的时候，其他列跟着变化。最后的结果是： [1 5 7； 2 3 4； 6 2 5] 10补充题：电路分析 电路分析常常涉及到对方程组的求解，这些方程常常是利用描述电流进入和离开节点的电流方程，或者描述电路中网络回路上的电压的电压方程得出的。下图描述了3个回路电压的方程

方程式如下： 11121222132423432 532()0()()0()0V R I R I I R I I R I R I I R I I R I V -++-=??-++-=??-++=? 假设5个电阻值为已知，2个电压值也为已知，求3个电流值。 11自行产生一个5行5列的数组，用两种方法得到最中间的三行三行矩阵。 12根据a=reshape(-4:4,3,3)做一个矩阵，然后（1）取出所有大于0的元素构成一个向量（可推广到从一个矩阵里筛选出符合条件的元素组成一个向量）（2）将原矩阵中大于0的元素正常显示，而小于等于0的元素全部用0来表示(可推广到将原矩阵中不符合条件的全用0来表示，符合条件的值不变)。(思考：大于0的正常显示，小于等于0的用-10来表示) 13建立如下一个元胞数组，现在要求计算第一个元胞第4行第2列加上第二个元胞+第三个元胞里的第二个元素+最后一个元胞的第二个元素。 a={pascal(4),'hello';17.3500,7:2:100} 14建立一个结构体的数组，包括3个人，字段有姓名，年龄，分数，其中分数由随机函数产生一个3行10列的数据表示了有10门课程，每门课程有三个阶段的分数。问题是： 问题1，如何找到第2个人的分数并显示出来 问题2，如何找到第2个人的每门课程3个阶段的平均分数并显示出来 问题3，全班同学（指这3个学生）的10门课程的每门课程的平均分如何计算出来？要求放到一个数组里。 问题4，找到这个班所有同学的姓名放到一个元胞数组里 15给定一个图像文件，格式是jpg ，通过inportdata 引入这个文件，查看引入后数据保存是一个数组，是100*100*3的一个数组，这是一个三维的数组，表明有100*100个点，每个点有RGB 三个方面决定的。现在，要求将这些点的中间部分30个象素宽度的正方形用黑点表示。 16手动构造一个长宽页是1000*1000*3的图像，每个图像的点有三个0-255之间的随机值构成。然后使用image 函数显示这个图像，再用imwrite 函数保存这个图像，使这个图像用图片浏览器也能够查看。 17给定一个矩阵 a = 1 2 3 4 5 6 编写一个M 函数，要求输入是a ，输出有三个：平均数，标准差，秩。程序运行

河南城建学院07240712燃气燃烧与应用B卷

☆ ☆ 密 封 线 内 不 要 答 题 ☆ ☆ 姓 名 学 号 班 级 河南城建学院2009—2010学年第二学期期末考试 《燃气燃烧与应用》试题（B 卷） 本套试卷共 3 页 一、填空题（每空1分，共19分）： 1．燃气成分一定时，层流扩散火焰的长度主要取决于 。 2.燃气热水器分 和 。 3.层流火焰传播理论包括 、扩散理论和 。 4.按空气供给的方法分类，可将燃烧器分为 、 、 。 5.影响火焰传播速度的因素有温度、 、压力、 、湿度和惰性气体的影响。 6.自由射流的实质是 。 7．确定气流紊动性的三个基本量是 、 、 。 8．常用的减噪装置有 、 和 。 9．燃烧设备运行的强度通常可分为 和 。 二、名词解释（每题3分，共15分）： 1.回火： 2.热量计温度： 3.华白数： 4.支链着火： 5、脱火速度极限 三、简答题(每小题5分，共20分) 1.简述大气式燃烧器中喷嘴的作用及结构形式。

☆☆密封线内不要答题☆☆ 2.简述进行完全预混式燃烧的条件。 3.简述在实际工程上,减少氮氧化合物生成量的措施。 4.简述强化燃烧过程的途径。 四、问答题（每小题7分，共28分） 1.什么是旋流，产生旋流的方法有哪些？ 2.试论述压差式热水器的工作原理。

☆☆密封线内不要答题☆☆3.试论述大气式燃烧器的特点及应用范围。 4. 试论述紊流火焰特点及紊流火焰比层流火焰传播快的原因。 五、计算题（18分） 1．已知某焦炉煤气容积成分如下： %, 7.3 %, 5. 22 %, 8.6 %, 57 4 2 4 2 H C CH CO H% 7.7 %, 3.2 2 2 N CO。 (1)求燃烧3 5Nm该焦炉煤气所需要的理论空气量； (2)当过剩空气系数1.1 = α时,求燃烧3 5Nm该焦炉煤气所需的实际空气量.