So while playing with some code to model a resonating spring (F=-k*pos), I discovered that the model I was using produced a very good sine wave approximation. This isn't normally something terribly interesting, as you can approximate sine/cosines very easily, but it was actually quite low complexity -- an iterative approximation with only 2 multiplies and 3 adds per sample. It can also generate the cos() (well, a 90 degree shifted signal) for each point for just one additional multiply.

The error is pretty low for the first few cycles, though after a bit it does drift in relation to the correct wave. I'm not going to spend too much more time on this, but if anybody wants to see if there's some way to correct it, go for it (it may just be rounding error, though, of course).

Here is the code:

// setup:
double cos_sc = PI/period_samples; // PI/n for a period of n
double pos=0; // actually sin(initial_state), 1 for cos(), 0 for sin()
double vel=1/cos_sc; // actually cos(initial_state)/cos_sc, 0 for cos(),1/cos_sc for sin()
double tmp = cos_sc*cos_sc;
double tmp2 = 1.0/(1.0+tmp);
double mul1 = (1.0-tmp)*tmp2;
double mul2 = tmp*2.0*tmp2;

// per-sample calculation
double output_value=pos;
// double cosine_output = vel*cos_sc;


// iterate to next sample
double newpos = pos*mul1 + vel*mul2;
vel -= (pos+newpos);
pos = newpos;

3 Comments

Another tag, yay

Recordings:

kerosene

7 Comments

Time to request a hearing, as I do not have any "graffiti" but they seem to think I do....

7 Comments

The smiley face was me...

Recordings:

bren - 1 -- [23:55]
bren - 2 -- [25:42]
bren - 3 -- [32:27]

5 Comments

We'll see if this works-- both with the kids and with the city...

7 Comments

(note the car #)

8 Comments

The 3 Cs: caffeine, calories, and cocktails... Of course you need the right amount of all, and the right sequencing, but if you get that all is forgiven.

And days like this don't hurt either...

2 Comments

Found my laptop with ergonomic keyboard!

Recordings:

freeform jam with newtchr

1 Comment

Can I just say it? It's such a great tool. I know, I know, it's for web pages, you say.. But I've been using it for tasks people would use perl for, and it's great. Converting Windows .rc files into SWELL-able C++, converting at&t assembly syntax to intel syntax for Jesusonic, such great things. Hooray for PHP. Easy to install on Win32 (just unzip), comes on OS X, it's a good time. A+.

Thank you that is all.

Recordings:

freeform jam with chr

2 Comments

Still hot out, for now.. Time to shotgun a contraband Bud Lite in the bathroom...

3 Comments

Woohoo the tiny maudio ekeys 37, $49, arrived... It's small, and light.

3 Comments

(thanks to Matt @ Open Labs)

7 Comments

(courtesy the Hilton Sonoma Wine Country Inn in Santa Rosa)

1 Comment

(and it ruled)

Comment...

"My friend Will, he's wicked smart."

Comment...