MATLAB與C語言

snowflying

雷亂圈 發表於 2015-12-2 06:36 PM

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

傳進去的是動態配置而來的

1. #include <stdio.h>
   
2. #include <stdlib.h>
   
3. #include <string.h>
   
4. #define print printf
   
5. #define N 3
   
6. #define M 8
   
7. 
   
8. void mutation(double **parent , double pm , int n , int m , double ***child)
   
9. {
   
10.     int i , j;
    
11.     double rd;
    
12.    
    
13.     srand(time(NULL));
    
14.     *child = (double**) malloc(sizeof(double**) * n);
    
15.    
    
16.     for(i = 0 ; i < n ; ++i)
    
17.         (*child)[i] = (double*) malloc(sizeof(double) * m);
    
18.    
    
19.     rd = (double) rand() / RAND_MAX;
    
20.    
    
21.    
    
22.     if(rd < pm)
    
23.     {
    
24.         int *mpoint = (int*) malloc(sizeof(int) * 100);
    
25.         
    
26.         for(i = 0 ; i < 100 ; ++i)
    
27.             mpoint[i] = (double) rand() / RAND_MAX * (m - 1) + 0.5;
    
28. 
    
29.         for(i = 0 ; i < n ; ++i)
    
30.             memcpy((*child)[i] , parent[i] , sizeof(double) * m);
    
31.             
    
32.         for(i = 0 ; i < 100 ; ++i)
    
33.             for(j = 0 ; j < n ; ++j)
    
34.                 (*child)[j][mpoint[i]] = abs(parent[j][mpoint[i]] - 1);
    
35. 
    
36.         free(mpoint);
    
37.     }
    
38.     else
    
39.     {
    
40.         for(i = 0 ; i < n ; ++i)
    
41.             memcpy((*child)[i] , parent[i] , sizeof(double) * m);
    
42.     }
    
43. }
    
44. int main()
    
45. {
    
46.     double **parent;
    
47.     double **child = NULL;
    
48.     double val = 0;
    
49.     int i , j;
    
50.    
    
51.     parent = (double**) malloc(sizeof(double*) * N);
    
52.    
    
53.     for(i = 0 ; i < N ; ++i)
    
54.         parent[i] = (double*) malloc(sizeof(double) * M);
    
55. 
    
56.     for(i = 0 ; i < N ; ++i)
    
57.         for(j = 0 ; j < M ; ++j)
    
58.             parent[i][j] = val++;
    
59. 
    
60.     mutation(parent , 0.3 , N , M , &child);
    
61.    
    
62.     for(i = 0 ; i < N ; ++i)
    
63.     {
    
64.         for(j = 0 ; j < M ; ++j)
    
65.         {
    
66.             print("%f " , child[i][j]);
    
67.         }
    
68.         putchar('\n');
    
69.     }
    
70.    
    
71.     for(i = 0 ; i < N ; ++i)
    
72.         free(parent[i]);
    
73.         
    
74.     free(parent);
    
75. 
    
76.     for(i = 0 ; i < N ; ++i)
    
77.         free(child[i]);
    
78.         
    
79.     free(child);
    
80.         
    
81.     return 0;
    
82. }

複製代碼

...

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

雷亂圈

snowflying 發表於 2015-12-3 02:56 AM

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

版大我將您所打的程式進行運行形成執行檔
但說這一段有錯
srand(time(NULL));

雷亂圈

雷亂圈 發表於 2015-12-3 11:17 AM

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

版大我將您所打的程式進行運行形成執行檔
但說這一段有錯
srand(time(NULL));

版大想多問一下當作學習用
假如直接是個二維陣列的話作法是如何...

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

雷亂圈

function [child1,child2] = crossover(parent1,parent2,pc)
stringlength=length(parent1(1,: ));
    if(rand < pc)
        cpoint =  round(rand*(stringlength-2))+1;
        child1=[parent1(:,1:cpoint) parent2(:,cpoint+1:stringlength)];
        child2=[parent2(:,1:cpoint) parent1(:,cpoint+1:stringlength)];
    else
    child1=parent1;
...

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

snowflying

雷亂圈 發表於 2015-12-3 10:57 PM

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

function [child1,child2] = crossover(parent1,parent2,pc)
stringlength=length(parent1(1,: ));
    if( ...

如果  parent1 和 parent2 維度一樣
...

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

雷亂圈

snowflying 發表於 2015-12-4 12:52 AM

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

如果  parent1 和 parent2 維度一樣
那麼共用同樣的 n、m 就行了

版大請問一下這兩段該如何編寫
cpoint =  round(rand*(stringlength-2))+1;
...

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

雷亂圈

版大
這兩段我翻書找資料幾天了但還是無法編寫出來請版大在幫一次忙吧
cpoint =  round(rand*(stringlength-2))+1;
child1=[parent1(:,1:cpoint) parent2(:,cpoint+1:stringlength)];

雷亂圈

謝謝版大的回答我在研究一下