|
From: Ilehi S. <ile...@ya...> - 2012-02-10 14:29:41
|
I use that expression in physics research, I want to fit this expression to determine the values of n, S, T, but this expression is calculated by the maple software and can not find a solution to loading by Gnuplot (conversion maple-gnuplot),in other hand Gnuplot does not understand the complex 'i', please help me to find a solution, thank you very much this is the expression f(x)=abs((.9756567752*exp(-(.1628102823*(1.+i))*sqrt(x))*(-5.084444444*10^(-21)*(5.208333335*10^21*exp(-0.450e-3*sqrt((6.28*i*x+1/t)/n)-67.50)*(S+sqrt((6.28*i*x+1/t)/n)*n)/(n*(2.25*10^10-(1.*(6.28*i*x+1/t))/n)*(sqrt((6.28*i*x+1/t)/n)*n-1.*S))-5.208333335*10^21*(S-1.5*10^5*n)/(n*(2.25*10^10-(1.*(6.28*i*x+1/t))/n)*(sqrt((6.28*i*x+1/t)/n)*n-1.*S)))*(1.+0.2763953196e-2*sqrt((6.28*i*x+1/t)/n)/((1.+i)*sqrt(x)))/(t*((6.28*i*x+1/t)/n-1.308996939*10^5*(1.+i)^2*x))+n*(1.-0.2763953196e-2*sqrt((6.28*i*x+1/t)/n)/((1.+i)*sqrt(x)))+(26.48148149/(t*n*(2.25*10^10-(1.*(6.28*i*x+1/t))/n)*(2.25*10^10-1.308996939*10^5*(1.+i)^2*x))+1973.518519/(2.25*10^10-1.308996939*10^5*(1.+i)^2*x))*(1.+414.5929793/((1.+i)*sqrt(x)))-1.405316647*10^(-23)*S*(5.208333335*10^21*exp(-0.450e-3*sqrt((6.28*i*x+1/t)/n)-67.50)*(S+sqrt((6.28*i*x+1/t)/n)*n)/(n*(2.25*10^10-(1.*(6.28*i*x+1/t))/n)*(sqrt((6.28*i*x+1/t)/n)*n-1.*S))-5.208333335*10^21*(S-1.5*10^5*n)/(n*(2.25*10^10-(1.*(6.28*i*x+1/t))/n)*(sqr t((6.28*i*x+1/t)/n)*n-1.*S))+5.208333335*10^21*exp(-0.450e-3*sqrt((6.28*i*x+1/t)/n)-67.50)/(n*(2.25*10^10-(1.*(6.28*i*x+1/t))/n))-5.208333335*10^21/(n*(2.25*10^10-(1.*(6.28*i*x+1/t))/n)))/((1.+i)*sqrt(x)))-1.024343225*exp((.1628102823*(1.+i))*sqrt (x))*(-5.084444444*10^(-21)*(5.208333335*10^21*exp(-0.450e-3*sqrt((6.28*i*x+1/t)/n)-67.50)*(S+sqrt((6.28*i*x+1/t)/n)*n)/(n*(2.25*10^10-(1.*(6.28*i*x+1/t))/n)*(sqrt((6.28*i*x+1/t)/n)*n-1.*S))-5.208333335*10^21*(S-1.5*10^5*n)/(n*(2.25*10^10-(1.*(6.28*i*x+1/t))/n)*(sqrt((6.28*i*x+1/t)/n)*n-1.*S)))*(1.-0.2763953196e-2*sqrt((6.28*i*x+1/t)/n)/((1.+i)*sqrt(x)))/(t*((6.28*i*x+1/t)/n-1.308996939*10^5*(1.+i)^2*x))+n*(1.+0.2763953196e-2*sqrt((6.28*i*x+1/t)/n)/((1.+i)*sqrt(x)))+(26.48148149/(t*n*(2.25*10^10-(1.*(6.28*i*x+1/t))/n)*(2.25*10^10-1.308996939*10^5*(1.+i)^2*x))+1973.518519/(2.25*10^10-1.308996939*10^5*(1.+i)^2*x))*(1.-414.5929793/((1.+i)*sqrt(x)))+1.405316647*10^(-23)*S*(5.208333335*10^21*exp(-0.450e-3*sqrt((6.28*i*x+1/t)/n)-67.50)*(S+sqrt((6.28*i*x+1/t)/n)*n)/(n*(2.25*10^10-(1.*(6.28*i*x+1/t))/n)*(sqrt((6.28*i*x+1/t)/n)*n-1.*S))-5.208333335*10^21*(S-1.5*10^5*n)/(n*(2.25*10^10-(1.*(6.28*i*x+1/t))/n)*(sqrt((6.28*i*x+1/t)/n)*n-1.*S))+5.208333335*10^21*exp(- 0.450e-3*sqrt((6.28*i*x+1/t)/n)-67.50)/(n*(2.25*10^10-(1.*(6.28*i*x+1/t))/n))-5.208333335*10^21/(n*(2.25*10^10-(1.*(6.28*i*x+1/t))/n)))/((1.+i)*sqrt(x)))-1.016888889*10^(-20)*(5.208333335*10^21*exp(-0.450e-3*sqrt((6.28*i*x+1/t)/n)-67.50)*(S+sqrt((6.28*i*x+1/t)/n)*n)/(n*(2.25*10^10-(1.*(6.28*i*x+1/t))/n)*(sqrt((6.28*i*x+1/t)/n)*n-1.*S))-5.208333335*10^21*(S-1.5*10^5*n)/(n*(2.25*10^10-(1.*(6.28*i*x+1/t))/n)*(sqrt((6.28*i*x+1/t)/n)*n-1.*S)))*(0.2434322478e-1-0.2763953196e-2*sqrt((6.28*i*x+1/t)/n)/((1.+i)*sqrt(x)))*exp(-0.450e-3*sqrt((6.28*i*x+1/t)/n))/(t*((6.28*i*x+1/t)/n-1.308996939*10^5*(1.+i)^2*x))+2.*n*(0.2434322478e-1+0.2763953196e-2*sqrt((6.28*i*x+1/t)/n)/((1.+i)*sqrt(x)))*exp(0.450e-3*sqrt((6.28*i*x+1/t)/n))+9.686178480*10^(-30)*(26.48148149/(t*n*(2.25*10^10-(1.*(6.28*i*x+1/t))/n)*(2.25*10^10-1.308996939*10^5*(1.+i)^2*x))+1973.518519/(2.25*10^10-1.308996939*10^5*(1.+i)^2*x))*(0.2434322478e-1-414.5929793/((1.+i)*sqrt(x))))/(.9751817626*exp(-(.1628102 823*(1.+i))*sqrt(x))-1.024841941*exp((.1628102823*(1.+i))*sqrt(x)))) experimental data Amplitude (f(x)) frequency(x) 3.3166247903554 1 3.54964786985977 0.914285714285714 3.75499667110372 0.828571428571429 3.92428337406972 0.8 4.06201920231798 0.771428571428571 4.23083916026124 0.714285714285714 4.38178046004133 0.685714285714286 4.56070170039655 0.714285714285714 4.7116875957559 0.685714285714286 4.9295030175465 0.6 4.9295030175465 0.628571428571428 5.05964425626941 0.6 5.23450093132096 0.628571428571428 5.46808924579693 0.571428571428571 5.81377674149945 0.542857142857143 6.04979338490167 0.542857142857143 6.27694193059009 0.514285714285714 6.64830805543787 0.485714285714286 6.92098258919931 0.457142857142857 7.3348483283569 0.428571428571428 7.61577310586391 0.4 7.90569415042095 0.4 8.14248119432891 0.4 8.40832920383116 0.4 8.73498712076898 0.371428571428571 9.04986187739901 0.371428571428571 9.42867965305853 0.371428571428571 9.91463564635635 0.4 10.4163333279998 0.371428571428571 10.7191417566893 0.342857142857143 11.0589330407594 0.371428571428571 11.3401940018679 0.371428571428571 11.8406080924925 0.371428571428571 12.169634341261 0.4 12.5019998400256 0.371428571428571 12.8062484748657 0.4 13.1529464379659 0.371428571428571 13.4126805672841 0.4 13.5683455144686 0.4 13.8130373198656 0.4 14.0712472794703 0.4 14.1067359796659 0.4 14.3527000944073 0.314285714285714 14.5945195193264 0.371428571428571 14.7105404387466 0.371428571428571 14.8963082674869 0.342857142857143 15.0864177325169 0.371428571428571 15.1986841535707 0.371428571428571 15.4175224987674 0.4 15.6204993518133 0.371428571428571 |
|
From: Hans-Bernhard B. <HBB...@t-...> - 2012-02-10 21:45:45
|
On 10.02.2012 15:29, Ilehi Soufiene wrote:
> I use that expression in physics research, I want to fit this
> expression to determine the values of n, S, T, but this expression
> is calculated by the maple software and can not find a solution to
> loading by Gnuplot (conversion maple-gnuplot),in other hand Gnuplot
> does not understand the complex 'i', please help me to find a
> solution,
That part is trivially fixed by defining 'i':
i = {0,1}
But there's more you'll have to change. For starters, '^' isn't a power
operator in gnuplot. You'll have to replace that by '**' --- or just
write those numbers in plain floating point notation, i.e. 5.08444E-21.
|