pow(readfile.t[k],j)请教

aaaaa5aa

最好将这个N宏定义

light511

可是我用这对参数构造,还是不行
XSJZA xsjzA(36,5);
单步调试第一次赋值时就报那个错了

[ 本帖最后由 light511 于 2009-7-8 19:04 编辑 ]

openspace

m=readfile.m

有readfile.p[k]
for(k=0;k<m;k++)
看看这个有没有越界

light511

这里没有越界,p[100],t[100],m=33

1. for(k=0;k<m;k++)
   
2.     A[i][j]+=pow(readfile.t[k],j);

复制代码

第一次赋值时i=0,j=0,k=0,就报那个错了.

[ 本帖最后由 light511 于 2009-7-8 19:15 编辑 ]

openspace

晕了

原型：extern float pow(float x, float y);
　　用法：#include <math.h>
　　功能：计算x的y次幂。

5^5  >> 36

light511

汗,这个地方是错了,我该了(n+1)*(n+1),但是那个错还是
第一次赋值时i=0,j=0,k=0,就报那个错了.
是不是我这个二维数组A [j]定义的有问题,好象这个数组不能被赋值

[ 本帖最后由 light511 于 2009-7-8 19:30 编辑 ]

openspace

一会半会儿看不出来

你把所有代码都贴出来
我跑一下看看

light511

代码如下

1. 
   
2. #include<stdio.h>
   
3. #include<string.h>
   
4. #include<math.h>
   
5. class ReadFile{
   
6.     public:
   
7.               int m,mP;
   
8.               const int sz;
   
9.               float *p;
   
10.               float *t;
    
11.                       int *P;
    
12.               ReadFile(int a):sz(a){
    
13.                       m=0;
    
14.                       mP=0;
    
15.                       p=new float[sz];
    
16.                       t=new float[sz];
    
17.                       P=new int[sz];
    
18.               }
    
19.               ~ReadFile(){
    
20.                       delete []p;
    
21.                       delete []t;
    
22.                       delete []P;
    
23.                       };
    
24.         bool readfile(char *);
    
25.         bool readfileP(char *);
    
26.         float func_16to10(char *);
    
27.       };
    
28. class XSJZA{
    
29.         public:
    
30.                 const int sz;
    
31.                 const int n;
    
32.                 double ** A;
    
33.                 bool createA(ReadFile& readfile);
    
34.                 XSJZA(int a,int b):sz(a),n(b){
    
35.                         A= new double*[sz];
    
36.                 }
    
37.                 ~XSJZA(){
    
38.                         delete []A;
    
39.                 }
    
40.         };
    
41. bool ReadFile::readfile(char * pathname)
    
42. {
    
43.         int i=0;
    
44.         char tmpp[7],tmpt[7];
    
45.         FILE * cfPtr;
    
46.     printf("pathname=%s\n",pathname);
    
47.         if((cfPtr=fopen(pathname,"r"))==NULL)
    
48.         {
    
49.                 printf("file can not be opened!\n");
    
50.                         return false;
    
51.         }
    
52.         else while(!feof(cfPtr))
    
53.         {
    
54.                 fscanf(cfPtr,"%s%s",&tmpp,&tmpt);
    
55.                 p[m]=func_16to10(tmpp);
    
56.                 t[m]=func_16to10(tmpt);
    
57.                 m++;
    
58.         }
    
59.         fclose(cfPtr);
    
60.         return true;
    
61. }
    
62. float ReadFile::func_16to10(char * str1)
    
63. {
    
64.         int i=0;
    
65.         int len=strlen(str1);
    
66.         float pt=0.0;
    
67. 
    
68.         for(i=len-1;i>=0;i--)
    
69.             if(*str1>=48 && *str1<=57)
    
70.                      pt+=(*str1++-48)*pow(16,i);
    
71.             else if(*str1>=65 && *str1<=70)
    
72.                      pt+=(*str1++-55)*pow(16,i);
    
73.                  else printf("你的数据16进制表示错误");
    
74.         return pt;
    
75. }
    
76. bool ReadFile::readfileP(char * pathname)
    
77. {
    
78.         int tmpP;
    
79.         FILE * cfPtr;
    
80.    
    
81.         if((cfPtr=fopen(pathname,"r"))==NULL)
    
82.         {
    
83.         //        MessageBoxA(NULL,"File can not be opened!","Filenotopen",MB_OK);
    
84.                         return false;
    
85.         }
    
86.         else while(!feof(cfPtr))
    
87.         {
    
88.                 fscanf(cfPtr,"%d",&tmpP);
    
89.                 P[mP]=tmpP;
    
90.                 mP++;
    
91.         }
    
92.         fclose(cfPtr);
    
93.         return true;
    
94. }
    
95. bool XSJZA::createA(ReadFile& readfile)
    
