00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055
00056
00057
00058
00059
00060
00061
00062
00063
00064
00065
00066
00067
00068
#ifndef __NOISEGENERATOR_H
00069
#define __NOISEGENERATOR_H
00070
00071
#include <stdlib.h>
00072
#include <time.h>
00073
#include <math.h>
00074
00075
00076
00077
00078 class NoiseGenerator{
00079
public:
00080 NoiseGenerator() {
00081
dimension=0;
00082 };
00083
00084 virtual ~NoiseGenerator(){}
00085
00086
00087
00088
00089 virtual void init(
unsigned int dimension) {
00090 this->dimension = dimension;
00091 };
00092
00093
00094
00095
00096
00097
virtual double generate(
double min,
double max) = 0;
00098
00099
00100
00101
00102
00103
00104
00105
00106 virtual void add(
double *value,
double p1,
double p2){
00107
for (
unsigned int i = 0; i <
dimension; i++){
00108 value[i]+=
generate(p1,p2);
00109 }
00110 }
00111
00112
protected:
00113
00114 double uniform(
double min=-0.1,
double max=0.1){
00115
return( (double(rand())/RAND_MAX)*(max-min)+min);
00116 }
00117 unsigned int dimension;
00118 };
00119
00120
00121 class WhiteUniformNoise :
public NoiseGenerator{
00122
public:
00123 WhiteUniformNoise(){}
00124 virtual ~WhiteUniformNoise(){}
00125 virtual double generate(
double min,
double max) {
00126
return uniform(min,max);
00127 };
00128 };
00129
00130
00131 class WhiteNormalNoise :
public NoiseGenerator{
00132
public:
00133 WhiteNormalNoise(){}
00134 virtual ~WhiteNormalNoise(){}
00135 virtual double generate(
double mean,
double stddev) {
00136
double x1=
uniform(0, 1);
00137
double x2=
uniform(0, 1);
00138
return( (sqrt(-2*log(x1)) *cos(2*M_PI*x2)) * stddev + mean) ;
00139 };
00140 };
00141
00142
00143 class ColorUniformNoise :
public NoiseGenerator{
00144
public:
00145
00146
00147
00148 ColorUniformNoise(
double tau=0.05)
00149 : tau(tau){
00150
sqrttau = sqrt(tau);
00151
mean1channel=0.0;
00152
mean=0;
00153 }
00154 virtual ~ColorUniformNoise(){
if(
mean) free(
mean);}
00155 virtual void init(
unsigned int dimension){
00156
NoiseGenerator::init(dimension);
00157
mean = (
double*)malloc(
sizeof(
double) * dimension);
00158
for (
unsigned int i=0; i<dimension; i++){
00159
mean[i]=0.0;
00160 }
00161 }
00162
00163 virtual double generate(
double min,
double max) {
00164
mean1channel+=
sqrttau *
uniform(min, max) -
tau *
mean1channel;
00165
return(mean1channel);
00166 }
00167
00168
00169
00170
00171
00172
00173 virtual void add(
double *value,
double min,
double max){
00174
for (
unsigned int i = 0; i < dimension; i++){
00175
mean[i]+=
sqrttau *
uniform(min, max) -
tau *
mean[i];
00176 value[i]+=mean[i];
00177 }
00178 }
00179
00180
protected:
00181 double tau;
00182 double sqrttau;
00183 double*
mean;
00184 double mean1channel;
00185 };
00186
00187
00188 class ColorNormalNoise :
public WhiteNormalNoise{
00189
public:
00190 ColorNormalNoise(
double tau=0.05)
00191 : tau(tau){
00192
sqrttau = sqrt(tau);
00193
mean = 0;
00194
mean1channel=0.0;
00195 }
00196
00197 virtual ~ColorNormalNoise(){
if(
mean) free(
mean);}
00198
00199 virtual void init(
unsigned int dimension){
00200
WhiteNormalNoise::init(dimension);
00201
mean = (
double*)malloc(
sizeof(
double) * dimension);
00202
for (
unsigned int i=0; i<dimension; i++){
00203
mean[i]=0.0;
00204 }
00205 }
00206
00207 virtual double generate(
double stddev,
double meanvalue) {
00208
mean1channel +=
sqrttau *
WhiteNormalNoise::generate(stddev, meanvalue) -
tau*
mean1channel;
00209
return(mean1channel);
00210 }
00211
00212
00213
00214
00215
00216
00217 virtual void add(
double *value,
double stddev,
double meanvalue){
00218
for (
unsigned int i = 0; i < dimension; i++){
00219
mean[i]+=
sqrttau *
WhiteNormalNoise::generate(stddev*
factor, meanvalue) -
tau*
mean[i];
00220 value[i]+=mean[i];
00221 }
00222 }
00223
00224
protected:
00225 double tau;
00226 double sqrttau;
00227 double*
mean;
00228 double mean1channel;
00229 double factor;
00230 };
00231
00232
00233 class SineWhiteNoise :
public NoiseGenerator{
00234
public:
00235
00236
00237
00238
00239
00240 SineWhiteNoise(
double omega,
double amplitude,
double phaseShift = M_PI/2,
00241
unsigned int channels = 0xFFFF)
00242 : omega(omega), amplitude(amplitude)
00243 ,
channels(
channels),
phaseShift(
phaseShift){
00244
t=0;
00245 }
00246
00247 virtual ~SineWhiteNoise(){}
00248
00249 virtual double generate(
double min,
double max) {
00250
t ++;
00251
return uniform(min, max) + sin(
t*
omega)*
amplitude*(max-min);
00252 }
00253
00254
00255
00256
00257
00258
00259 virtual void add(
double *value,
double min,
double max){
00260
00261
for (
unsigned int i = 0; i < dimension; i++){
00262
if(i <
channels)
00263 value[i]+=sin(
t*
omega+i*
phaseShift)*
amplitude*(max-min);
00264 value[i]+=
uniform(min, max);
00265 }
00266
t ++;
00267 }
00268 void setOmega(
double omega){
00269 this->omega=omega;
00270 }
00271 void setPhaseShift(
double phaseShift){
00272 this->phaseShift=phaseShift;
00273 }
00274
00275
protected:
00276 long int t;
00277 double omega;
00278 double amplitude;
00279 unsigned int channels;
00280 double phaseShift;
00281
00282 };
00283
00284
00285
#endif
00286