Beispiel zur Schnellen Fourier-Transformationx := Array([-1,-1,-1,-1,1,1,1,1]); # real parts of data y := Array([0,0,0,0,0,0,0,0]); # imaginary parts of dataFFT(3,x,y): # transform datax;y; # imaginary parts of transformed data zip((a,b)->a+b*I, x, y): convert(%, list); iFFT(3,x,y): # check results x; y; y := map(fnormal, y);Gl\344ttung einer verrauschten harmonischen Schwingung durch Faltung mit einer GaussverteilungDefinition einer St\366rfunktion (Rauschen)noise := stats[random, normald]: Harmonische Schwingung mit \374berlagertem Rauschenre_data := Array([seq(sin(0.0625*k)+0.1*noise(),
k=1..2^8)]);im_data := Array([seq(0, k=1..2^8)]):xcoords := Array([seq(0.0625*k, k=1..2^8)]): Grafische Darstellungplotdata := convert(zip((a,b)->[a,b], xcoords,
re_data), list):plot(plotdata, style=POINT); Definition einer Gaussfunktion zur Faltung (Kern,kernel)re_kernel := Array([seq(exp(-100.0*(k/2^8)^2),
k=1..2^8)]):im_kernel := Array([seq(0, k=1..2^8)]): FFT in den Frequezbereich von verrauschten Daten und KernFFT(8, re_data, im_data):FFT(8, re_kernel, im_kernel):Punktweise Multiplikation der transformierten Datendata := zip((a,b)->(a+b*I), re_data, im_data):kernel := zip((a,b)->(a+b*I), re_kernel, im_kernel):newdata := zip((a,b)->a*b, data, kernel):new_re_data := map(Re, newdata):new_im_data := map(Im, newdata): R\374cktransformation mittels inverser FFTiFFT(8, new_re_data, new_im_data):Grafische Drastellung der gegl\344tteten Datenplotdata := convert(zip((a,b)->[a,b], xcoords,
new_re_data), list):plot(plotdata, style=POINT); TTdSMApJM1JUQUJMRV9TQVZFLzI0MDU1MlgqJSlhbnl0aGluZ0c2IjYiW2dsISElISEhIikiKSQhKysrKys1ISIqJCErISoqKioqKioqKiEjNSQhCismKioqKioqKioqRiwkISsmKSoqKioqKioqRiwkIisrKysrNUYpJCIrISoqKioqKioqKkYsJCIrJioqKioqKioqKkYsJCIrJikqKioqKioKKipGLEYmCg==TTdSMApJM1JUQUJMRV9TQVZFLzI0MjA3MlgqJSlhbnl0aGluZ0c2IjYiW2dsISElISEhIikiKSQiIiFGKCQiKysrKytEISM+RickISsrKysrREYrCkYnRixGJ0YpRiYKTTdSMApJM1JUQUJMRV9TQVZFLzI0MDU1MlgqJSlhbnl0aGluZ0c2IjYiW2dsISElISEhIikiKSQhKysrKys1ISIqJCErISoqKioqKioqKiEjNSQhCismKioqKioqKioqRiwkISsmKSoqKioqKioqRiwkIisrKysrNUYpJCIrISoqKioqKioqKkYsJCIrJioqKioqKioqKkYsJCIrJikqKioqKioKKipGLEYmCg==TTdSMApJM1JUQUJMRV9TQVZFLzI0MjA3MlgqJSlhbnl0aGluZ0c2IjYiW2dsISElISEhIikiKSQiIiFGKCQiKysrKytEISM+RickISsrKysrREYrCkYnRixGJ0YpRiYKTTdSMApJM1JUQUJMRV9TQVZFLzI0MDU1MlgqJSlhbnl0aGluZ0c2IjYiW2dsISElISEhIikiKSQhKysrKys1ISIqJCErISoqKioqKioqKiEjNSQhCismKioqKioqKioqRiwkISsmKSoqKioqKioqRiwkIisrKysrNUYpJCIrISoqKioqKioqKkYsJCIrJioqKioqKioqKkYsJCIrJikqKioqKioKKipGLEYmCg==TTdSMApJM1JUQUJMRV9TQVZFLzI0MjA3MlgqJSlhbnl0aGluZ0c2IjYiW2dsISElISEhIikiKSQiIiFGKCQiKysrKytEISM+RickISsrKysrREYrCkYnRixGJ0YpRiYKTTdSMApJNFJUQUJMRV9TQVZFLzQ4MzExODhYKiUpYW55dGhpbmdHNiI2IltnbCEhJSEhISIpIikkIiIhRihGJ0YnJCEiIUYoRidGKUYnRidGJgo=TTdSMApJNFJUQUJMRV9TQVZFLzUzNjEzNjRYKiUpYW55dGhpbmdHNiI2IltnbCEhJSEhISJbW2wiW1tsJCIrZSYzZy0kISIpJCIrbyYzaWAkRikkCiIrYFs6JjMpRikkISsiKT5sNyIpRikkIStULGRXQEYpJCErMGcuTjVGKSQhK09faig+KCEiKiQhKylvJHpUT0Y2JCErRy1uIlIkRjYkIikKOW46YEY2JCErITM3XUQiRjYkIStMOm9eOEY2JCErJ0dlKDMpKiEjNSQhKjR0cnQlRjYkISoxIlFFb0Y2JCErOz4kPiRwRkMkIStZNms5NwpGNiQhK151ImV3IkY2JCErNVtMbGFGQyQhKiZbJylwd0ZDJCIrQSdwZzYiRjYkIitxXHR2JilGQyQiK2A9Jyo+PUY2JCEraCd6Z18oRkMkCiIrZi8kWyYpKUZDJCErNWZ4SXhGQyQiKyYpUikpKWYnRkMkIiteKjRacyJGNiQiK3dMIz47J0ZDJCIrOEpHWzVGNiQhKjlbKlxwRjYkIisKWFhYJmYiRjYkISopW3FWO0ZDJCIrJUgxOD4iRkMkIiorTjFpIkZDJCIrdktuXCQqRkMkISgpR1AhKkY2JCErNVtiOz9GQyQhKzs6Syk+JApGQyQiK0F0MUI9RjYkISolUicqR3lGNiQhK2pYUi83RjYkIiotYVVTJEZDJCIrOnkpKWY2RjYkISskSF94dCJGQyQiKz4/OzY5RjYkIiskCipcb188RjYkISk7PVs1RkMkIitebkJEIilGQyQhK1UnPl4vIkY2JCIrUSkqSFI9RjYkIiteSSkqUV5GQyQiK2lfQmBERjYkIitcJnlbVCIKRjYkIStCYEg/PkY2JCIrKVx2Jlw7RjYkIStoYzpeQEZDJCIqSWU9RCdGQyQiK0YzUkNARjYkISpvZVkhKSlGNiQiKnB3JmV2RjYkISskPgoxY20iRjYkIip0KG8nKUdGNiQiKnR0aihvRjYkISstX1lIOUY2JCIrPjwoPmsiRjYkIiolZSlSZiVGNiQhKz52THI6RjYkIitsbHExNUY2CiQhKyhbd3QtJEY2JCIpJjRZKydGNiQiK1NTQjBvRkMkIisvWCw2OkY2JCErdlwxKlsoRkMkIitOcm4kPSJGNiQiK3FYXjc7RjYkISs4KHoKTU0jRjYkISskeVFnJkdGQyQiKmFSQS4iRjYkIithX2AmKSoqRkMkIStiNjkwOUY2JCErWDYhW00nRkMkIit0cUZZZEZDJCErdl8hXGYiRgpDJCErLEtLVzZGNiQiK1xeV2VHRkMkIisvKikqNCsiRjYkIis7K18uNUY2JCErTHonUnkjRjYkIipQKCo0JSopRkMkIiszSEZ5NkY2JCIrCjIvRVM1RjYkIitgUjwneiNGNiQhKnQzMzEkRjYkIioqeXJzcUY2JCIqLls8LCJGNiQhK3NCc15vRkMkIit1bG5vTUZDJCIrKlJCXmQpRkMKJCErKFJAYmsiRkMkISsiPm9BKyNGQyQhKz1FVi5CRjYkIitHbCNcTCJGNiQiK28ubyVII0ZDJCIrI3B5WV0iRjYkISsuQHUqPiJGNiQiKwpBeHgjKkdGQyQiK3A3cSE0IkY2JCEraGhzUThGNiQiKyQ0WVsrIkY2JCIrXSJHTGEiRkMkISssKCpwXnRGQyQhKjlZUT8oRkMkIitxMy0mCjMqRkMkIStMUVdtNkY2JCErIlxWQmsmRkMkIis6aSY9aCJGNiQhKzxQZy5DRkMkISppMEFaJ0ZDJCIrb2pBYz9GNiQiKkZYNkonRjYkISoKQnM8aSRGNiQhKWQnUVglRjYkIStLWypwYiJGNiQiKyh5dENFI0Y2JCIqWiJ5SzVGKSQiKXZRKFsjRikkIik4TT4lKUYpJCIpIzQ4USdGKQokIik3TT4lKUYpJCIpdFEoWyNGKUZeW2wkIisienRDRSNGNiQhKz1bKnBiIkY2JCEpZydRWCVGNiQhKklzPGkkRjYkIipHWDZKJ0Y2JCIrCnFqQWM/RjYkISprL0FaJ0ZDJCErVU9nLkNGQyQiKzZpJj1oIkY2JCEqUlZCayZGNiQhK1NRV202RjYkIitDMy0mMypGQyQhKilmJVE/KEYKQyQhK1onKnBedEZDRmB5JCIrMGglWysiRjYkIStjaHNROEY2JCIrcTdxITQiRjYkIis0eHgjKkdGQ0ZmeCQiKylveVldIkY2JCIrYS5vJQpII0ZDRmB4JCErQEVWLkJGNiQhK2sib0ErI0ZDJCEqXEBiayJGNiQiKlVCXmQpRjYkIipidydvTUY2RmR3JCIrLyFbPCwiRkMkIitleXJzCnFGQyQhK04mMzMxJEZDJCIrZVI8J3ojRjYkIis4L0VTNUY2JCIrK0hGeTZGNiQiKjMoKjQlKilGQyQhK016J1J5I0Y2JCIrQCtfLjVGNiQKIisoKikpKjQrIkY2JCIrd11XZUdGQyQhKy1LS1c2RjYkISteXyFcZiJGQyQiK3FxRllkRkMkIStRNiFbTSdGQyQhK2A2OTA5RjYkIis+YApgJikqKkZDJCIrdCVSQS4iRkMkIStEKFFnJkdGQyQhKzkoek1NI0Y2JCIrImU5RGgiRjYkIitScm4kPSJGNiQhK2dcMSpbKEZDRmR0JCIrCnFTQjBvRkMkIikyaC9nRjYkISslW3d0LSRGNiQiK3BscTE1RjYkIStGdkxyOkY2JCIqbiYpUmYlRjYkIis8PCg+ayJGNkZkcyQiKzNQUHcKb0ZDJCIqYyhvJylHRjYkISsiPTFjbSJGNiQiKnR3JmV2RjYkISpwZVkhKSlGNiQiK0wzUkNARjYkIio/ZT1EJ0ZDJCErXmM6XkBGQyQiKwo7YmRcO0Y2RmByJCIrYyZ5W1QiRjYkIituX0JgREY2JCIrLkopKlFeRkMkIitPKSpIUj1GNiQhK1YnPl4vIkY2JCIrdG5CRCIpRkMkISlwCjxbNUZDJCIrJilcb188RjYkIitHPzs2OUY2JCErO0J2UDxGQyQiKzx5KSlmNkY2JCIqXWBVUyRGQyQhK2dYUi83RjYkISooUicqR3lGNiQKIitHdDFCPUY2JCErKWVAJCk+JEZDJCEreVpiOz9GQyQhKCRHUCEqRjYkIipcdCdcJCpGNiQiKUJqPztGNiQiKkQxOD4iRjZGZG8kIitMWApYJmYiRjYkIStkJFsqXHBGQyQiK1lKR1s1RjYkIiteTSM+OydGQyQiK3kqNFpzIkY2JCIreFIpKSlmJ0ZDJCErY2V4SXhGQyQiKyVbSVsmCikpRkMkIStoJnpnXyhGQyQiK1c9Jyo+PUY2JCIrPVx0diYpRkMkIitEJ3BnNiJGNiQhKjVoKXB3RkMkIStBXUxsYUZDJCErX3UiZXciRjYKRkokIStTPCQ+JHBGQyQhKzY3UUVvRkMkIStWRTxQWkZDJCErNCNlKDMpKkZDJCErKVwib144RjYkISsyQCxiN0Y2JCIpTW46YEY2JCErQAotbiJSJEY2JCErWk96VE9GNiQhK2Bfaig+KEY2JCErKSpmLk41RikkIStSLGRXQEYpJCErJCk+bDciKUYpJCIrXls6JjMpRikkIitsJjNpCmAkRilGJgo=