也学C语言，写了个程序有点问题

yjmwxwx

原书里面的一段 ch04.html.zip (19.34 KB, 下载次数: 73)

2017-12-27 04:54 上传

点击文件名下载附件

输入









结果




网页版计算器


程序写的有问题和网页版对应不起来，请问各位老师哪里有问题呢？

1. #include <stdio.h>
   
2. void main() {
   
3.   float shurur [] = {0.3535,0.3535,0.6464,1.0607,0.3535,-1.0607,-1.3535,-0.3535};
   
4.   float shurui [] = {0,0,0,0,0,0,0,0};
   
5.   float shuchur[8];
   
6.   float shuchui[8];
   
7.   float b0r,b0i,b1r,b1i,b2r,b2i,b3r,b3i,b4r,b4i,b5r,b5i,b6r,b6i,b7r,b7i;
   
8.   float d0r,d0i,d1r,d1i,d2r,d2i,d3r,d3i,d4r,d4i,d5r,d5i,d6r,d6i,d7r,d7i;
   
9.   float c5r,c5i,c7r,c7i;
   
10. 
    
11.       /*一*/
    
12.        b0r=shurur[0]+shurur[4];
    
13.        b0i=shurui[0]+shurui[4];
    
14.        b4r=shurur[0]-shurur[4];
    
15.        b4i=shurui[0]-shurui[4];
    
16.        b2r=shurur[2]+shurur[6];
    
17.        b2i=shurui[2]+shurui[6];
    
18.        b6r=shurur[2]-shurur[6];
    
19.        b6i=shurui[2]-shurui[6];
    
20.        b1r=shurur[1]+shurur[5];
    
21.        b1i=shurui[1]+shurui[5];
    
22.        b5r=shurur[1]-shurur[5];
    
23.        b5i=shurui[1]-shurui[5];
    
24.        b3r=shurur[3]+shurur[7];
    
25.        b3i=shurui[3]+shurui[7];
    
26.        b7r=shurur[3]-shurur[7];
    
27.        b7i=shurui[3]-shurui[7];
    
28. 
    
29.        /*二*/
    
30.        d0r=b0r+b2r;
    
31.        d0i=b0i+b2i;
    
32.        d2r=b0r-b2r;
    
33.        d2i=b0i-b2i;
    
34.       
    
35.        d6i=b6r;
    
36.        d4i=b4i+d6i;
    
37.        d4r=b4r;
    
38.        d6r=b4r;
    
39. 
    
40.        d1r=b1r+b3r;
    
41.        d1i=b1i+b3i;
    
42.        d3r=b1r-b3r;
    
43.        d3i=b1i-b3i;
    
44. 
    
45.        d5r=b5r;
    
46.        d5i=b7r;
    
47.        d7r=b5r;
    
48.        d7i=d7r;
    
49.    
    
50.       
    
51.        /*三*/
    
52.        shuchur[0]=d0r+d1r;
    
53.        shuchui[0]=d0i-d1i;
    
54. 
    
55.        shuchur[1]=d4r;
    
56.        d5i=1.414*d5r;
    
57.        shuchui[1]=d4i+d5i;
    
58. 
    
59.        shuchur[2]=d2r;
    
60.        shuchui[2]=d2i+d3r;
    
61. 
    
62. 
    
63.        shuchur[3]=d6r;
    
64.        d7i=-1.414*d7r;
    
65.        shuchui[3]=d6i+d7i;
    
66.        shuchur[3]=d6r;
    
67.       
    
68.        d5i=d5r*1.414;
    
69.        shuchur[5]=d4r;
    
70.        shuchui[5]=d4i-d5i;
    
71. 
    
72.        shuchur[4]=d0r-d1r;
    
73.        shuchui[4]=d0i-d1i;
    
74. 
    
75. 
    
76.        shuchur[6]=d2r;
    
77.        shuchui[6]=d3r;
    
78. 
    
79.        shuchur[7]=d6r;
    
80.        shuchui[7]=d6i-d7i;
    
81.       
    
82.        printf("%f   %f\n",shuchur[0],shuchui[0]);
    
83.        printf("%f   %f\n",shuchur[1],shuchui[1]);
    
84.        printf("%f   %f\n",shuchur[2],shuchui[2]);
    
85.        printf("%f   %f\n",shuchur[3],shuchui[3]);
    
86.        printf("%f   %f\n",shuchur[4],shuchui[4]);
    
87.        printf("%f   %f\n",shuchur[5],shuchui[5]);
    
88.        printf("%f   %f\n",shuchur[6],shuchui[6]);
    
89.        printf("%f   %f\n",shuchur[7],shuchui[7]);
    
90. 
    
91. 
    
92. }
    

