Gaussian white noise code in C

Hello all,

I'm in need of a C program which implements the Gaussian white noise
algorithm in C. The program is to accept a binary file as input and
alter it according to the algorithm.
I'm affraid this project must run on Windoze :( 
Thank you in advance

Jess
0
jessex (1)
12/7/2003 9:48:52 AM
comp.dsp 20141 articles. 0 followers. allnor (8507) is leader. Post Follow

6 Replies
1123 Views

Similar Articles

[PageSpeed] 49
"Jess" <jessex@kmfms.com> wrote in message
news:63bbea8b.0312070148.4a07e00@posting.google.com...
> Hello all,
>
> I'm in need of a C program which implements the Gaussian white noise
> algorithm in C. The program is to accept a binary file as input and
> alter it according to the algorithm.
> I'm affraid this project must run on Windoze :(

Jess,

What does Windows have to do with your question?  Presumably you have an
ANSI C compiler that's suitable for generating code that will run on the
usual platforms, no?

What does the "binary file as input" have to do with the question?  Are you
meaning to add the noise to the input data or.....?  If you intend to add
it, then you'll be concerned with scaling the noise to match the intended
application.  This includes the number represenation you'll be using (fixed,
floating, etc.)

Check out:
http://www.dspguru.com/howto/tech/wgn.htm
http://www.mathtools.net/Java/Random_numbers/C_and_C++/

Fred


0
Fred
12/7/2003 7:51:10 PM
On Sun, 07 Dec 2003 09:48:52 +0000, Jess wrote:

> Hello all,
> 
> I'm in need of a C program which implements the Gaussian white noise
> algorithm in C. The program is to accept a binary file as input and
> alter it according to the algorithm.
> I'm affraid this project must run on Windoze :( Thank you in advance
> 


Although I have not tried it, white noise should have a flat spectrum.
You could try setting all values to 1 in a complex array and performing
an inverse FFT on it.

Possibly a tweak or two is in order.

This should be guaranteed white noise !
0
Robin48gx
12/8/2003 12:05:35 AM
Robin48gx <robin48gx@hotmail.com> wrote in message > 
> Although I have not tried it, white noise should have a flat spectrum.
> You could try setting all values to 1 in a complex array and performing
> an inverse FFT on it.
> 
> Possibly a tweak or two is in order.
> 
> This should be guaranteed white noise !

Well, if you do this, you'll get an impulse - autocorrelation of white noise!

Kiran
0
kiranatiisc
12/8/2003 4:21:37 AM
Robin48gx wrote:

> Although I have not tried it, white noise should have a flat spectrum.
> You could try setting all values to 1 in a complex array and performing
> an inverse FFT on it.
> 
> Possibly a tweak or two is in order.
> 
> This should be guaranteed white noise !

This can not be less white noise!!  This is not even noise!!
This is a deterministic signal!!

Maybe (and only *maybe* -- as in, I'm almost sure that it
won't work, but there's a slight chance that it might) if
you set the magnitudes of the FFT to 1, and set the phase
to uniformly distributed random values between -pi and pi...
Maybe you'll get something similar to white noise...

And BTW, how do you expect this to produce *Gaussian* white
noise?  How do you know that it's not going to be uniformly
distributed white noise?

Carlos
--

0
Carlos
12/8/2003 5:49:44 AM
Robin48gx wrote:

> On Sun, 07 Dec 2003 09:48:52 +0000, Jess wrote:
> 
> 
>>Hello all,
>>
>>I'm in need of a C program which implements the Gaussian white noise
>>algorithm in C. The program is to accept a binary file as input and
>>alter it according to the algorithm.
>>I'm affraid this project must run on Windoze :( Thank you in advance
>>
> 
> 
> 
> Although I have not tried it, white noise should have a flat spectrum.
> You could try setting all values to 1 in a complex array and performing
> an inverse FFT on it.
> 
> Possibly a tweak or two is in order.
> 
> This should be guaranteed white noise !

Don't offer guarantees too freely. With one set of phases, it could be
the spectrum of an impulse. With another, of a chirp. Anyhow, summing 12 
to 20 random numbers is easier than a broadband IFFT.

Jerry
-- 
Engineering is the art of making what you want from things you can get.
�����������������������������������������������������������������������

0
Jerry
12/8/2003 4:52:41 PM
What you need to search for is the Box-Muller algorithm. So simpe that even
the meanest programmer should be able to implement it in C.

Nick Felgate

