测试2，测试2

info = fitsinfo('spec-0446-51899-0276.fits')

rowend = info.BinaryTable.Rows;

data1 = fitsread('spec-0446-51899-0276.fits','binarytable',1)

flux = data1{1}

wave = data1{2}

syms x y

x=wave

y=flux

f=fit(y,x,'poly1')

plot(f,x,y)

x1=wave(2000:3000);

y1=flux(2000:3000);

plot(x1,y1)

flux1='a*x^b'

startpoints=[-1.5 5]

f1=fit(wave,flux,flux1,'start',startpoints,'exclude',wave<3.7&wave>3.85)

plot(f1,wave,flux)

title('Fit with data points')

 

 

3.1415926