RE: Fit and InverseLaplaceTransform. Improvement
- To: mathgroup at smc.vnet.net
- Subject: [mg23565] RE: Fit and InverseLaplaceTransform. Improvement
- From: "SANCHEZ DE LEON, Guillermo" <gsl at fab.enusa.es>
- Date: Sat, 20 May 2000 03:10:38 -0400 (EDT)
- Sender: owner-wri-mathgroup at wolfram.com
[Sorry I sent this question a few minutes ago, but I have yet found a a solution. I will apreciate any suggestion to improve my method] > Dear friend, I have the LaplaceTransform: > > In[1]:= Cout[s_] = (3*(7876.7668/(10.41635 + s) - 7891.21/(12.09768 + s))* > (off*(PS + s*Tu) + s*(PS + (on + s)*Tu)))/ > (off*(3*PS + 3*s*Tu + PS*s*Tu + p*s*(PS + s*Tu)) + > s*(3*PS + 3*on*Tu + on*PS*Tu + 3*s*Tu + PS*s*Tu + p*s*(PS + (on + > s)*Tu))); > > Where "off", "Tu", "PS", "on", "PS" are coeffs. to be calculated. I have > the experimental "data" of Cout[t]. For this reason I need to fit data to > a model given by " InverseLaplaceTransform[Cout[s], s, t] ". The > experimental Cout[t] data are: > > In[2]:=data = {{0.03, 27.5}, {0.08, 48.61}, {0.28, 133.39}, > {0.33, 106.18}, {0.38, 92.11}, {0.48, 60.47}, > {0.55, 56.54}, {0.65, 27.73}, {0.75, 23.25}, {0.85, 15.83}, > {0.95, 13.82}, {1.05, 10.58}, {1.15, 7.83}, > {1.25, 7.35}, {1.35, 6.08}, {1.45, 5.21}, > {1.6, 4.18}, {1.8, 3.48}, {2., 2.86}, > {2.2, 2.42}, {2.40, 2.11}, > {2.6, 2.00}, {2.78, 1.7}, {3., 1.58}, > {3.40, 1.58}, {3.75, 1.44}, {4.25, 1.35}, > {4.75, 1.3}, {5.25, 1.29}, {6.75, 0.81}, {7.25, 0.78}, > {7.75, 0.73}, {8.25, 0.75}, {11., 0.68}}; > > I have found the following method. > > In[3]:=sse2[(Tu_)?NumberQ, (off_)?NumberQ, (on_)?NumberQ, (p_)?NumberQ, (PS_)?NumberQ] := Block[{sol, q2}, q2[t_] = InverseLaplaceTransform[ Cout[s], s, t]; Plus @@ Apply[(q2[#1] - #2)^2 & , data, {1}]] > In[4]:=FindMinimum[sse2[Tu, off, on, p, PS], {Tu, 10, 20}, > {off, 0.1, 0.3}, {on, 0.4, 0.7}, {p, 1, 2}, {PS, 4, 5}, MaxIterations > -> 50] > out[4]:0{605.6289546197199, {Tu -> 6.467180826261866, off -> 0.0968608641750308, on -> 0.6629560337902439, p -> 0.971307914603113, PS -> 2.709414690357603}} Any suggestion to improve it > Guillermo Sanchez > guillerm at gugu.usal.es > gsl at fab.enusa.es > > > >