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

35353593535359姓名：方健學號652081701073問題：使用matlab的文件計算乙醇精餾過程的物料平衡。并且會使用for循環(huán)語句對于含零元素多的矩陣進行構造。工業(yè)生產乙醇的二級精餾過程如圖所示過程的處理量為10000kg/h,組成為的水，乙醇和機物（重量餾塔的回流比為3，第一精餾塔塔頂產物含乙醇，而第二精餾塔頂含乙醇。第一精餾塔塔底含進料有機物，剩余的在第二塔頂，兩塔塔底物料中都不含有乙醇。求解每股物料的量。分析，本問題是一個物料衡算問題，由于涉及未知數(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調用函數(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總結通本次求解精餾過程中物料衡算們可以知道對于構造矩陣中含有零元