我有一個Visual C++ 版本的程式碼 要轉成Dev-C++版本

liforlove

chevylin0802 發表於 2013-6-23 12:17 AM [url=forum.php?mod=redirect&goto=findpost&pid=241374516&ptid=8794453][/url]
將三種計算方式改寫成物件，標頭檔改成c++專屬的，存成c++專案，原先的副檔名改為cpp ...

改寫成物件什麼意思??
先改標頭檔 這樣下面的計算方式是不是連double的部分也要改
...

瀏覽完整內容，請先 註冊 或 登入會員

chevylin0802

liforlove 發表於 2013-6-23 01:23 AM

下載: 訪客無法瀏覽下載點，請先 註冊 或 登入會員

改寫成物件什麼意思??
先改標頭檔 這樣下面的計算方式是不是連double的部分也要改
然後do改成for

...

瀏覽完整內容，請先 註冊 或 登入會員

liforlove

chevylin0802 發表於 2013-6-23 12:09 PM

下載: 訪客無法瀏覽下載點，請先 註冊 或 登入會員

暫時的版本
試過應該可以

這一樣有用到副程式嗎?
如果沒有我從int main(int argc, char *argv[])這行下面開始
...

瀏覽完整內容，請先 註冊 或 登入會員

liforlove

liforlove 發表於 2013-6-23 12:54 PM [url=forum.php?mod=redirect&goto=findpost&pid=241404008&ptid=8794453][/url]
這一樣有用到副程式嗎?
如果沒有我從int main(int argc, char *argv[])這行下面開始
但一開始的class就說 ...

那你的using namespace std;
...

瀏覽完整內容，請先 註冊 或 登入會員

chevylin0802

liforlove 發表於 2013-6-23 01:29 PM

下載: 訪客無法瀏覽下載點，請先 註冊 或 登入會員

那你的using namespace std;

跟int main(int argc, char *argv[])這兩行呢?

不一定需要using namespace std

int main() 這段我是照你原本的code 沿用的
...

瀏覽完整內容，請先 註冊 或 登入會員

liforlove

chevylin0802 發表於 2013-6-23 03:06 PM

下載: 訪客無法瀏覽下載點，請先 註冊 或 登入會員

不一定需要using namespace std

int main() 這段我是照你原本的code 沿用的

還是不行欸.....
整段如下

...

瀏覽完整內容，請先 註冊 或 登入會員

chevylin0802

liforlove 發表於 2013-6-23 04:57 PM

下載: 訪客無法瀏覽下載點，請先 註冊 或 登入會員

還是不行欸.....
整段如下

1. 打開新專案, 設為C++專案及Console Application
2. 清除原本的程式碼
3. 將我的程式碼複製並貼上去
...

瀏覽完整內容，請先 註冊 或 登入會員

liforlove

chevylin0802 發表於 2013-6-23 05:00 PM

下載: 訪客無法瀏覽下載點，請先 註冊 或 登入會員

1. 打開新專案, 設為C++專案及Console Application
2. 清除原本的程式碼
3. 將我的程式碼複製並貼上去

痾......不行就是不行
...

瀏覽完整內容，請先 註冊 或 登入會員

liforlove

可以是可以了
只不過答案跟原本的似乎不太一樣...左邊為C++版本  右邊為原版

chevylin0802

liforlove 發表於 2013-6-23 06:19 PM

下載: 訪客無法瀏覽下載點，請先 註冊 或 登入會員

可以是可以了
只不過答案跟原本的似乎不太一樣...左邊為C++版本  右邊為原版

1. #include <cstdlib>
   
2. #include <iostream>
   
3. #include <cmath>
   
4. 
   
5. class NewtonRoot {
   
6.       private:
   
7.               double x0, x1, eps, max_itera;
   
8.       public:
   
9.              NewtonRoot(double x0, double x1, double eps, int max_itera)
   
10.                         { this->x0 = x0; this->x1 = x1; this->eps = eps; this->max_itera = max_itera; };
    
11.              ~NewtonRoot() { };
    
12.              double DiffResult();
    
13.              double DiffmResult();
    
14.              double XnResult();
    
15.              double func(double x);
    
16. };
    
17. 
    
18. double NewtonRoot::func(double x)
    
19. {
    
20.      double t = x-1.235;
    
21.      return pow(t, 5.0);
    
22. }
    
23. 
    
24. // -------------------------------------------------------
    
25. double NewtonRoot::DiffResult(void)
    
