|
|
|
- あきたけ おおふさ
- 6 years ago
- Views:
Transcription
1 B
2 RGB RGB HSV RGB HSV HSV RGB RGB RGB
3 OpenGL OpenGL 3D OpenGL 3D D A 57 B 60 C 64 D OpenGL (3Ddoppler.c) 69 E 3Ddoppler.c (3DdopplerFunction.c) 84 3
4 D 3D OpenGL 1.3 [1] 1 (IR) (UV ) 4
![1 1.4 X γ 1 X 1: ([1] ) 1.](/docs-images/66/54419044/images/5-0.jpg "4 2 400nm 700nm 400nm X γ 700nm 700nm 546nm")
5 1 1.4 X γ 1 X 1: ([1] ) nm 700nm 400nm X γ 700nm 700nm 546nm 435nm RGB RGB 3 5
6 nm 750nm 405nm 420nm 605nm 490nm 580nm RGB 1nm 2 2: 3 0 π 1 2 π 3: 6
7 km/s [6, p.3]
8 : ([6, p.3] ) RGB
9 2 RGB 2 RGB RGB 2.1 RGB HSV RGB R G B RGB RGB RGB RGB HSV HSV H S V % 0% 100% 0% 5: HSV ([2] ) 9
10 2.2 RGB HSV 2 RGB 2.2 RGB HSV RGB HSV 0% 100% RGB HSV [2] RGB 3 M = max(r, G, B) m = min(r, G, B) V V = M (2.1) V = 0 S H 0 S S = 255 M m M H H RGB (2.2) H = 60 G B M m { H = B R M m { H = R G M m (2.3) (2.4) (2.5) R max (2.3) G max (2.4) B max (2.5) V RGB RGB H H int while(h < 0) H += 360; H = H % 360; 10
11 2 RGB 2.3 HSV RGB 2.3 HSV RGB HSV RGB S 0 (2.2) M m RGB R = G = B = V i = H 60 (2.6) F = H 60 i (2.7) { M = V 1 S 255 { N = V 1 S 255 F { K = V 1 S 255 (1 F ) (2.8) (2.9) (2.10) (2.6) H 60 (2.6) i R G B R = V, G = K, B = M (i = 0) R = N, G = V, B = M (i = 1) R = M, G = V, B = K (i = 2) R = M, G = N, B = V (i = 3) R = K, G = M, B = V (i = 4) R = V, G = M, B = N (i = 5) (2.11) RGB HSV RGB RGB 2.4 RGB RGB RGB RGB RGB RGB 11
12 2.5 RGB 2 RGB 2.5 RGB RGB RGB RGB 700nm 546nm 435nm (CIE) 380nm 420nm 570nm 500nm 405nm 380nm 750nm RGB RGB 700nm R = nm 570nm RG 255 RGB R LIMIT P URP LE = 380 LIMIT RED = 750 RED = 700 Y ELLOW = 570 GREEN = 546 BLUE = 436 P URP LE = 405 Wavelength if(limit_red > Wavelength && Wavelength > RED){ res = LIMIT_RED - RED; R = * ((LIMIT_RED - Wavelength)/reS); else{ if(red >= Wavelength && Wavelength >= YELLOW){ R = 255; else{ if(yellow > Wavelength && Wavelength > GREEN){ res = YELLOW - GREEN; R = (255.0 * ((YELLOW - Wavelength)/reS)); else{ if(blue > Wavelength && Wavelength >= PURPLE){ res = BLUE - PURPLE; R = * ((BLUE - Wavelength)/reS); else{ if(purple > Wavelength && Wavelength > LIMIT_PURPLE){ res = PURPLE - LIMIT_PURPLE; R = (255.0 * ((PURPLE - Wavelength)/reS)); else{ R = 0; GB 12
13 v 0 x 0 τ 0 x 0 + vt v t 1 = x 0 /c t 2 = t + (x 0 vt)/c n t = (n 1)τ 0 γ { 1 1 β 2 τ = t 2 t 1 = τ 0 γ(1 β) = τ β (3.1) β γ β = v c γ = (1 β2 ) 1 2 6: 6 0 θ τ 0 L 6 t 0 = 0 x = x 0 h x v x 1 n x 2 n + 1 x 2 x 1 x 2 x 1 = vτ 0 γ (3.2) 13
14 3.1 3 n n + 1 t 1 t 2 t 1 = x 1 v + r 1 c t 2 = x 2 v + r 2 c t 2 t 1 τ (3.3) (3.4) t 2 t 1 = x 2 x 1 { x2 v + r { 2 x1 c v + r 1 c + (r 2 r 1 ) = (x 2 x 1 ) v = c(vτ 0γ) cv (r 2 r 1 ) lim x 2 x 1 (x 2 x 1 ) + vτ 0γ c = dr dx c (r 2 r 1 ) (x 2 x 1 ) (3.5) = x r = cos(π θ) = cos θ (3.6) (3.5) (3.6) τ τ = { τ 0 γ 1 + v (r 2 r 1 ) c (x 2 x 1 ) = τ 0 γ(1 β cos θ) 1 β cos θ = τ 0 (1 β 2 ) 1 2 λ f λ = λ 0 1 β cos θ (1 β 2 ) 1 2 (3.7) (3.8) f = f 0 (1 β 2 ) β cos θ (3.9) θ = 0 (3.7) (3.1) β θ (3.8) (3.9) 14
15 (3.8) 100nm 7 8 % 7 β = 0.2 β = β = 0.2 β = 0.5 β = β = 0.2 β 80% 7 β = 0.99 β 1 θ = 90 7: 8: (3.8) 9 580nm
16 : 16
17 x y 10 x v cos θ y v sin θ S x V S 10: x x (4.1) t = γ(t + V c 2 x ) x = γ(x + V t ) (4.1) 10 x v x = v cos θ v y = v sin θ (4.2) v x = v cos θ v y = v sin θ (4.3) v x = x t = v x + V 1 + V v x c 2 (4.4) 17
18 4.1 4 v y = y t = v y γ(1 + V v x c 2 ) (4.5) (4.4) (4.5) (4.2) (4.3) θ v = c v = c tan θ = v sin θ v cos θ = sin θ 1 β 2 cos θ + β (4.6) (4.6) θ β θ (4.6) t = γ(t V c 2 x) x = γ(x V t) (4.7) tan θ = sin θ 1 β 2 cos θ + β (4.8) (4.8) θ β θ 18
19 (4.6) 11 β = 0 θ = 30 θ = 60 θ = 90 θ = 90 β : (4.6) (4.6) RGB HSV V 19
20 4.3 4 V = V θ 2 π θ 2 π (4.9) θ > θ RGB 9 12: 20
21 : 14: β β = 0.5 β = 0.9 β = 0.5 PC 21
22 : 15 A A A c v c v v A 16 (a) (b) (c) (d) James Terrel 1959 [7] 22
23 : : 23
24 B A A 1 A 2 O A 18 v 1 c v v 1 18 v 2 c v < v 2 v 1 < v 2 A 1 B v 1 A 1 B O v 2 18 A 1 B 18 A 1 B B O A 1 O B O A B 18 18: v c v ( ) ( ) 24
25 A C A A 1 C A 2 A 1 A 2 c v v c v = v c v < v A C C A 3 A 1 19: A 1 A 2 25
26 (a) c (a) A 1 (a) (a) (b) (a) A 2 (a) (c) A 3 c c c c c v (a) (c) θ 1 θ 3 20: 26
27 B B A B [3, p.15] θ (5.1) a a v c v c a b β = v c b 1 β 2 θ θ tan θ = β 1 β 2 (5.1) (5.1) 21: [3, p.15] (5.1) 27
28 [5, p ] NO
29 A A A A 22 z C A z z v z L 22 C A C z C z z C C L a b C T v c 22: A B C z z C A z z A = z C + b A z z A = z C + b + v t C T A C t A t C t B 29
30 t c 22 c(t t) = (L + a) 2 + z A 2 ct = L 2 + zc 2 z A = z C + b + v t (5.2) z A ct c t = (L + a) 2 + z A 2 L 2 + zc 2 c t = (L + a) 2 + (z C + b + v t) 2 ( L 2 + zc 2 c t)2 = (L + a) 2 + (z C + b + v t) 2 (5.3) (L 2 + zc 2 ) ((L + a)2 + (z C + b ) 2 ) = X L 2 + zc 2 + β(z C + b ) = Y β = v (5.3) v t c 0 = (c 2 t 2 ) t 2 (2c L 2 + zc 2 + 2v(z C + b )) t + X 2c L 2 + zc 2 + 2v(z C + b ) ± (2c L 2 + zc 2 + 2v(z C + b )) 2 4(c 2 v 2 )X t = 2(c 2 v 2 ) t = Y ± Y 2 (1 β 2 )X c(1 β 2 ) { β v t = (1 β 2 Y ± Y ) 2 (1 β 2 )X (5.4) (5.4) v t z A = z C + b + v t z x A C z b 30
31 : 24: 23 31
32 : BC c 26: β = : β = : β =
33 β : 30: 24 31: 25 33
34 AB Terrell [7, p.3] 32: Terrell 0.8c π 1 2π π AB Terrell
35 (5.1) (5.1) B 33: (5.4) a v t (5.4) 0 (5.4) L = 0 y = 0 z C = 0 a = a h = 0 b = 0 X = a 2 Y = 0 v t v t = aβ tan θ = β (5.5) v t v t tan θ (5.1) 1 β 2 (5.1) 0 tan θ = β 35
36 z 3 2 z 2 3 y h 2 c(t t) = (L + a) 2 + (y + h) 2 + z A 2 ct = L 2 + y 2 + zc 2 z A = z C + b + v t (5.6) (5.6) 2 v t ct c t = (L + a) 2 + (y + h) 2 + z A 2 L 2 + y 2 + zc 2 c t = (L + a) 2 + (y + h) 2 + (z C + b + v t) 2 ( L 2 + y 2 + zc 2 c t)2 = (L + a) 2 + (y + h) 2 + (z C + b + v t) 2 (5.7) (L 2 +y 2 +zc 2 ) ((L+a)2 +(y+h) 2 +(z C +b ) 2 ) = X L 2 + y 2 + zc 2 + β(z C + b ) = Y β = v c t = v t = 2c L 2 + y 2 + zc 2 + 2v(z C + b ) ± (2c L 2 + y 2 + zc 2 + 2v(z C + b )) 2 4(c 2 v 2 )X { β (1 β 2 Y ± Y ) 2 (1 β 2 )X 2(c 2 v 2 ) (5.8) (5.8) v t 2 z A = z C + b + v t y, h (5.4) 36
37 (5.8) z x y 3 z x y z x x, y, z y, z x 34: 3 35: 34 36: 34 37
38 c 3 y c 0.36c 37: 0 38: : c 2 40: 0.8c 38
39 : 42: RGB 39
40 6 6 3D 3D OpenGL 6.1 OpenGL OpenGL(Open Graphics Library ) 2D/3D OS API( ) Windows UNIX MacOS C java OpenGL / C java C++ OpenGL GLUT(OpenGL Utility Toolkit) GLEW(OpenGL Extension Wrangler Library) OpenGL 2D/3D 1 3D OpenGL OpenGL GLUT 43: [9] 10 44: 40
41 6 6.2 OpenGL 3D 6.2 OpenGL 3D MAP C MAP 5 x< x y< y z< z // #define MAX_LENGTH 30 // (x ) #define MAX_DEPTH 220 // (z ) #define MAX_HEIGHT 5 // (y ) 41
42 6.2 OpenGL 3D 6 OpenGL 45: 46: : 45 48:
43 6 6.2 OpenGL 3D 1: PageUp PageDown a 15 d 15 w x Home r f g 1 end ( 50 ) 8 Z 9 Z Esc 43
44 6.3 OpenGL 3D OpenGL 3D 3D (5.8) c z : 0.8c 50:
45 6 6.3 OpenGL 3D 51: 0.8c 52: : 0.8c 54: c
46 6.3 OpenGL 3D 6 55: 49 56: 0.8c : 0.8c 58: c 46
47 D 6.4 3D c 59: 0.2c 60: 0.2c c c θ = θ = 0 47
48 6.4 3D 6 61: 0.2c 0.2c θ = : 0.2c 45 48
49 D θ 1 θ : 49
50 6.4 3D 6 64: 0.2c 59 65: : 0.2c 59 67: 0 66 θ = 90 θ = c 59 50
51 c 68: 0.2c : 70: ( ) 51
52 : 0.8c : 0.8c 73: 0.8c c
53 : 75: c 76: 53
54 7 7 OpenGL 3D 1 3D 1 1 [8] 20km/h 1 54
55 7 web HSV 55
56 [1] wikipedia 95%B0 ( ) [2] wikipedia HSV BA%E9%96%93 ( ) [3] 1996 [4] 2011 [5] / (2001) [6] U Kraus First-person visualizations of the special and general theory of relativity EUROPEAN JOURNAL PHYSICS [p.1-p.6] 2007 [7] James Terrell Invisibility of the Lorentz Contraction PHYSICAL REVIEW 116 [p.1041-p.1045] 1959 [8] George Gamow( )/Russell Stannard( ) (2001) [9] MICC OpenGL (2010) 56
57 A A : : 77 79: c 57
58 A 80: 0 81: 0 82: 1 83: 2 84: 3 85: 0.8c 86: 0.8c 87: 0.8c 88: 0.8c 89: 0.8c 58
59 A 90: 0 91: 0.8c 59
60 B 92: 93: 92 B ( ) // #define BETA 0.95 // /LengthPer //100 1/ #define LengthPer 10 typedef int boolean; #define true 1 #define false 0 #define OB false // #define FIXATION true // ( ) #define DOPPLER false // // #define MAX_LENGTH 30 // (x ) #define MAX_DEPTH 70 // (z ) #define MAX_HEIGHT 5 // (y ) //1 #define MAX_Lambda 2 // #define CAMERAx 15 #define CAMERAy 1 #define CAMERAz 20 // #define CGlength 15.0 // #define GAZEx CAMERAx #define GAZEy CAMERAy 60
61 B #define GAZEz CAMERAz-CGlength // #define PointSize 100.f // #define CRITERIA_X 5.0 #define CRITERIA_Y 0.0 #define CRITERIA_Z -2.0 // #define ROTATION 15.0 // #define M_PI #define TurnoverRate M_PI*ROTATION/180.0 // #define Stride 1.0 // #define LigPosX 0.4f #define LigPosY 0.2f #define LigPosZ 0.3f #define LigPosW 0.f // #define AMB_R 0.5f #define AMB_G 0.5f #define AMB_B 0.5f #define AMB_A 1.f // #define DIFF_R 1.f #define DIFF_G 1.f #define DIFF_B 1.f #define DIFF_A 1.f // #define SPEC_R 1.f #define SPEC_G 1.f #define SPEC_B 1.f #define SPEC_A 1.f // #define glshin 100.2f // // #define GROUND_SPEC_R 0.8f #define GROUND_SPEC_G 0.8f #define GROUND_SPEC_B 0.8f 61
62 B #define GROUND_SPEC_A 1.f // #define GROUND_AandD_R 0.8f #define GROUND_AandD_G 0.8f #define GROUND_AandD_B 0.8f #define GROUND_AandD_A 1.f // // #define LIMIT_PURPLE 380 #define LIMIT_RED 750 // (CIE) 1931 CIE // #define RED 700 #define GREEN 546 #define BLUE 436 // #define ORANGE 625 #define YELLOW 570 #define CYAN 500 #define INDIGO 420 #define PURPLE 405 /* nm nm nm nm nm nm */ /* #define ORANGE 605 #define YELLOW 580 #define CYAN 490 #define INDIGO 420 #define PURPLE 405 */ typedef struct{ boolean object; int WaveSum; int lambda[max_lambda]; double Yi[MAX_Lambda]; XYZobject; // (MAX_Lambda) 62
63 B typedef struct{ boolean object; double z; int WaveSum; int lambda[max_lambda]; double Yi[MAX_Lambda]; XYobject; typedef struct{ boolean object; boolean visibility; int WaveSum; int lambda[max_lambda]; double Yi[MAX_Lambda]; boolean left; //X boolean right; //X boolean buttom; //Y boolean top; //Y boolean front; //Z boolean back; //Z XYZvisibility; typedef struct{ double R; double G; double B; RGB; typedef struct{ double V; double S; double H; HSV; typedef struct{ boolean object; double R; double G; double B; XYZ; int round1(double str); void Ground(); HSV RGBtoHSV(RGB rgb); RGB HSVtoRGB(HSV hsv); int redjudg(int H0); int greenjudg(int H0); int bluejudg(int H0); RGB DopplerSynthesis(double X, double Y, double Z, int lambda[], double Yi[], int sum); double Lorentz(double l, double vc); 63
64 C double Deformation(double LL, double Y, double zc, double a, double h, double b, double vc); C #include <stdio.h> #include <stdlib.h> #include <my/myhead.h> #pragma warning(disable : 4996) /* 1 :XYZ (1 XYZ ) 2 : x, y, z, n,, Yi, sum ( 7 ) */ int main(int argc, char *argv[]){ int x1, y1, z1, x2, y2, z2; int MAXx=0, MAXy=0, MAXz=0; boolean ob; char name[64]; FILE *fp; XYZobject xyz[max_length][max_height][max_depth]; int x, y, z, tmpx=-1, tmpy=-1, tmpz=-1; int lambda, sum, count; double Yi; // if(argc > 2){ printf(" \n"); scanf("%d", &x); exit(1); /*argc=1 */ if(argc==1){ do{ printf("3d \n"); printf(" x(%d ), y(%d ), z(%d ) 3 \n", MAX_LENGTH, MAX_HEIGHT, MAX_DEPTH); scanf("%d %d %d", &MAXx, &MAXy, &MAXz); while(!(1<=maxx && MAXx<=MAX_LENGTH && 1<=MAXy && MAXy<=MAX_HEIGHT && 1<=MAXz && MAXz<=MAX_DEPTH) for(x=0; x<maxx; x++){ for(y=0; y<maxy; y++){ for(z=0; z<maxz; z++){ xyz[x][y][z].object = false; xyz[x][y][z].wavesum = 0; 64
65 C xyz[x][y][z].lambda[0] = 0; xyz[x][y][z].yi[0] = 0.0; printf(" \n\n"); else{ if ((fp = fopen(argv[1], "r")) == NULL){ fprintf(stderr, " %s \n", argv[1]); // scanf("%d", &x); exit(1); else{ /* */ fscanf(fp, "%d %d %d", &MAXx, &MAXy, &MAXz); while(fscanf(fp, "%d %d %d %d %d %lf %d", &x, &y, &z, &ob, &lambda, &Yi, &sum)!= EOF){ // tmp if(!(tmpx==x && tmpy==y && tmpz==z)){ // count = 0; xyz[x][y][z].object = ob; xyz[x][y][z].wavesum = sum; xyz[x][y][z].lambda[count] = lambda; xyz[x][y][z].yi[count] = Yi; else{ // count++; xyz[x][y][z].lambda[count] = lambda; xyz[x][y][z].yi[count] = Yi; tmpx = x; tmpy = y; tmpz = z; /* */ fclose(fp); 65
66 C printf(" \n\n"); /* */ /* RGB 4 */ /* */ while(1){ printf(" 0 x=%d y=%d z=%d \n", MAXx, MAXy, MAXz); printf(" ( )\n"); printf(" < \n"); scanf("%d %d %d", &x1, &y1, &z1); /* */ if((0<=x1 && x1<maxx) && (0<=y1 && y1<maxy) && (0<=z1 && z1<maxz)){ /* or */ while(1){ printf(" x:%d %d y:%d %d z:%d %d \n", x1, MAXx, y1, MAXy, z1, MAXz); printf(" ( )\n"); scanf("%d %d %d", &x2, &y2, &z2); if(((0<=x2 && x2<maxx && x1<=x2) && (0<=y2 && y2<maxy && y1<=y2) && (0<=z2 && z2<maxz && z1<=z2)) (x2==-1 && y2==-1 && z2==-1)) break; printf(" \n"); else{ /* */ if(x1==-1 && y1==-1 && z1==-1) break; /* */ if(x2!=-1 && y2!=-1 && z2!=-1){ while(1){ printf(" %d %d \n", true, false); scanf("%d", &ob); if(ob==false ob==true) break; if(ob == 1){ // count = 0; /* */ 66
67 C while(max_lambda > count){ while(1){ printf(" (nm) (0 1.0) \n"); scanf("%d %lf", &lambda, &Yi); // if(0.0 <= Yi && Yi <= 1.0) break; else printf(" \n \n"); /* */ for(x=x1; x<=x2; x++){ for(y=y1; y<=y2; y++){ for(z=z1; z<=z2; z++){ xyz[x][y][z].object = true; xyz[x][y][z].wavesum = count + 1; xyz[x][y][z].lambda[count] = lambda; xyz[x][y][z].yi[count] = Yi; count++; if(max_lambda > count){ // printf(" \n"); printf(" 1 0 \n"); scanf("%d", &ob); if(ob==0) break; printf(" %d \n", MAX_Lambda - count); printf(" \n"); else{ /* */ for(x=x1; x<=x2; x++){ for(y=y1; y<=y2; y++){ for(z=z1; z<=z2; z++){ xyz[x][y][z].object = false; xyz[x][y][z].wavesum = 0; 67
68 C printf(" \n"); /* */ while(1){ printf(" \n"); scanf("%s", name); if ((fp = fopen(name, "w")) == NULL) { fprintf(stderr, " %s \n", name); else break; /* 1 xyz */ fprintf(fp, "%d\t %d\t %d\n", MAXx, MAXy, MAXz); for(x=0; x<maxx; x++){ for(y=0; y<maxy; y++){ for(z=0; z<maxz; z++){ count = 0; fprintf(fp, "%d\t %d\t %d\t %d\t %d\t %f\t %d\n", x, y, z, xyz[x][y][z].object, xyz[x][y][z].lambda[count], xyz[x][y][z].yi[count], xyz[x][y][z].wavesum); if(xyz[x][y][z].wavesum > 1){ /* */ count = 1; while(xyz[x][y][z].wavesum > count){ fprintf(fp, "%d\t %d\t %d\t %d\t %d\t %f\t %d\n", x, y, z, xyz[x][y][z].object, xyz[x][y][z].lambda[count], xyz[x][y][z].yi[count], xyz[x][y][z].wavesum); count++; fclose(fp); return 0; 68
69 D OPENGL (3DDOPPLER.C) D OpenGL (3Ddoppler.c) #include <stdio.h> #include <stdlib.h> #include <string.h> #include <math.h> #include <my/myhead.h> #include <gl/glut.h> #pragma warning(disable : 4996) #define DRAW_LENGTH 1.0/LengthPer /* Esc PageUp PageDown Home g 1 adwx a: d: w: x: f ( ) r */ // XYZvisibility xyz[max_length][max_height][max_depth]; XYobject xyzsave1[max_length*lengthper][max_height*lengthper]; XYobject xyzsave2[max_length*lengthper][max_height*lengthper]; // int MAXx=0, MAXy=0, MAXz=0; // double CameraX=CAMERAx, CameraY=CAMERAy, CameraZ=Lorentz(CAMERAz, BETA); // double GazeX=GAZEx, GazeY=GAZEy, GazeZ=GAZEz; // ( z 0 ) double thetaxz=acos((gazez-cameraz)/cglength); // // // double observerx=camerax, observery=cameray, observerz=lorentz(cameraz, BETA); double CriteriaX=observerX+CRITERIA_X, CriteriaY=observerY+CRITERIA_Y, CriteriaZ=observerZ+CRITERI 69
70 D OPENGL (3DDOPPLER.C) // ( ) GLfloat lightamb[] = {AMB_R, AMB_G, AMB_B, AMB_A; GLfloat lightdiff[] = {DIFF_R, DIFF_G, DIFF_B, DIFF_A; GLfloat lightspec[] = {SPEC_R, SPEC_G, SPEC_B, SPEC_A; // GLfloat lightpos[] = {LigPosX, LigPosY, LigPosZ, LigPosW; // ( ) GLfloat Spec[] = { f, f, f, 1.f; GLfloat Shin = glshin; // void Architecture(FILE *fp){ int x, y, z, ob, lambda, sum; int count=0, tmpsum, tmpx, tmpy, tmpz; double Yi; fscanf(fp, "%d %d %d", &MAXx, &MAXy, &MAXz); while(fscanf(fp, "%d %d %d %d %d %lf %d", &x, &y, &z, &ob, &lambda, &Yi, &sum)!= EOF){ count = 1; tmpsum = sum; tmpx = x; tmpy = y; tmpz = z; xyz[x][y][z].object = ob; xyz[x][y][z].lambda[0] = lambda; xyz[x][y][z].yi[0] = Yi; xyz[x][y][z].wavesum = sum; // while((tmpsum>count) && (tmpx==x && tmpy==y && tmpz==z) && MAX_Lambda>count){ fscanf(fp, "%d %d %d %d %d %lf %d", &x, &y, &z, &ob, &lambda, &Yi, &sum); xyz[x][y][z].lambda[count] = lambda; xyz[x][y][z].yi[count] = Yi; count++; void cube(double x, double y, double z, RGB rgb){ float Rf, Gf, Bf; // glcolormaterial(gl_front_and_back, GL_AMBIENT); 70
71 D OPENGL (3DDOPPLER.C) Rf = (float)(rgb.r*0.9); Gf = (float)(rgb.g*0.9); Bf = (float)(rgb.b*0.9); glcolor4f(rf, Gf, Bf, 1.f); // glcolormaterial(gl_front_and_back, GL_DIFFUSE); Rf = (float)rgb.r; Gf = (float)rgb.g; Bf = (float)rgb.b; glcolor4f(rf, Gf, Bf, 1.f); // glpushmatrix(); gltranslatef(x, y, z); glutsolidcube(draw_length); glpopmatrix(); void clearsave1(){ int x, y; for(x=0; x<max_length*lengthper; x++){ for(y=0; y<max_height*lengthper; y++){ xyzsave1[x][y].object = false; void clearsave2(){ int x, y; for(x=0; x<max_length*lengthper; x++){ for(y=0; y<max_height*lengthper; y++){ xyzsave2[x][y].object = false; void Save1(int x, int y, int z, double X, double Y, double Z){ int arrayx, arrayy; // 2 double xdif, ydif, zdif; // double CoordinateZ1, oldz; RGB rgb; // arrayx = round1(x/draw_length); arrayy = round1(y/draw_length); // xdif = X - CriteriaX; ydif = Y - CriteriaY; 71
72 D OPENGL (3DDOPPLER.C) zdif = Z - Lorentz(CriteriaZ, BETA); // z XYZ, xyz beta xyzsave1[arrayx][arrayy].z = Deformation(CRITERIA_X, CRITERIA_Y, CRITERIA_Z, xdif, ydif, zdif, BET // XYZ 7 rgb = DopplerSynthesis(observerX-X, observery-y, observerz-z, xyz[x][y][z].lambda, xyz[x][y][z].yi // z // z if(z>0){ if(xyzsave2[arrayx][arrayy].object==true){ zdif = (Z-DRAW_LENGTH) - Lorentz(CriteriaZ, BETA); oldz = Deformation(CRITERIA_X, CRITERIA_Y, CRITERIA_Z, xdif, ydif, zdif, BETA); for(coordinatez1=oldz; CoordinateZ1 < xyzsave1[arrayx][arrayy].z; CoordinateZ1 += DRAW_LENGTH){ cube(x, Y, CoordinateZ1, rgb); cube(x, Y, xyzsave1[arrayx][arrayy].z, rgb); // xyzsave1[arrayx][arrayy].object = true; xyzsave1[arrayx][arrayy].wavesum = xyz[x][y][z].wavesum; memcpy(xyzsave1[arrayx][arrayy].lambda, xyz[x][y][z].lambda, MAX_Lambda); memcpy(xyzsave1[arrayx][arrayy].yi, xyz[x][y][z].yi, MAX_Lambda); void Save2(int x, int y, int z, double X, double Y, double Z){ int arrayx, arrayy; // 2 double xdif, ydif, zdif; // double CoordinateZ2, oldz; RGB rgb; // arrayx = round1(x/draw_length); arrayy = round1(y/draw_length); // xdif = X - CriteriaX; ydif = Y - CriteriaY; zdif = Z - Lorentz(CriteriaZ, BETA); // z XYZ, xyz beta xyzsave2[arrayx][arrayy].z = Deformation(CRITERIA_X, CRITERIA_Y, CRITERIA_Z, xdif, ydif, zdif, BET // XYZ 7 rgb = DopplerSynthesis(observerX-X, observery-y, observerz-z, xyz[x][y][z].lambda, xyz[x][y][z].yi // z // z 72
73 D OPENGL (3DDOPPLER.C) if(z>0){ if(xyzsave1[arrayx][arrayy].object==true){ zdif = (Z-DRAW_LENGTH) - Lorentz(CriteriaZ, BETA); oldz = Deformation(CRITERIA_X, CRITERIA_Y, CRITERIA_Z, xdif, ydif, zdif, BETA); for(coordinatez2=oldz; CoordinateZ2 < xyzsave2[arrayx][arrayy].z; CoordinateZ2 += DRAW_LENGTH){ cube(x, Y, CoordinateZ2, rgb); cube(x, Y, xyzsave2[arrayx][arrayy].z, rgb); // xyzsave2[arrayx][arrayy].object = true; xyzsave2[arrayx][arrayy].wavesum = xyz[x][y][z].wavesum; memcpy(xyzsave2[arrayx][arrayy].lambda, xyz[x][y][z].lambda, MAX_Lambda); memcpy(xyzsave2[arrayx][arrayy].yi, xyz[x][y][z].yi, MAX_Lambda); void XYdrawSave1(int z, double Z, int count){ int x, y; double Y, X; double sidetop, sidebuttom; // for(y=0; y<=maxy-1; y++){ for(x=0; x<=maxx-1; x++){ if(xyz[x][y][z].visibility == true){ if((count%lengthper==0 && xyz[x][y][z].front==true) ((count+1)%lengthper==0 && xyz[x][y][z].bac if(xyz[x][y][z].front==true && count%lengthper==0){ //front for(y=y; Y<y+1-DRAW_LENGTH/10.0; Y+=DRAW_LENGTH){ for(x=x; X<x+1-DRAW_LENGTH/10.0; X+=DRAW_LENGTH){ Save1(x,y, z, X, Y, Z); else{ if(xyz[x][y][z].back==true && (count+1)%lengthper==0){ //back for(y=y; Y<y+1-DRAW_LENGTH/10.0; Y+=DRAW_LENGTH){ for(x=x; X<x+1-DRAW_LENGTH/10.0; X+=DRAW_LENGTH){ Save1(x,y, z, X, Y, Z); else{ // sidetop = 1.0; sidebuttom = 0.0; //top if(xyz[x][y][z].top==true){ sidetop -= DRAW_LENGTH; Y=y+1.0-DRAW_LENGTH; for(x=x; X<x+1-DRAW_LENGTH/100.0; X+=DRAW_LENGTH){ Save1(x,y, z, X, Y, Z); 73
74 D OPENGL (3DDOPPLER.C) //buttom if(xyz[x][y][z].buttom==true){ sidebuttom += DRAW_LENGTH; Y=y; for(x=x; X<x+1-DRAW_LENGTH/10.0; X+=DRAW_LENGTH){ Save1(x,y, z, X, Y, Z); //right if(xyz[x][y][z].right==true){ X=x+1-DRAW_LENGTH; for(y=y+sidebuttom; Y<y+sideTop-DRAW_LENGTH/10.0; Y+=DRAW_LENGTH){ Save1(x,y, z, X, Y, Z); //left if(xyz[x][y][z].left==true){ X=x; for(y=y+sidebuttom; Y<y+sideTop-DRAW_LENGTH/10.0; Y+=DRAW_LENGTH){ Save1(x,y, z, X, Y, Z); void XYdrawSave2(int z, double Z, int count){ int x, y; double Y, X; double sidetop, sidebuttom; // for(y=0; y<=maxy-1; y++){ for(x=0; x<=maxx-1; x++){ if(xyz[x][y][z].visibility == true){ if((count%lengthper==0 && xyz[x][y][z].front==true) ((count+1)%lengthper==0 && xyz[x][y][z].bac if(xyz[x][y][z].front==true && count%lengthper==0){ //front for(y=y; Y<y+1-DRAW_LENGTH/10.0; Y+=DRAW_LENGTH){ for(x=x; X<x+1-DRAW_LENGTH/10.0; X+=DRAW_LENGTH){ Save2(x,y, z, X, Y, Z); else{ if(xyz[x][y][z].back==true && (count+1)%lengthper==0){ //back for(y=y; Y<y+1-DRAW_LENGTH/10.0; Y+=DRAW_LENGTH){ for(x=x; X<x+1-DRAW_LENGTH/10.0; X+=DRAW_LENGTH){ Save2(x,y, z, X, Y, Z); else{ 74
75 D OPENGL (3DDOPPLER.C) // sidetop = 1.0; sidebuttom = 0.0; //top if(xyz[x][y][z].top==true){ sidetop -= DRAW_LENGTH; Y=y+1.0-DRAW_LENGTH; for(x=x; X<x+1-DRAW_LENGTH/100.0; X+=DRAW_LENGTH){ Save2(x,y, z, X, Y, Z); //buttom if(xyz[x][y][z].buttom==true){ sidebuttom += DRAW_LENGTH; Y=y; for(x=x; X<x+1-DRAW_LENGTH/10.0; X+=DRAW_LENGTH){ Save2(x,y, z, X, Y, Z); //right if(xyz[x][y][z].right==true){ X=x+1-DRAW_LENGTH; for(y=y+sidebuttom; Y<y+sideTop-DRAW_LENGTH/10.0; Y+=DRAW_LENGTH){ Save2(x,y, z, X, Y, Z); //left if(xyz[x][y][z].left==true){ X=x; for(y=y+sidebuttom; Y<y+sideTop-DRAW_LENGTH/10.0; Y+=DRAW_LENGTH){ Save2(x,y, z, X, Y, Z); // void testsave2(int x, int y, int z, double X, double Y, double Z){ double xdif, ydif, zdif; // double zz; RGB rgb; // xdif = X - CriteriaX; ydif = Y - CriteriaY; 75
76 D OPENGL (3DDOPPLER.C) zdif = Z - Lorentz(CriteriaZ, BETA); // z XYZ, xyz beta zz = Deformation(CRITERIA_X, CRITERIA_Y, CRITERIA_Z, xdif, ydif, zdif, BETA); // XYZ 7 rgb = DopplerSynthesis(observerX-X, observery-y, observerz-z, xyz[x][y][z].lambda, xyz[x][y][z].yi // z cube(x, Y, zz, rgb); // // void TESTfunc(){ double X, Y, Z=0.0; double startx, endx; startx = CAMERAx ; endx = CAMERAx ; for(z=0.0; Z<1.0; Z+=DRAW_LENGTH){ for(x=startx; X<endX; X+=DRAW_LENGTH){ for(y=0.0; Y<10.0; Y+=DRAW_LENGTH){ testsave2(11, 0, 3, X, Y, Z); // // // void display(void){ int z; int count=0, i; double Z, lorentz, lorentzdrawlength; //float Rf, Gf, Bf; //RGB TMP; glclear(gl_color_buffer_bit GL_DEPTH_BUFFER_BIT); glloadidentity(); glulookat(camerax, CameraY, CameraZ, GazeX, GazeY, GazeZ, 0.0, 1.0, 0.0); // gllightfv(gl_light0, GL_POSITION, lightpos); glenable(gl_color_material); glmaterialf(gl_front_and_back, GL_SHININESS, Shin); // Ground(); 76
77 D OPENGL (3DDOPPLER.C) //test // if(ob==true){ glcolor4f(1.f, 0.f, 0.f, 1.f); glpushmatrix(); gltranslatef(camerax, CAMERAy, Lorentz(CAMERAz, BETA)); glutsolidcube(0.5); glpopmatrix(); glcolor4f(1.f, 1.f, 0.f, 1.f); glpushmatrix(); gltranslatef(camerax, CAMERAy, Deformation(CRITERIA_X, CRITERIA_Y, CRITERIA_Z, -CRITERIA_X, -CRITE //glutsolidcube(0.5); glpopmatrix(); //TESTfunc(); if(fixation==false){ observerx=camerax; observery=cameray; observerz=cameraz; else{ observerx=camerax; observery=cameray; observerz=cameraz; CriteriaX=observerX+CRITERIA_X; CriteriaY=observerY+CRITERIA_Y; CriteriaZ=observerZ+CRITERIA_Z; count = 0; //Z lorentz = Lorentz(1.0, BETA); lorentzdrawlength = lorentz * DRAW_LENGTH; //Z Vector //count Save Save2 for(z=0; z<=maxz-1; z++){ for(i=0; i<lengthper; i++){ Z = count * lorentzdrawlength; if(count%2==0){ clearsave2(); XYdrawSave2(z, Z, count); else{ clearsave1(); XYdrawSave1(z, Z, count); count++; 77
78 D OPENGL (3DDOPPLER.C) gldisable(gl_color_material); glutswapbuffers(); void reshape(int w, int h){ glviewport(0,0,w,h); glmatrixmode(gl_projection); glloadidentity(); gluperspective(30.0, double(w)/h, 1.0, ); glmatrixmode(gl_modelview); void idle(void){ glutpostredisplay(); // // void otherinit(void){ glclearcolor(1.f, 1.f, 1.f, 1.f); glcleardepth(1.f); glenable(gl_depth_test); // ( ) gllightfv(gl_light0, GL_AMBIENT, lightamb); gllightfv(gl_light0, GL_DIFFUSE, lightdiff); gllightfv(gl_light0, GL_SPECULAR, lightspec); // glenable(gl_light0); glenable(gl_lighting); // glenable(gl_normalize); // ( ) void keyboard(unsigned char key, int x, int y){ double height; //[ESC] if(key == 27){ exit(0); switch(key){ 78
79 D OPENGL (3DDOPPLER.C) case a : // thetaxz = thetaxz + TurnoverRate; if(thetaxz>m_pi) thetaxz = thetaxz - M_PI*2; GazeX = CameraX + (CGlength * sin(thetaxz)); GazeZ = CameraZ + (CGlength * cos(thetaxz)); break; case d : // thetaxz = thetaxz - TurnoverRate; if(thetaxz<-1*m_pi) thetaxz = thetaxz + M_PI*2; GazeX = CameraX + (CGlength * sin(thetaxz)); GazeZ = CameraZ + (CGlength * cos(thetaxz)); break; case w : // GazeY++; break; case x : // GazeY--; break; case 4 : //90 thetaxz = thetaxz + M_PI/2.0; if(thetaxz>m_pi) thetaxz = thetaxz - M_PI*2; GazeX = CameraX + (CGlength * sin(thetaxz)); GazeZ = CameraZ + (CGlength * cos(thetaxz)); break; case 6 : //90 thetaxz = thetaxz - M_PI/2.0; if(thetaxz<-1*m_pi) thetaxz = thetaxz + M_PI*2; GazeX = CameraX + (CGlength * sin(thetaxz)); GazeZ = CameraZ + (CGlength * cos(thetaxz)); break; case 2 : // thetaxz = thetaxz - M_PI; if(thetaxz<-1*m_pi) thetaxz = thetaxz + M_PI*2; GazeX = CameraX + (CGlength * sin(thetaxz)); GazeZ = CameraZ + (CGlength * cos(thetaxz)); break; case f : // GazeY = CameraY; break; case r : // GazeX=CameraX; GazeY=CameraY; GazeZ=CameraZ-CGlength; thetaxz=acos((gazez-cameraz)/cglength); break; case g : // 1 height = CameraY; CameraY = 1.0; GazeY += height; break; 79
80 D OPENGL (3DDOPPLER.C) case 8 : // CameraZ--; GazeZ--; break; case 9 : // CameraZ++; GazeZ++; break; // void specialkey(int key, int x, int y){ double theta; switch(key){ case GLUT_KEY_UP: // theta = thetaxz; CameraX += Stride * sin(theta); GazeX += Stride * sin(theta); CameraZ += Stride * cos(theta); GazeZ += Stride * cos(theta); break; case GLUT_KEY_DOWN: // theta = thetaxz - M_PI; if(theta>-m_pi) theta = theta + M_PI*2; CameraX += Stride * sin(theta); GazeX += Stride * sin(theta); CameraZ += Stride * cos(theta); GazeZ += Stride * cos(theta); break; case GLUT_KEY_LEFT: // theta = thetaxz + (M_PI/2.0); if(theta>m_pi) theta = theta - M_PI*2; CameraX += Stride * sin(theta); GazeX += Stride * sin(theta); CameraZ += Stride * cos(theta); GazeZ += Stride * cos(theta); break; case GLUT_KEY_RIGHT: // theta = thetaxz - (M_PI/2.0); if(theta>-m_pi) theta = theta + M_PI*2; CameraX += Stride * sin(theta); GazeX += Stride * sin(theta); CameraZ += Stride * cos(theta); GazeZ += Stride * cos(theta); break; case GLUT_KEY_PAGE_UP: // CameraY++; GazeY++; 80
81 D OPENGL (3DDOPPLER.C) break; case GLUT_KEY_PAGE_DOWN: // CameraY--; GazeY--; break; case GLUT_KEY_HOME: // CameraX=CAMERAx; CameraY=CAMERAy; GazeX = CameraX; GazeY = CameraY; CameraZ=CAMERAz; GazeZ = CameraZ-CGlength; break; case GLUT_KEY_END: //test 50 theta = thetaxz; CameraX += Stride*50.0 * sin(theta); GazeX += Stride*50.0 * sin(theta); CameraZ += Stride*50.0 * cos(theta); GazeZ += Stride*50.0 * cos(theta); break; // void surfaceoptimization(int x, int y, int z){ //X if(x==0){ xyz[x][y][z].left = true; else{ if(xyz[x-1][y][z].object==false) xyz[x][y][z].left = true; else xyz[x][y][z].left = false; //X if(x==maxx-1){ xyz[x][y][z].right = true; else{ if(xyz[x+1][y][z].object==false) xyz[x][y][z].right = true; else xyz[x][y][z].right = false; //Y if(y==0){ xyz[x][y][z].buttom = true; else{ if(xyz[x][y-1][z].object==false) 81
82 D OPENGL (3DDOPPLER.C) xyz[x][y][z].buttom = true; else xyz[x][y][z].buttom = false; //Y if(y==maxy-1){ xyz[x][y][z].top = true; else{ if(xyz[x][y+1][z].object==false) xyz[x][y][z].top = true; else xyz[x][y][z].top = false; //Z if(z==0){ xyz[x][y][z].front = true; else{ if(xyz[x][y][z-1].object==false) xyz[x][y][z].front = true; else xyz[x][y][z].front = false; //Z if(z==maxz-1){ xyz[x][y][z].back = true; else{ if(xyz[x][y][z+1].object==false) xyz[x][y][z].back = true; else xyz[x][y][z].back = false; // void Optimization(){ int x, y, z; for(x=0; x<maxx; x++){ for(y=0; y<maxy; y++){ for(z=0; z<maxz; z++){ if(xyz[x][y][z].object == true){ xyz[x][y][z].visibility = false; if(x==0 x==maxx-1 y==0 y==maxy-1 z==0 z==maxz-1) xyz[x][y][z].visibility = true; else{if(!(xyz[x-1][y-1][z-1].object==true && xyz[x-1][y-1][z].object==true && xyz[x-1][y-1][z+1].o xyz[x-1][y][z-1].object ==true && xyz[x-1][y][z].object ==true && xyz[x-1][y][z+1].object ==t 82
83 D OPENGL (3DDOPPLER.C) xyz[x-1][y+1][z-1].object==true && xyz[x-1][y+1][z].object==true && xyz[x-1][y+1][z+1].object==t xyz[x+1][y-1][z-1].object==true && xyz[x+1][y-1][z].object==true && xyz[x+1][y-1][z+1].object==t xyz[x+1][y+1][z-1].object==true && xyz[x+1][y+1][z].object==true && xyz[x+1][y+1][z+1].object==t xyz[x][y-1][z-1].object ==true && xyz[x][y-1][z].object ==true && xyz[x][y-1][z+1].object ==t xyz[x][y][z-1].object ==true && xyz[x][y][z+1].object ==true && xyz[x][y+1][z-1].object ==true && xyz[x][y+1][z].object ==true && xyz[x][y+1][z+1].object ==t xyz[x+1][y][z-1].object ==true && xyz[x+1][y][z].object ==true && xyz[x+1][y][z+1].object ==t )){ // xyz[x][y][z].visibility = true; if(xyz[x][y][z].visibility == true){ surfaceoptimization(x, y, z); // // int main(int argc, char *argv[]){ FILE *file; int end; //RGB rgb; if ((file = fopen(argv[1], "r")) == NULL){ fprintf(stderr, " %s \n", argv[1]); scanf("%d", &end); exit(1); else{ Architecture(file); /* */ fclose(file); printf(" \n\n"); // Optimization(); glutinit(&argc,argv); glutinitdisplaymode(glut_rgba GLUT_DEPTH GLUT_DOUBLE); glutinitwindowsize(640, 480); glutcreatewindow("3d Pseudo Shooting"); // glutdisplayfunc(display); glutreshapefunc(reshape); glutidlefunc(idle); 83
84 E 3DDOPPLER.C (3DDOPPLERFUNCTION.C) // otherinit(); // glutkeyboardfunc(keyboard); // glutspecialfunc(specialkey); // glutignorekeyrepeat(gl_true); glutmainloop(); return 0; E 3Ddoppler.c (3DdopplerFunction.c) #include <stdio.h> #include <stdlib.h> #include <string.h> #include <math.h> #include <my/myhead.h> #include <gl/glut.h> #pragma warning(disable : 4996) int round1(double str){ int intr; intr = str+0.5; return intr; // void Ground(){ glcolormaterial(gl_front_and_back,gl_specular); glcolor4f(ground_spec_r, GROUND_SPEC_G, GROUND_SPEC_B, GROUND_SPEC_A); glcolormaterial(gl_front_and_back,gl_ambient_and_diffuse); glcolor4f(ground_aandd_r, GROUND_AandD_G, GROUND_AandD_B, GROUND_AandD_A); glnormal3f(0,1,0); glbegin(gl_quads); glvertex3f((float)max_length*(-100.f),-0.5f,(float)max_depth*(-100.f)); glvertex3f((float)max_length*100.f,-0.5f,(float)max_depth*(-100.f)); glvertex3f((float)max_length*100.f,-0.5f,(float)max_depth*100.f); 84
85 E 3DDOPPLER.C (3DDOPPLERFUNCTION.C) glvertex3f((float)max_length*(-100.f),-0.5f,(float)max_depth*100.f); glend(); //RGB HSV HSV HSV RGBtoHSV(RGB rgb){ int min, max; HSV hsv; if(rgb.r >= rgb.g){ if(rgb.r >= rgb.b) max = rgb.r; else max = rgb.b; if(rgb.g <= rgb.b) min = rgb.g; else min = rgb.b; else{ if(rgb.g >= rgb.b) max = rgb.g; else max = rgb.b; if(rgb.r <= rgb.b) min = rgb.r; else min = rgb.b; hsv.v = max; //max 0 if(hsv.v == 0.0){ hsv.s=0.0; hsv.h=0.0; return hsv; hsv.s = * (max-min)/max; //S=0 H if(hsv.s == 0){ hsv.h = 0; return hsv; 85
86 E 3DDOPPLER.C (3DDOPPLERFUNCTION.C) // if(rgb.r >= rgb.g){ if(rgb.r >= rgb.b) hsv.h = 60.0 * (double)(rgb.g-rgb.b)/(max-min); else hsv.h = 60.0 * (4.0+(double)(rgb.R-rgb.G)/(max-min)); else{ if(rgb.g >= rgb.b) hsv.h = 60.0 * (2.0+(double)(rgb.B-rgb.R)/(max-min)); else hsv.h = 60.0 * (4.0+(double)(rgb.R-rgb.G)/(max-min)); // while(!(0.0 <= hsv.h && hsv.h < 360.0)){ hsv.h = hsv.h - (360.0 * (hsv.h/abs(hsv.h))); return hsv; //HSV RGB RGB RGB HSVtoRGB(HSV hsv){ RGB rgb; int i; double F, M, N, K, V; if(hsv.s==0){ rgb.r = hsv.v; rgb.g = hsv.v; rgb.b = hsv.v; return rgb; i = (int)floor(hsv.h/60.0); F = hsv.h/60.0-i; M = (hsv.v * (1.0 - hsv.s/255.0)); N = (hsv.v * (1.0 - F * hsv.s/255.0)); K = (hsv.v * (1.0 - (1.0-F) * hsv.s/255.0)); V = hsv.v; //printf("%d \n", i); switch(i){ case 0: rgb.r = V; 86
87 E 3DDOPPLER.C (3DDOPPLERFUNCTION.C) rgb.g = K; rgb.b = M; break; case 1: rgb.r = N; rgb.g = V; rgb.b = M; break; case 2: rgb.r = M; rgb.g = V; rgb.b = K; break; case 3: rgb.r = M; rgb.g = N; rgb.b = V; break; case 4: rgb.r = K; rgb.g = M; rgb.b = V; break; case 5: rgb.r = V; rgb.g = M; rgb.b = N; break; return rgb; // RGB // red H0 H int redjudg(int H0){ double res=0; double H = (double)h0; int x; if(limit_red>h && H>RED){ res = (double)(limit_red - RED); x = (int)(255.0 * ((-H + (double)limit_red)/res)); else{ if(red>=h && H>=YELLOW){ x = 255; else{ if(yellow>h0 && H0>GREEN){ 87
88 E 3DDOPPLER.C (3DDOPPLERFUNCTION.C) res = (double)(yellow - GREEN); x = (int)( (255.0 * ((-H + (double)yellow)/res))); else{ if(blue>h0 && H0>=PURPLE){ res = (double)(blue - PURPLE); x = (int)(255.0 * ((-H + (double)blue)/res)); else{ if(purple>h0 && H0>LIMIT_PURPLE){ res = (double)(purple - LIMIT_PURPLE); x = (int)( (255.0 * ((-H + (double)purple)/res))); else{ x = 0; return x; // green int greenjudg(int H0){ double grs = 0; double H = (double)h0; int x; if(red>=h0 && H0>YELLOW){ grs = (double)(red - YELLOW); x = (int)(255 * ((-H + (double)red)/grs)); else{ if(yellow>=h0 && H0>CYAN){ x = 255; else{ if(cyan>=h0 && H0>BLUE){ grs = (double)(cyan - BLUE); x = (int)(255 - (255 * ((-H +(double)cyan)/grs))); else{ x = 0; return x; 88
89 E 3DDOPPLER.C (3DDOPPLERFUNCTION.C) // blue int bluejudg(int H0){ double bls = 0; double H = (double)h0; int x; if(green>=h0 && H0>=CYAN){ bls = (double)(green - CYAN); x = (int)(255 * ((-H + (double)green)/bls)); else{ if(cyan>h0 && H0>=PURPLE){ x = 255; else{ if(purple>h0 && H0>LIMIT_PURPLE){ bls = (double)(purple - LIMIT_PURPLE); x = (int)(255 - (255 * ((-H +(double)purple)/bls))); else{ x = 0; return x; // RGB // // 1 RGB RGB DopplerSynthesis(double X, double Y, double Z, int lambda[], double Yi[], int sum){ RGB rgb; HSV hsv; int newwavelength; double theta; // rgb.r = sum; theta = acos(z / sqrt(x*x + Y*Y + Z*Z)); if(doppler==true) newwavelength = (double)lambda[0] * (1.0-(BETA*cos(theta))) / sqrt(1.0-(beta*beta)); else newwavelength = (double)lambda[0]; 89
90 E 3DDOPPLER.C (3DDOPPLERFUNCTION.C) // RGB rgb.r = redjudg(newwavelength); rgb.g = greenjudg(newwavelength); rgb.b = bluejudg(newwavelength); hsv = RGBtoHSV(rgb); hsv.v = hsv.v * Yi[0]; rgb = HSVtoRGB(hsv); rgb.r = rgb.r/255.0; rgb.g = rgb.g/255.0; rgb.b = rgb.b/255.0; return rgb; double Lorentz(double l, double vc){ l = l * sqrt(1.0-(vc*vc)); return l; // Deformation z double Deformation(double LL, double Y, double zc, double a, double h, double b, double vc){ double zc = zc, L=LL, y=y; double vt; vt = ( sqrt(l*l+y*y+zc*zc)+vc*(zc+b) - sqrt( pow(sqrt(l*l+y*y+zc*zc)+vc*(zc+b), 2.0) - (1.0-vc*vc)*((L*L+y*y+zc*zc)-((L+a)*(L+a)+(y+h)*(y+h)+(zc+b)*(zc+b))) ) ) * (vc/(1.0-vc*vc)); return zc+b+vt; 90
#include <stdio.h> 2 #include <stdlib.h> 3 #include <GL/glut.h> 4 Program 1 (OpenGL GameSample001) 5 // 6 static bool KeyUpON = false; // 7 sta
1 1. 1 #include 2 #include 3 #include 4 Program 1 (OpenGL GameSample001) 5 // 6 static bool KeyUpON = false; // 7 static bool KeyDownON = false; // 8 static bool KeyLeftON
1 1. Program 1 OpenCV (OpenCV Sample001) 1 /* 2 - > - > - >VC++ 3 ( ) 4 C:\opencv\build\include 5 ( ) 6 C:\opencv\build\x86\vc10\lib 7 - > - > - > - >
1 1. Program 1 OpenCV (OpenCV Sample001) 1 /* 2 - > - > - >VC++ 3 ( ) 4 C:\opencv\build\include 5 ( ) 6 C:\opencv\build\x86\vc10\lib 7 - > - > - > - > 8 (240 O p e n C V ) 9 opencv_core240d.lib 10 opencv_imgproc240d.lib
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
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
/*p7-1-1*/
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 /* e8.c Copyright (c) 2003 by T. HAYASHI and
Microsoft PowerPoint - info_eng3_05ppt.pptx
インタラクティブシステム構築法 第 5 回 OpenGL と GLUT の使い方 (3) 埼玉大学情報システム工学科小林貴訓 シェーディング 光源の設定を有効にする glenable(gl_lighting); // 光源の設定を有効にする glenable(gl_light0); //0 番目の光源を有効にする (8 個まで設定可能 ) 光源の位置 GLfloat light0pos[] = {
Microsoft PowerPoint - 04.pptx
初期化 コールバック関数の登録 glutmainloop() 描画関数 マウス処理関数 キーボード処理関数などの関数ポインタを登録する イベント待ちの無限ループ 再描画? no マウス入力? no キーボード入力? no yes yes yes 描画関数の呼び出し マウス処理関数の呼び出し キーボード処理関数の呼び出し void keyboard(unsigned char key, int x,
$ ls -l $ ls -l -a $ ls -la $ ls -F $ ls <dirname> <dirname> $ cd <dirname> <dirname> $ cd $ pwd $ cat <filename> <filename> $ less <filename> <filena
$ pwd /home1/t0/t0903 / / /home1/t0/t0903 / / /home1/t0/t0903 / /... ~ $ ls $ ls -a $ ls -l $ ls -l -a $ ls -la $ ls -F $ ls $ cd $ cd $ pwd $ cat
3D グラフィックス処理の一般過程 1. 3D グラフィックス処理の一般過程
3. 3D ビューイング 1. 3Dグラフィックス処理の一般過程 2. 射影と射影変換 3. ビューボリュームとクリッピング 4. 陰面処理とデプスバッファ 5. ビューポート変換 6. 3Dグラフィックスを描く 7. モデルビュー変換 3D グラフィックス処理の一般過程 1. 3D グラフィックス処理の一般過程 3D グラフィックス処理の一般過程 1. モデリング変換 座標系の異なる複数のオブジェクトを仮想世界に配置し,
libaux.dvi
AUX OpenGL 1 OpenGL (AUX libaux.a) OpenGL Programming Guide () OpenGL 1 OpenGL OS (API) OS OS OS OpenGL Windows Windows X X OpenGL Programming Guide AUX toolkit AUX OS OpenGL SGI OpenGL OS OpenGL AUX Windows
2 T ax 2 + 2bxy + cy 2 + dx + ey + f = 0 a + b + c > 0 a, b, c A xy ( ) ( ) ( ) ( ) u = u 0 + a cos θ, v = v 0 + b sin θ 0 θ 2π u = u 0 ± a
2 T140073 1 2 ax 2 + 2bxy + cy 2 + dx + ey + f = 0 a + b + c > 0 a, b, c A xy u = u 0 + a cos θ, v = v 0 + b sin θ 0 θ 2π u = u 0 ± a cos θ, v = v 0 + b tan θ π 2 < θ < π 2 u = u 0 + 2pt, v = v 0 + pt
第3章 OpenGL の基礎
3 OpenGL April 11, 2017 1 / 28 3.1 ( ) OpenGL OpenGL 2 / 28 3.2 OpenGL OpenGL OpenGL (Open Graphics Library) Silicon Graphics, Inc. 2 3 API (Application Program Interface) [4] UNIX OS Windows Macintosh
Original : Hello World! (0x0xbfab85e0) Copy : Hello World! (0x0x804a050) fgets mstrcpy malloc mstrcpy (main ) mstrcpy malloc free fgets stream 1 ( \n
2008 3 10 1 mstrcpy char *mstrcpy(const char *src); mstrcpy malloc (main free ) stdio.h fgets char *fgets(char *s, int size, FILE *stream); s size ( ) stream FILE ( man ) 40 ( ) %./a.out String : test
新版明解Ｃ言語 実践編
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,
卒 業 研 究 報 告.PDF
C 13 2 9 1 1-1. 1-2. 2 2-1. 2-2. 2-3. 2-4. 3 3-1. 3-2. 3-3. 3-4. 3-5. 3-5-1. 3-5-2. 3-6. 3-6-1. 3-6-2. 4 5 6 7-1 - 1 1 1-1. 1-2. ++ Lisp Pascal Java Purl HTML Windows - 2-2 2 2-1. 1972 D.M. (Dennis M Ritchie)
第3章 OpenGL の基礎
3 OpenGL April 20, 2012 1 / 23 31 ( ) OpenGL OpenGL 2 / 23 32 OpenGL OpenGL OpenGL (Open Graphics Library) Silicon Graphics, Inc 2 3 API (Application Program Interface) [4] UNIX OS Windows Macintosh CAD
新・明解Ｃ言語 ポインタ完全攻略
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)
£Ã¥×¥í¥°¥é¥ß¥ó¥°ÆþÌç (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
2 : 2008/12/ /01/ G :
2 : 2008/12/08 2008/01/16 075730G : 1 project draw main.cpp 1 draw main.cpp /* * main.cpp * draw * * Created by C-T on 08/12/08. * Copyright 2008 MyCompanyName. All rights reserved. * */ #include
新・明解Ｃ言語 実践編
第 1 章 見 21 1-1 見えないエラー 見 List 1-1 "max2x1.h" a, b max2 List 1-1 chap01/max2x1.h max2 "max2x1.h" #define max2(a, b) ((a) > (b)? (a) : (b)) max2 List 1-2 List 1-2 chap01/max2x1test.c max2 #include
II ( ) prog8-1.c s1542h017%./prog8-1 1 => 35 Hiroshi 2 => 23 Koji 3 => 67 Satoshi 4 => 87 Junko 5 => 64 Ichiro 6 => 89 Mari 7 => 73 D
II 8 2003 11 12 1 6 ( ) prog8-1.c s1542h017%./prog8-1 1 => 35 Hiroshi 2 => 23 Koji 3 => 67 Satoshi 4 => 87 Junko 5 => 64 Ichiro 6 => 89 Mari 7 => 73 Daisuke 8 =>. 73 Daisuke 35 Hiroshi 64 Ichiro 87 Junko
PowerPoint 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
( ) 1 1: 1 #include <s t d i o. h> 2 #include <GL/ g l u t. h> 3 #include <math. h> 4 #include <s t d l i b. h> 5 #include <time. h>
2007 12 5 1 2 2.1 ( ) 1 1: 1 #include 2 #include 3 #include 4 #include 5 #include 6 7 #define H WIN 400 // 8 #define W WIN 300 // 9
演算増幅器
コンピュータグラフィックス 2 前回は GLUT を使った簡単な 2 次元グラフィックスについて習った 今週は以下の項目について 補足していく イベント駆動型プログラムの動作について コンピュータグラフィックスの座標系 イベント駆動型プログラム従来のプログラムとの違いこれまでに学習してきたプログラムは上から下に順次実行され 条件分岐や繰り返し処理によって プログラムの流れ (flow: フロー )
1 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)
r07.dvi
19 7 ( ) 2019.4.20 1 1.1 (data structure ( (dynamic data structure 1 malloc C free C (garbage collection GC C GC(conservative GC 2 1.2 data next p 3 5 7 9 p 3 5 7 9 p 3 5 7 9 1 1: (single linked list 1
ohp07.dvi
19 7 ( ) 2019.4.20 1 (data structure) ( ) (dynamic data structure) 1 malloc C free 1 (static data structure) 2 (2) C (garbage collection GC) C GC(conservative GC) 2 2 conservative GC 3 data next p 3 5
£Ã¥×¥í¥°¥é¥ß¥ó¥°(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
第7章 レンダリング
7 April 11, 2017 1 / 59 7.1 ( ) CG 3 ( ) 2 / 59 7.2 7.2.1 ( ) 3 (rendering) 1 / (hidden line/surface calculation) a (outer normal algorithm) b Z (Z-buffer algorithm) c (scan-line algorithm) 2 (shading)
第7章 レンダリング
7 May 18, 2012 1 / 60 71 ( ) CG 3 ( ) 2 / 60 72 71 ( ) 3 (rendering) 1 / (hidden line/surface calculation) a (outer normal algorithm) b Z (Z-buffer algorithm) c (scan-line algorithm) 2 (shading) a (flat
6 6.1 sound_wav_files flu00.wav.wav 44.1 khz 1/44100 spwave Text with Time spwave t T = N t N 44.1 khz t = 1 sec j t f j {f 0, f 1, f 2,, f N 1
6 6.1 sound_wav_files flu00.wav.wav 44.1 khz 1/44100 spwave Text with Time spwave t T = t 44.1 khz t = 1 sec 44100 j t f j {f 0, f 1, f 2,, f 1 6.2 T {f 0, f 1, f 2,, f 1 T ft) f j = fj t) j = 0, 1, 2,,
2 2 GLUI 2 GLUI 2.1 GLUI GLUI OpenGL OpenGL glut OpenGL glut C++ Linux, Windows (Visual C++, gcc), Macintosh glut glut GUI glut GUI CG glmultmatrix()
1 20 (2) OpenGL+GUI(GLUI) 3DCG TA 2008 10 27 1 OpenGL OpenGL GUI GLUI 1.1 http://www.cyber.t.u-tokyo.ac.jp/~kuni/enshu2008/ 1.2 TA 1.3 2008/12/4( ) EyeToy 2 2 GLUI 2 GLUI 2.1 GLUI GLUI OpenGL OpenGL glut
C
C 1 2 1.1........................... 2 1.2........................ 2 1.3 make................................................ 3 1.4....................................... 5 1.4.1 strip................................................
p = 1, 2, cos 2n + p)πj = cos 2nπj 2n + p)πj, sin = sin 2nπj 7.1) f j = a ) 0 + a p + a n+p cos 2nπj p=1 p=0 1 + ) b n+p p=0 sin 2nπj 1 2 a 0 +
7 7.1 sound_wav_files flu00.wav.wav 44.1 khz 1/44100 spwave Text with Time spwave T > 0 t 44.1 khz t = 1 44100 j t f j {f 0, f 1, f 2,, f 1 = T t 7.2 T {f 0, f 1, f 2,, f 1 T ft) f j = fj t) j = 0, 1,
C による数値計算法入門 ( 第 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
comment.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) /* -----------------------------------------------------------------
謗域･ｭ逕ｨppt
情報工学 2017 年度後期第 5 回 [11 月 1 日 ] 静岡大学 工学研究科機械工学専攻ロボット 計測情報講座創造科学技術大学院情報科学専攻 三浦憲二郎 講義日程 第 6 回 11 月 8 日画像処理パート第 1 回 第 7 回 11 月 15 日 CGパート第 6 回 第 8 回 11 月 22 日 CGパート第 7 回 第 9 回 11 月 29 日 CGパート試験 講義アウトライン [11
double rx[natom], ry[natom], rz[natom]; 原子の座標 速度 力 ポテンシャルエ double vx[natom], vy[natom], vz[natom]; ネルギーを受ける配列を準備 double fx[natom], fy[natom], fz[natom
![double rx[natom], ry[natom], rz[natom]; 原子の座標 速度 力 ポテンシャルエ double vx[natom], vy[natom], vz[natom]; ネルギーを受ける配列を準備 double fx[natom], fy[natom], fz[natom](/thumbs/88/115544059.jpg "double rx[natom], ry[natom], rz[natom]; 原子の座標 速度 力 ポテンシャルエ double vx[natom], vy[natom], vz[natom]; ネルギーを受ける配列を準備 double fx[natom], fy[natom], fz[natom")
GLUI による MD の GUI 化 前提条件 :GLUI のプログラミング環境が整っていること 3 原子の MD コード ( 下図 ) viewer ウィンドウ内のマウス左クリックで MD 開始 右クリックで MD 停止 control パネルは solid/wireframe を切り替えるチェックボタン 球の滑らかさと半径を決める窓 ( スピナー ) オブジェクトを回転 移動 拡大縮小させるコントローラ
: CR (0x0d) LF (0x0a) line separator CR Mac LF UNIX CR+LF MS-DOS WINDOWS Japan Advanced Institute of Science and Technology
I117 8 1 School of Information Science, Japan Advanced Institute of Science and Technology : CR (0x0d) LF (0x0a) line separator CR Mac LF UNIX CR+LF MS-DOS WINDOWS Japan Advanced Institute of Science and
ohp08.dvi
19 8 ( ) 2019.4.20 1 (linked list) ( ) next ( 1) (head) (tail) ( ) top head tail head data next 1: 2 (2) NULL nil ( ) NULL ( NULL ) ( 1 ) (double linked list ) ( 2) 3 (3) head cur tail head cur prev data
double float
2015 3 13 1 2 2 3 2.1.......................... 3 2.2............................. 3 3 4 3.1............................... 4 3.2 double float......................... 5 3.3 main.......................
£Ã¥×¥í¥°¥é¥ß¥ó¥°ÆþÌç (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)(),
r08.dvi
19 8 ( ) 019.4.0 1 1.1 (linked list) ( ) next ( 1) (head) (tail) ( ) top head tail head data next 1: NULL nil ( ) NULL ( NULL ) ( 1 ) (double linked list ) ( ) 1 next 1 prev 1 head cur tail head cur prev
‚æ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
PC Windows 95, Windows 98, Windows NT, Windows 2000, MS-DOS, UNIX CPU
1. 1.1. 1.2. 1 PC Windows 95, Windows 98, Windows NT, Windows 2000, MS-DOS, UNIX CPU 2. 2.1. 2 1 2 C a b N: PC BC c 3C ac b 3 4 a F7 b Y c 6 5 a ctrl+f5) 4 2.2. main 2.3. main 2.4. 3 4 5 6 7 printf printf
実際の株価データを用いたオプション料の計算
2002 2 20 1 1 3 2 3 2.1 : : : : : : : : : : : : : : : : : : : : : : : : : : : : 5 2.1.1 : : : : : : : : : : : : : : : : : : : : 5 2.1.2 : : : : : : : : : : : : : : : : : : : : 6 2.2 : : : : : : : : : :
2 2 2 OpenGL (R,G,B,A) 2.1 OpenGL (x y) width height pixels void glreadpixels(glint x, GLint y, GLsizei width, GLsizei height, GLenum format, GLenum t
1 22 (2) OpenGL+GUI(GLUI) 3DCG TA 2010 10 18 1 OpenGL OpenGL GUI GLUI 1.1 http://www.cyber.t.u-tokyo.ac.jp/~kuni/enshu2010/ 1.2 TA 1.3 2010/12/6( ) USB 2 2 2 OpenGL (R,G,B,A) 2.1 OpenGL (x y) width height
untitled
Q 8 1 8.1 (C++) C++ cin cout 5 C++ 16 6 p.63 8.3 #include 7 showbase noshowbase showpoint noshowpoint 8.3 uppercase 16 nouppercase 16 setfill(int) setprecision(int) setw(int) setbase(int) dec
Microsoft Word - C.....u.K...doc
C uwêííôöðöõ Ð C ÔÖÐÖÕ ÐÊÉÌÊ C ÔÖÐÖÕÊ C ÔÖÐÖÕÊ Ç Ê Æ ~ if eíè ~ for ÒÑÒ ÌÆÊÉÉÊ ~ switch ÉeÍÈ ~ while ÒÑÒ ÊÍÍÔÖÐÖÕÊ ~ 1 C ÔÖÐÖÕ ÐÊÉÌÊ uê~ ÏÒÏÑ Ð ÓÏÖ CUI Ô ÑÊ ÏÒÏÑ ÔÖÐÖÕÎ d ÈÍÉÇÊ ÆÒ Ö ÒÐÑÒ ÊÔÎÏÖÎ d ÉÇÍÊ
[1] #include<stdio.h> main() { printf("hello, world."); return 0; } (G1) int long int float ± ±
![[1] #include<stdio.h> main() { printf(hello, world.); return 0; } (G1) int long int float ± ±](/thumbs/99/143405735.jpg "[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,
tuat1.dvi
( 1 ) http://ist.ksc.kwansei.ac.jp/ tutimura/ 2012 6 23 ( 1 ) 1 / 58 C ( 1 ) 2 / 58 2008 9 2002 2005 T E X ptetex3, ptexlive pt E X UTF-8 xdvi-jp 3 ( 1 ) 3 / 58 ( 1 ) 4 / 58 C,... ( 1 ) 5 / 58 6/23( )
[ 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:
![[ 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:](/thumbs/97/132863283.jpg "[ 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 ) ;
超初心者用
3 1999 10 13 1. 2. hello.c printf( Hello, world! n ); cc hello.c a.out./a.out Hello, world printf( Hello, world! n ); 2 Hello, world printf n printf 3. ( ) int num; num = 100; num 100 100 num int num num
P06.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,
main
14 1. 12 5 main 1.23 3 1.230000 3 1.860867 1 2. 1988 1925 1911 1867 void JPcalendar(int x) 1987 1 64 1 1 1 while(1) Ctrl C void JPcalendar(int x){ if (x > 1988) printf(" %d %d \n", x, x-1988); else if(x
BW BW
Induced Sorting BW 11T2042B 2015 3 23 1 1 1.1................................ 1 1.2................................... 1 2 BW 1 2.1..................................... 2 2.2 BW.................................
1) OOP 2) ( ) 3.2) printf Number3-2.cpp #include <stdio.h> class Number Number(); // ~Number(); // void setnumber(float n); float getnumber();
: : :0757230G :2008/07/18 2008/08/17 1) OOP 2) ( ) 3.2) printf Number3-2.cpp #include class Number Number(); // ~Number(); // void setnumber(float n); float getnumber(); private: float num; ;
ex14.dvi
1,, 0, b (b b 2 b ) n k n = n j b j, (0 n j b 1), n =(n k n k 1...n 1 n 0 ) b, n j j j +1, 0,...,b 1 (digit). b b, n b 1 ñ, ñ = k (b 1 n j )b j b N, n b n, n = b N n, n =ñ+1 b N, n m n + m (mod b N ),
file"a" file"b" fp = fopen("a", "r"); while(fgets(line, BUFSIZ, fp)) {... fclose(fp); fp = fopen("b", "r"); while(fgets(line, BUFSIZ, fp)) {... fclose
I117 9 2 School of Information Science, Japan Advanced Institute of Science and Technology file"a" file"b" fp = fopen("a", "r"); while(fgets(line, BUFSIZ, fp)) {... fclose(fp); fp = fopen("b", "r"); while(fgets(line,
/* drawing function */ function_graph();// drawing graph glflush(); int main(int argc, char **argv ) glutinit( &argc, argv ); glutinitdisplaymode( GLU
OpenGL ( ) #include #include #define PI 3.14159265 void function_graph() int j; float x, y; glbegin( GL_LINE_STRIP );// sine curve by line glcolor3f( 0.0f, 1.0f, 1.0f );// line color
II 3 yacc (2) 2005 : Yacc 0 ~nakai/ipp2 1 C main main 1 NULL NULL for 2 (a) Yacc 2 (b) 2 3 y
II 3 yacc (2) 2005 : Yacc 0 ~nakai/ipp2 1 C 1 6 9 1 main main 1 NULL NULL 1 15 23 25 48 26 30 32 36 38 43 45 47 50 52 for 2 (a) 2 2 1 Yacc 2 (b) 2 3 yytext tmp2 ("") tmp2->next->word tmp2 yytext tmp2->next->word
ohp03.dvi
19 3 ( ) 2019.4.20 CS 1 (comand line arguments) Unix./a.out aa bbb ccc ( ) C main void int main(int argc, char *argv[]) {... 2 (2) argc argv argc ( ) argv (C char ) ( 1) argc 4 argv NULL. / a. o u t \0
Microsoft Word - no206.docx
3.2 双方向リスト 今までのリストは 前から順にたどることしかできませんでした 今度は逆にもたどることができる 双方向リストを扱います この場合は 構造体には次を表すポインタの他に前を表すポインタを持つ ことになります 今回は最初と最後をポインタを使うと取り扱いが面倒になるので 最初 (start) と最後 (end) を どちらとも構造体 ( 値は意味を持たない ) を使うことにします こうすることによって
/* do-while */ #include <stdio.h> #include <math.h> int main(void) double val1, val2, arith_mean, geo_mean; printf( \n ); do printf( ); scanf( %lf, &v
1 http://www7.bpe.es.osaka-u.ac.jp/~kota/classes/jse.html kota@fbs.osaka-u.ac.jp /* do-while */ #include #include int main(void) double val1, val2, arith_mean, geo_mean; printf( \n );
1 level Level swtich ButtonPress ButtonRelease Expose Level
UNIX 4 2D/3D Grahpics,GUI :2-3 - 045708G 045726E 045730C 045735D 045759B 045762B 1 level1 1 11 2 12 4 13 6 14 6 2 Level2 6 21 6 211 swtich 11 212 ButtonPress 11 213 ButtonRelease 12 214 Expose 12 22 12
A/B (2018/10/19) Ver kurino/2018/soft/soft.html A/B
A/B (2018/10/19) Ver. 1.0 kurino@math.cst.nihon-u.ac.jp http://edu-gw2.math.cst.nihon-u.ac.jp/ kurino/2018/soft/soft.html 2018 10 19 A/B 1 2018 10 19 2 1 1 1.1 OHP.................................... 1
準備 計算結果を可視化するために OpenGL を 利用する. 2
2. 2 次元粒子法シミュレーション (+ 少しだけ OpenGL) 茨城大学工学部 教授乾正知 準備 計算結果を可視化するために OpenGL を 利用する. 2 OpenGL 3 次元コンピュータグラフィックス用の標準的なライブラリ. 特に CAD やアート, アニメーション分野 ( ゲーム以外の分野 ) で広く利用されている. OpenGL は仕様がオープンに決められており, 企業から独立した団体が仕様を管理している.
untitled
2004/12/21 2/2 (11/16) DT-MRI (11/30) /OpenGL 12/7 12/14 (12/21) 1/11 (1/18) OpenGL ~ ~ OpenGL Silicon Graphics, OpenGL ~ ~ OpenGL OpenGL Utility Library (GLU) OpenGL. OpenGL. OpenGL Utility Toolkit (GLUT)
mstrcpy char *mstrcpy(const char *src); mstrcpy malloc (main free ) stdio.h fgets char *fgets(char *s, int size, FILE *stream); s size ( )
2008 3 10 1 mstrcpy char *mstrcpy(const char *src); mstrcpy malloc (main free ) stdio.h fgets char *fgets(char *s, int size, FILE *stream); s size ( ) stream FILE ( man ) 40 ( ) %./a.out String : test
Excel ではじめる数値解析 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. このサンプルページの内容は, 初版 1 刷発行時のものです.
Excel ではじめる数値解析 サンプルページ この本の定価 判型などは, 以下の URL からご覧いただけます. http://www.morikita.co.jp/books/mid/009631 このサンプルページの内容は, 初版 1 刷発行時のものです. Excel URL http://www.morikita.co.jp/books/mid/009631 i Microsoft Windows
Microsoft Word - Cプログラミング演習(12)
第 12 回 (7/9) 4. いくつかのトピック (5)main 関数の引数を利用したファイル処理 main 関数は, 起動する環境から引数を受け取ることができる 例えば 次に示すように,main 関数に引数を用いたプログラムを作成する 01 /* sample */ 02 /* main 関数の引数 */ 03 #include 04 05 main(int argc, char
プログラミング方法論 II 第 14,15 回 ( 担当 : 鈴木伸夫 ) 問題 17. x 座標と y 座標をメンバに持つ構造体 Point を作成せよ 但し座標 は double 型とする typedef struct{ (a) x; (b) y; } Point; 問題 18. 問題 17 の
プログラミング方法論 II 第 14,15 回 ( 担当 : 鈴木伸夫 ) 問題 17. x 座標と y 座標をメンバに持つ構造体 Point を作成せよ 但し座標 は double 型とする typedef struct{ (a) x; (b) y; Point; 問題 18. 問題 17 の Point を用いて 2 点の座標を入力するとその 2 点間の距 離を表示するプログラムを作成せよ 平方根は
2009 T
T060061 Wii Visual C++ 2008 Express Edition Visual C++ 2008 Express Edition ++ ++ Paint.net ++ 1 2009 T060061 2 1 4 2 4 2.1 Visual C++ 2008 Express Edition.......................... 4 2.2.....................................
C B
C 095707B 2010 6 8 1 LEVE1 2 1.1 LEVEL 1.1................................................ 2 1.1.1 1................................................ 2 1.1.2 1.2..............................................
プログラミング基礎
C プログラミング 演習 アルゴリズム基礎論 演習 第 10 回 今後の予定 12/22( 月 ) 期末試験 (60 分間 ) 場所 :A1611 時間 :16:20~17:20 課題の最終提出締切 :12/19( 金 ) これ以降の新規提出は評価されない 12/22までに最終状況を提示するので, 提出したのに や になってる人は自分の提出内容や提出先を再確認した上で12/26までに問い合わせること
1 OpenGL OpenGL OpenGL OpenGL
2008 OpenGL 2009 2 27 1 OpenGL 4 1.1 OpenGL.............................. 4 1.2 OpenGL............... 4 1.2.1............... 4 1.2.2............................. 5 2 OpenGL 6 2.1.......................
$ pwd /home1/t0/t0903 / / /home1/t0/t0903 / / /home1/t0/t0903 / /... ~ $ ls $ ls -a
$ pwd /home1/t0/t0903 / / /home1/t0/t0903 / / /home1/t0/t0903 / /... ~ $ ls $ ls -a $ ls -l $ ls -l -a $ ls -la $ ls -F $ ls $ cd $ cd $ pwd $ cat
Ｃ言語によるアルゴリズムとデータ構造
Algorithms and Data Structures in C 4 algorithm List - /* */ #include List - int main(void) { int a, b, c; int max; /* */ Ÿ 3Ÿ 2Ÿ 3 printf(""); printf(""); printf(""); scanf("%d", &a); scanf("%d",
XMPによる並列化実装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
C 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 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](/thumbs/88/114913715.jpg "C 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..............................................
( ) ( ) 30 ( ) 27 [1] p LIFO(last in first out, ) (push) (pup) 1
![( ) ( ) 30 ( ) 27 [1] p LIFO(last in first out, ) (push) (pup) 1](/thumbs/96/128219471.jpg "( ) ( ) 30 ( ) 27 [1] p LIFO(last in first out, ) (push) (pup) 1")
() 2006 2 27 1 10 23 () 30 () 27 [1] p.97252 7 2 2.1 2.1.1 1 LIFO(last in first out, ) (push) (pup) 1 1: 2.1.2 1 List 4-1(p.100) stack[] stack top 1 2 (push) (pop) 1 2 void stack push(double val) val stack
コンピュータ概論
4.1 For Check Point 1. For 2. 4.1.1 For (For) For = To Step (Next) 4.1.1 Next 4.1.1 4.1.2 1 i 10 For Next Cells(i,1) Cells(1, 1) Cells(2, 1) Cells(10, 1) 4.1.2 50 1. 2 1 10 3. 0 360 10 sin() 4.1.2 For
3.1 stdio.h iostream List.2 using namespace std C printf ( ) %d %f %s %d C++ cout cout List.2 Hello World! cout << float a = 1.2f; int b = 3; cout <<
C++ C C++ 1 C++ C++ C C++ C C++? C C++ C *.c *.cpp C cpp VC C++ 2 C++ C++ C++ [1], C++,,1999 [2],,,2001 [3], ( )( ),,2001 [4] B.W. /D.M.,, C,,1989 C Web [5], http://kumei.ne.jp/c_lang/ 3 Hello World Hello
£Ã¥×¥í¥°¥é¥ß¥ó¥°ÆþÌç (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 =
program.dvi
2001.06.19 1 programming semi ver.1.0 2001.06.19 1 GA SA 2 A 2.1 valuename = value value name = valuename # ; Fig. 1 #-----GA parameter popsize = 200 mutation rate = 0.01 crossover rate = 1.0 generation
memo
数理情報工学演習第一 C ( 第 12 回 ) 2016/07/11 DEPARTMENT OF MATHEMATICAL INFORMATICS 1 今日の内容 : ファイルの入出力 コマンドライン引数 2 分探索 クイックソート ( ライブラリ ) 文字列検索 2 ファイル操作の手続き : ファイル操作 ファイルからのデータ読み込み ファイルへのデータ書き出し 基本的な手順 読みこむ / 書き出すファイルを開く
I. Backus-Naur BNF S + S S * S S x S +, *, x BNF S (parse tree) : * x + x x S * S x + S S S x x (1) * x x * x (2) * + x x x (3) + x * x + x x (4) * *
2015 2015 07 30 10:30 12:00 I. I VI II. III. IV. a d V. VI. 80 100 60 1 I. Backus-Naur BNF S + S S * S S x S +, *, x BNF S (parse tree) : * x + x x S * S x + S S S x x (1) * x x * x (2) * + x x x (3) +
x h = (b a)/n [x i, x i+1 ] = [a+i h, a+ (i + 1) h] A(x i ) A(x i ) = h 2 {f(x i) + f(x i+1 ) = h {f(a + i h) + f(a + (i + 1) h), (2) 2 a b n A(x i )
![x h = (b a)/n [x i, x i+1 ] = [a+i h, a+ (i + 1) h] A(x i ) A(x i ) = h 2 {f(x i) + f(x i+1 ) = h {f(a + i h) + f(a + (i + 1) h), (2) 2 a b n A(x i )](/thumbs/92/107866678.jpg "x h = (b a)/n [x i, x i+1 ] = [a+i h, a+ (i + 1) h] A(x i ) A(x i ) = h 2 {f(x i) + f(x i+1 ) = h {f(a + i h) + f(a + (i + 1) h), (2) 2 a b n A(x i )")
1 f(x) a b f(x)dx = n A(x i ) (1) ix [a, b] n i A(x i ) x i 1 f(x) [a, b] n h = (b a)/n y h = (b-a)/n y = f (x) h h a a+h a+2h a+(n-1)h b x 1: 1 x h = (b a)/n [x i, x i+1 ] = [a+i h, a+ (i + 1) h] A(x
r03.dvi
19 ( ) 019.4.0 CS 1 (comand line arguments) Unix./a.out aa bbb ccc ( ) C main void... argc argv argc ( ) argv (C char ) ( 1) argc 4 argv NULL. / a. o u t \0 a a \0 b b b \0 c c c \0 1: // argdemo1.c ---
1 return main() { main main C 1 戻り値の型 関数名 引数 関数ブロックをあらわす中括弧 main() 関数の定義 int main(void){ printf("hello World!!\n"); return 0; 戻り値 1: main() 2.2 C main
C 2007 5 29 C 1 11 2 2.1 main() 1 FORTRAN C main() main main() main() 1 return 1 1 return main() { main main C 1 戻り値の型 関数名 引数 関数ブロックをあらわす中括弧 main() 関数の定義 int main(void){ printf("hello World!!\n"); return
:30 12:00 I. I VI II. III. IV. a d V. VI
2017 2017 08 03 10:30 12:00 I. I VI II. III. IV. a d V. VI. 80 100 60 1 I. Backus-Naur BNF X [ S ] a S S ; X X X, S [, a, ], ; BNF X (parse tree) (1) [a;a] (2) [[a]] (3) [a;[a]] (4) [[a];a] : [a] X 2 222222
untitled
II 4 Yacc Lex 2005 : 0 1 Yacc 20 Lex 1 20 traverse 1 %% 2 [0-9]+ { yylval.val = atoi((char*)yytext); return NUM; 3 "+" { return + ; 4 "*" { return * ; 5 "-" { return - ; 6 "/" { return / ; 7 [ \t] { /*
演算増幅器
スペースインベーダーもどき 1000 行プログラムの参考として スペースインベーダーもどきのプログラムを配布する いくつか習って いないものもあるので 補足の説明を加えていく 文字列の描画 文字の描画は glutbitmapcharacter() を用いる これは以下のようにして利用する int i; char *str = "Display String"; glcolor3f(0.0, 0.0,
x(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)
lexex.dvi
(2018, c ) http://istksckwanseiacjp/ ishiura/cpl/ 4 41 1 mini-c lexc,, 2 testlexc, lexc mini-c 1 ( ) mini-c ( ) (int, char, if, else, while, return 6 ) ( ) (+, -, *, /, %, &, =, ==,!=, >, >=,
P02.ppt
int If 2 1 ,,, 3 a + b ab a - b ab a * b ab a / b ab a % b ab a + b 4 2 list0201.c /, % /*/ int vx, vy; puts(""); printf("vx"); scanf("%d", &vx); printf("vy"); scanf("%d", &vy); printf("vx + vy = %d\n",
2 1 list0415.c forfor #include int i, j; for (i = 1; i
A04-133 21 2 9 1 3 1.1.......................................... 3 1.2..................................... 3 1.3................................... 4 2 5 2.1.................................... 5 2.2...............................
2004 2005 2 2 1G01P038-0 1 2 1.1.............................. 2 1.2......................... 2 1.3......................... 3 2 4 2.1............................ 4 2.2....................... 4 2.3.......................
ex01.dvi
,. 0. 0.0. C () /******************************* * $Id: ex_0_0.c,v.2 2006-04-0 3:37:00+09 naito Exp $ * * 0. 0.0 *******************************/ #include int main(int argc, char **argv) double