复制代码

scu319hy

是对得上的啊，只是精度不够，两者的计算误差不一样而已。你改用double精度可以提高一倍。如果要更高精度的结果，某些编译器可以支持long double。能支持80位或更长的浮点数。

leonliu_2811

建议你边看C语言书边看别人写的代码, 看懂了再自己写.

glff9

C语言太麻烦，不好学呀

zystbbs

感觉好多人研究fft

yjmwxwx

scu319hy 发表于 2017-12-27 10:45
是对得上的啊，只是精度不够，两者的计算误差不一样而已。你改用double精度可以提高一倍。如果要更高精度的 ...

错误找到了，后面两步分写错了，没乘-1

yjmwxwx

这书里怎么不给个范例，自己写这么简单的几行代码都写不对，我又改了改，各位老师麻烦再给看看

1. #include <stdio.h>
   
2. void main() {
   
3.   float shurur[] = {0.3535,0.3535,0.6464,1.0607,0.3535,-1.0607,-1.3535,-0.3535};
   
4.   float b0r,b0i,b1r,b1i,b2r,b2i,b3r,b3i,b4r,b4i,b5r,b5i,b6r,b6i,b7r,b7i;
   
5.   float d0r,d1r,d2r,d3r,d4r,d5r,d6r,d7r;
   
6.   float d0i,d1i,d2i,d3i,d4i,d5i,d6i,d7i;
   
7.   float c0r,c1r,c2r,c3r,c4r,c5r,c6r,c7r,c0i,c1i,c2i,c3i,c4i,c5i,c6i,c7i;
   
8.       /*一*/
   
9.        b0r=shurur[0]+shurur[4];
   
10.        b4r=shurur[0]-shurur[4];
    
11.        b2r=shurur[2]+shurur[6];
    
12.        b6r=shurur[2]-shurur[6];
    
13.        b1r=shurur[1]+shurur[5];
    
14.        b5r=shurur[1]-shurur[5];
    
15.        b3r=shurur[3]+shurur[7];
    
16.        b7r=shurur[3]-shurur[7];
    
17. 
    
18.        /*二*/
    
19.           d0r=b0r+b2r;
    
20.           d0i=0;
    
21.           d2r=b0r-b2r;
    
22.           d2i=0;
    
23. 
    
24.           b6i=b6r*(-1);
    
25.           b4i=0;
    
26.           d4i=b4i+b6i;
    
27.           d4r=b4r;
    
28.           d6r=b4r;
    
29.           d6i=b4i-b6i;
    
30. 
    
31.           d1r=b1r+b3r;
    
32.           d1i=0;
    
33.           d3r=b1r-b3r;
    
34.           d3i=0;
    
35. 
    
36.           b7i=b7r*(-1);
    
37.           d5i=b4i+b7i;
    
38.           d5r=b7r;
    
39.           d7r=b5r;
    
40.           d7i=b4i-b7i;
    
41.       
    
42.        /*三*/
    
43.           c0r=d0r+d1r;
    
44.           c0i=d0i;
    
45.           c4r=d0r-d1r;
    
46.           c4i=d1i;
    
47. 
    
48.           d5i=d5r*1.414;
    
49.           c1r=d4r;
    
50.           c5i=d4i+d5i;
    
51.           c5r=d4r;
    
52.           c1i=d4i-d5i;
    
53. 
    
54.           d3i=d3r*(-1);
    
55.           c2r=d2r;
    
56.           c2i=d2i+d3i;
    
57.           c6i=d2i-d3i;
    
58.           c6r=d2r;
    
59. 
    
60.           d7i=d7r*(-1.4147);
    
61.           c3r=d6r;
    
62.           c3i=d6i+d7i;
    
63.           c7r=d6r;
    
64.           c7i=d6i-d7i;
    
65.           
    
66. 
    
67.               
    
68.        printf("%f         %f \n",c0r,c0i);
    
69.        printf("%f         %f \n",c1r,c1i);
    
70.        printf("%f         %f \n",c2r,c2i);
    
71.        printf("%f         %f \n",c3r,c3i);
    
72.        printf("%f         %f \n",c4r,c4i);
    
73.        printf("%f         %f \n",c5r,c5i);
    
74.        printf("%f         %f \n",c6r,c6i);
    
75.        printf("%f         %f \n",c7r,c7i);
    
76. 
    
77. 
    
78. }
    

