计算机化工应用超越方程求解.pptxVIP

计算机化工应用第二讲;问题1—线性方程组求解策略;方法一;Private Sub Command1_Click() Dim i, j, m, n As Integer Dim a(), z(), x(), w, aa(), s, t, k, l n = InputBox(n) ReDim a(n + 2, 2 * n), z(n + 2, 2 * n), x(n + 1), aa(n + 2, 2 * n) For i = 1 To n For j = 1 To n + 1 a(i, j) = InputBox(输入系数矩阵A( i , j )) Next j Next I For i = 1 To n If i = n Then GoTo 200 For t = i + 1 To n If Abs(a(i, i)) Abs(a(t, i)) Then For s = i To n + 1 aa(t, s) = a(i, s) a(i, s) = a(t, s) a(t, s) = aa(t, s) Next s Else End If Next t ;具体应用; 方法2- excel 规划求解;规划求解器求解线性方程组;方法3-matlab求解;问题2 非线性方程求解;方法1—二分法;实际问题求解;Dim a2, a3, a4, a5, b2, b3, b5, bv As Double Dim c2, c3, c5, tc, t, p As Double Private Sub Command1_Click() Dim a, b, x, x1, x2, y, k, y1, y2 As Double a2 = -4391473.1 a3 = 233734790 a4 = -8196792900 a5 = 113229830000 b2 = 4501.7239 b3 = -102972.05 b5 =bv = 20.101853 c2 = c3 = 5081973600 c5 = -3229376000000 tc = 304.2 温度可修改 t = InputBox(TEMPRESURE, ℃) t = 273.15 + t 压力可修改 p = InputBox(PRESURE) a = bv + 1 %二分法初始左边起点，保证解的有效性，不同方程，有不同设置 b = a ;计算结果（程序见第二章）;方法2-利用excel求解;方程设置;宏编程;方法3-matlab求解;Matlab求解摩尔体积核心代码;计算结果;问题3-非线性方程组求解;实际问题-反应平衡计算;实际问题-反应平衡计算;Matlab程序;Matlab程序;计算结果;Excel计算;Excel计算;Excel计算