APPLIED NUMERICAL METHODS USING MATLAB phần 9 docx

... tmpg*abs(dS)’; %Eqs. (9. 4.8), (9. 4 .9) end d(i-N_b) = tmpf*abs(dS)’; %Eq. (9. 4.10) end d=d-A12(1:N_i,1:N_b)*c(1:N_b)’; %Eq. (9. 4.10) c(N_b + 1:N_n) = A12(1:N_i,N_b+1:N_n)\d; %Eq. (9. 4.7) for s = 1:N_s for ... parameters dialog window y 2 0 0 1 2 3 x Figure 9. 14 Procedure and results of using PDEtool for Example 9. 3 /9. 8. 444 PARTIAL DIFFERENTIAL EQUATIONS 9. Noting that the solu...

Ngày tải lên: 09/08/2014, 12:22

... interpolation/extrapolation. Applied Numerical Methods Using MATLAB  , by Yang, Cao, Chung, and Morris Copyr ight  2005 John Wiley & Sons, I nc., ISBN 0-471- 698 33-4 117 126 INTERPOLATION ... -0.8552 0. 499 8 -0.6456 -0.5774 0 0.2403 0.8552 0.5184 0.8 090 0.1100 0.5774 0 0 >>err = U*S*V’-A %to check if the result is right err = 1.0e-015* -0.2220 -0.2220 00 0.4441 0 2.5...

Ngày tải lên: 09/08/2014, 12:22

... 4.6212e-2 2 .98 22e-2 8.4103e-2 0.0001 9. 4278e-3 9. 4277e-3 0.00001 2.1853e-1 2 .98 58e-3 8. 493 7e-2 (P5.12.1,2) 0.001 1.2 393 e-5 1.3545e-5 0.0001 8.3626e-3 5.0315e-6 6.4849e-6 0.00001 1.3846e -9 8.8255e-7 (P5.13.1) ... solving ODEs, it gives us a Applied Numerical Methods Using MATLAB  , by Yang, Cao, Chung, and Morris Copyr ight  2005 John Wiley & Sons, I nc., ISBN 0-...

Ngày tải lên: 09/08/2014, 12:22

... ||A i x − b i || x ||A i x − b i || inconsistent 1.5000 4.4409e-15 A 1 x = b 1 0 (1 .98 60e-15) 1.5000 0.3143 1.78 89 A 2 x = b 2 0.6286 0 .94 29 A 3 x = b 3 A 4 x = b 4 2.5000 1.2247 0.0000 A 5 x = b 5 A 6 x ... Why? %nm1p 19: Quantization Error x = 2-2^-50; for n = 1:2^3 x = x+2^-52; fprintf(’%20.18E\n’,x) end 1.20 Avoiding Large Errors/Overflow/Underflow (a) For x = 9. 8 201 and y = 10.2...

Ngày tải lên: 09/08/2014, 12:22

... 3*x2-2.5 = 0’) x1 = [ 2.] x2 = [ 0.500000] [ -1.2064 59] [ 0 .94 1336] [0.6032 29 -0. 392 630*i] [-1. 095 668 -0.540415e-1*i] [0.6032 29 +0. 392 630*i] [-1. 095 668 +0.540415e-1*i] >>S = solve(’x^3 - y^3 ... 1 49 (cf) Note that the objective of the WLS scheme is to put greater emphasis on more reliable data. %do_wlse1 for Ex.3.7 clear, clf x=[135 792 46810]; %input data y = [0.0831 0...

Ngày tải lên: 09/08/2014, 12:22

... cos(π/4) h 1 = 0.1000000000 0.67060 297 29 −0.03650380828 h 2 = 0.0100000000 0.7035 594 917 0.03 295 65188 −0.0035472 895 0 h 3 = 0.0010000000 0.7067531100 0.003 193 6183 −0.00035367121 h 4 = 0.0001000000 ... 10 −k D 2k|x=π/4 D 2k − D 2(k−1) D 2k|x=π/4 − cos(π/4) h 1 = 0.1000000000 0.70 592 88 590 −0.00117 792 2 19 h 2 = 0.0100000000 0.707 094 996 1 0.0011661371 −0.00001178505 h 3 = 0.00100...

Ngày tải lên: 09/08/2014, 12:22

... fminunc() sim anl() [0, 0] [2 .90 35, 2 .90 35] [2 .90 35, 2 .90 36] [2 .90 36, 2 .90 36] [2. 896 6, 2 .90 36] (f o =−156.66) (f o =−156.66) (f o =−156.66) (f o =−156.66) [−0.5,−1.0] [2 .90 35, −2.7468] [−2.7468, −2.7468] ... margin xo_s = 1.2118 0.5765 fo_s = 0.5322 %min value xo_steep = 1.2768 0. 598 9 fo_steep = 0.2 899 %not a minimum co_steep = -1.2768 -0. 598 9 -1.5386 0.1525 %violating...

Ngày tải lên: 09/08/2014, 12:22

... ∂ 2 f/∂x 2 i (1) [0. 696 5 −0.1423] −, − M (6) [−1. 692 6 −0.1183] (2) [2.5463 −0.1 896 ] (7) [−2.6573 −2.82 19] +, + m (3) [2.52 09 2 .90 27] +, + G (8) [−0.3227 −2.4257] (4) [−0.3865 2 .90 49] (9) [2.5216 −2. 894 6] +, ... cos(x 2 − 2 .9) whose gradient vector function is g(x) =∇f(x) =  4x 3 1 − 24x 1 − 4 + 20 sin(x 1 − 2.5) cos(x 2 − 2 .9) 4x 3 2 − 32x 2 − 5 + 20 cos(x 1 − 2.5) sin(x 2...

Ngày tải lên: 09/08/2014, 12:22

... Heine- mann, Boston, 199 9. [R-2] Recktenwald, G. W., Numerical Methods with MATLAB, Prentice-Hall, Upper Sad- dle River, NJ, 2000. [S-1] Schilling, R. J., and Harris, S. L., Applied Numerical Methods for ... Company, Reading, MA, 198 4. [M-1] Mathews, J. H., and Fink, K. D., Numerical Methods Using MATLAB, Prentice- Hall, Upper Saddle River, NJ, 199 9. [M-2] Maron, Melvin J...

Ngày tải lên: 09/08/2014, 12:22