#include "stdio.h"
#include "stdlib.h"
#include "math.h"
#include "iostream.h"
void Gauss();
void main()
{
Gauss();
}
void Gauss()
{
int h_size;
float siz,sigma;
int i, j;
printf("Please input size of gaussian core/n");
scanf("%d",&h_size);
printf("Please input sigma:/n");
scanf("%f",&sigma);
siz=(float)(h_size-1)/2;
float **a,**b; a=new float*[h_size];
for(i=0;i<h_size;i++) a[i]=new float[h_size];
b=new float*[h_size];
for( i=0;i<h_size;i++) b[i]=new float[h_size];
for(i=0;i<h_size;i++)
{
for(j=0;j<h_size;j++)
{
a[i][j]=-siz+j;
printf("%4.2f ",a[i][j]);
}
printf("/n");
}
printf("/n");
for( i=0;i<h_size;i++)
{
for(j=0;j<h_size;j++)
{
b[i][j]=a[j][i];
printf("%4.2f ",b[i][j]);
}
printf("/n");
}
printf("/n");
float h_sum=0;
for( i=0;i<h_size;i++)
{
for(j=0;j<h_size;j++)
{
a[i][j]=a[i][j]*a[i][j];
b[i][j]=b[i][j]*b[i][j];
a[i][j]=-(a[i][j]+b[i][j])/(2*sigma*sigma);
a[i][j]=(float)exp(a[i][j]);
if(a[i][j]<0.0001) a[i][j]=0;
h_sum=h_sum+a[i][j];
}
}
for(i=0;i<h_size;i++)
{
for(j=0;j<h_size;j++)
{
a[i][j]=a[i][j]/h_sum;
}
}
for(i=0;i<h_size;i++)
{
for(j=0;j<h_size;j++)
{
printf("%4.4f ",a[i][j]);
}
printf("/n");
}
}
gs1.rar （29个子文件） 
gs1 
gs1.dsp 4KB gs1.plg 240B gs1.suo 7KB 1.cpp 1KB 1.dsp 3KB gs1.ncb 33KB 1.dsw 527B 1.opt 48KB gs1.vcproj.CONSEN-WXL.WDD.user 1KB gs1.opt 53KB gs1.sln 871B 1.plg 236B gs1.vcproj 5KB Debug 1.ilk 202KB BuildLog.htm 4KB vc60.idb 41KB vc90.pdb 44KB vc90.idb 27KB 1.obj 6KB gs1.exe.embed.manifest 2KB gs1.ilk 202KB 1.pdb 361KB vc60.pdb 60KB gs1.pch 286KB
1.exe 196KB gs1.pdb 361KB gs1.exe 196KB 1.ncb 33KB gs1.dsw 531B资源评论
