工程數(shù)學(xué)作業(yè)第五次方健 (10.29)

35353593535359姓名：方健學(xué)號652081701073問題：使用matlab的文件計算乙醇精餾過程的物料平衡。并且會使用for循環(huán)語句對于含零元素多的矩陣進行構(gòu)造。工業(yè)生產(chǎn)乙醇的二級精餾過程如圖所示過程的處理量為10000kg/h,組成為的水，乙醇和機物（重量餾塔的回流比為3，第一精餾塔塔頂產(chǎn)物含乙醇，而第二精餾塔頂含乙醇。第一精餾塔塔底含進料有機物，剩余的在第二塔頂，兩塔塔底物料中都不含有乙醇。求解每股物料的量。分析，本問題是一個物料衡算問題，由于涉及未知數(shù)較多，故是線性方程組的求解，對水用W示，A表示乙醇，而有機物料用R表示，包括再沸器和冷凝器的第一蒸餾塔有如下物料衡算方程。蒸餾塔：262RR26

塔底物料不含有乙醇回流比為3W5A精餾塔的回流比為5R5434AA343蒸餾塔：712A7RR7WW11A1113W11

回流比為39

11三個獨立的線性約束

A77AR0.08R2故求解Ax=BX=A\B

A

W

A

W

W

A

W

A

W

A

W

W

A

故編寫m文如下=%Ai=1:19A，B

在的commandwindow調(diào)用函數(shù)solve_distillationsolve_distillationAColumnsthrough101.000000-1.333300000001.000000-1.333300000

0

1.0000

0

0-1.3333

0

0000

0

0

1.0000

0

0

00-3.0000

00

0

0

0

1.0000

0

000

-3.00000

0

0

0

0

1.0000

000

00

0

0

0

0

0

1.000001.0000

00

0

0

0

0

0

01.0000

0

00

0

0

0

0

0

001.0000

-0.66670

0

0

0

0

0

000

1.00000

0

0

0

0

0

01.0000

0

00

0

0

0

0

0

000

00

0

0

0

0

0

000

00

0

0

0

0

0

000

00

0

0

0

0

0

000

00

0

0

0

0

0

00-1.0000

00

0

0

0

0

0

000

00

0

0

0

0

0

000

00

0

0

0

0

0

000

-1.0000

Columns11through19000000000000000000000

0

0

0

000-3.0000

0

0

0

0000

0

0

0

0000

0

0

0

0001.0000

0

00

0

0000

0

0

0

0001.0000

0

00

0

00001.0000

0-1.3333

00

0

0000

1.00000

-1.33330

0

0000

01.0000

00

0

-3.0000000

0

01.0000

0

00-3.000000

0

0

1.0000

01.00000-1.00000

00

00

1.0000

000

0

0

0

01.0000-0.05260

0

0

0

001.0000B0000

00800080001000100000000000x=1.0e+03*1.86674.00000.80001.40003.00000.60007.53330.80000.46671.00000.20000.21054.00000.15793.00000.41400.20000.05261.0000故求解如下

A

W

A

W

W

A

W

A

W

A

W

W

A

=1.0e+03*1.86674.00000.8000

1.40003.00000.60007.53330.80000.46671.00000.20000.21054.00000.15793.00000.41400.20000.05261.0000（單位為kg/h總結(jié)通本次求解精餾過程中物料衡算們可以知道對于構(gòu)造矩陣中含有零元