matlab拟合正弦曲线的问题

Posted 2023-02-25

我有一系列的随时间轴变化的测量数据，这些数据有噪声，但基本是正弦曲线，我使用下列程序进行拟合，对于单组数据效果倒是不错。但是这样的问题是sin函数中的四个参数都要先根据测量点的大致位置，人工进行初步计算得出
%对timep拟合
timepd=double(timep);
plength=1:length(timepd);
f=inline('fitp(1)*sin(fitp(2)*x+fitp(3))+fitp(4)','fitp','x');
fitp=nlinfit(plength,timepd,f,[4 2*pi/19.625 -12.5*19.625/2/pi 61]);
f_sinp=fitp(1)*sin(fitp(2)*plength+fitp(3))+fitp(4);
如何使程序完全自动化呢，不用人工介入的。因为我同时要进行多组这样的数据要处理，所以必须自动化才行，也没法使用cftool，因为这样也没法进行n组数据处理。
请高手解答。答案中请不要给出使用cftool的方法，我这边使用过cftool，拟合效果就很不好，振幅和周期都对不上。谢谢了。

conn.asp的部分字段如下：
dim conn,connstr,db,rs,rs_s,rs_s1
db="data/$database$169852.mdb"
set conn = server.createobject("adodb.connection")
connstr="provider=microsoft.jet.oledb.4.0;data source=" & server.mappath(""&db&"")
conn.open connstr
if error then
response.write "connect database error!"
response.end
end if 参考技术A 不好意思，可能是我说的不是很清楚，我用的是MDI框架，我想在主窗体(MainFrm)上做一个刷新的按钮(RefreshButton),在我点RefreshButton的时候，自动抓取当前活动Tabsheet上的DBGrid.DataSource.DataSet。