How to extend the fitting curve in command
6 views (last 30 days)
Show older comments
I want to fit a set of data with custom function. I wrote:
"" Wavelength=VarName3
Wavelength =
1.0e+03 *
0.6176
0.6273
0.6371
0.6482
0.6572
0.6684
0.6801
0.6902
0.7014
0.7103
0.7218
0.7333
0.7431
0.7524
0.7608
0.7721
0.7825
0.7934
0.8044
0.8179
0.8267
0.8365
0.8469
0.8565
0.8658
0.8771
0.8907
0.9028
0.9110
0.9215
0.9316
0.9407
0.9503
0.9605
0.9707
0.9797
0.9945
1.0004
1.0115
>> Intensity=VarName4
Intensity =
795.5000
788.5000
773.5000
727.0000
716.5000
694.5000
656.5000
626.5000
600.0000
595.5000
593.0000
589.5000
545.0000
540.5000
516.5000
509.0000
475.5000
470.0000
432.5000
394.0000
369.0000
357.5000
355.5000
329.0000
315.5000
307.0000
302.5000
300.0000
293.0000
291.0000
279.0000
275.0000
269.0000
265.0000
242.0000
240.0000
234.5000
233.0000
224.0000
>> f = fit(Wavelength , Intensity, 'a*100000000000*(x.^-5)*(exp(1.4E07/(T*x))-1).^-1', 'StartPoint', [10000,0.0014] )
f =
General model:
f(x) = a*100000000000*(x.^-5)*(exp(1.4E07/(T*x))-1).^-1
Coefficients (with 95% confidence bounds):
T = 8917 (8406, 9428)
a = 8.682e+06 (7.45e+06, 9.913e+06)
>> plot( f,Wavelength,Intensity,'o' )
""
And I got the plot as follow
Now I want to extend the fitting curve further like this figure:
What should I do??? I need your help
0 Comments
Accepted Answer
Ahmet Cecen
on 10 Aug 2014
All you need to do is calculate Y axis for that entire range using the fit you calculated. Pass x=0:1:1400 through f(x), then plot x vs f(x).
5 Comments
Ahmet Cecen
on 10 Aug 2014
Edited: Ahmet Cecen
on 10 Aug 2014
OOOH. Use the hold command to overlay the two figures. Say:
x=0:1:1400;
yhat=f(x);
plot(x,yhat); %As before
hold on;
plot(X,Y,'o');%Check http://www.mathworks.com/help/matlab/ref/linespec.html for marker types.
X and Y would be your actual training data, so I think X is wavelength and Y is intensity. Use the GUI afterwards to polish the figure.
More Answers (2)
Sagar bendre
on 12 Aug 2016
HI, I am plotting a PSD with pwelch. I want to increase in max range of wavelength. Ufortunatelly the range by software is holding till 3.14 cyc/m. I want to see max wavelength of my data..
pl let me know...
0 Comments
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!