|
|
- あやか うづき
- 4 years ago
- Views:
Transcription
1 LU LU LU SOR ( ) SOR n x 1,, x n a 11 x 1 + a 1n x n = b 1 a n1 x 1 + a nn x n = b n
2 4 2 x 1 x 1 = b 1 a 12 x 2 a 1n x n a 11 n 1 x 2,, x n x n x n x n 1 x n x n n n Ax = b a 11 a 12 a 1n b 1 a 21 a 22 a 2n b 2 a n1 a n2 a nn b n 2 1 a 21 /a 11 a 2i = a 2i a 1i a 21 a 11, b 2 = b 2 b 1 a 21 a 11 i = 1,, n a 11 a 12 a 1n b 1 0 a 22 a 2n b a 32/a 22 0 a n2 a nn b n a i3 = a 3i a 2i b 3 = b 3 b 2 a 32 a 22 a 32 a 22, i = 2,, n a 11 a 12 a 1n b 1 0 a 22 a 2n b a 3n b a nn b n
3 4 3 A 0 a 11 a 12 a 1n b 1 0 a 22 a 2n b a (n 1) nn b (n 1) n 412 x n = b (n 1) n /a (n 1) n,n x n 1 = (b (n 2) n 1 a (n 2) n 1,n x n)/a (n 2) n 1,n 1 i n a (n 1) in /a (n 1) nn n 1 i ā b1 0 ā 22 0 b ā nn bn x i = b i /ā ii, i = 1, n 413 a (k 1) kk 0 0 a (k 1) lk 0 l 0 0 a (k 1) lk, l = k + 1,, n 0 A 414
4 4 4 1 // gauss elimination partial pivotting 2 #include "pchh" 3 #include <iostream> 4 #include <mathh> 5 #define N 4 6 int 7 main(int ac, char av[]) 8 { 9 double a[n][n+1] = {{10, 10, 30, 40, 10}, 10 {20, 10, 50, 10, 50}, 11 {30, 60, 20, 10, 80}, 12 {20, 20, 20, 30, 20}}; 13 int i, j, k, l, pivot; 14 double x[n]; 15 double p, q, t; for(i=0; i < N; i++) { 18 t = 0; 19 pivot = i; 20 for(k = i; k < N; k++) { 21 if(fabs(a[k][i]) > t) { 22 t = fabs(a[k][i]); 23 pivot = k; 24 } 25 } 26 if(pivot!= i) { 27 for(j=i; j < N+1; j++) { 28 t = a[i][j]; 29 a[i][j] = a[pivot][j]; 30 a[pivot][j] = t; 31 } 32 } 33 p = a[i][i]; 34 for(j = i+1; j < N; j++) { 35 q = a[j][i]/p; 36 for(k = i+1; k < N+1; k++) 37 a[j][k] = a[j][k] q a[i][k]; 38 a[j][i] = 0; 39 } 40 } 41 for(i = N 1; i >=0; i ) { 42 x[i] = a[i][n]; 43 for(j = N 1; j>i; j ) { 44 x[i] = x[i] a[i][j] x[j]; 45 } 46 x[i] = x[i]/a[i][i]; 47 } 48 printf("array a is\n"); 49 for(i=0; i < N; i++) { 50 for(j=0; j < N+1; j++) { 51 printf("%105f, ", a[i][j]); 52 } 53 printf("\n"); 54 } 55 printf("solution x is \n"); 56 for(i=0; i < N; i++) { 57 printf("x%d = %2517e\n", i+1, x[i]); 58 } 59 return 0; 60 }
5 x + y 3z 4w = 1 2x + y + 5z + w = 5 2x + y + (5 + ϵ)z + w = 8 2x + 2y + 2z 3w = 2 x = 2 y = ϵ 51 5ϵ ϵ = z = 3 ϵ w = 4ϵ 24 5ϵ x = 20, y = , z = 30000, w = ϵ = x = 20, y = , z = , w = ϵ = array a is , , , , , , , , , , , , , , , , , , , , solution x is x1 = e+00 x2 = e+05 x3 = e+04 x4 = e+04 ϵ = array a is , , , , , , , , , , , , , , , , , , , ,
6 4 6 solution x is x1 = e+00 x2 = e+05 x3 = e+04 x4 = e+04 ϵ = array a is , , , , , , , , , , , , , , , , , , , , solution x is x1 = e+00 x2 = e+07 x3 = e+06 x4 = e+06 ϵ = array a is , , , , , , , , , , , , , , , , , , , , solution x is x1 = e+00 x2 = e+07 x3 = e+06 x4 = e+06 ϵ = array a is , , , , , , , , , , , , , , , , , , , , index is 2, 3, 0, 1, solution x is x1 = e+00 x2 = e+07 x3 = e+06 x4 = e+06
7 4 7 ϵ i i i 1 1 p = a[i][i]; 2 for(j = i+1; j < N+1; j++) 3 a[i][j] = a[i][j]/p; // a {i,i} 1 4 a[i][i] = 1; 5 for(j = i+1; j < N; j++) { 6 q = a[j][i]; 7 for(k = i+1; k < N+1; k++) 8 a[j][k] = a[j][k] q a[i][k]; 9 a[j][i] = 0; 10 } 1 p = a[i][i]; 2 for(j = i+1; j < N+1; j++) 3 a[i][j] = a[i][j]/p; // a {i,i} 1 4 a[i][i] = 1; 5 for(j = 0; j < N; j++) { // i+1n 1 6 if (j == i) continue; // i 7 q = a[j][i]; 8 for(k = i+1; k < N+1; k++) 9 a[j][k] = a[j][k] q a[i][k]; 10 a[j][i] = 0; 11 } A a[i][n] i i + 1
8 LU Ax = b j i a j,i /a i,i A 1 a ji /a ii 1 i i i M i = a i+1,i /a i,i 1 a n,i /a i,i 1 i j A = a 11 a 12 a 1n a 21 a 22 a 2n a n1 a n2 a nn 2 a 11 a 12 a 1n 0 a A 22 a 2n = 0 a n2 a nn M 1 = m m n1 0 1 A A = M 1 A m i1 = a i1 /a 11 A (k) = M k M 2 M 1 A M k = m k+1,k 1 0 m nk 0 1
9 4 9 m jk = a (k 1) jk /a (k 1) kk n 1 A (n 1) = â 11 â 12 â 1n 0 â 22 â 2n 0 0 â 3n 0 0 â nn = U U L = M 1 n 1 M 1 2 M 1 1 A = M 1 n 1 M 1 n 2 M 1 1 U = LU L M k M k = m k+1,k 1 0 m n,k 0 1 L = 1 m k,1 1 m k+1,1 m k+1,k 1 m n,1 m n,k m n,k+1 1 U D 1 U â m â 22 0 A = LDU = m n1 m n â nn 1 â 12 â 1n 0 1 â 2n LDU â ij = â ij/â ii A (a ij = a ji ) L T = U A = LDL T A LU Ax = LUx = b Ux = y
10 4 10 Ly = b Ux = y L U y 1 = b 1 y 2 = b 2 m 2,1 y 1 y n = b n (m n,1 y m n,n 1 y n 1 ) x n = y n /â n,n x n 1 = (y n 1 â n 1,n y n )/â n 1,n 1 x 1 = (y 1 â 1,n y n â 1,n 1 y n 1 â 1,2 y 2 )/â 1,1 431 LU 1 // lu decomposition 2 #include "pchh" 3 #include <iostream> 4 #include <mathh> 5 #define N 4 6 int 7 main(int ac, char av[]) 8 { 9 double a[n][n] = { 10 {10, 10, 30, 40}, 11 {20, 10, 50, 10}, 12 {30, 60, 20, 10}, 13 {20, 20, 20, 30} }; 14 int i, j, k, l; 15 double b[n][n]; 16 double c[n][n]; 17 double p, q, t; 18 double s; 19 printf("array A is\n"); 20 for (i = 0; i < N; i++) { 21 for (j = 0; j < N; j++) { 22 printf("%105f, ", a[i][j]); 23 } 24 printf("\n"); 25 } 26 for (i = 0; i < N; i++) { 27 for (j = 0; j < N; j++) { 28 b[i][j] = 0; 29 }
11 b[i][i] = 1; 31 } for (i = 0; i < N; i++) { 34 for (j = i + 1; j < N; j++) { 35 q = a[j][i] / a[i][i]; 36 for (k = i + 1; k < N; k++) 37 a[j][k] = a[j][k] q a[i][k]; 38 a[j][i] = 0; 39 b[j][i] = q; 40 } 41 } 42 printf("array L is\n"); 43 for (i = 0; i < N; i++) { 44 for (j = 0; j < N; j++) { 45 printf("%105f, ", b[i][j]); 46 } 47 printf("\n"); 48 } 49 printf("array U is\n"); 50 for (i = 0; i < N; i++) { 51 for (j = 0; j < N; j++) { 52 printf("%105f, ", a[i][j]); 53 } 54 printf("\n"); 55 } 56 for (i = 0; i < N; i++) { 57 for (j = 0; j < N; j++) { 58 s = 0; 59 for (k = 0; k < N; k++) { 60 s = s + b[i][k] a[k][j]; 61 } 62 c[i][j] = s; 63 } 64 } printf("array L*U is\n"); 67 for (i = 0; i < N; i++) { 68 for (j = 0; j < N; j++) { 69 printf("%105f, ", c[i][j]); 70 } 71 printf("\n"); 72 } 73 return 0; 74 } 432 L b a b L 1 U b a
12 LU i j P ij P ij = i j A i j P ij P ij A (i j P ij AP ij ) j i a j,i /a i,i A 1 i a ji /a ii 1 j i i i a i+1,i /a i,i 1 a n,i /a i,i 1 LU U L A A A A = LU U = , L = A A = = U
13 A 0 0 n 0 0 Ax (k) b 0, k {x (k), k = 1, 2, } A F E D D = a a a nn, E = A = D + E + F a a n1 a n2 0, F = 0 a 12 a 1n 0 0 a 2n = a 1 11 {b 1 (a 12 x (k) 2 + a 13 x (k) a 1n x (k) 2 = a 1 22 {b 2 (a 21 x (k) 1 + a 23 x (k) a 2n x (k) j = a 1 jj {b j (a j1 x (k) a jj 1 x (k) j 1 + a jj+1x (k) j a jnx (k) n = a 1 nn{b n (a n1 x (k) 1 + a n2 x (k) a nn 1 x (k) n 1 )} = D 1 {b (E + F)x (k) } = D 1 (E + F)x (k) + D 1 b
14 x j x j 1 = a 1 11 {b 1 (a 12 x (k) 2 + a 13 x (k) a 1n x (k) 2 = a 1 22 {b 2 (a a 23 x (k) a 2n x (k) j = a 1 jj {b j (a j a jj 1 j 1 + a jj+1 x (k) j a jnx (k) n = a 1 nn{b n (a n1 1 + a n a nn 1 n 1 )} = D 1 {b E Fx (k) } = (D + E) 1 Fx (k) + (D + E) 1 b 443 SOR ( ) j j x (k) j 1 ω SOR ξ (k+1) 1 = a 1 11 {b 1 (a 12 x (k) 2 + a 13 x (k) a 1n x (k) 1 = x (k) 1 + ω(ξ (k+1) 1 x (k) 1 ) ξ (k+1) 2 = a 1 22 {b 2 (a a 23 x (k) a 2n x (k) 2 = x (k) 2 + ω(ξ (k+1) 2 x (k) 2 ) ξ (k+1) j = a 1 jj {b j (a j a jj 1 j 1 + a jj+1 x (k) j a jnx (k) j = x (k) j + ω(ξ (k+1) j x (k) j ) ξ (k+1) n = a 1 nn{b n (a n1 1 + a n a nn 1 n 1 )} n = x (k) n + ω(ξ n (k+1) x (k) n ) ξ (k+1) = D 1 {b E Fx (k) } = x (k) + ω(ξ (k+1) x (k) ) ξ = (I + ωd 1 E) 1 {(1 ω)i ωd 1 F}x (k) + ω(d + ωe) 1 b
15 = Mx (k) + b M < 1 x (k) () M M Mx M = sup x 0 x M // Jacobi method 2 #include "pchh" 3 #include <iostream> 4 #include <mathh> 5 #define N 4 6 int 7 main(int ac, char av[]) 8 { 9 double a[n][n] = { 10 {50, 10, 30, 40}, 11 {10, 50, 10, 10}, 12 {30, 60, 80, 10}, 13 {20, 20, 20, 90} }; 14 double b[n] = { 10, 50, 80, 20 }; int i, j, k, n; 17 double x[n] = { 10, 10, 10, 10 }, y[n]; 18 double e, p; 19 e = 00; 20 for (i = 0; i < N; i++) { 21 p = 0; 22 for (j = 0; j < N; j++) 23 p = p + a[i][j] x[j]; 24 p = p b[i]; 25 e = e + fabs(p); 26 } 27 n = 0; 28 while (e > 10e 5) { 29 n++; 30 printf("%d: err = %f\n", n, e); 31 for (i = 0; i < N; i++) { 32 y[i] = b[i]; 33 for (j = 0; j < i; j++) 34 y[i] = y[i] a[i][j] x[j]; 35 for (j = i + 1; j < N; j++) 36 y[i] = y[i] a[i][j] x[j]; 37 y[i] = y[i] / a[i][i]; 38 } 39 for (i = 0; i < N; i++)
16 x[i] = y[i]; e = 00; 43 for (i = 0; i < N; i++) { 44 p = 0; 45 for (j = 0; j < N; j++) 46 p = p + a[i][j] x[j]; 47 p = p b[i]; 48 e = e + fabs(p); 49 } 50 } 51 for (i = 0; i < N; i++) { 52 printf("x%d = %2517e\n", i + 1, x[i]); 53 } 54 return 0; 55 } 446 while 1 while (e > 10e 5) { 2 n++; 3 printf("%d: err = %f\n", n, e); 4 for (i = 0; i < N; i++) { 5 x[i] = b[i]; 6 for (j = 0; j < i; j++) 7 x[i] = x[i] a[i][j] x[j]; 8 for (j = i + 1; j < N; j++) 9 x[i] = x[i] a[i][j] x[j]; 10 x[i] = x[i] / a[i][i]; 11 } 447 SOR while
17 #define omega while (e > 10e 5) { 3 n++; 4 printf("%d: err = %f\n", n, e); 5 for (i = 0; i < N; i++) { 6 xi = b[i]; 7 for (j = 0; j < i; j++) 8 xi = xi a[i][j] x[j]; 9 for (j = i + 1; j < N; j++) 10 xi = xi a[i][j] x[j]; 11 xi = xi / a[i][i]; 12 x[i] = x[i] + omega (xi x[i]); 13 } SOR SOR ω = SOR ω x + 3y + 4z = 10 2x + y + 5z = 7 6x + 5y + z = A = LU 43 3x + y + z = 0 x + 3y + z = 4 x + y + 3z = 6 SOR omega SOR omega = 1
Gauss
15 1 LU LDL T 6 : 1g00p013-5 1 6 1.1....................................... 7 1.2.................................. 8 1.3.................................. 8 2 Gauss 9 2.1.....................................
More informationQR
1 7 16 13 1 13.1 QR...................................... 2 13.1.1............................................ 2 13.1.2..................................... 3 13.1.3 QR........................................
More information[1] #include<stdio.h> main() { printf("hello, world."); return 0; } (G1) int long int float ± ±
[1] #include printf("hello, world."); (G1) int -32768 32767 long int -2147483648 2147483647 float ±3.4 10 38 ±3.4 10 38 double ±1.7 10 308 ±1.7 10 308 char [2] #include int a, b, c, d,
More informationsim98-8.dvi
8 12 12.1 12.2 @u @t = @2 u (1) @x 2 u(x; 0) = (x) u(0;t)=u(1;t)=0fort 0 1x, 1t N1x =1 x j = j1x, t n = n1t u(x j ;t n ) Uj n U n+1 j 1t 0 U n j =1t=(1x) 2 = U n j+1 0 2U n j + U n j01 (1x) 2 (2) U n+1
More information2000 7 17 iii,,.,,,,.,.,,. =,,.,,.,.,,,,,,, fortran,.,,,,.,,, iv (i),., 10,, 10 1. (ii),, fortran, fortran, C,.. (i),., 10 23.. (ii),.,, fortran Pasal,, for, if then else. C., fortran C.,.,,.,.,,.,,.,.,,,,.
More informationXMPによる並列化実装2
2 3 C Fortran Exercise 1 Exercise 2 Serial init.c init.f90 XMP xmp_init.c xmp_init.f90 Serial laplace.c laplace.f90 XMP xmp_laplace.c xmp_laplace.f90 #include int a[10]; program init integer
More informationŠéŒØ‘÷†u…x…C…W…A…fi…l…b…g…‘†[…NfiüŒå†v(fl|ŁŠ−Ù) 4. −mŠ¦fiI’—Ÿ_ 4.1 −mŠ¦ŁªŁz‡Ì„v”Z
( ) 4. 4.1 2009 1 14 ( ) ( ) 4. 2009 14.1 14 ( ) 1 / 41 1 2 3 4 5 4.1 ( ) 4. 2009 14.1 14 ( ) 2 / 41 X i (Ω)
More informationC による数値計算法入門 ( 第 2 版 ) 新装版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 新装版 1 刷発行時のものです.
C による数値計算法入門 ( 第 2 版 ) 新装版 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/009383 このサンプルページの内容は, 新装版 1 刷発行時のものです. i 2 22 2 13 ( ) 2 (1) ANSI (2) 2 (3) Web http://www.morikita.co.jp/books/mid/009383
More information#define N1 N+1 double x[n1] =.5, 1., 2.; double hokan[n1] = 1.65, 2.72, 7.39 ; double xx[]=.2,.4,.6,.8,1.2,1.4,1.6,1.8; double lagrng(double xx); main
=1= (.5, 1.65), (1., 2.72), (2., 7.39).2,.4,.6,.8, 1., 1.2, 1.4, 1.6 1 1: x.2 1.4128.4 1.5372.6 1.796533.8 2.198 1.2 3.384133 1.4 4.1832 1.6 5.1172 8 7 6 5 y 4 3 2 1.5 1 1.5 2 x 1: /* */ #include
More informationkiso2-09.key
座席指定はありません 計算機基礎実習II 2018 のウェブページか 第9回 ら 以下の課題に自力で取り組んで下さい 計算機基礎実習II 第7回の復習課題(rev07) 第9回の基本課題(base09) 第8回試験の結果 中間試験に関するコメント コンパイルできない不完全なプログラムなど プログラミングに慣れていない あるいは複雑な問題は 要件 をバラして段階的にプログラムを作成する exam08-2.c
More informationLINEAR ALGEBRA I Hiroshi SUZUKI Department of Mathematics International Christian University
LINEAR ALGEBRA I Hiroshi SUZUKI Department of Mathematics International Christian University 2002 2 2 2 2 22 2 3 3 3 3 3 4 4 5 5 6 6 7 7 8 8 9 Cramer 9 0 0 E-mail:hsuzuki@icuacjp 0 3x + y + 2z 4 x + y
More informationGauss Strassen LU LU LU LU 22 5 Gauss LU
I, 208 3 23 3 2 4 3 6 3 6 32 6 32 6 322 7 323 Gauss 8 324 Strassen 9 325 0 326 0 327 LU 0 33 4 LU 2 4 2 4 2 42 3 43 4 42 8 43 0 LU 20 44 LU 22 5 Gauss LU 23 5 23 52 Gauss 23 53 LU 24 53 24 532 25 533 LU
More information第7章 有限要素法のプログラミング
April 3, 2019 1 / 34 7.1 ( ) 2 Poisson 2 / 34 7.2 femfp.c [1] main( ) input( ) assem( ) ecm( ) f( ) solve( ) gs { solve( ) output( ) 3 / 34 7.3 fopen() #include FILE *fopen(char *fname, char
More information2 [1] 7 5 C 2.1 (kikuchi-fem-mac ) input.dat (cat input.dat type input.dat ))
2.3 2007 8, 2011 6 21, 2012 5 29 1 (1) (2) (3) 2 Ω Poisson u = f (in Ω), u = 0 (on Γ 1 ), u n = 0 (on Γ 2) f Γ 1, Γ 2 Γ := Ω n Γ Ω (1980) Fortran : kikuchi-fem-mac.tar.gz 1 ( fem, 6701) tar xzf kikuchi-fem-mac.tar.gz
More informationjoho12.ppt
n φ 1 (x),φ 2 (x),,φ n (x) (x i, f i ) Q n c 1,,c n (,f k ) q n = c i φ i (x) x Q n i=1 c 1 = 0 c 1 n ( c 1 ) = q n f k 2 Q n ( c 1 ) = q 2 2 n ( ) 2 f k q n ( ) + f k Q n c 1 = c 1 q n 2 q n( ) q n q
More informationC言語による数値計算プログラミング演習
5. 行列の固有値問題 n n 正方行列 A に対する n 個の固有値 λ i (i=1,,,n) と対応する固有ベクトル u i は次式を満たす Au = λ u i i i a11 a1 L a1 n u1i a1 a a n u i A =, ui = M O M M an 1 an L ann uni これらはまとめて, つぎのように書ける 5.1 ヤコビ法 = Λ, = [ u1 u u
More information(Basic Theory of Information Processing) 1
(Basic Theory of Information Processing) 1 10 (p.178) Java a[0] = 1; 1 a[4] = 7; i = 2; j = 8; a[i] = j; b[0][0] = 1; 2 b[2][3] = 10; b[i][j] = a[2] * 3; x = a[2]; a[2] = b[i][3] * x; 2 public class Array0
More informationKrylov (b) x k+1 := x k + α k p k (c) r k+1 := r k α k Ap k ( := b Ax k+1 ) (d) β k := r k r k 2 2 (e) : r k 2 / r 0 2 < ε R (f) p k+1 :=
127 10 Krylov Krylov (Conjugate-Gradient (CG ), Krylov ) MPIBNCpack 10.1 CG (Conjugate-Gradient CG ) A R n n a 11 a 12 a 1n a 21 a 22 a 2n A T = =... a n1 a n2 a nn n a 11 a 21 a n1 a 12 a 22 a n2 = A...
More information(2018 2Q C) [ ] 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
(2018 2Q C) [ ] 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 R y 2 a p =, c q = b d p + a + c q = b + d q p P q a p = c R c b
More information1. 1 BASIC PC BASIC BASIC BASIC Fortran WS PC (1.3) 1 + x 1 x = x = (1.1) 1 + x = (1.2) 1 + x 1 = (1.
Section Title Pages Id 1 3 7239 2 4 7239 3 10 7239 4 8 7244 5 13 7276 6 14 7338 7 8 7338 8 7 7445 9 11 7580 10 10 7590 11 8 7580 12 6 7395 13 z 11 7746 14 13 7753 15 7 7859 16 8 7942 17 8 Id URL http://km.int.oyo.co.jp/showdocumentdetailspage.jsp?documentid=
More informationII 2 3.,, A(B + C) = AB + AC, (A + B)C = AC + BC. 4. m m A, m m B,, m m B, AB = BA, A,, I. 5. m m A, m n B, AB = B, A I E, 4 4 I, J, K
II. () 7 F 7 = { 0,, 2, 3, 4, 5, 6 }., F 7 a, b F 7, a b, F 7,. (a) a, b,,. (b) 7., 4 5 = 20 = 2 7 + 6, 4 5 = 6 F 7., F 7,., 0 a F 7, ab = F 7 b F 7. (2) 7, 6 F 6 = { 0,, 2, 3, 4, 5 },,., F 6., 0 0 a F
More informationスライド 1
数値解析 平成 29 年度前期第 14 週 [7 月 10 日 ] 静岡大学工学研究科機械工学専攻ロボット 計測情報分野創造科学技術大学院情報科学専攻 三浦憲二郎 期末試験 7 月 31 日 ( 月 ) 9 10 時限 A : 佐鳴会議室 B : 佐鳴ホール 講義アウトライン [7 月 10 日 ] 関数近似と補間 最小 2 乗近似による関数近似 ( 復習 ) ラグランジュ補間 形状処理工学の基礎
More informationC 2 / 21 1 y = x 1.1 lagrange.c 1 / Laglange / 2 #include <stdio.h> 3 #include <math.h> 4 int main() 5 { 6 float x[10], y[10]; 7 float xx, pn, p; 8 in
C 1 / 21 C 2005 A * 1 2 1.1......................................... 2 1.2 *.......................................... 3 2 4 2.1.............................................. 4 2.2..............................................
More informationスライド 1
数値解析 平成 24 年度前期第 13 週 [7 月 11 日 ] 静岡大学創造科学技術大学院情報科学専攻工学部機械工学科計測情報講座 三浦憲二郎 講義アウトライン [7 月 11 日 ] 関数近似と補間 最小 2 乗近似による関数近似 ラグランジュ補間 形状処理工学の基礎 点列からの曲線の生成 T.Kanai, U.Tokyo 関数近似 p.116 複雑な関数を簡単な関数で近似する関数近似 閉区間
More informationスライド 1
数値解析 2019 年度前期第 13 週 [7 月 11 日 ] 静岡大学創造科学技術大学院情報科学専攻工学部機械工学科計測情報講座 三浦憲二郎 講義アウトライン [7 月 11 日 ] 関数近似と補間 最小 2 乗近似による関数近似 ラグランジュ補間 T.Kanai, U.Tokyo 関数近似 p.116 複雑な関数を簡単な関数で近似する 関数近似 閉区間 [a,b] で定義された関数 f(x)
More informationX G P G (X) G BG [X, BG] S 2 2 2 S 2 2 S 2 = { (x 1, x 2, x 3 ) R 3 x 2 1 + x 2 2 + x 2 3 = 1 } R 3 S 2 S 2 v x S 2 x x v(x) T x S 2 T x S 2 S 2 x T x S 2 = { ξ R 3 x ξ } R 3 T x S 2 S 2 x x T x S 2
More informationAHPを用いた大相撲の新しい番付編成
5304050 2008/2/15 1 2008/2/15 2 42 2008/2/15 3 2008/2/15 4 195 2008/2/15 5 2008/2/15 6 i j ij >1 ij ij1/>1 i j i 1 ji 1/ j ij 2008/2/15 7 1 =2.01/=0.5 =1.51/=0.67 2008/2/15 8 1 2008/2/15 9 () u ) i i i
More information: : : : ) ) 1. d ij f i e i x i v j m a ij m f ij n x i =
1 1980 1) 1 2 3 19721960 1965 2) 1999 1 69 1980 1972: 55 1999: 179 2041999: 210 211 1999: 211 3 2003 1987 92 97 3) 1960 1965 1970 1985 1990 1995 4) 1. d ij f i e i x i v j m a ij m f ij n x i = n d ij
More information(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
(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 R y 2 a p =, c q = b d p + a + c q = b + d q p P q a p = c R c b
More information数値計算法
12.1 電気回路網に関するキルヒホッフの法則による解法 1 工学的諸問題を多元連立 1 次方程式で表現することができる. 例えば, 荷物を最短の時間と最低のコストで輸送するためにはどのようなルートで物流を行うか という問題, 工場の部品の在庫の状況からいかに最小のコストで製品をつくるか という問題, 機械要素の運動の問題, 電気回路の解析の問題など, いくつか挙げられる. つまり, 計算機で多元連立方程式を解くことができれば,
More information1 28 6 12 7 1 7.1...................................... 2 7.1.1............................... 2 7.1.2........................... 2 7.2...................................... 3 7.3...................................
More information(Jacobi Gauss-Seidel SOR ) 1. (Theory of Iteration Method) Jacobi Gauss-Seidel SOR 2. Jacobi (Jacobi s Iteration Method) Jacobi 3. Gauss-Seide
03 9 (Jacobi Gauss-Seidel SOR (Theory of Iteration Method Jacobi Gauss-Seidel SOR Jacobi (Jacobi s Iteration Method Jacobi 3 Gauss-Seidel (Gauss-Seidel Method Gauss-Seidel 4 SOR (SOR Method SOR 9 Ax =
More information行列代数2010A
(,) A (,) B C = AB a 11 a 1 a 1 b 11 b 1 b 1 c 11 c 1 c a A = 1 a a, B = b 1 b b, C = AB = c 1 c c a 1 a a b 1 b b c 1 c c i j ij a i1 a i a i b 1j b j b j c ij = a ik b kj b 1j b j AB = a i1 a i a ik
More informationMicrosoft Word - COMP-MATH-2016-FULLTEXT.doc
#13 2017/1/17 16 連立一次方程式を解くその 2 1) 対角線上に 0 が現れる場合を回避する 0x + 3y = 9 2x + 4y = 16 これを kadai-26 で解いてみなさい 以下のようにエラーになる 0.00000000 3.00000000 9.00000000 2.00000000 4.00000000 16.00000000-1.#IND0000-1.#IND0000-1.#IND0000-1.#IND0000-1.#IND0000-1.#IND0000
More informationOHP.dvi
0 7 4 0000 5.. 3. 4. 5. 0 0 00 Gauss PC 0 Gauss 3 Gauss Gauss 3 4 4 4 4 3 4 4 4 4 3 4 4 4 4 3 4 4 4 4 u [] u [3] u [4] u [4] P 0 = P 0 (),3,4 (,), (3,), (4,) 0,,,3,4 3 3 3 3 4 4 4 4 0 3 6 6 0 6 3 6 0 6
More information資料
PC PC C VMwareをインストールする Tips: VmwareFusion *.vmx vhv.enable = TRUE Tips: Windows Hyper-V -rwxr-xr-x 1 masakazu staff 8552 7 29 13:18 a.out* -rw------- 1 masakazu staff 8552 7 29
More information+ + + + n S (n) = + + + + n S (n) S (n) S 0 (n) S (n) 6 4 S (n) S (n) 7 S (n) S 4 (n) 8 6 S k (n) 0 7 (k + )S k (n) 8 S 6 (n), S 7 (n), S 8 (n), S 9 (
k k + k + k + + n k 006.7. + + + + n S (n) = + + + + n S (n) S (n) S 0 (n) S (n) 6 4 S (n) S (n) 7 S (n) S 4 (n) 8 6 S k (n) 0 7 (k + )S k (n) 8 S 6 (n), S 7 (n), S 8 (n), S 9 (n), S 0 (n) 9 S (n) S 4
More informationθ (t) ω cos θ(t) = ( : θ, θ. ( ) ( ) ( 5) l () θ (t) = ω sin θ(t). ω := g l.. () θ (t) θ (t)θ (t) + ω θ (t) sin θ(t) =. [ ] d dt θ (t) ω cos θ(t
7 8, /3/, 5// http://nalab.mind.meiji.ac.jp/~mk/labo/text/furiko/ l (, simple pendulum) m g mlθ (t) = mg sin θ(t) () θ (t) + ω sin θ(t) =, ω := ( m ) ( θ ) sin θ θ θ (t) + ω θ(t) = ( ) ( ) g l θ(t) = C
More informationii 3.,. 4. F. ( ), ,,. 8.,. 1. (75% ) (25% ) =7 24, =7 25, =7 26 (. ). 1.,, ( ). 3.,...,.,.,.,.,. ( ) (1 2 )., ( ), 0., 1., 0,.
(1 C205) 4 10 (2 C206) 4 11 (2 B202) 4 12 25(2013) http://www.math.is.tohoku.ac.jp/~obata,.,,,..,,. 1. 2. 3. 4. 5. 6. 7. 8. 1., 2007 ( ).,. 2. P. G., 1995. 3. J. C., 1988. 1... 2.,,. ii 3.,. 4. F. ( ),..
More information実際の株価データを用いたオプション料の計算
2002 2 20 1 1 3 2 3 2.1 : : : : : : : : : : : : : : : : : : : : : : : : : : : : 5 2.1.1 : : : : : : : : : : : : : : : : : : : : 5 2.1.2 : : : : : : : : : : : : : : : : : : : : 6 2.2 : : : : : : : : : :
More information第9回 配列(array)型の変数
第 12 回 配列型の変数 情報処理演習 ( テキスト : 第 4 章, 第 8 章 ) 今日の内容 1. 配列の必要性 2. 配列の宣言 3. 配列変数のイメージ 4. 配列変数を使用した例 5. 範囲を超えた添字を使うと? 6. 多次元配列変数 7. 多次元配列変数を使用した例 8. データのソーティング 9. 今日の練習問題 多数のデータ処理 1. 配列の必要性 ( テキスト 31 ページ )
More information新版明解C言語 実践編
2 List - "max.h" a, b max List - max "max.h" #define max(a, b) ((a) > (b)? (a) : (b)) max List -2 List -2 max #include "max.h" int x, y; printf("x"); printf("y"); scanf("%d", &x); scanf("%d", &y); printf("max(x,
More informationmain.dvi
y () 5 C Fortran () Fortran 32bit 64bit 2 0 1 2 1 1bit bit 3 0 0 2 1 3 0 1 2 1 bit bit byte 8bit 1byte 3 0 10010011 2 1 3 0 01001011 2 1 byte Fortran A A 8byte double presicion y ( REAL*8) A 64bit 4byte
More informationA11 (1993,1994) 29 A12 (1994) 29 A13 Trefethen and Bau Numerical Linear Algebra (1997) 29 A14 (1999) 30 A15 (2003) 30 A16 (2004) 30 A17 (2007) 30 A18
2013 8 29y, 2016 10 29 1 2 2 Jordan 3 21 3 3 Jordan (1) 3 31 Jordan 4 32 Jordan 4 33 Jordan 6 34 Jordan 8 35 9 4 Jordan (2) 10 41 x 11 42 x 12 43 16 44 19 441 19 442 20 443 25 45 25 5 Jordan 26 A 26 A1
More information数学Ⅱ演習(足助・09夏)
II I 9/4/4 9/4/2 z C z z z z, z 2 z, w C zw z w 3 z, w C z + w z + w 4 t R t C t t t t t z z z 2 z C re z z + z z z, im z 2 2 3 z C e z + z + 2 z2 + 3! z3 + z!, I 4 x R e x cos x + sin x 2 z, w C e z+w
More information数値計算
数値計算 垣谷公徳 17 号館 3 階電子メール : kimi@ee.ous.ac.jp プログラミング言語の一般論 データ型 ( 定数と変数 配列 ) 代入 基本演算 ( 四則演算 ) 入出力 分岐 繰返処理 関数 外部手続き 1 2 入力関数 入出力 getchar, getc, fgetc ; 一文字入力 gets, fgets, fread ; 文字列 ( データ列 ) 入力 scanf,
More information離散最適化基礎論 第 11回 組合せ最適化と半正定値計画法
11 okamotoy@uec.ac.jp 2019 1 25 2019 1 25 10:59 ( ) (11) 2019 1 25 1 / 38 1 (10/5) 2 (1) (10/12) 3 (2) (10/19) 4 (3) (10/26) (11/2) 5 (1) (11/9) 6 (11/16) 7 (11/23) (11/30) (12/7) ( ) (11) 2019 1 25 2
More information1 4 2 EP) (EP) (EP)
2003 2004 2 27 1 1 4 2 EP) 5 3 6 3.1.............................. 6 3.2.............................. 6 3.3 (EP)............... 7 4 8 4.1 (EP).................... 8 4.1.1.................... 18 5 (EP)
More information1 1.1 C 2 1 double a[ ][ ]; 1 3x x3 ( ) malloc() malloc 2 #include <stdio.h> #include
1 1.1 C 2 1 double a[ ][ ]; 1 3x3 0 1 3x3 ( ) 0.240 0.143 0.339 0.191 0.341 0.477 0.412 0.003 0.921 1.2 malloc() malloc 2 #include #include #include enum LENGTH = 10 ; int
More information第1章
...1... 1... 1...3... 3... 5... 6...9... 9... 9... 9... 10... 12... 14 Windows...16... 16... 16...21... 21... 25...30...31...32...33 3 1 2 3 4 5 6 7 8 n x, y i 0,1,, n 1 i i y P x ) i 0, 1,, n 1 n 1 i
More information£Ã¥×¥í¥°¥é¥ß¥ó¥°ÆþÌç (2018) - Â裶²ó ¨¡ À©¸æ¹½Â¤¡§·«¤êÊÖ¤· ¨¡
(2018) 2018 5 24 ( ) while ( ) do while ( ); for ( ; ; ) while int i = 0; while (i < 100) { printf("i = %3d\n", i); i++; while int i = 0; i while (i < 100) { printf("i = %3d\n", i); i++; while int i =
More information9 8 7 (x-1.0)*(x-1.0) *(x-1.0) (a) f(a) (b) f(a) Figure 1: f(a) a =1.0 (1) a 1.0 f(1.0)
E-mail: takio-kurita@aist.go.jp 1 ( ) CPU ( ) 2 1. a f(a) =(a 1.0) 2 (1) a ( ) 1(a) f(a) a (1) a f(a) a =2(a 1.0) (2) 2 0 a f(a) a =2(a 1.0) = 0 (3) 1 9 8 7 (x-1.0)*(x-1.0) 6 4 2.0*(x-1.0) 6 2 5 4 0 3-2
More information1 1.1 C 2 1 double a[ ][ ]; 1 3x x3 ( ) malloc() 2 double *a[ ]; double 1 malloc() dou
1 1.1 C 2 1 double a[ ][ ]; 1 3x3 0 1 3x3 ( ) 0.240 0.143 0.339 0.191 0.341 0.477 0.412 0.003 0.921 1.2 malloc() 2 double *a[ ]; double 1 malloc() double 1 malloc() free() 3 #include #include
More information(1) (2) (1) (2) 2 3 {a n } a 2 + a 4 + a a n S n S n = n = S n
. 99 () 0 0 0 () 0 00 0 350 300 () 5 0 () 3 {a n } a + a 4 + a 6 + + a 40 30 53 47 77 95 30 83 4 n S n S n = n = S n 303 9 k d 9 45 k =, d = 99 a d n a n d n a n = a + (n )d a n a n S n S n = n(a + a n
More informationnumb.dvi
11 Poisson kanenko@mbkniftycom alexeikanenko@docomonejp http://wwwkanenkocom/ , u = f, ( u = u+f u t, u = f t ) 1 D R 2 L 2 (D) := {f(x,y) f(x,y) 2 dxdy < )} D D f,g L 2 (D) (f,g) := f(x,y)g(x,y)dxdy (L
More information£Ã¥×¥í¥°¥é¥ß¥ó¥°(2018) - Âè11²ó – ½ÉÂꣲ¤Î²òÀ⡤±é½¬£² –
(2018) 11 2018 12 13 2 g v dv x dt = bv x, dv y dt = g bv y (1) b v 0 θ x(t) = v 0 cos θ ( 1 e bt) (2) b y(t) = 1 ( v 0 sin θ + g ) ( 1 e bt) g b b b t (3) 11 ( ) p14 2 1 y 4 t m y > 0 y < 0 t m1 h = 0001
More informationj x j j j + 1 l j l j = x j+1 x j, n x n x 1 = n 1 l j j=1 H j j + 1 l j l j E
8 9 7 6 4 2 3 5 1 j x j j j + 1 l j l j = x j+1 x j, n x n x 1 = n 1 l j j=1 H j j + 1 l j l j E a n 1 H = ae l j, j=1 l j = x j+1 x j, x n x 1 = n 1 j=1 l j, l j = ±l l > 0) n 1 H = ϵ l j, j=1 ϵ e x x
More information£Ã¥×¥í¥°¥é¥ß¥ó¥°ÆþÌç (2018) - Â裱£²²ó ¡Ý½ÉÂꣲ¤Î²òÀ⡤±é½¬£²¡Ý
(2018) 2018 7 5 f(x) [ 1, 1] 3 3 1 3 f(x) dx c i f(x i ) 1 0 i=1 = 5 ) ( ) 3 ( 9 f + 8 5 9 f(0) + 5 3 9 f 5 1 1 + sin(x) θ ( 1 θ dx = tan 1 + sin x 2 π ) + 1 4 1 3 [a, b] f a, b double G3(double (*f)(),
More information12.2 電気回路網に関するキルヒホッフの法則による解法 2 多元連立 1 次方程式の工学的応用についての例を 2 つ示す.1 つはブリッジ T 型回路, もう 1 つはホーイストンブリッジ回路である. 示された回路図と与えられた回路定数からキルヒホッフの法則を使って多元連立 1 次方程式を導出する
12.2 電気回路網に関するキルヒホッフの法則による解法 2 多元連立 1 次方程式の工学的応用についての例を 2 つ示す.1 つはブリッジ T 型回路, もう 1 つはホーイストンブリッジ回路である. 示された回路図と与えられた回路定数からキルヒホッフの法則を使って多元連立 1 次方程式を導出する. その式のパラメータを前述のプログラムに反映させてシミュレーションを行う. 12.2.0 キルヒホッフの法則
More information1 4 1 ( ) ( ) ( ) ( ) () 1 4 2
7 1995, 2017 7 21 1 2 2 3 3 4 4 6 (1).................................... 6 (2)..................................... 6 (3) t................. 9 5 11 (1)......................................... 11 (2)
More informationcomment.dvi
( ) (sample1.c) (sample1.c) 2 2 Nearest Neighbor 1 (2D-class1.dat) 2 (2D-class2.dat) (2D-test.dat) 3 Nearest Neighbor Nearest Neighbor ( 1) 2 1: NN 1 (sample1.c) /* -----------------------------------------------------------------
More informationc-all.dvi
III(994) (994) from PSL (9947) & (9922) c (99,992,994,996) () () 2 3 4 (2) 2 Euler 22 23 Euler 24 (3) 3 32 33 34 35 Poisson (4) 4 (5) 5 52 ( ) 2 Turbo 2 d 2 y=dx 2 = y y = a sin x + b cos x x = y = Fortran
More information(I) GotoBALS, http://www-is.amp.i.kyoto-u.ac.jp/ kkimur/charpoly.html 2
sdmp Maple - (Ver.2) ( ) September 27, 2011 1 (I) GotoBALS, http://www-is.amp.i.kyoto-u.ac.jp/ kkimur/charpoly.html 2 (II) Nehalem CPU GotoBLAS Intel CPU Nehalem CPU, GotoBLAS, Hyper-Thread technology
More informationp _
Index ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap. ap.
More information1 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 =
1 8, : 8.1 1, z = ax + by + c ax by + z c = a b +1 x y z c = 0, (0, 0, c), n = ( a, b, 1). f = a ii x i + i
More information第10章 アイソパラメトリック要素
June 5, 2019 1 / 26 10.1 ( ) 2 / 26 10.2 8 2 3 4 3 4 6 10.1 4 2 3 4 3 (a) 4 (b) 2 3 (c) 2 4 10.1: 3 / 26 8.3 3 5.1 4 10.4 Gauss 10.1 Ω i 2 3 4 Ξ 3 4 6 Ξ ( ) Ξ 5.1 Gauss ˆx : Ξ Ω i ˆx h u 4 / 26 10.2.1
More information08 p Boltzmann I P ( ) principle of equal probability P ( ) g ( )g ( 0 ) (4 89) (4 88) eq II 0 g ( 0 ) 0 eq Taylor eq (4 90) g P ( ) g ( ) g ( 0
08 p. 8 4 k B log g() S() k B : Boltzmann T T S k B g g heat bath, thermal reservoir... 4. I II II System I System II II I I 0 + 0 const. (4 85) g( 0 ) g ( )g ( ) g ( )g ( 0 ) (4 86) g ( )g ( 0 ) 0 (4
More informationx(t) + t f(t, x) = x(t) + x (t) t x t Tayler x(t + t) = x(t) + x (t) t + 1 2! x (t) t ! x (t) t 3 + (15) Eular x t Teyler 1 Eular 2 Runge-Kutta
6 Runge-KuttaEular Runge-Kutta Runge-Kutta A( ) f(t, x) dx dt = lim x(t + t) x(t) t 0 t = f(t, x) (14) t x x(t) t + dt x x(t + dt) Euler 7 t 1 f(t, x(t)) x(t) + f(t + dt, x(t + dt))dt t + dt x(t + dt)
More information£Ã¥×¥í¥°¥é¥ß¥ó¥°ÆþÌç (2018) - Â裵²ó ¨¡ À©¸æ¹½Â¤¡§¾ò·ïʬ´ô ¨¡
(2018) 2018 5 17 0 0 if switch if if ( ) if ( 0) if ( ) if ( 0) if ( ) (0) if ( 0) if ( ) (0) ( ) ; if else if ( ) 1 else 2 if else ( 0) 1 if ( ) 1 else 2 if else ( 0) 1 if ( ) 1 else 2 (0) 2 if else
More informationPowerPoint Presentation
p.130 p.198 p.208 2 double weight[num]; double min, max; min = max = weight[0]; for( i= 1; i i < NUM; i++ ) ) if if ( weight[i] > max ) max = weight[i]: if if ( weight[i] < min ) min = weight[i]: weight
More information‚æ4›ñ
( ) ( ) ( ) A B C D E F G H I J K L M N O P Q R S T U V W X Y Z a b c d e f g h i j k l m n o p q r s t u v w x y z 0 1 2 3 4 5 6 7 8 9 (OUS) 9 26 1 / 28 ( ) ( ) ( ) A B C D Z a b c d z 0 1 2 9 (OUS) 9
More information新・明解C言語 ポインタ完全攻略
2 1-1 1-1 /* 1-1 */ 1 int n = 100; int *p = &n; printf(" n %d\n", n); /* n int */ printf("*&n %d\n", *&n); /* *&n int */ printf(" p %p\n", p); /* p int * */ printf("&*p %p\n", &*p); /* &*p int * */ printf("sizeof(n)
More informationA
A 2563 15 4 21 1 3 1.1................................................ 3 1.2............................................. 3 2 3 2.1......................................... 3 2.2............................................
More informationInformatics 2014
C 計算機の歴史 手回し計算機 新旧のソロバン バベッジの階差機関 スパコン ENIAC (1946) パソコン 大型汎用計算機 電卓 現在のコンピュータ Input Output Device Central Processing Unit I/O CPU Memory OS (Operating System) OS Windows 78, Vista, XP Windows Mac OS X
More information( ) ( ) 1729 (, 2016:17) = = (1) 1 1
1729 1 2016 10 28 1 1729 1111 1111 1729 (1887 1920) (1877 1947) 1729 (, 2016:17) 12 3 1728 9 3 729 1729 = 12 3 + 1 3 = 10 3 + 9 3 (1) 1 1 2 1729 1729 19 13 7 = 1729 = 12 3 + 1 3 = 10 3 + 9 3 13 7 = 91
More information2014 3 10 5 1 5 1.1..................................... 5 2 6 2.1.................................... 6 2.2 Z........................................ 6 2.3.................................. 6 2.3.1..................
More informationuntitled
1 1 Ax = b A R m m A b R m x R m A shift-and invert Lanczos - LU CG A = LU LU Ly = b Ux = y A LU A A = LL T 1 LU b,, Vol. 11, No. 4, pp. 14 18 (2006). x * x (0), x (1), x (2), A Ap A # x (n+1) = Cx (n)
More informationii
ii iii 1 1 1.1..................................... 1 1.2................................... 3 1.3........................... 4 2 9 2.1.................................. 9 2.2...............................
More informationPart () () Γ Part ,
Contents a 6 6 6 6 6 6 6 7 7. 8.. 8.. 8.3. 8 Part. 9. 9.. 9.. 3. 3.. 3.. 3 4. 5 4.. 5 4.. 9 4.3. 3 Part. 6 5. () 6 5.. () 7 5.. 9 5.3. Γ 3 6. 3 6.. 3 6.. 3 6.3. 33 Part 3. 34 7. 34 7.. 34 7.. 34 8. 35
More information2008 ( 13 ) C LAPACK 2008 ( 13 )C LAPACK p. 1
2008 ( 13 ) C LAPACK LAPACK p. 1 Q & A Euler http://phase.hpcc.jp/phase/mppack/long.pdf KNOPPIX MT (Mersenne Twister) SFMT., ( ) ( ) ( ) ( ). LAPACK p. 2 C C, main Asir ( Asir ) ( ) (,,...), LAPACK p.
More informationFAX780CL_chap-first.fm
FAX-780CL ABCDEFGHIα 01041115:10 :01 FAX-780CL α 1 1 2 3 1 2 f k b a FAX-780CL α n p q 09,. v m t w FAX-780CL A BC B C D E F G H I c i c s s i 0 9 V X Q ( < N > O P Z R Q: W Y M S T U V i c i k
More information115 9 MPIBNCpack 9.1 BNCpack 1CPU X = , B =
115 9 MPIBNCpack 9.1 BNCpack 1CPU 1 2 3 4 5 25 24 23 22 21 6 7 8 9 10 20 19 18 17 16 X = 11 12 13 14 15, B = 15 14 13 12 11 16 17 18 19 20 10 9 8 7 6 21 22 23 24 25 5 4 3 2 1 C = XB X dmat1 B dmat2 C dmat
More information[ 1] 1 Hello World!! 1 #include <s t d i o. h> 2 3 int main ( ) { 4 5 p r i n t f ( H e l l o World!! \ n ) ; 6 7 return 0 ; 8 } 1:
005 9 7 1 1.1 1 Hello World!! 5 p r i n t f ( H e l l o World!! \ n ) ; 7 return 0 ; 8 } 1: 1 [ ] Hello World!! from Akita National College of Technology. 1 : 5 p r i n t f ( H e l l o World!! \ n ) ;
More informationuntitled
20010916 22;1017;23;20020108;15;20; 1 N = {1, 2, } Z + = {0, 1, 2, } Z = {0, ±1, ±2, } Q = { p p Z, q N} R = { lim a q n n a n Q, n N; sup a n < } R + = {x R x 0} n = {a + b 1 a, b R} u, v 1 R 2 2 R 3
More information20 9 19 1 3 11 1 3 111 3 112 1 4 12 6 121 6 122 7 13 7 131 8 132 10 133 10 134 12 14 13 141 13 142 13 143 15 144 16 145 17 15 19 151 1 19 152 20 2 21 21 21 211 21 212 1 23 213 1 23 214 25 215 31 22 33
More information1 R n (x (k) = (x (k) 1,, x(k) n )) k 1 lim k,l x(k) x (l) = 0 (x (k) ) 1.1. (i) R n U U, r > 0, r () U (ii) R n F F F (iii) R n S S S = { R n ; r > 0
III 2018 11 7 1 2 2 3 3 6 4 8 5 10 ϵ-δ http://www.mth.ngoy-u.c.jp/ ymgmi/teching/set2018.pdf http://www.mth.ngoy-u.c.jp/ ymgmi/teching/rel2018.pdf n x = (x 1,, x n ) n R n x 0 = (0,, 0) x = (x 1 ) 2 +
More information( ) 1.1 Polychoric Correlation Polyserial Correlation Graded Response Model Partial Credit Model Tetrachoric Correlation ( ) 2 x y x y s r 1 x 2
1 (,2007) SPSSver8 1997 (2002) 1. 2. polychoric correlation coefficient (polyserial correlation coefficient) 3. (1999) M-plus R 1 ( ) 1.1 Polychoric Correlation Polyserial Correlation Graded Response Model
More informationlinearal1.dvi
19 4 30 I 1 1 11 1 12 2 13 3 131 3 132 4 133 5 134 6 14 7 2 9 21 9 211 9 212 10 213 13 214 14 22 15 221 15 222 16 223 17 224 20 3 21 31 21 32 21 33 22 34 23 341 23 342 24 343 27 344 29 35 31 351 31 352
More information( )
7..-8..8.......................................................................... 4.................................... 3...................................... 3..3.................................. 4.3....................................
More information2 (2016 3Q N) c = o (11) Ax = b A x = c A n I n n n 2n (A I n ) (I n X) A A X A n A A A (1) (2) c 0 c (3) c A A i j n 1 ( 1) i+j A (i, j) A (i, j) ã i
[ ] (2016 3Q N) a 11 a 1n m n A A = a m1 a mn A a 1 A A = a n (1) A (a i a j, i j ) (2) A (a i ca i, c 0, i ) (3) A (a i a i + ca j, j i, i ) A 1 A 11 0 A 12 0 0 A 1k 0 1 A 22 0 0 A 2k 0 1 0 A 3k 1 A rk
More information1 Introduction 1 (1) (2) (3) () {f n (x)} n=1 [a, b] K > 0 n, x f n (x) K < ( ) x [a, b] lim f n (x) f(x) (1) f(x)? (2) () f(x)? b lim a f n (x)dx = b
1 Introduction 2 2.1 2.2 2.3 3 3.1 3.2 σ- 4 4.1 4.2 5 5.1 5.2 5.3 6 7 8. Fubini,,. 1 1 Introduction 1 (1) (2) (3) () {f n (x)} n=1 [a, b] K > 0 n, x f n (x) K < ( ) x [a, b] lim f n (x) f(x) (1) f(x)?
More information1 Edward Waring Lagrange n {(x i, y i )} n i=1 x i p i p i (x j ) = δ ij P (x) = p i p i (x) = n y i p i (x) (1) i=1 n j=1 j i x x j x i x j (2) Runge
Edwrd Wring Lgrnge n {(x i, y i )} n x i p i p i (x j ) = δ ij P (x) = p i p i (x) = n y i p i (x) () n j= x x j x i x j (2) Runge [] [2] = ξ 0 < ξ < ξ n = b [, b] [ξ i, ξ i ] y = /( + 25x 2 ) 5 2,, 0,,
More information+ 1 ( ) I IA i i i 1 n m a 11 a 1j a 1m A = a i1 a ij a im a n1 a nj a nm.....
+ http://krishnathphysaitama-uacjp/joe/matrix/matrixpdf 1 ( ) I IA i i i 1 n m a 11 a 1j a 1m A = a i1 a ij a im a n1 a nj a nm (1) n m () (n, m) ( ) n m B = ( ) 3 2 4 1 (2) 2 2 ( ) (2, 2) ( ) C = ( 46
More informationn ξ n,i, i = 1,, n S n ξ n,i n 0 R 1,.. σ 1 σ i .10.14.15 0 1 0 1 1 3.14 3.18 3.19 3.14 3.14,. ii 1 1 1.1..................................... 1 1............................... 3 1.3.........................
More information) ] [ h m x + y + + V x) φ = Eφ 1) z E = i h t 13) x << 1) N n n= = N N + 1) 14) N n n= = N N + 1)N + 1) 6 15) N n 3 n= = 1 4 N N + 1) 16) N n 4
1. k λ ν ω T v p v g k = π λ ω = πν = π T v p = λν = ω k v g = dω dk 1) ) 3) 4). p = hk = h λ 5) E = hν = hω 6) h = h π 7) h =6.6618 1 34 J sec) hc=197.3 MeV fm = 197.3 kev pm= 197.3 ev nm = 1.97 1 3 ev
More informationP06.ppt
p.130 p.198 p.208 2 1 double weight[num]; double min, max; min = max = weight[0]; for( i= 1; i < NUM; i++ ) if ( weight[i] > max ) max = weight[i]: if ( weight[i] < min ) min = weight[i]: weight 3 maxof(a,
More informationgengo1-11
関数の再帰定義 自然数 n の階乗 n! を計算する関数を定義してみる 引数は整数 返却値も整数 n! = 1*2*3*... * (n 1)*n である ただし 0! = 1 とする int factorial(int n) int i, tmp=1; if( n>0 ) for(i=1; i
More informationkiso2-06.key
座席指定があります Linux を起動して下さい 第6回 計算機基礎実習II 計算機基礎実習II 2018 のウェブページか ら 以下の課題に自力で取り組んで下さい 第5回の復習課題(rev05) 第6回の基本課題(base06) 第5回課題の回答例 ex05-2.c 1. キーボードから整数値 a を入力すると a*a*a の値を出力することを繰り返すプログラムを作成しなさい 2. ただし 入力された
More information