File: convolve.inp

package info (click to toggle)
gretl 2025a-1
  • links: PTS, VCS
  • area: main
  • in suites: trixie
  • size: 64,792 kB
  • sloc: ansic: 425,462; sh: 4,916; makefile: 3,251; cpp: 2,777; xml: 606; perl: 364
file content (15 lines) | stat: -rw-r--r-- 284 bytes parent folder | download | duplicates (8) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
 
# define two polynomials
a = { 1, 0.5,    0, 0 }'
b = { 1, 0.3, -0.8, 0 }'
# perform the transforms
f = fft(a ~ b)
# complex-multiply the two transforms 
fc = cmult(f[,1:2], f[,3:4]) 
print fc
# compute the coefficients of c = a * b via the inverse transform
c = ffti(fc)
print c