I want to bring forward
1) f1[x_]:= Exp[-(x-y)^2/(2*(s^2))];
2) f2[x_]:= .1/(x + .1);
3)ff[d_]:=NIntegrate[f1[d-x] f2[x],{x,0,d}];
4) Plot[{f1[x],f2[x],ff[x]},{x,0,10}]
1) bis 4) is mathematica language
in Stata-language.
I'd like to have
integral(f1(d-x)*f2(x))dx from 0 to 0.1 /d = 0.1
integral(f1(d-x)*f2(x))dx from 0 to 0.2 /d = 0.2
.
.
.......................................................... /d = 10
in order to plot my convolution-function.
next step will be deconvolution,I hope .........
any advice will be appreciate very much.
andreas aschbacher
Dear fellows !
>
> I have two functions :
>
> i) f1[x_]:= Exp[-(x-y)^2/(2*(s^2))];
> ii) f2[x_]:= .1/(x + .1);
> with:
> s = 1;
> y = 1;
> ->
> now I make convolution-integral:
> ff[d_]:=NIntegrate[f1[d-x] f2[x],{x,0,d}];
> ->
> now I plot the three functions with:
> Plot[{f1[x],f2[x],ff[x]},{x,0,10}]
>
> // all the syntax is from mathematica //
>
> this is a new way to "compose" two functions,not additive as I always did
> with LevenbergMarquardt,now I have composition in the way of convolution!
> -I try to change parameters s and y to "fit" measurement-curve.
> my question:
> is there a deconvolution-possibility in Stata ??
> with -denormix- I don't have the possibility of de-convolution
> in the sense of above,despite it was useful to get to know it.
>
> andreas aschbacher
>
>
>
>
> --
> +++ NEU bei GMX und erstmalig in Deutschland: T�V-gepr�fter Virenschutz
> +++
> 100% Virenerkennung nach Wildlist. Infos: http://www.gmx.net/virenschutz
>
> *
> * For searches and help try:
> * http://www.stata.com/support/faqs/res/findit.html
> * http://www.stata.com/support/statalist/faq
> * http://www.ats.ucla.edu/stat/stata/
>
--
+++ NEU bei GMX und erstmalig in Deutschland: T�V-gepr�fter Virenschutz +++
100% Virenerkennung nach Wildlist. Infos: http://www.gmx.net/virenschutz
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