26. {
    
27.      double y0= this->func(this->x0), y1;
    
28.      double diff, x2;
    
29.      double x1t, x0t;
    
30.      int i=0;
    
31.      
    
32.      x0t=this->x0;
    
33.      x1t=this->x1;
    
34.      
    
35.      do{
    
36.            y1=this->func(x1t);
    
37.            diff = (y1 - y0) / (x1t - x0t);
    
38.            x2 = x1t - y1 / diff;
    
39.            x0t = x1t;
    
40.            y0 = y1;
    
41.            x1t = x2;
    
42.            ++i;
    
43.      }while(fabs(x1t-x0t)>this->eps && i<this->max_itera);
    
44.      print_f("> times: %4d , ", i);
    
45.      return x2;
    
46. }
    
47. 
    
48. // -------------------------------------------------------
    
49. double NewtonRoot::DiffmResult(void)
    
50. {
    
51.      const int m = 5;
    
52.      double y0 = this->func(this->x0), y1;
    
53.      double diff, x2;
    
54.      double x0t, x1t;     
    
55.      int i=0;
    
56.      
    
57.      x0t = this->x0;
    
58.      x1t = this->x1;
    
59.      
    
60.      do{
    
61.            y1=this->func(x1t);
    
62.            diff = (y1 - y0) / (x1t - x0t);
    
63.            x2 = x1t - m * y1 / diff;
    
64.            x0t = x1t;
    
65.            y0 = y1;
    
66.            x1t = x2;
    
67.            ++i;
    
68.      }while(fabs(x1t-x0t)>this->eps && i<this->max_itera);
    
69.      print_f("> times: %4d , ", i);
    
70.      return x2;
    
71. }
    
72. 
    
73. // -------------------------------------------------------
    
74. double NewtonRoot::XnResult(void)
    
75. {
    
76.      double y0 = this->func(this->x0), y1;
    
77.      double diff, x2, delta;
    
78.      double rate=1.0;
    
79.      double x0t, x1t;
    
80.      int i=0;
    
81. 
    
82.      x0t = this->x0;
    
83.      x1t = this->x1;
    
84.      srand(time(NULL));
    
85. 
    
86.      do{
    
87.            delta = y1 = this->func(x1t);
    
88.            // if(fabs(y1) < fabs(y0)) rate=1.0;
    
89.            // else if(rate > 5E-5) rate*=0.5;
    
90.            rate = 5.0 * rand() / RAND_MAX; // more faster
    
91. 
    
92.            diff = (y1 - y0) / (x1t - x0t);
    
93.            x2 = x1t - rate * y1 / diff;         
    
94.            x0t = x1t; y0 = y1; x1t = x2;
    
95.            ++i;
    
96.      }while(fabs(delta)>this->eps && i<this->max_itera);
    
97.      print_f("> times: %4d , ", i);
    
98.      return x2;
    
99. }
    
100. 
     
101. int main()
     
102. {
     
103.      const double eps=1E-9;
     
104.      const int max_iterator=1000;
     
105.      double x, low, up;
     
106. 
     
107.      low = -10.0, up = low + 1.0;
     
108.      class NewtonRoot *n = new NewtonRoot(low, up, eps, max_iterator);
     
109.      x = n->DiffResult();
     
110.      print_f("NewtonRoot.Diff  func(%+.15e) = %+.15e\n", x, n->func(x));
     
111.      x = n->XnResult();
     
112.      print_f("NewtonRoot.Xn    func(%+.15e) = %+.15e\n", x, n->func(x));
     
113.      x = n->DiffmResult();
     
114.      print_f("NewtonRoot.Diffm func(%+.15e) = %+.15e\n", x, n->func(x));
     
115.      
     
116.      system("PAUSE");
     
117.      
     
118.      delete n;
     
119.      
     
120.      return EXIT_SUCCESS;
     
121. }

複製代碼

我用這個版本試
都很正常耶...

瀏覽完整內容，請先 註冊 或 登入會員

snowflying

liforlove 發表於 2013-6-23 06:19 PM

下載: 訪客無法瀏覽下載點，請先 註冊 或 登入會員

可以是可以了
只不過答案跟原本的似乎不太一樣...左邊為C++版本  右邊為原版

chevylin0802 前輩的 code 我這邊跑也很正常
...

瀏覽完整內容，請先 註冊 或 登入會員

liforlove

snowflying 發表於 2013-6-24 12:36 AM

下載: 訪客無法瀏覽下載點，請先 註冊 或 登入會員

chevylin0802 前輩的 code 我這邊跑也很正常
剛甚至試過砍掉 system("pause"); 並拿到 FreeBSD 上跑，也  ...

是為了數學寫程式!!! 我們只有上課交而且才半學期......
...

瀏覽完整內容，請先 註冊 或 登入會員