建筑工程管理工程数学作业精编Word文档下载推荐.docx

建筑工程管理工程数学作业精编Word文档下载推荐.docx

《建筑工程管理工程数学作业精编Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《建筑工程管理工程数学作业精编Word文档下载推荐.docx（12页珍藏版）》请在冰豆网上搜索。

0.54380.53490.53440.5370

0.55460.54610.54560.5481

0.56490.55670.55620.5587

0.57470.56680.56630.5688

0.58400.57650.57600.5784

0.59290.58570.58530.5876

0.60140.59460.59410.5965

0.60960.60310.60270.6049

0.61750.61130.61090.6131

0.62510.61920.61880.6210

0.63250.62690.62650.6286

0.63960.63430.63390.6360

0.64660.64150.64110.6432

0.65330.64850.64820.6501

0.65990.65540.65500.6569

0.66640.66200.66170.6636

0.67270.66860.66830.6700

0.67880.67500.67470.6764

0.68490.68120.68100.6826

0.69090.68740.68710.6887

0.69670.69350.69320.6947

0.70250.69940.69920.7007

0.70820.70530.70510.7065

0.71390.71110.71090.7123

0.71950.71690.71670.7180

0.72500.72260.72240.7236

0.73050.72820.72800.7292

0.73590.73370.73360.7347

0.74130.73930.73910.7402

0.74670.74470.74460.7456

0.75200.75020.75000.7510

0.75730.75560.75540.7564

0.76250.76090.76080.7617

0.76780.76630.76620.7670

0.77300.77160.77150.7723

0.77820.77690.77680.7775

0.78330.78210.78200.7828

0.78850.78740.78730.7880

0.79370.79260.79250.7931

0.79880.79780.79770.7983

0.80390.80300.80290.8035

0.80900.80810.80810.8086

0.81410.81330.81320.8137

0.81920.81840.81840.8189

0.82430.82360.82350.8240

0.82930.82870.82860.8291

0.83440.83380.83370.8341

0.83950.83890.83890.8392

0.84450.84400.84400.8443

0.84960.84910.84900.8494

0.85460.85420.85410.8544

0.85960.85920.85920.8595

0.86470.86430.86430.8645

0.86970.86940.86930.8696

0.87470.87440.87440.8746

0.87980.87950.87950.8797

0.88480.88450.88450.8847

0.88980.88960.88950.8897

0.89480.89460.89460.8947

0.89990.89960.89960.8998

0.90490.90470.90470.9048

0.90990.90970.90970.9098

0.91490.91470.91470.9148

0.91990.91980.91980.9199

0.92490.92480.92480.9249

0.92990.92980.92980.9299

0.93490.93480.93480.9349

0.93990.93990.93990.9399

0.94500.94490.94490.9449

0.95000.94990.94990.9499

0.95500.95490.95490.9549

0.96000.95990.95990.9600

0.96500.96490.96490.9650

0.97000.96990.96990.9700

0.97500.97500.97500.9750

0.98000.98000.98000.9800

0.98500.98500.98500.9850

0.99000.99000.99000.9900

0.99500.99500.99500.9950

达到相同精度J迭代的迭代次数为：

4024

达到相同精度G-S迭代的迭代次数为：

2000

达到相同精度sor迭代的迭代次数为：

478

sor迭代最佳松弛因子：

1.7000

由结果可见对于此题达到相同精度迭代次数sor迭代<

G-S迭代<

J迭代

练习题2：

设，对于线性方程组建立迭代法

（1）求出的范围使迭代法收敛，

（2）求出最优使迭代法的渐进收敛速度最大。

迭代矩阵B=

记（=1,2,3,…,n）为A的特征值，为A的最大和最小的特征值

∵A对称正定，∴其特征值为正，>

又B=

所以B的特征值为

所以﹤1

解得0﹤﹤

又为A的最大特征值

所以使迭代法收敛的的范围是

0﹤﹤

（2）因为迭代矩阵的谱半径越小，迭代收敛越快.

所以

使迭代法的渐进收敛速度最大,最优松弛因子

练习题3：

对某电路的分析，能够归结为下面的线性方程组，其中

R（1,1）=31；

R（1,2）=-13；

R（1,6）=-10；

R（2,1）=-13；

R（2,2）=35；

R（2,3）=-9；

R（2,5）=-11；

R（3,2）=-9；

R（3,3）=31；

R（3,4）=-10；

R（4,3）=-10；

R（4,4）=79；

R（4,5）=-30；

R（4,9）=-9；

R（5,4）=-30；

R（5,5）=57；

R（5,6）=-7；

R（5,8）=-5；

R（6,5）=-7；

R（6,6）=47；

R（6,7）=-30；

R（7,6）=-30；

R（7,7）=41；

R（8,5）=-5；

R（8,8）=27；

R（8,9）=-2；

R（9,4）=-9；

