#include <stdio.h> #include <numeric.h> #define N 7 /* u array size */ #define M 2 /* response function dimension */ #define N2 N-M+1 int main() { int i,j,n,m; array double u[N] = {6,19,32,45,58,71,42}, v[M] = {6,7}; array double q[N2], r[N]; printf("u=%f\n",u); printf("v=%f\n",v); deconv(u,v,q); /* real v, not specify remaind result */ deconv(u,v,q,r); /* real v, specified remaind result */ printf("q=%f\n",q); printf("r=%f\n",r); }
u=6.000000 19.000000 32.000000 45.000000 58.000000 71.000000 42.000000 v=6.000000 7.000000 q=1.000000 2.000000 3.000000 4.000000 5.000000 6.000000 r=0.000000 5.000000 32.000000 45.000000 58.000000 71.000000 42.000000