1081阅读
3回复

[求助]请高手帮我看看这个牛顿迭代法求解数值相速度的程序

UID ：43237

注册：2009-10-07
登录：2014-01-15
发帖：55
等级：仿真一级

0楼发表于: 2010-08-21 12:05:44

这个程序时按照葛德彪书上第43页的公式编写的。但是不知道为什么，迭代结果不对。[table=100%][tr][td]//这个程序用牛顿迭代法求解数值色散方程。具体算法请参见葛书43页!
#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#define DELTA_LAMBDA 0.2
#define PI 3.1415926535
double function( const double, const double );
int main()
{
double beta=1.0, beta1=0.0;
double error;
double phi = 1.0;
FILE *fp;

fp = fopen("data.dat","w");
while ( phi <= 360.0 )
{
  error = 1.0;
  while( fabs(beta - beta1) > 1e-6 )
  {
   beta1 = beta - function(beta,phi);
   beta = beta1;
  }
  printf( "%lg %lg\n",phi,1.0/beta );
  fprintf( fp, "%lg %lg\n",phi,1.0/beta );
  phi = phi + 1.0;
  beta=3.0;
  beta1=0.0;
}
fclose(fp);
}
double function( const double beta, const double phi )
{
double result, tmp;
double A, B, C; //定义三个系数；
double a = 0.5;
C = (PI*DELTA_LAMBDA)*(PI*DELTA_LAMBDA);
A = PI*DELTA_LAMBDA*cos(phi);
B = PI*DELTA_LAMBDA*sin(phi);
result = (sin(A*beta)*si ..

未注册仅能浏览部分内容，查看全部内容及附件请先登录或注册