matlab中二维高斯模板的计算

[cpp]  view plain copy

1. #include <iostream>  
2. #include <math.h>  
3. using namespace std;  
4. #define PI 3.1415926  
5. #define n 7  
6. //#define center_value 4  
7. #define sigma  0.84089642  
8. void makegaussian()  
9. {  
10.     int center_x = n/2,center_y = n/2;  
11.       
12.     double result[n][n];  
13.     double distance;  
14.     double temp;  
15.     double h_sum = 0;  
16.     for(int i=0;i<n;i++)  
17.     {  
18.         for(int j=0;j<n;j++)  
19.         {  
20.             distance = (center_x - i) * (center_x - i)+ (center_y - j) * (center_y - j);  
21.               
22.             temp = exp( (0-distance)/(2*sigma*sigma) )/2*PI*sigma*sigma;  
23.             h_sum += temp;  
24.             result[i][j] = temp;  
25.         }  
26.     }  
27.       
28.     cout<<h_sum<<endl;  
29.   
30.     //归一化  
31.     for(int i=0;i<n;i++)  
32.     {  
33.         for(int j=0; j<n; j++)  
34.         {  
35.             result[i][j]/=h_sum;  
36.         }  
37.         cout<<endl;  
38.     }  
39.       
40.       
41.     for(int i=0;i<n;i++)  
42.     {  
43.         for(int j=0; j<n; j++)  
44.         {  
45.             printf("%0.8f  ",result[i][j]);  
46.         }  
47.         cout<<endl;  
48.     }  
49.       
50.   
51. }  
52.   
53. int main()  
54. {  
55.     makegaussian();  
56.     system("pause");  
57. }  

 本文转自 http://blog.csdn.net/tianguokaka/article/details/6318136