http://w.atwiki.jp/ferret_modoki/ 工房ferret ja 2009-06-04T12:35:42+09:00 1244086542 https://w.atwiki.jp/ferret_modoki/pages/81.html fb 2009-06-04T12:35:42+09:00 1244086542 https://w.atwiki.jp/ferret_modoki/pages/80.html i 2009-06-04T12:33:07+09:00 1244086387 https://w.atwiki.jp/ferret_modoki/pages/79.html t 2009-06-04T12:29:38+09:00 1244086178 https://w.atwiki.jp/ferret_modoki/pages/78.html i 2009-06-04T12:28:03+09:00 1244086083 https://w.atwiki.jp/ferret_modoki/pages/77.html #include <stdio.h> #include <stdlib.h> #include <math.h> main() { int i,j,k,a,b; double m,n,h,h2,o; static double L[700][701]; printf("分割数n="); scanf("%lf",&n); h = 1 / n; h2 = h*h-2; a = n+1; b = n+2; printf("n=%lf, h=%lf, h^2-2=%lf, a=%d, b=%d\n\n",n,h,h2,a,b); /*行列入力*/ L[0][0] = 1; L[n][n] = 1; for(i=1; i<a-1; i++){ L[i][i] = h2; if(i != 1){ L[i][i-1] = 1; } if(i != a-2){ L[i][i+1] = 1; } } for(i=1; i<n; i++){ L[i][a] = -1*h*h*h*i; } /*掃き出し実行*/ for(k=0; k<a; k++){ for(i=0; i<a; i++){ if(i != k){ m = L[i][k]; for(j=0; j<b; j++){ L[i][j] = L[i][j] - ((m / L[k][k]) * L[k][j]); } } } m = L[k][k]; for(j=k; j<b; j++){ L[k][j] = L[k][j] / m; } } /* 出力 */ for(i=0;i<a;i++){ m = (i*h); if( m==.2 || m==.4 || m==.6 || m==.8 || (.599999<m && m<.600001)) printf("u%lf = %lf \n",m,L[i][a]); } return 0; } /* 10 u0.200000 = 0.036132 u0.400000 = 0.062842 u0.600000 = 0.071084 u0.800000 = 0.052550 20 u0.200000 = 0.036106 u0.400000 = 0.062798 u0.600000 = 0.071035 u0.800000 = 0.052514 40 u0.200000 = 0.036100 u0.400000 = 0.062787 u0.600000 = 0.071022 u0.800000 = 0.052505 80 u0.200000 = 0.036098 u0.400000 = 0.062784 u0.600000 = 0.071019 u0.800000 = 0.052503 160 u0.200000 = 0.036098 u0.400000 = 0.062783 u0.600000 = 0.071019 u0.800000 = 0.052503 320 u0.200000 = 0.036098 u0.400000 = 0.062783 u0.600000 = 0.071018 u0.800000 = 0.052503 640 u0.200000 = 0.036098 u0.400000 = 0.062783 u0.600000 = 0.071018 u0.800000 = 0.052502 参考文献：　http://ufcpp.net/study/linear/determinant.html#sweeping */ 2009-05-29T14:55:45+09:00 1243576545 https://w.atwiki.jp/ferret_modoki/pages/76.html -[[掃き出しプログラム]] 2009-05-29T14:55:36+09:00 1243576536 https://w.atwiki.jp/ferret_modoki/pages/68.html 各授業に関するページ -[[画像情報処理]] -[[計算力学]] 2009-05-29T14:55:12+09:00 1243576512 https://w.atwiki.jp/ferret_modoki/pages/75.html /* ---------------------------------------- Image Processing pg1.c Image Ratio Processing / NDVI Gray Image Output (c)Ichinoseki National College Advanced Course ------------------------------------------ */ #include <stdio.h> #include <stdlib.h> #include <math.h> main(argc,argv) int argc; char *argv[]; { static int i,j,m; static double a1,a2,a3; static int pixel,line; FILE *fp1,*fp2,*fp3; /* Command parameter check */ if (argc!=6){ printf("Usage : command input_file1 input_file2 pixel line output_file\n"); return 1; } /* Image size parameter load */ sscanf(argv[3],"%d",&pixel); sscanf(argv[4],"%d",&line); /* Input image file1 check */ if (( fp1=fopen(argv[1],"r") )==NULL){ printf("Can't open input_file\n"); return 1; } /* Input image file1 check */ if (( fp2=fopen(argv[2],"r") )==NULL){ printf("Can't open input_file\n"); return 1; } /* Output image file check */ if (( fp3=fopen(argv[5],"w") )==NULL){ printf("Can't open output_file\n"); return 1; } /* Skip header data of Input Image */ fseek(fp1,24,0); fseek(fp2,24,0); /* ppm header output */ fprintf(fp3,"P6\n"); fprintf(fp3,"# INCT\n"); fprintf(fp3,"%s %s\n",argv[3],argv[4]); fprintf(fp3,"255\n"); /* NDVI processing */ for (j=0; j<line; j++) { for (i=0; i<pixel; i++) { a1=(double)getc(fp1); a2=(double)getc(fp2); a3=(a1-a2)/(a1+a2); m=(0.2+a3)*212.5; /*NDVI値が-0.2以下ならば黒を出力。それ以外ならPseudo colorを出力。*/ if(a3<-0.2){ putc(0,fp3); putc(0,fp3); putc(0,fp3); } else{ putc(127.0*(1.0+cos(m*9.5/255.0)),fp3); putc(255.0*cos(m*1.5/255.0),fp3); putc(255.0*(1-a1/255.0),fp3); } } } fclose(fp1); fclose(fp2); fclose(fp3); return 0; } 2009-05-28T12:00:11+09:00 1243479611 https://w.atwiki.jp/ferret_modoki/pages/74.html /* ---------------------------------------- Image Processing pg1.c Image Ratio Processing / NDVI Gray Image Output (c)Ichinoseki National College Advanced Course ------------------------------------------ */ #include <stdio.h> #include <stdlib.h> #include <math.h> main(argc,argv) int argc; char *argv[]; { static int i,j,m; static double a1,a2,a3; static int pixel,line; FILE *fp1,*fp2,*fp3; /* Command parameter check */ if (argc!=6){ printf("Usage : command input_file1 input_file2 pixel line output_file\n"); return 1; } /* Image size parameter load */ sscanf(argv[3],"%d",&pixel); sscanf(argv[4],"%d",&line); /* Input image file1 check */ if (( fp1=fopen(argv[1],"r") )==NULL){ printf("Can't open input_file\n"); return 1; } /* Input image file1 check */ if (( fp2=fopen(argv[2],"r") )==NULL){ printf("Can't open input_file\n"); return 1; } /* Output image file check */ if (( fp3=fopen(argv[5],"w") )==NULL){ printf("Can't open output_file\n"); return 1; } /* Skip header data of Input Image */ fseek(fp1,24,0); fseek(fp2,24,0); /* ppm header output */ fprintf(fp3,"P6\n"); fprintf(fp3,"# INCT\n"); fprintf(fp3,"%s %s\n",argv[3],argv[4]); fprintf(fp3,"255\n"); /* NDVI processing */ for (j=0; j<line; j++) { for (i=0; i<pixel; i++) { a1=(double)getc(fp1); a2=(double)getc(fp2); a3=(a1-a2)/(a1+a2); m=(1+a3)*128; /*NDVI値が0.3以上ならば赤を出力。それ以外ならグレー画像を出力。*/ if(a3>0.3){ putc(255,fp3); putc(0,fp3); putc(0,fp3); } else{ putc(m,fp3); putc(m,fp3); putc(m,fp3); } } } fclose(fp1); fclose(fp2); fclose(fp3); return 0; } 2009-05-28T11:59:41+09:00 1243479581 https://w.atwiki.jp/ferret_modoki/pages/73.html /* ---------------------------------------- Image Processing pg1.c Image Ratio Processing / NDVI Gray Image Output (c)Ichinoseki National College Advanced Course ------------------------------------------ */ #include <stdio.h> #include <stdlib.h> #include <math.h> main(argc,argv) int argc; char *argv[]; { static int i,j,m; static double a1,a2,a3; static int pixel,line; FILE *fp1,*fp2,*fp3; /* Command parameter check */ if (argc!=6){ printf("Usage : command input_file1 input_file2 pixel line output_file\n"); return 1; } /* Image size parameter load */ sscanf(argv[3],"%d",&pixel); sscanf(argv[4],"%d",&line); /* Input image file1 check */ if (( fp1=fopen(argv[1],"r") )==NULL){ printf("Can't open input_file\n"); return 1; } /* Input image file1 check */ if (( fp2=fopen(argv[2],"r") )==NULL){ printf("Can't open input_file\n"); return 1; } /* Output image file check */ if (( fp3=fopen(argv[5],"w") )==NULL){ printf("Can't open output_file\n"); return 1; } /* Skip header data of Input Image */ fseek(fp1,24,0); fseek(fp2,24,0); /* ppm header output */ fprintf(fp3,"P5\n"); fprintf(fp3,"# INCT\n"); fprintf(fp3,"%s %s\n",argv[3],argv[4]); fprintf(fp3,"255\n"); /* NDVI processing */ for (j=0; j<line; j++) { for (i=0; i<pixel; i++) { a1=(double)getc(fp1); a2=(double)getc(fp2); a3=(a1-a2)/(a1+a2); m=(1+a3)*128; putc(m,fp3); } } fclose(fp1); fclose(fp2); fclose(fp3); return 0; } 2009-05-28T11:59:07+09:00 1243479547