鱼C论坛

 找回密码
 立即注册
查看: 1407|回复: 4

报错板之相关问题,谢谢!!

[复制链接]
发表于 2020-6-10 21:30:25 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能^_^

您需要 登录 才可以下载或查看,没有账号?立即注册

x
error: stray '\246' in program|
error: stray '\302' in program|
请问一下,报错板出现这个是哪里出错了呢(翻译我查了,但是我还是不懂)
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复

使用道具 举报

发表于 2020-6-10 21:31:21 | 显示全部楼层
不发代码就想让别人帮你?做梦
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

发表于 2020-6-10 23:52:52 | 显示全部楼层
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复

使用道具 举报

 楼主| 发表于 2020-6-13 01:59:24 | 显示全部楼层

#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#define N 200
#define P 4
int MatrixInvSymmetric(double *a,int n)//求逆
{
    int i,j,k,m;
    double w,g,*b;
    b = (double*)malloc(n*sizeof(double));
    for (k=0; k<n; k++)
    {
        w=a[0];
        if (fabs(w)<1.0e-8)
        {
            free(b);
            return(-2);
        }
        m=n-k-1;
        for (i=1; i<n; i++)
        {
            g=a[i*n];
            b[i]=g/w;
            if (i<=m)
                b[i]=-b[i];
            for (j=1; j<=i; j++)
                a[(i-1)*n+j-1]=a[i*n+j]+g*b[j];
        }
        a[n*n-1]=1.0/w;
        for (i=1; i<n; i++)
            a[(n-1)*n+i-1]=b[i];
    }
    for (i=0; i<n-1; i++)
        for (j=i+1; j<n; j++)
            a[i*n+j]=a[j*n+i];
    free(b);
    return(2);
}

void readBinaryFile(char *filenames, double *arr, int n)
{
    //数据的读取
    FILE *fp;
    if((fp=fopen(filenames,"rb"))==NULL)
    {
        printf("cannot open file!\n");
        exit(0);
    }
    fread(arr,sizeof(double),n,fp);
    fclose(fp);
}

void EstimateBeta(double *β,double *Y,double *X)
{
    double M[N*N]= {0},Q[P*1]= {0};
    int i,j,k;
    for(i=0; i<N; i++)
    {
        for(j=0; j<N; j++)
        {
            for(k=0; k<N; k++)
            {
                M[N*i+j]+=X[i*N+k]*X[j*N+k];
            }
        }
    }
    MatrixInvSymmetric(M,N);
    for(i=0; i<P; i++)
    {
        for(j=0; j<N; j++)
        {
            Q[i]+=X[j*N+i]*Y[j];
        }
    }
    for(i=0; i<P; i++)
    {
        for(j=0; j<P; j++)
        {
            β[i]+=M[i*P+j]*Q[j];
        }
    }
}


int main()
{
    double Y[N], X[N*P],β[P]= {0};
    readBinaryFile("yb.dat",Y,N);
    readBinaryFile("xb.dat",X,N*P);
    EstimateBeta(β,Y,X);
    return 0;
}
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

 楼主| 发表于 2020-6-13 02:21:47 | 显示全部楼层
啊,我懂啦,谢谢大家
想知道小甲鱼最近在做啥?请访问 -> ilovefishc.com
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

小黑屋|手机版|Archiver|鱼C工作室 ( 粤ICP备18085999号-1 | 粤公网安备 44051102000585号)

GMT+8, 2025-1-13 15:49

Powered by Discuz! X3.4

© 2001-2023 Discuz! Team.

快速回复 返回顶部 返回列表