"Jess" <jessex@kmfms.com> wrote in message
news:63bbea8b.0312070148.4a07e00@posting.google.com...
> Hello all,
>
> I'm in need of a C program which implements the Gaussian white noise
> algorithm in C. The program is to accept a binary file as input and
> alter it according to the algorithm.
> I'm affraid this project must run on Windoze :(
> Thank you in advance
>
> Jess


0
Nick
12/12/2003 1:28:42 AM
Reply:
Similar Artilces:

US-TX-Austin: AutoCAD Design Tech, 3-5yrs exp drafting, AutoCAD2000 or ver.14; c (45341157609)
US-TX-Austin: AutoCAD Design Tech, 3-5yrs exp drafting, AutoCAD2000 or ver.14; c (45341157609) ============================================================================================== Position: AutoCAD Design Tech Reference: SMC01967 Location: Austin TX Duration: contract Skills: Associates degree in Design, AutoCAD or equilevant 3-5yrs exp in drafting, using document management systems, AutoCAD and computer-aided design software. AutoCAD 2000 or version 14 direct exp. Manual Drafting ...

morton code
Hello, I would like to use morton code to index gis data in a db. Could you give me some pointers to relevant litterature ? I have already googled a lot, I know how to compute the code, but how do I code the windows queries ? Any help highly appreciated. Noel Frankinet ...

Google Mathematica code search
Google code search came out today. It has an operator that allows users to restrict the domain to Mathematica code. It allows regular expressions in the search string. I was able to locate code in .m files in my Google code project. I was able to locate code in .nb files on the Internet. the relevant operator is: lang:mathematica here is the website: http://www.google.com/codesearch This would be a good search to add to the SetFaceAndFont2 pallette. -- http://chris.chiasson.name/ This should help me out a lot figuring out the mistakes I'm making trying to g...

Divine Divinity Cheat Codes
Are where? Ken Murphy <kmurph@worldnet.att.net> once tried to test me with: > Are where? > I don't think there are any. -- Knight37 U.N. Repressentative: So, Mr. Evil - Dr. Evil: It's Dr. Evil, I didn't spend six years in Evil Medical School to be called "mister," thank you very much. -- "Austin Powers" "Ken Murphy" <kmurph@worldnet.att.net> wrote in message news:jjgni0p7knlp1kt2grkoaitf4gp90kcl9a@4ax.com... > Are where? Try the Larian message boards for some links and sites with DD info. Pretty cool stuff once you fin...

reducing high frequency noise
Hello There, I have a problem with noise on my data' s. The problem I get data's from a laser vibrometer. A Matrix with one row and 2000 columns. If I plot the function I get a decreasing sinusoidal wave with noise on the curve. With which algorithm I get a fuction without these noise. I have to try the curve fitting (Polynomial) but I think it is rubbish. Is there somebody who can tell me a solution for my problem. Thank you for you help and advanced. Best regards Sebastian Hi Sebastian, You'll have to filter your datas with a highpass filter. 1. Use the FIR1 function to ...

additive white gaussian noise
Hi Does anyone know the algorithm or have a source code for applying Additive white gaussian noise? thanx very much -M- Mesvak wrote: > Does anyone know the algorithm or have a source code for applying > Additive white gaussian noise? For algorithms also try 'comp.programming'. Another newsgroup to visit for math stuff is 'sci.math.num-analysis'. Here we mostly discuss C++ language issues, and dispensation of algorithms and sources is usually done elsewhere. Like the Web. FYI. V -- Please remove capital 'A's when replying by e-mail I do not respond t...

Audio jamming using white gaussian noise (matlab code)
Hi anybdy have any idea how can a baseband audio sound can be supressed or jammed using a white gaussian noise and how the result can be shown in Matlab. any help will be appreciated. -Shruti ...

