|
|
Регистрация | << Правила форума >> | FAQ | Пользователи | Календарь | Поиск | Сообщения за сегодня | Все разделы прочитаны |
|
Опции темы | Поиск в этой теме | Опции просмотра |
#1
|
|||
|
|||
перевод с языка си в делфи
пожалоста помагите
Код:
//Cubic spline interpolation program //when we have two columns of data x and y in input file: // //x0 y0 //x1 y1 //... //xn yn // //and we want to find such function f(x) //where f(xi) = yi //and f(x) is cubic function on every [x_k-1, x_k] segment //and f(x), f'(x), f''(x) are continual //the result is four columns of cubic polinom coefficients #include <math.h> #include <stdio.h> #include <process.h> float *x, *y, *h, *l, *delta, *lambda, *c, *d, *b; int N; char filename[256]; FILE* InFile=NULL; void count_num_lines(){ //count number of lines in input file - number of equations int nelf=0; //non empty line flag do{ nelf = 0; while(fgetc(InFile)!='\n' && !feof(InFile)) nelf=1; if(nelf) N++; }while(!feof(InFile)); N--; } void readmatrix(){ int i=0; //read matrixes a and b from input file for(i=0; i<N+1; i++){ fscanf(InFile, "%f", &x[i]); fscanf(InFile, "%f", &y[i]); } } void allocmatrix(){ //allocate memory for matrixes x = new float[N+1]; y = new float[N+1]; h = new float[N+1]; l = new float[N+1]; delta = new float[N+1]; lambda = new float[N+1]; c = new float[N+1]; d = new float[N+1]; b = new float[N+1]; } void freematrix(){ delete [] x; delete [] y; delete [] h; delete [] l; delete [] delta; delete [] lambda; delete [] c; delete [] d; delete [] b; } void printresult(){ int k=0; printf("\nA[k]\tB[k]\tC[k]\tD[k]\n"); for(k=1; k<=N; k++){ printf("%f\t%f\t%f\t%f\n", y[k], b[k], c[k], d[k]); } } void testresult(){ float start = x[0]; float end = x[N]; float step = (end - start)/20; FILE* OutFile = fopen("test.txt", "wt"); for(float s = start; s<=end; s+= step){ //find k, where s in [x_k-1; x_k] for(int k=1; k<=N; k++){ if(s>=x[k-1] && s<=x[k]){ break; } } float F = y[k] + b[k]*(s-x[k]) + c[k]*pow(s-x[k], 2) + d[k]*pow(s-x[k], 3); fprintf(OutFile, "%f\t%f\n", s, F); } fclose(OutFile); } void cls(){ for(int i=0; i<25; i++) printf("\n"); } void main(){ int k=0; cls(); do{ printf("\nInput filename: "); scanf("%s", filename); InFile = fopen(filename, "rt"); }while(InFile==NULL); count_num_lines(); rewind(InFile); allocmatrix(); readmatrix(); for(k=1; k<=N; k++){ h[k] = x[k] - x[k-1]; if(h[k]==0){ printf("\nError, x[%d]=x[%d]\n", k, k-1); return; } l[k] = (y[k] - y[k-1])/h[k]; } delta[1] = - h[2]/(2*(h[1]+h[2])); lambda[1] = 1.5*(l[2] - l[1])/(h[1]+h[2]); for(k=3; k<=N; k++){ delta[k-1] = - h[k]/(2*h[k-1] + 2*h[k] + h[k-1]*delta[k-2]); lambda[k-1] = (3*l[k] - 3*l[k-1] - h[k-1]*lambda[k-2]) / (2*h[k-1] + 2*h[k] + h[k-1]*delta[k-2]); } c[0] = 0; c[N] = 0; for(k=N; k>=2; k--){ c[k-1] = delta[k-1]*c[k] + lambda[k-1]; } for(k=1; k<=N; k++){ d[k] = (c[k] - c[k-1])/(3*h[k]); b[k] = l[k] + (2*c[k]*h[k] + h[k]*c[k-1])/3; } printresult(); testresult(); freematrix(); } lmikle: пользуемся тегами!!! Последний раз редактировалось lmikle, 09.01.2013 в 04:49. |