#include <stdio.h>
#include <math.h>
struct xuanzhi
{
  double x,y,c,r;
}
xz[4]={{3,10,1,1},{11.5,12.5,1,1},{3,3,1,1},{10,6,1,1};

main()
{
    struct xuanzhi;
    int i;
    double d[4];
    double x0=0.0,y0=0.0,min=0.0,TC=0.0,t=0.0,m=0.0,w=0.0,x1=0.0,y1=0.0,m1=0.0,m2=0.0,w2=0.0,t1=0.0,m3=0.0,w3=0.0;
     for(i=0;i<4;i++)
       {
           m+=xz[i].c*xz[i].r*xz[i].x;
           w+=xz[i].c*xz[i].r;
           m1+=xz[i].c*xz[i].r*xz[i].y;
       }
           x0=m/w;
           y0=m1/w;
     for(i=0;i<4;i++)
       {
            d[i]=sqrt((xz[i].x-x0)*(xz[i].x-x0)+(xz[i].y-y0)*(xz[i].y-y0));
            m3+=xz[i].c*xz[i].r*xz[i].x/d[i];
            t1+=xz[i].c*xz[i].r*xz[i].y/d[i];
            w3+=xz[i].c*xz[i].r/d[i];
        }
            x1=m3/w3;
            y1=t1/w3;
    do
    {   
            x0=x1;
            y0=y1;
       for(i=0;i<4;i++)
       {
            d[i]=sqrt((xz[i].x-x0)*(xz[i].x-x0)+(xz[i].y-y0)*(xz[i].y-y0));
            m2+=xz[i].c*xz[i].r*xz[i].x/d[i];
            t+=xz[i].c*xz[i].r*xz[i].y/d[i];
            w2+=xz[i].c*xz[i].w/d[i];
        }
        x1=m2/w2;
        y1=t/w2;
        m2=0.0;
        t=0.0;
        w2=0.0;
    }while(fabs(x0-x1)-0.1>0&&fabs(y0-y1)-0.1>0);
    printf("x0=%.2lf,y0=%.2lf\n",x0,y0);
    for(i=0;i<4;i++)
    {
        d[i]=sqrt((xz[i].x-x0)*(xz[i].x-x0)+(xz[i].y-y0)*(xz[i].y-y0));
        min=xz[i].c*xz[i].r*d[i];
        TC+=min;
    }
    printf("x1=%.2lf\n  y1=%.2lf\n", x1, y1);
    getch();
}

xz[4]={{3,10,1,1},{11.5,12.5,1,1},{3,3,1,1},{10,6,1,1};
这个是什么东西啊？！
看起来有点怪怪的！

这段代码太乱，结构体后没有分号，xz数组没有类型，并且少了“}”，用得着采用全局变量.

是不是最后的 getch()应该改成putch（）？

拜托楼主给个注释啊……这看的很纠结啊……非要在纸上拿笔算……

其实拿笔纸上算算就好了，正计算中……彻底纠结中……

#include <stdio.h>
#include <math.h>
struct xuanzhi
{
  double x,y,c,r;
}
xz[4]= {{3,10,1,1},{11.5,12.5,1,1},{3,3,1,1},{10,6,1,1}};//////////////少个大括号
void main()                // 最好加反回类型，没有就用void，不加也没问题但vc6会有警告
{
    struct xuanzhi;
    int i;
    double d[4];
    double x0=0.0,y0=0.0,min=0.0,TC=0.0,t=0.0,m=0.0,w=0.0,x1=0.0,y1=0.0,m1=0.0,m2=0.0,w2=0.0,t1=0.0,m3=0.0,w3=0.0;
     for(i=0;i<4;i++)
       {
           m+=xz[i].c*xz[i].r*xz[i].x;
           w+=xz[i].c*xz[i].r;
           m1+=xz[i].c*xz[i].r*xz[i].y;
       }
           x0=m/w;
           y0=m1/w;
     for(i=0;i<4;i++)
       {
            d[i]=sqrt((xz[i].x-x0)*(xz[i].x-x0)+(xz[i].y-y0)*(xz[i].y-y0));
            m3+=xz[i].c*xz[i].r*xz[i].x/d[i];
            t1+=xz[i].c*xz[i].r*xz[i].y/d[i];
            w3+=xz[i].c*xz[i].r/d[i];
        }
            x1=m3/w3;
            y1=t1/w3;
    do
    {  
            x0=x1;
            y0=y1;
       for(i=0;i<4;i++)
       {
            d[i]=sqrt((xz[i].x-x0)*(xz[i].x-x0)+(xz[i].y-y0)*(xz[i].y-y0));
            m2+=xz[i].c*xz[i].r*xz[i].x/d[i];
            t+=xz[i].c*xz[i].r*xz[i].y/d[i];
            w2+=xz[i].c*xz[i].w/d[i];                        //////////// 结构里没有.w成员,不知道你想用的是哪个成员，自己改吧
        }
        x1=m2/w2;
        y1=t/w2;
        m2=0.0;
        t=0.0;
        w2=0.0;
    }while(fabs(x0-x1)-0.1>0&&fabs(y0-y1)-0.1>0);
    printf("x0=%.2lf,y0=%.2lf\n",x0,y0);
    for(i=0;i<4;i++)
    {
        d[i]=sqrt((xz[i].x-x0)*(xz[i].x-x0)+(xz[i].y-y0)*(xz[i].y-y0));
        min=xz[i].c*xz[i].r*d[i];
        TC+=min;
    }
    printf("x1=%.2lf\n  y1=%.2lf\n", x1, y1);
    getch();                                            ////////// 要用这个函数要包函特定的头文件，具体哪个我忘了
}


[ 本帖最后由 bccn250 于 2010-5-28 19:15 编辑 ]

谢谢，我试试

编译器上显示哪里错了的

#include<stdio.h>
#include<math.h>
struct xuanzhi
{
  double x,y,c,r;
}
xz[4]={{3,10,1,1},{11.5,12.5,1,1},{3,3,1,1},{10,6,1,1}};

void main()
{
    struct xuanzhi;
    int i;
    double d[4];
    double x0=0.0,y0=0.0,min=0.0,TC=0.0,t=0.0,m=0.0,w=0.0,x1=0.0,y1=0.0,m1=0.0,m2=0.0,w2=0.0,t1=0.0,m3=0.0,w3=0.0;
     for(i=0;i<4;i++)
       {
           m+=xz[i].c*xz[i].r*xz[i].x;
           w+=xz[i].c*xz[i].r;
           m1+=xz[i].c*xz[i].r*xz[i].y;
       }
           x0=m/w;
           y0=m1/w;
     for(i=0;i<4;i++)
       {
            d[i]=sqrt((xz[i].x-x0)*(xz[i].x-x0)+(xz[i].y-y0)*(xz[i].y-y0));
            m3+=xz[i].c*xz[i].r*xz[i].x/d[i];
            t1+=xz[i].c*xz[i].r*xz[i].y/d[i];
            w3+=xz[i].c*xz[i].r/d[i];
        }
            x1=m3/w3;
            y1=t1/w3;
    do
    {   
            x0=x1;
            y0=y1;
       for(i=0;i<4;i++)
       {
            d[i]=sqrt((xz[i].x-x0)*(xz[i].x-x0)+(xz[i].y-y0)*(xz[i].y-y0));
            m2+=xz[i].c*xz[i].r*xz[i].x/d[i];
            t+=xz[i].c*xz[i].r*xz[i].y/d[i];
            w2+=xz[i].c*xz[i].y/d[i];
        }
        x1=m2/w2;
        y1=t/w2;
        m2=0.0;
        t=0.0;
        w2=0.0;
    }while(fabs(x0-x1)-0.1>0&&fabs(y0-y1)-0.1>0);
    printf("x0=%.2lf,y0=%.2lf\n",x0,y0);
    for(i=0;i<4;i++)
    {
        d[i]=sqrt((xz[i].x-x0)*(xz[i].x-x0)+(xz[i].y-y0)*(xz[i].y-y0));
        min=xz[i].c*xz[i].r*d[i];
        TC+=min;
    }
    printf("x1=%.2lf\n  y1=%.2lf\n", x1, y1);

}