|
Last change
on this file since 8717 was
7276,
checked in by ole, 16 years ago
|
|
Merged numpy branch back into the trunk.
In ~/sandpit/anuga/anuga_core/source
svn merge -r 6246:HEAD ../../branches/numpy .
In ~/sandpit/anuga/anuga_validation
svn merge -r 6417:HEAD ../branches/numpy_anuga_validation .
In ~/sandpit/anuga/misc
svn merge -r 6809:HEAD ../branches/numpy_misc .
For all merges, I used numpy version where conflicts existed
The suites test_all.py (in source/anuga) and validate_all.py passed using Python2.5 with numpy on my Ubuntu Linux box.
|
|
File size:
417 bytes
|
| Line | |
|---|
| 1 | import numpy as num |
|---|
| 2 | from pylab import randn, conv, plot, psd, subplot, show |
|---|
| 3 | |
|---|
| 4 | |
|---|
| 5 | # Create signal |
|---|
| 6 | dt = 0.01 |
|---|
| 7 | t = num.arange(0, 10, dt) # Time vector |
|---|
| 8 | signal = 0.1*num.sin(2*num.pi*t) |
|---|
| 9 | |
|---|
| 10 | # Create noise |
|---|
| 11 | noise = conv(randn(len(t)), num.exp(-t/0.05))*dt |
|---|
| 12 | noise = noise[:len(t)] |
|---|
| 13 | |
|---|
| 14 | # Create noisy signal and plot |
|---|
| 15 | noisy_signal = signal + noise |
|---|
| 16 | |
|---|
| 17 | subplot(211) |
|---|
| 18 | plot(t, noisy_signal) |
|---|
| 19 | subplot(212) |
|---|
| 20 | psd(noisy_signal, 512, 1/dt) |
|---|
| 21 | |
|---|
| 22 | show() |
|---|
| 23 | |
|---|
Note: See
TracBrowser
for help on using the repository browser.
