ケミカルエンジニアのためのExcelを用いた化学工学計算法



Similar documents
untitled

コンピュータ概論

Excel Excel Excel = Excel III 7 ( ) 1

2X Y Y X θ 1, θ 2,... Y = f(x, θ 1, θ 2,...) θ k III 8 ( ) 1 / 39

( )

untitled

II Time-stamp: <05/09/30 17:14:06 waki> ii

コンピュータ概論

( a 3 = 3 = 3 a a > 0(a a a a < 0(a a a

x, y x 3 y xy 3 x 2 y + xy 2 x 3 + y 3 = x 3 y xy 3 x 2 y + xy 2 x 3 + y 3 = 15 xy (x y) (x + y) xy (x y) (x y) ( x 2 + xy + y 2) = 15 (x y)

関数のグラフを描こう

() n C + n C + n C + + n C n n (3) n C + n C + n C 4 + n C + n C 3 + n C 5 + (5) (6 ) n C + nc + 3 nc n nc n (7 ) n C + nc + 3 nc n nc n (

D0020.PDF

Excel Excel Excel = Excel ( ) 1

B 5 (2) VBA R / B 5 ( ) / 34

D0050.PDF

Excel ではじめる数値解析 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです.

1

II R n k +1 v 0,, v k k v 1 v 0,, v k v v 0,, v k R n 1 a 0,, a k a 0 v 0 + a k v k v 0 v k k k v 0,, v k σ k σ dimσ = k 1.3. k

i

改訂版 :基本的な文字化の原則(Basic Transcription System for Japanese: BTSJ)

2009 IA I 22, 23, 24, 25, 26, a h f(x) x x a h

My関数の作成演習問題集

PSCHG000.PS

N88 BASIC 0.3 C: My Documents 0.6: 0.3: (R) (G) : enterreturn : (F) BA- SIC.bas 0.8: (V) 0.9: 0.5:

Microsoft Word - VBA基礎(3).docx


●70974_100_AC009160_KAPヘ<3099>ーシス自動車約款(11.10).indb

超簡単、売上入力画面作成

方程式の解法

untitled

Part () () Γ Part ,

2 X Y Y X θ 1,θ 2,... Y = f (X,θ 1,θ 2,...) θ k III 8 ( ) 1 / 39

(1) (2) (1) (2) 2 3 {a n } a 2 + a 4 + a a n S n S n = n = S n

?

all.dvi

45 VBA Fortran, Pascal, C Windows OS Excel VBA Visual Basic Excel VBA VBA Visual Basic For Application Microsoft Office Office Excel VBA VBA Excel Acc

6. Euler x

OABC OA OC 4, OB, AOB BOC COA 60 OA a OB b OC c () AB AC () ABC D OD ABC OD OA + p AB + q AC p q () OABC 4 f(x) + x ( ), () y f(x) P l 4 () y f(x) l P

Lecture on

REALV5_A4…p_Ł\1_4A_OCF

untitled

「都市から地方への人材誘致・移住促進に関する調査」

<91498EE88CA D815B2E786C73>

〔 大 会 役 員 〕

橡本体資料+参考条文.PDF

5.. z = f(x, y) y y = b f x x g(x) f(x, b) g x ( ) A = lim h g(a + h) g(a) h g(x) a A = g (a) = f x (a, b)

~~濱田のジイサンとの出会い~~

x () g(x) = f(t) dt f(x), F (x) 3x () g(x) g (x) f(x), F (x) (3) h(x) = x 3x tf(t) dt.9 = {(x, y) ; x, y, x + y } f(x, y) = xy( x y). h (x) f(x), F (x


1990 IMO 1990/1/15 1:00-4:00 1 N N N 1, N 1 N 2, N 2 N 3 N 3 2 x x + 52 = 3 x x , A, B, C 3,, A B, C 2,,,, 7, A, B, C

untitled

< >


4 4 4 a b c d a b A c d A a da ad bce O E O n A n O ad bc a d n A n O 5 {a n } S n a k n a n + k S n a a n+ S n n S n n log x x {xy } x, y x + y 7 fx

1. A0 A B A0 A : A1,...,A5 B : B1,...,B

untitled

6.1 工 程 管 理 (ネットワーク 工 程 表 ) 159 VBA( 実 行 結 果 例 ) 出 力 結 果 シート 出 力 結 果 を 図 に 描 くと 下 図 のようになる. 図 6.3 ネットワーク 工 程 表 ( 出 力 結 果 より 作 図 )


Appendix A BASIC BASIC Beginner s All-purpose Symbolic Instruction Code FORTRAN COBOL C JAVA PASCAL (NEC N88-BASIC Windows BASIC (1) (2) ( ) BASIC BAS

(1) (2) (3) (4) (1) (a) (b) (c) (d) kg 9.8 N 5.0 kg 19.6 m/s kg m/s 8.0kg (2) 1 r=1.0m ABC QA =1

Microsoft PowerPoint - VBA解説1.ppt [互換モード]

1 1 [1] ( 2,625 [2] ( 2, ( ) /

04年度LS民法Ⅰ教材改訂版.PDF

Microsoft Word - 課題1_2.doc

AccessVBA−‹ŠpŁÒ-flO“Z

() x + y + y + x dy dx = 0 () dy + xy = x dx y + x y ( 5) ( s55906) 0.7. (). 5 (). ( 6) ( s6590) 0.8 m n. 0.9 n n A. ( 6) ( s6590) f A (λ) = det(a λi)

Microsoft Word - 触ってみよう、Maximaに2.doc

a n a n ( ) (1) a m a n = a m+n (2) (a m ) n = a mn (3) (ab) n = a n b n (4) a m a n = a m n ( m > n ) m n 4 ( ) 552


(2016 2Q H) [ ] R 2 2 P = (a, b), Q = (c, d) Q P QP = ( ) a c b d (a c, b d) P = (a, b) O P ( ) a p = b P = (a, b) p = ( ) a b R 2 {( ) } R 2 x = x, y

2/ 土 :30 11:20 似通った科目名がありますので注意してください. 受験許可されていない科目を解答した場合は無効 整理番号と科目コードは受験許可証とよく照合し正確に記入

Microsoft PowerPoint - Visualプログラミング

II

csv csv

431 a s a s a s d a s a 10 d s 11 f a 12 g s 13 a 14 a 15


, 1. x 2 1 = (x 1)(x + 1) x 3 1 = (x 1)(x 2 + x + 1). a 2 b 2 = (a b)(a + b) a 3 b 3 = (a b)(a 2 + ab + b 2 ) 2 2, 2.. x a b b 2. b {( 2 a } b )2 1 =

1 1 ( ) ( % mm % A B A B A 1

5.. z = f(x, y) y y = b f x x g(x) f(x, b) g x ( ) A = lim h 0 g(a + h) g(a) h g(x) a A = g (a) = f x (a, b)

R R 16 ( 3 )

( [1]) (1) ( ) 1: ( ) 2 2.1,,, X Y f X Y (a mapping, a map) X ( ) x Y f(x) X Y, f X Y f : X Y, X f Y f : X Y X Y f f 1 : X 1 Y 1 f 2 : X 2 Y 2 2 (X 1

a, b a bc c b a a b a a a a p > p p p 2, 3, 5, 7,, 3, 7, 9, 23, 29, 3, a > p a p [ ] a bp, b p p cq, c, q, < q < p a bp bcq q a <

270万回再生レポート

4 (induction) (mathematical induction) P P(0) P(x) P(x+1) n P(n) 4.1 (inductive definition) A A (basis ) ( ) A (induction step ) A A (closure ) A clos

( )/2 hara/lectures/lectures-j.html 2, {H} {T } S = {H, T } {(H, H), (H, T )} {(H, T ), (T, T )} {(H, H), (T, T )} {1

Microsoft PowerPoint - vp演習課題

untitled

sec13.dvi

/ 2 n n n n x 1,..., x n 1 n 2 n R n n ndimensional Euclidean space R n vector point R n set space R n R n x = x 1 x n y = y 1 y n distance dx,

08‘͆i”O“Z†j

2 p T, Q

1/1 lim f(x, y) (x,y) (a,b) ( ) ( ) lim limf(x, y) lim lim f(x, y) x a y b y b x a ( ) ( ) xy x lim lim lim lim x y x y x + y y x x + y x x lim x x 1

(1) + b = b +, (2) b = b, (3) + 0 =, (4) 1 =, (5) ( + b) + c = + (b + c), (6) ( b) c = (b c), (7) (b + c) = b + c, (8) ( + b)c = c + bc (9

PowerPoint プレゼンテーション

PROSTAGE[プロステージ]

DVIOUT

1 8, : 8.1 1, 2 z = ax + by + c ax by + z c = a b +1 x y z c = 0, (0, 0, c), n = ( a, b, 1). f = n i=1 a ii x 2 i + i<j 2a ij x i x j = ( x, A x), f =

2 log 3 9 log 0 0 a log log 3 9 y 3 y = = 9 y = 2 0 y = 0 a log 0 0 a = a 9 2 = 3 log 9 3 = 2 a 0 = a = a log a a = log a = 0 log a a =. l

P ZR.indd

Transcription:

VBA 7.1 f ()= ( f ( )) y = f ()(1) y = f ( )( ) + f ( ) (1) 1 1 1 = f ( ) f ( ) (2) 1 n = = y = f() y = () 1 n+1 = n (f( n )f( n ))

1.7 1 2 log()2 145

7.2 f ( ) 1-1 -1 f ( ) (3) (4) 34 5 1 1) ( ) ( ) '( = f f f ) ' 1 f = + ) ( f ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) '( ) 1 1 1 1 1 1 + = = = f f f f f f f f f ( 146

X 1 X1 X2: 1 FXf( 1 ) FX1f( ) NX1: DX: 147

Sub () '-------------------- Dim X As Double, X1 As Double, X2 As Double 'X,X1 Dim N1 As Integer ' N1 EPS =.1 ' X = Cells(3, 3) 'Ecel X X1 = Cells(4, 3) 'Ecel X1 1: FX = X ^ 3-2 * X ^ 2 - X + 2 ' FX FX1 = X1 ^ 3-2 * X1 ^ 2 - X1 + 2 ' FX1 X2 = (X * FX1 - X1 * FX) / (FX1 - FX) ' X2 DX = Abs((X2 - X1) / X1) 'X1 X2 DX N1 = N1 + 1 ' N1 '-------------------- Cells(12 + N1, 2) = N1 ' N1 Cells(12 + N1, 3) = X 'X Cells(12 + N1, 4) = X1 'X1 Cells(12 + N1, 5) = X2 ' X2 Cells(12 + N1, 6) = FX ' FX Cells(12 + N1, 7) = FX1 ' FX1 '********************************** '----------DX EPS 2 ---------- If DX < EPS Then GoTo 2 '---------- X1 X 1 ---------- Else X = X1 ' X1 X X1 = X2 ' X2 X1 GoTo 1 End If '*************************************************** 2: Cells(8, 3) = X1 ' X1 Cells(9, 3) = N1 ' N1 End Sub 7.3 148

f()= f( )f( 1 )f( )f( 1 )< 2 1 y=f()[, 1 ] [, 1 ] 1 1 [, 1 ] c1 c1 ( + 1 )/2 2f( c1 )f( )f( c1 )f( )> = c1 f( c1 ) f( )f( c1 )f( )< 1 = c1 31,2 = 1-1 /.1 f( ) 1 f( c2 ) f( c3 ) f( c1 ) c = c1 c3 c2 + 2 1 f( c )f( )>c f( c )f( )<c1 149

Sub () Dim tin As Range Dim tout As Range Set tin = Range("B1:B1") Set tout = Range("C1:C1") '------------- = tin.cells(1, 1) 1 = tin.cells(2, 1) f = ^ 3-2 * ^ 2 - + 2 f1 = 1 ^ 3-2 * 1 ^ 2-1 + 2 If f * f1 > Then MsgBo "" Eit Sub End If Do c = ( + 1) / 2 fc = c ^ 3-2 * c ^ 2 - c + 2 y = f * fc If y = Then Eit Do ElseIf y > Then = c Else 1 = c End If n = n + 1 If n = 5 Then c 1 MsgBo "5" Eit Sub End If Loop Until Abs((1 - ) / 1) <.1 tout.cells(1, 1) = "=" & c tout.cells(2, 1) = "=" & n End Sub 15

7.4 i y i y i f ( i ) f ( i ) e i = y i f ( i ) y=a 1 +a S y e i S y=a 1 +a a a 1 y=a 1 +a S=e i ={y i -(a 1 i +a )} S (S a )= (S a 1 )= a,a S 1 5 2 -yf () = a + a 1 a a 1 X 1 2 3 4 5 6 7 Y 1.5 2.4 3.3 4. 4.8 5.8 6. 151

y=a 1 +a y e i Ecel 152

-y -y X X Y Y 153

X/Y/ (S) (F) 1 R 154

L E OK 155

3 Sub SAISYO1() Dim NP As Integer Dim XX(1), YY(1) As Double Dim targetin, targetout As Range NP = Cells(2, 2) For i = 1 To NP XX(i) = Cells(3 + i, 2) YY(i) = Cells(3 + i, 3) Net i SX2 = SX = SSX = SY = SXSY = For i = 1 To NP SX2 = SX2 + XX(i) ^ 2 SX = SX + XX(i) 156

SY = SY + YY(i) SXSY = SXSY + XX(i) * YY(i) Net i SSX = SX ^ 2 A = (SX2 * SY - SX * SXSY) / (NP * SX2 - SSX) A1 = (-SX * SY + SXSY * NP) / (NP * SX2 - SSX) Cells(2, 5) = "Y=A+A1*X" Cells(3, 5) = "A=" Cells(4, 5) = "A1=" Cells(3, 6) = A Cells(4, 6) = A1 End Sub 157

7.5 158

1 12wt.% 25wt.% 1g F[g]V[g]L[g] w 1 w 2 y 1 y 2 1 2 F=V+L (1) w 1 Fy 1 V 1 L (2) w 2 Fy 2 V 2 L (3) y 1 w 1.12w 2 1.12.88y 1 y 2 1 1.25 2 1.25.75 F1 (1)(2) V+L1 (1).121V+.25L.25LV12 (2) (1) (2) 159

1 1 V 1 =.25 L 12 (4) 1 1 1 (5).25 12 (5) V 52 [g], L48 [g] 2 98.5wt. SO 3 (SO 3 22wt.%78wt.%) 1g/sSO 3 97wt.% H 2 SO 4 98.5wt.%H 2 SO 4 H 2 SO 4 95wt.% (98.5wt.% H 2 SO 4 ) H 2 SO 4 98SO 3 8S32 (a) 95wt.% H 2 SO 4 (b), y, z, w [g/s] i [g/s] A 16

A 1 + = i + y (1) (.22)( 1) 32 32 32 +.95 =.985y (2) 8 98 98 (.78)(1)=i (3) B + w = z (4).95+.985.97z (5) (1)(5) 1 1 1 1 A B y 8.8 1 1 1 z = 1 i 78.95.97.985 w AB AB A.312 B.3216 1.312 1.95 1.3216 1.97 1 1 1 y 8.8 1 z = i 78.985 w = 151, y= 173 z= 353., i= 78 w= 21 a) 95wtH 2 SO 4 151g/sb 98.5wt H SO 4 173g/s 161

7.6 y = f() [a, b]n, 1, 2, n y, y 1, y 2, y n P, P 1, P 2, P n n S = S 1 + S 2 + S 3 + S n (1) h=1/n S 1 =(h/2)(y +y 1 ), S 2 =(h/2)(y 1 +y 2 ),, S n =(h/2)(y n1 +y n ) (1) S=(h/2)(y +y 1 )+(h/2)(y 1 +y 2 )+(h/2)(y 2 +y 3 )++(h/2)(y n 1+y n ) =(h/2)y + y n + 2(y 1 + y 2 + y n 1) y y = f() P P 1 P 2 P -1 P Pn S 1 S 2 S Sn y y 1 y 2 y 1 y y n a 1 2 n = y 1 b = S 1 1 h y 162

Sub fore1() Dim a, b, h, V As Single a = Cells(6, 3) b = Cells(6, 4) h = Cells(6, 6) M = (b- a) / h SS = N = = For N = 1 To M - 1 X = a + N * h SS = SS + 3 * X ^ 2 Net N S =.5 * h * (3 * a ^ 2 + 3 * b ^ 2 + 2 * SS) Cells(12, 3) = S End Sub y=3 2 =~2 163

. 8.1 adsorption 4mg/ 4m 3 6g mg/g/g nw (1) + V[] C [mg/] : W[g] (C -C)V (2) nw(c -C)V (3) C C V n = V = ( C C ) (4) W W V (3) 1 W 1g : n[mol/g] V n g / g) = ( C C W 2 C 4 3 ( 4 ) = ( C 4) 6 164

-2/3 Langmuir C n KC n = 1+ KC 1 165

1 2 3 3g 4mg/ 4m 3 C[mg/] [g/g] 3g 3 2 2 2 4 5 15g 4mg/ 4m 3 C[mg/] [g/g] 15g 4 5 4 166

8.2 (a) Henry nkc (1) K 1 (a) C n C p (a) (b) (c) (d)bet 1 167

(b)langmuir 2 1 [-]r [mols -1 ]a[mols -1 ] 2 Langmuir ra (2) r[mols -1 ](1-)[-]C[mol m -3 ] r b(1-)c (3) b[m 3 s -1 ] ab(1-)c (4) b C θ = (5) a + b C 168

n n [molg -1 ] b/a K[m 3 mol -1 ](K)θ = n n KC n = 1+ KC (6) (c)(freundlich) 1(c) nc 1/ (7) n[molg -1 ][-] (d)bet(brunauer-emmett-teller) BET 1(d) 1 BET 3 q = q K c ( 1 c) ( 1 + Kc c) 169

3 8.3 Langmuir 12 1 Kn n K 1 n, C, C/n C [molm -3 ] n [molg -1 ] C/n [gm -3 ] 46 225 19 52.9 21.4 2.68 1.99 1.49 1.5.6 172 128 73.2 5.4 35.7 (6) 1 K C = 1+ n n KC (9) C C n (1 + KC ) C n KC 1 1 = n K n = C = C (1) 17

=1n =1(n K)(1)C/nC =1n =1n K 1 C/nC 1 4 =1/n =.36 =1/(n K)=35.7 n K n =3.27 molg -1 K =.86 m 3 mol -1 n KCn(Cn) 5 = 1 / n 1 n K 4 171

n [ molg -1 ] Langmuir n : C [ molm -3 ] 5 12 C n 8.4 6 7 V[m 3 ] C [molm -3 ]W[g] n[molg -1 ]WnW (C -C b )VC b [molm -3 ] nw(c -C b )V (11) n=-v/w(c b -C ) (12) 1 172

(12)-V/W 2 3 4 C C 1 C 2 C 3 6 n [ molg -1 ] 3. 2.5 2. 1.5 1..5 q 1 V W c 1 c 1 2 3 4 5 C [ molm -3 ] 7 173

-V/W-V/W=-1/1=-.1 m 3 g -1 2 mol m -3 (4-2)1=38 mol 3. 2.5 2. 1.5 1..5 q 1 V? - W c 1 c 1 2 3 4 5 174