复制代码

里面还有错误，但是不知道为什么结果对了

zhangxiyi5277

这是算法呀，，前辈直接在c上搞呀，，感觉好痛苦。

平凡8236

没有拿来就能用的程序！！

yjmwxwx

谢谢各位老师，我找到问题啦，原来是虚数计算的问题！

我小学文化没学过虚数计算

yjmwxwx

我说怎么不能理解（0.707+0.707i）×（1.414+1.414i）=1.999i呢，原来根本没学过。

早知道也好好上上学啦，我上初中就没听过课课桌里全是垃圾没书。

yjmwxwx

又改了下，这次没问题了，谢谢各位老师指导，终于把书里的计算方法写成了程序，虽然写的很烂，但是零基础学个东西不容易。

1. #include<stdio.h>
   
2. void main() {
   
3.   float shurur[] = {0.3535,0.3535,0.6464,1.0607,0.3535,-1.0607,-1.3535,-0.3535};
   
4.   float sr, si,ac,bd,ad,bc;
   
5.   float b0r,b0i,b1r,b1i,b2r,b2i,b3r,b3i,b4r,b4i,b5r,b5i,b6r,b6i,b7r,b7i;
   
6.   float d0r,d1r,d2r,d3r,d4r,d5r,d6r,d7r;
   
7.   float d0i,d1i,d2i,d3i,d4i,d5i,d6i,d7i;
   
8.   float c0r,c1r,c2r,c3r,c4r,c5r,c6r,c7r,c0i,c1i,c2i,c3i,c4i,c5i,c6i,c7i;
   
9.       /*一*/
   
10.        b0r=shurur[0]+shurur[4];
    
11.        b4r=shurur[0]-shurur[4];
    
12.        b2r=shurur[2]+shurur[6];
    
13.        b6r=shurur[2]-shurur[6];
    
14.        b1r=shurur[1]+shurur[5];
    
15.        b5r=shurur[1]-shurur[5];
    
16.        b3r=shurur[3]+shurur[7];
    
17.        b7r=shurur[3]-shurur[7];
    
18. 
    
19.        /*二*/
    
20.        d0r=b0r+b2r;
    
21.        d0i=0;
    
22.        d2r=b0r-b2r;
    
23.        d2i=0;
    
24. 
    
25. 
    
26.        b4i=0;
    
27.        b6i=0;
    
28.        sr=0;
    
29.        si=-1;
    
30.        ac=b6r*sr;
    
31.        bd=b6i*si;
    
32.        ad=b6i*sr;
    
33.        bc=b6r*si;
    
34.        b6r=ac-bd;
    
35.        b6i=ad+bc;
    
36.        d4r=b4r+b6r;
    
37.        d4i=b4i+b6i;
    
38.        d6r=b4r-b6r;
    
39.        d6i=b4i-b6i;
    
40. 
    
41.        d1r=b1r+b3r;
    
42.        d1i=0;
    
43.        d3r=b1r-b3r;
    
44.        d3i=0;
    
45. 
    
46. 
    
47.        b6i=0;
    
48.        b7i=0;
    
49.        sr=0;
    
50.        si=-1;
    
51.        ac=b7r*sr;
    
52.        bd=b7i*si;
    
53.        ad=b7i*sr;
    
54.        bc=b7r*si;
    
55.        b7r=ac-bd;
    
56.        b7i=ad+bc;
    
57.        d5r=b5r+b7r;
    
58.        d5i=b5i+b7i;
    
59.        d7r=b5r-b7r;
    
60.        d7i=b5i-b7i;
    
61. 
    
62.       
    
63.        /*三*/
    
64.          
    
65.          
    
66.           c0r=d0r+d1r;
    
67.           c0i=d0i+d1i;
    
68.           c4r=d0r-d1r;
    
69.           c4i=d0i-d1i;
    
70. 
    
71.          
    
72.           sr=0.7071;
    
73.           si=-0.7071;
    
74.           ac=d5r*sr;
    
75.           bd=d5i*si;
    
76.           ad=d5i*sr;
    
77.           bc=d5r*si;
    
78.           d5r=ac-bd;
    
79.           d5i=ad+bc;
    
80.           c1r=d4r+d5r;
    
81.           c1i=d4i+d5i;
    
82.           c5r=d4r-d5r;
    
83.           c5i=d4i-d5i;
    
84. 
    
85.           sr=0;
    
86.           si=-1;
    
87.           ac=d3r*sr;
    
88.           bd=d3i*si;
    
89.           ad=d3i*sr;
    
90.           bc=d3r*si;
    
91.           d3r=ac-bd;
    
92.           d3i=ad+bc;
    
93.           c2r=d2r+d3r;
    
94.           c2i=d2i+d3i;
    
95.           c6r=d2r-d3r;
    
96.           c6i=d2i-d3i;
    
97. 
    
98.           sr=-0.7071;
    
99.           si=-0.7071;
    
100.           ac=d7r*sr;
     
101.           bd=d7i*si;
     
102.           ad=d7i*sr;
     
103.           bc=d7r*si;
     
104.           d7r=ac-bd;
     
105.           d7i=ad+bc;
     
106.           c3r=d6r+d7r;
     
107.           c3i=d6i+d7i;
     
108.           c7r=d6r-d7r;
     
109.           c7i=d6i-d7i;
     
110.               
     
111.        printf("%f         %f \n",c0r,c0i);
     
112.        printf("%f         %f \n",c1r,c1i);
     
113.        printf("%f         %f \n",c2r,c2i);
     
114.        printf("%f         %f \n",c3r,c3i);
     
115.        printf("%f         %f \n",c4r,c4i);
     
116.        printf("%f         %f \n",c5r,c5i);
     
117.        printf("%f         %f \n",c6r,c6i);
     
118.        printf("%f         %f \n",c7r,c7i);
     
119. 
     
120. 
     
121. }
     

复制代码

cuiningkun

看着头都大了，找到就好了。。。

onetty

慢慢琢磨吧，这个得循序渐进。

袁国钟

C语言太麻烦，不好学呀


呵呵，如果用汇编，可能想都不要想了。C语言只是一些计算公式，只要公式对了，数据类型对了基本就差不多了。我这里说句多余的，有些人把C言看的太那个了，我也说不清了，就像楼主就是一个，我看过他好多帖子，只是在装没知识，其实我也不知道，呵呵，楼主回答一下吧
。