算法及编程语言 - 声振论坛 - 振动,动力学,声学,信号处理,故障诊断

yflll 发表于 2005-12-28 13:56

高手，高手，高高手！

#include<STDIO.H> #include<CONIO.H> #include<MATH.H> main() { int i=1,n=4; double xm; float t=0.0; float f(float x,float t); float Simpson(float(*f)(float,float),float,float,int); for(t=0.0;t<=2.0;i++) { xm=Simpson(f,0.0,t,n); printf("x(%f)=%lf\n",t,xm); xm=0.0; t=t+0.025; n=n*(i+1); }
}
float Simpson(float (*f)(float x,float t),float a,float b,int n) { int k; float s,s1,s2=0.0; float h=(b-a)/n; s1=f(a+h/2,b); for(k=1;k<=n-1;k++) { s1=s1+f(a+k*h,b); s2=s2+f(a+k*h,b); } s=h/6*(f(a,b)+4*s1+2*s2+f(b,b)); return s; }
 float f(float x,float t) { return sin(10*x)*exp(1.25(x-t))*sin(24.968730*(t-x))*(25000/24968.730444); } 各位高手，看看我的程序问题出在哪？ 我是利用辛普森求积公式求解，sin(10*x)*exp(1.25(x-t))*sin(24.968730*(t-x))*(25000/24968.730444);在（0，t）上的积分，t的取值只能计算十五次，就提示我：floating point error;Divide by 0 这是怎么回事？请高手指点！！！

yflll 发表于 2005-12-28 15:34

没人呢？

Silence 发表于 2005-12-28 15:58

不是高手，斗胆一试，建议搂主以后提问用务实一点的方法，呵呵，把想帮忙的同志都吓跑了 #include<STDIO.H> #include<CONIO.H> #include<MATH.H> void main() { int i=1,n=50; double xm; double t=0.0; double f(double x,double t); double Simpson(double(*f)(double,double),double,double,int); for(t=0.0;t<=2.0;i++) { n = 50; xm=Simpson(f,0.0,t,n); printf("%d\n",i); printf("x(%f)=%lf\n",t,xm); xm=0.0; t=t+double(0.025); n=n*(i+1); } } double Simpson(double (*f)(double x,double t),double a,double b,int n) { int k; double s,s1,s2=0.0; double h=(b-a)/n; s1=f(a+h/2,b); for(k=1;k<=n-1;k++) { s1=s1+f(a+k*h,b); s2=s2+f(a+k*h,b); } s=h/6*(f(a,b)+4*s1+2*s2+f(b,b)); return s; }
 double f(double x,double t) { double temp; temp = sin(10*x)*exp(1.25*(x-t)); temp = temp*sin(24.968730*(t-x))*(25000/24968.730444); return temp; } 现在都能算过去了，结果自己找别的数学软件比较吧。

lang473036 发表于 2006-4-30 11:03

我不知道你到底想干什么??但是我把你的程序用VC++运行了一下修正一错误可以运行了啊 #include<STDIO.H> #include<CONIO.H> #include<MATH.H> main() { int i=1,n=4; double xm; float t=0.0; float f(float x,float t); float Simpson(float(*f)(float,float),float,float,int); for(t=0.0;t<=2.0;i++) { xm=Simpson(f,0.0,t,n); printf("x(%f)=%lf\n",t,xm); xm=0.0; t=t+0.025; n=n*(i+1); } return 0; }
float Simpson(float (*f)(float x,float t),float a,float b,int n) { int k; float s,s1,s2=0.0; float h=(b-a)/n; s1=f(a+h/2,b); for(k=1;k<=n-1;k++) { s1=s1+f(a+k*h,b); s2=s2+f(a+k*h,b); } s=h/6*(f(a,b)+4*s1+2*s2+f(b,b)); return s; }
 float f(float x,float t) { return sin(10*x)*exp(1.25*(x-t))*sin(24.968730*(t-x))*(25000/24968.730444); }

页: [1]

声振论坛's Archiver

高手，高手，高高手！