96. {
    
97.         int i=0,j=0,k=0,m=readfile.m;
    
98. 
    
99.     if(0==i)
    
100.         {
     
101.                 for(j=0;j<(n+1)*(n+1);j++)
     
102.                         if(j>=0 && j<=n)
     
103.                                 for(k=0;k<m;k++)
     
104.                                 //        A[i][j]+=pow(readfile.t[k],j);
     
105.                                 A[i][j]+=1.0;
     
106.                                 else if((j>=n+1) && (j<=2*n+1))
     
107.                                         for(k=0;k<m;k++)
     
108.                                                 A[i][j]+=pow(readfile.t[k],(j-(n+1)))*readfile.p[k];
     
109.                                         else if((j>=2*(n+1)) && (j<=3*n+2))
     
110.                                                 for(k=0;k<m;k++)
     
111.                                                         A[i][j]+=pow(readfile.t[k],(j-2*(n+1)))*readfile.p[k]*readfile.p[k];
     
112.                                                 else if((j>=3*(n+1)) && (j<=4*n+3))
     
113.                                                         for(k=0;k<m;k++)
     
114.                                                                 A[i][j]+=pow(readfile.t[k],(j-3*(n+1)))*readfile.p[k]*readfile.p[k]*readfile.p[k];
     
115.                                                         else if((j>=4*(n+1)) && (j<=5*n+4))
     
116.                                                                 for(k=0;k<m;k++)
     
117.                                                                         A[i][j]+=pow(readfile.t[k],(j-4*(n+1)))*readfile.p[k]*readfile.p[k]*readfile.p[k]*readfile.p[k];
     
118.                                                                 else if((j>=5*(n+1)) && (j<=6*n+5))
     
119.                                                                         for(k=0;k<m;k++)
     
120.                                                                                 A[i][j]+=pow(readfile.t[k],(j-5*(n+1)))*readfile.p[k]*readfile.p[k]*readfile.p[k]*readfile.p[k]*readfile.p[k];
     
121.         }
     
122.         for(i=1;i<(n+1)*(n+1);i++)
     
123.         {       
     
124.                 for(j=0;j<(n+1)*(n+1);j++)
     
125.                         if(j>=0 && j<=n)
     
126.                                 for(k=0;k<m;k++)
     
127.                                         A[i][j]+=pow(readfile.t[k],j)*A[0][i];
     
128.                 else if((j>=n+1) && (j<=2*n+1))
     
129.                     for(k=0;k<m;k++)
     
130.                         A[i][j]+=pow(readfile.t[k],(j-(n+1)))*readfile.p[k]*A[0][i];
     
131.                                         else if((j>=2*(n+1)) && (j<=3*n+2))
     
132.                                                 for(k=0;k<m;k++)
     
133.                                                         A[i][j]+=pow(readfile.t[k],(j-2*(n+1)))*readfile.p[k]*readfile.p[k]*A[0][i];
     
134.                                                 else if((j>=3*(n+1)) && (j<=4*n+3))
     
135.                                                         for(k=0;k<m;k++)
     
136.                                                                 A[i][j]+=pow(readfile.t[k],(j-3*(n+1)))*readfile.p[k]*readfile.p[k]*readfile.p[k]*A[0][i];
     
137.                                                         else if((j>=4*(n+1)) && (j<=5*n+4))
     
138.                                                                 for(k=0;k<m;k++)
     
139.                                                                         A[i][j]+=pow(readfile.t[k],(j-4*(n+1)))*readfile.p[k]*readfile.p[k]*readfile.p[k]*readfile.p[k]*A[0][i];
     
140.                                                                 else if((j>=5*(n+1)) && (j<=6*n+5))
     
141.                                                                         for(k=0;k<m;k++)
     
142.                                                                                 A[i][j]+=pow(readfile.t[k],(j-5*(n+1)))*readfile.p[k]*readfile.p[k]*readfile.p[k]*readfile.p[k]*readfile.p[k]*A[0][i];
     
143.         }
     
144.        
     
145.         for(i=0;i<(n+1)*(n+1);i++)
     
146.         {
     
147.                 for(j=0;j<(n+1)*(n+1);j++)
     
148.                         printf("%f\t",A[i][j]);
     
149.                 printf("\n");
     
150.         }
     
151.         return true;       
     
152. }
     
153. int main()
     
154. {
     
155.         ReadFile readfile(100);
     
156.         ReadFile & Readfile=readfile;
     
157.         printf("readfile.sz=%d\n",readfile.sz);
     
158.         readfile.readfile("F:\\test\\PRODUCT\\p-t.txt");
     
159.         readfile.readfileP("F:\\test\\PRODUCT\\P.txt");
     
160.         for(int i=0;i<readfile.m;i++)
     
161.                 printf("%f\t%f\t%d\n",readfile.p[i],readfile.t[i],readfile.P[i]);
     
162.         printf("readfile.m=%d \t readfile.mP=%d\n",readfile.m,readfile.mP);
     
163.        
     
164.         XSJZA xsjzA(36,5);
     
165.         xsjzA.createA(Readfile);
     
166.         return 0;
     
167. }
     
168. 
     

复制代码

需要这两个文件,我再给你
F:\\test\\PRODUCT\\p-t.txt
F:\\test\\PRODUCT\\P.txt

light511

F:\\test\\PRODUCT\\p-t.txt内容如下:
919A16 1018
919A1C 1017
919A04 1018
919A55 1017
919A64 1018
919A48 1017
919A29 1018
919A33 1017
919A71 1018
919A5C 1017
919A69 1018
919A7B 1017
919A43 1018
919A61 1017
9199EF 1018
9199F5 1017
919A01 1017
9199E3 1017
9199C4 1018
9199B7 1017
91998D 1018
9199A7 1017
9199A2 1018
91997C 1017
919953 1018
919951 1017
91996E 1018
919962 1017
91995D 1018
919960 1017
919949 1018
919964 1018
919956 1018
F:\\test\\PRODUCT\\P.txt内容如下:
918
917
918
917
918
917
918
917
918
917
918
917
918
917
918
917
917
917
918
917
918
917
918
917
918
917
918
917
918
917
918
918
918

openspace

double ** A;
      

                XSJZA(int a,int b):sz(a),n(b){
                        A= new double*[sz];         // 这里只创建了一维，没有二维
                }

应该这样
A= new double*[sz];
for (int i = 0; i < sz; i++)
   A = new double[sz];

另外，new好像不初始化基本类型的数组元素，具体忘了