[求助]各位大侠帮我看看我的程序有什么问题啊

妍妍惠子 · 发表于 2006-6-6 11:14

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

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

x

我在做卡尔曼滤波器的基本软件。有个求矩阵逆的地方，我用了全选主元高斯约当法。可４阶方阵还能出个结果（但结果不对）；到了５阶，编译时依然没错，就是运行时老出“can't run target cpu”之类到话。我用到是DSP的软件CCS，c6701 simulater map1各位高手帮我看看吧,先谢谢了！！！ #include "stdio.h" #define N 5//矩阵维数 extern int dcinv(double a[],int n);
void main() { int i,j; /*static double x[N*N]={1,0,0,0,0,0, 0,1,0,0,0,0, 0,0,1,0,0,0, 0,0,0,1,0,0, 0,0,0,0,1,0, 0,0,0,0,0,1};*/ static double x[N*N]={0,0,0,1,2, 0,0,0,2,3, 1,1,0,0,0, 0,1,1,0,0, 0,0,1,0,0,}; /*inverse of x=0 0 1 -1 1 0 0 0 1 -1 0 0 0 0 1 -3 2 0 0 0 2 -1 0 0 0*/ /*double x[N*N]={3,-2,0,0, 5,-3,0,0, 0, 0,3,4, 0, 0,1,1};*/ /*inverse of x=-3 2 0 0 -5 3 0 0 0 0 -1 4 0 0 1 -3*/ double y[N*N]; for(i=0;i<N;i++) for(j=0;j<N;j++) y[i*N+j]=x[i*N+j]; printf("begin\n"); i=dcinv(x,N); printf("end\n"); if(i!=0) { printf("mat x is:\n"); for(i=0;i<N;i++) { for(j=0;j<N;j++) printf("%e ",y[i*N+j]); printf("\n"); } printf("\n"); printf("mat x-1 is:\n"); for(i=0;i<N;i++) { for(j=0;j<N;j++) printf("%e ",x[i*N+j]); printf("\n"); } printf("\n"); } }

妍妍惠子 · 发表于 2006-6-6 11:15

这是求逆子函数 #include "stdio.h" #include "math.h" #include "stdlib.h"
int dcinv(double a[],int n) { int *is,*js,i,j,k,l,u,v; double d,p; is=malloc(n*sizeof(int)); js=malloc(n*sizeof(int)); for(k=0;k<=n-1;k++) { d=0.0; for(i=k;i<=n-1;i++) for(j=k;j<=n-1;j++) {l=i*n+j; p=fabs(a[l]); if(p>d){d=p;is[k]=i;js[k]=j;} } /*选出第k行第k列的右下角子阵中绝对值最大 的元素，放入d中，并将此元素的行列号存入 is[k]、js[k]中*/ if(d+1.0==1.0) {free(is);free(js);printf("err**not inv\n"); return(0); } /*若A为奇异矩阵，则求逆失败，返回错误信息*/ if(is[k]!=k) for(j=0;j<=n-1;j++) {u=k*n+j;v=is[k]*n+j; p=a;a=a[v];a[v]=p; }//行交换 if(js[k]!=k) for(i=0;i<=n-1;i++) {u=i*n+k;v=i*n+js[k]; p=a;a=a[v];a[v]=p; }//列交换 /*通过行交换列交换将最大绝对值元素换到主元素位置上*/ l=k*n+k; a[l]=1.0/a[l]; /*a[k][k]=1/a[k][k]*/ for(j=0;j<=n-1;j++) if(j!=k) {u=k*n+j;a=a*a[v];} /*a[k][j]=a[k][j]*a[k][k]*/ for(i=0;i<=n-1;i++) if(i!=k) for(j=0;j<=n-1;j++) if(j!=k) {u=i*n+j; a=a-a[i*n+k]*a[k*n+j]; } /*a[j]=a[j]-a[k]*a[k][j]*/ for(i=0;i<=n-1;i++) if(i!=k) {u=i*n+k;a=-a*a[l]; } /*a[k]=-a[k]*a[k][k]*/ /*for(i=0;i<=n-1;i++) for(j=0;j<=n-1;j++) printf("a[%n][%n]=%e ",i,j,a[i*n+j]); printf("\n");*/ } for(k=n-1;k>=0;k--) {if(js[k]!=k) for(j=0;j<=n-1;j++) {u=k*n+j;v=js[k]*n+j; p=a;a=a[v];a[v]=p; } if(is[k]!=k) for(i=0;i<=n-1;i++) {u=i*n+k;v=i*n+is[k]; p=a;a=a[v];a[v]=p; } } free(is); free(js); return(1); }

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

[C/C++] [求助]各位大侠帮我看看我的程序有什么问题啊

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