Stata The Stata listserver
[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

Re: st: convolution


From   "Andreas Aschbacher" <aa_surf@gmx.at>
To   statalist@hsphsun2.harvard.edu
Subject   Re: st: convolution
Date   Mon, 22 Mar 2004 15:09:28 +0100 (MET)

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/



© Copyright 1996–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index