Gaussian smoothing using rlft3 (numerical recipes)
Hello, I need some help in using the rlft3 (Numerical Recipes in c++ book, Chapter 12) to apply a Gaussian smoothing to a 2D image. How do I create a suitable filter function (page 535)? Thanks, Joao j.f.c.neves@googlemail.com wrote: > I need some help in using the rlft3 (Numerical Recipes in c++ book, > Chapter 12) to apply a Gaussian smoothing to a 2D image. How do I > create a suitable filter function (page 535)? Wrong newsgroup. This one discusses the C++ language proper, not algorithms that happen to be implemented in C++ (see http://www.parashift.com/c++-faq-lite/how-to-pos...

Gaussian RV
could you please help me solving this matlab problem >>generate a sequence of N random numbers that correspond to a Gaussian RV with mean 65 and variance 225. Call this sequence seq1. Plot seq1, and plot its histogram. eeahmad@gmail.com wrote in message <15a9c762-ae01-469d-82dd- dde08bbd1981@59g2000hsb.googlegroups.com>... > could you please help me solving this matlab problem > > >>generate a sequence of N random numbers that correspond to a Gaussian RV with mean 65 and variance 225. Call this sequence seq1. Plot seq1, and plot its histogram. What h...

Memoryless gaussian source
I want to generate memoryless unit-variance Gaussian vector of length k. It is equivalent to generating one value for each of k independent normal random variables (a single realization of length k of a Gaussian process ). I suppose randn([1,k]) in Matlab will give me k values which would be realizations of a single Gaussian random variable. In order to acheive my goal, should I use randn() in a loop (k times) with a different state each time i.e. inp = []; for i = 1:k inp = [inp,randn('state',sum(100*clock))]; end ...

Another rich white ape beats his wife
http://gawker.com/bill-o-reilly-accused-of-domestic-violence-in-custody-b-1705006992 What the fuck is wrong with these animals? Shoot every last one and we can have a real civilization. tmelmosfire <tmelmosfire@gmail.com> wrote: > http://gawker.com/bill-o-reilly-accused-of-domestic-violence-in-custody-b-1705006992 > What the fuck is wrong with these animals? Shoot every last one and we can have a real civilization. OJ FTW On Monday, May 18, 2015 at 5:24:10 PM UTC-4, tmelmosfire wrote: > http://gawker.com/bill-o-reilly-accused-of-domestic-violence-in-custody-b-1...

pgsql-server: Code review for recent dbsize changes.
Log Message: ----------- Code review for recent dbsize changes. Fix some thinkos, enforce coding style and message style standards, improve documentation. Modified Files: -------------- pgsql-server/contrib/dbsize: README.dbsize (r1.3 -> r1.4) (http://developer.postgresql.org/cvsweb.cgi/pgsql-server/contrib/dbsize/README.dbsize.diff?r1=1.3&r2=1.4) dbsize.c (r1.14 -> r1.15) (http://developer.postgresql.org/cvsweb.cgi/pgsql-server/contrib/dbsize/dbsize.c.diff?r1=1.14&r2=1.15) dbsize.sql.in (r1.3 -> r1.4) (http://...

Anybody know something about an Aztec C Compiler for the C128...
.... which was able to generate code for all three modes of the Commodore 128 (= CP/M with the Z80, C128 native mode, C64 mode) ? It would be very, very nice to get any further information about it, and of course, any download link. The Aztec company itself was closed, so the software could be called 'abandonware', and it should not harm anybody to offer downloads with an Aztec compiler. Regards Peter -- * Try http://www.z80.eu for CP/M computer and software infos. "Peter Dassow (remove 'NOSPAM.')" <peter.dassow@NOSPAM.z80.eu> wrote in message news:foe...

Higher Order Statistics for Gaussian Distribution
Hi Everyone, I'm calculating higher order Cumulants (order>2) for Gaussian Distribution. Theoretically these Cumulants should be equal to zero. But I'm not getting these Cumulants equal to zero. The methodology is as under: Data is generated by randn(1,4096) built in function in Matlab. function [rcum] = cumulants(data,ncumulant) %ncumulant represent total number of cumulants to be calculated. %data is Gaussian in my case, generated using randn(1,4096). rcum = 0; n = 1; cum(1,n) = 0; n = 2; while(n <= ncumulant) kn = moment(data,...

Google Icon Vase Speaker
Many computer speakers have a rather grim look and too many wires, so consider this sound alternative: a single 2-way surround sound speaker with innovative audio lens technology, a tweeter tuned for clarity. T for more info http://www.googleiloveu.blogspot.com/ ...

temporary break in a code
Good morning I have a stupid question. I have a program, running for days to do simulations. This slows down my computer, so from time to time, I would like to interrupt my simulations and then to resume them when I am done doing other stuff I am sure there is a way, but I could not find it sorry for the stupid question!! Francesco > I have a stupid question. I have a program, running for days to do > simulations. This slows down my computer, so from time to time, I > would like to interrupt my simulations and then to resume them when > I > am done doing other stuff > &g...

error in loading data into the code
hello , I have started python last week,so maybe my question is a bit stupid. I got this error in my simple python code while trying to load data into the code. I have added this libaraies from pylab import * from scipy import * import matplotlib.mlab as mlab from numpy import * from scipy import optimize but get this following error : /usr/lib/pymodules/python2.6/matplotlib/mlab.py:1267: DeprecationWarning: use numpy.loadtxt warnings.warn("use numpy.loadtxt", DeprecationWarning) # 2009/06/13 any idea ? thanks,Nakisa nakisa <nakisa.nooraee@gmail.com> ...

C++/G++ Benchmark.
Hi Looking for G++/C++ Benchmark source codes links/resources readily available for free download for linux os. BR Mukesh K S. ...

US-TX-Austin: RF Mixed Signal Design Eng, Analog design exp., wireless design; C (45315357603)
US-TX-Austin: RF Mixed Signal Design Eng, Analog design exp., wireless design; C (45315357603) ============================================================================================== Position: RF Mixed Signal Design Eng. Reference: SMC01602 Location: Austin TX Duration: C-P Skills: MSEE or MSCE or BSEE plus 5 years of experience. Exp designing mixed-signal integrated circuit design at the transistor level. Strong emphasis on analog design required. Wireless design experience at GHz operating frequencies ...

A Code Point's Tale: There and Back Again
This is probably obvious in the docs and I'm just missing it, but here goes: So, I see there is str.each_codepoint, which I want to use in a function to convert Unicode Strings to a list of Unicode code points. But what can I do if I have a list of Unicode code points and want to convert them back into a String? -- Posted via http://www.ruby-forum.com/. I hope this is what u r looking for http://ruby-unicode.rubyforge.org/doc/ -- Posted via http://www.ruby-forum.com/. Hi, On 30.04.2011 06:12, Terry Michaels wrote: > This is probably obvious in the docs and I'm just missing ...

Verilog Netlest Reader Code, ATPG Code
I am looking for example code in C/C++ or Java that reads net lists, and formats them in the form for a ATPG algorithm. An example that worked with the D-Algorithm would be a good start. Second I am also looking for descriptions of one of the more advanced(faster) ATPG algorithms that is good enough to actual produce code from it. Anyone know of a good source? ...

modify txt file in C
hi everyone, i met some problem in this. i used fseek to locate position prepared to be written. for example, fseek(stream,16,seek_set) then fputs(time,stream) while time will be users input. the content of the txt file to be modified aabb, ccddeeff, 30gghhiijjkk obvious the string 30(at byte position 16) is going to be changed. but if users input were 3000 then gg(position 18) will be overwrote. i will get the result aabb, ccddeeff, 3000hhiijjkk but all i need is aabb, ccddeeff, 3000gghhiijjkk or EVEN users input were 300000 i still need the result to be aabb, ccddeeff, 300000gghhiijjk...

Gaussian complex random vector with Matlab
Hello, I would like to know how to generate a gaussian complex vector ( ~N(0,1) ) with Matlab. Thanks you for your help. Vincent here's the code, clear all; close all; a = randn(1,1000); b = randn(1,1000); c = a + j*b; c = c - mean(c); c = c/std(c); mean(c) std(c) rgds bharat pathak Arithos Designs www.Arithos.com DSP design consultancy and Training company. >Hello, > >I would like to know how to generate a gaussian complex vector >( ~N(0,1) ) with Matlab. >Thanks you for your help. > >Vincent > ...

Objective-C FAQ #7
Archive-name: computer-lang/Objective-C/faq Posting-Frequency: monthly Frequently Asked Questions - comp.lang.objective-c compiled by David Stes (stes@pandora.be) February 14 2004 Contents * Contents * 1. About this FAQ + 1.1 Where can I find the latest version of the FAQ ? * 2. Objective-C Compiler Commands + 2.1 What's the file suffix for Objective-C source ? + 2.2 How do I compile .m files with the Stepstone com...

this code is memory overflow or not?
#include char *code; void main() { char buf[8] = "book"; strcpy(code, buf); } blackrosezy@gmail.com wrote: > #include > > char *code; > > void main() > > { > > char buf[8] = "book"; > > strcpy(code, buf); > > } See FAQ 7.1 in comp.lang.c: http://www.faqs.org/faqs/C-faq/faq/index.html Cheers! --M blackrosezy@gmail.com wrote: > #include It seems your newsreader ate your include statement. I will assume it is <string.h> (or even <cstring> with appropriate using declarations). > char *code; You never init...