那个大师帮我改点程序啊，不长，应该挺简单的但我不懂

qiaoruijie · 发表于 2006-5-22 18:54

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？我要加入

x

原来的程序如下 /*e.g:y'=y-2x/y,x∈[0,0.6]
/*　　y(0)=1 /*使用经典四阶龙格-库塔算法进行高精度求解　 /* y(i＋1)＝yi＋( K1+ 2*K2 +2*K3+ K4)/6 /* K1=h*f(xi,yi) /* K2=h*f(xi+h/2,yi+K1/2) /* K3=h*f(xi+h/2,yi+K2/2) /* K4=h*f(xi+h,yi+K3) */ #include <stdio.h> #include <math.h>
float f(float den,float p0) //要求解的微分方程的右部的函数 e.g: y-2x/y { float rus; // den=w/(W0+sl); // rus=k*A*f/Wp*sqrt(RTp)-(k-1)*..... rus=p0-2*den/p0; return(rus); }
//float fden() //{
//}
 void main() { float x0; //范围上限 int x1; //范围下限 float h; //步长 int n; //计算出的点的个数 float k1,k2,k3,k4; float y[3]; //用于存放计算出的常微分方程数值解 int i=0; int j;
//以下为函数的接口 printf("intput x0:"); scanf("%f",&x0);
 printf("input x1:"); scanf("%f",&x1);
 printf("input y[0]:"); scanf("%f",&y[0]);
 printf("input h:"); scanf("%f",&h);
// 以下为核心程序 n=((x1-x0)/h); n=3; for(j=0;j<n;j++) { 
 k1=h*f(x0,y); //求K1 k2=h*f((x0+h/2),(y+k1/2)); //求K2 k3=h*f((x0+h/2),(y+k2/2)); //求K3 k4=h*f((x0+h),(y+k3)); //求K4 y[i+1]=y+((k1+2*k2+2*k3+k4)/6); //求yi+1 x0+=float(0.2); printf("y[%f]=%f\n",x0,y[i+1]); ++i; } }
 y'是微分的意思，要把原来的方程改成y'=-(400+4000*cos(x))+(3200000-8000000*sin(x))/y,,,x∈[-180°,180°],,,中间要取100个点，也就是h＝3.6°，，，y(0)=16000，度数也可以用实数表示。 那个大师会帮帮我，感激不尽！加我QQ也行157544051

gghhjj · 发表于 2006-8-21 07:09

参考C语言算法程序集吧，四阶龙格库塔法程序是现成的

账号		自动登录	找回密码
密码			我要加入

[应用数学] 那个大师帮我改点程序啊，不长，应该挺简单的但我不懂

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。