R（9,8）=-2；

R（9,9）=29；

V=（-15,27,-23,0,-20,12,-7,7,10）T

其余元素为零。

要求：

（1）用高斯列主元消去法求解该方程组；

（2）用SOR方法迭代求解该方程组，误差，近似最佳松弛因子由试算法确定，设

SOR方法迭代法近似最佳松弛因子为w=1.18,迭代次数为n=12

高斯列主元消去法和SOR方法迭代比较:

高斯法值SOR法值

-0.2892-0.2892

0.34540.3454

-0.7128-0.7128

-0.2206-0.2206

-0.4304-0.4304

0.15430.1543

-0.0578-0.0578

0.20110.2011

0.29020.2902

第四章作业

分别用不动点迭代法和牛顿迭代法求解方程

其中初值，计算精度为。

不动点迭代法使用进行迭代

f=inline（'

-0.9*x^2+1.7*x+2.5'

'

x'

）;

df=inline（'

-1.8*x+1.7'

x0=5;

e=1e-6;

n1=0;

x1=17/9+25/9/x0;

n1=n1+1;

while（norm（x1-x0）>

=e）&

（n1<

=1000）

x0=x1;

disp（[n1,x1]）;

end

x_b=x1

n=0;

x1=x0-feval（f,x0）/feval（df,x0）;

n=n+1;

（n<

%输出结果

clc;

fprintf（'

不动点迭代次数为n=%g'

n1）;

结果为x=%g\n'

x_b）;

牛顿迭代迭代次数为n=%g，结果为'

n）;

x=%g\n'

x1）;

结果：

不动点迭代次数为n=15结果为x=2.8601

牛顿迭代迭代次数为n=5，结果为x=2.8601

分别用不动点迭代法、牛顿迭代法和逆Broyden秩1方法求解方程组

且对结果进行比较。

其中逆Broyden秩1方法的初始矩阵分别取I和，取初值，计算精度。

迭代次数

x1

x2

x3

不动点迭代

17.0000

0.9055

1.0852

0.6722

牛顿迭代

4.0000

逆Broyden秩1方法

（H0=I）

16.0000

（H0=F’（x0）’）

9.0000

理想和非理想气体定律

考虑固定容器中壹定量的气体的压强、体积和温度之间，有如下定律：

理想气体定律

其中p为压强，V为体积，n为摩尔数，R为普适气体常数，T为绝对温度；

非理想气体的vanderWaals方程

其中称为摩尔体积，a,b是和具体气体有关的经验常数。

在某化工设计中，需要精确估计二氧化碳和氧气在不同温度和压强下的摩尔体积，以选择适当的容器，其中各参数分别为

R=0.082054Latm/（molK）；

二氧化碳：

a=3.592，b=0.04267；

氧气：

a=1.360，b=0.03183；

设计压强分别为1，10，50，100atm；

温度分别为400，600，800K。

（1）用vanderWaals方程计算在不同条件下的摩尔体积；

（2）将上述计算结果和理想气体的计算结果进行比较。

根据题意，利用牛顿迭代法求解问题

（1），对于不同的气体，分别输出结果：

clear;

disp（'

结果输出4*3矩阵，1到4行压力分别为1、10、50、100atm，1至3列温度分别为400、600、800K'

）

用vanderWaals方程计算得到的摩尔体积v（单位为mol/L）结果为：

'

form=1:

2;

ifm==1

a=1.360;

b=0.03183;

else

a=3.592;

b=0.04267;

P=[1,10,50,100];

R=0.082054;

T=[400,600,800];

x0=1;

e=1e-8;

v=zeros（4,3）;

fori=1:

4

forj=1:

3

（p+a/x^2）*（x-b）-R*t'

a'

b'

p'

R'

t'

p-a/x^2+2*a*b/x^3'

p=P（i）;

a=a;

b=b;

t=T（j）;

R=R;

x1=x0-f（x0,a,b,p,R,t）/df（x0,a,b,p）;

r=sqrt（norm（f（x1,a,b,p,R,t）））;

v（i,j）=x1;

v_n=v;

气体为氧气'

气体为二氧化碳'

disp（v_n）

输出摩尔体积v（单位为mol/L）理想值：

T=[400600800];

p=[11050100];

v_i=zeros（4,3）;

4;

3;

v_i（i,j）=R*T（j）/p（i）;

disp（v_i）;

结果输出4*3矩阵，1到4行压力分别为1、10、50、100atm，1至3列温度分别为400、600、800K

气体为氧气

32.812049.236665.6543

3.27284.92766.5756

0.64840.98991.3247

0.32200.49860.6689

气体为二氧化碳

32.754749.202165.6312

3.21464.89316.5525

0.58530.95551.3022

0.25280.46480.6473

32.821649.232465.6432

3.28224.92326.5643

0.65640.98461.3129

0.32820.49230.6564