數(shù)學建模使用MATLAB進行數(shù)據(jù)擬合(共5頁)

1、精選優(yōu)質(zhì)文檔-傾情為你奉上1.線性最小二乘法x=19 25 31 38 44'y=19.0 32.3 49.0 73.3 97.8'r=ones(5,1),x.2;ab=ry % if AB=C then B=ACx0=19:0.1:44;y0=ab(1)+ab(2)*x0.2;plot(x,y,'o',x0,y0,'r')運行結果：2. 多項式擬合方法x0=1990 1991 1992 1993 1994 1995 1996;y0=70 122 144 152 174 196 202;a=polyfit(x0,y0,1)y97=polyval(

2、a,1997)x1=1990:0.1:1997;y1=a(1)*x1+a(2);plot(x1,y1)hold onplot(x0,y0,'*')plot(1997,y97,'o')3.最小二乘優(yōu)化3.1 lsqlin 函數(shù)例四：x=19 25 31 38 44'y=19.0 32.3 49.0 73.3 97.8'r=ones(5,1),x.2;ab=lsqlin(r,y)x0=19:0.1:44;y0=ab(1)+ab(2)*x0.2;plot(x,y,'o',x0,y0,'r')3.2 lsqcurvefit

3、 函數(shù)（1）定義函數(shù)function f=fun1(x,tdata);f=x(1)+x(2)*exp(-0.02*x(3)*tdata); %其中x(1)=a，x(2)=b，x(3)=k（2）td=100:100:1000;cd=4.54 4.99 5.35 5.65 5.90 6.10 6.26 6.39 6.50 6.59;x0=0.2 0.05 0.05;x=lsqcurvefit(fun1,x0,td,cd) % x(1)=a，x(2)=b，x(3)=kt=100:10:1000;c=x(1)+x(2)*exp(-0.02*x(3)*t);plot(t,c)hold onplot(td

4、,cd,'*')3.3 lsqnonlin 函數(shù)（1）定義函數(shù)function f=fun2(x);td=100:100:1000;cd=4.54 4.99 5.35 5.65 5.90 6.10 6.26 6.39 6.50 6.59;f=x(1)+x(2)*exp(-0.02*x(3)*td)-cd;（2）主函數(shù)x0=0.2 0.05 0.05; %初始值是任意的x=lsqnonlin(fun2,x0)3.4 lsqnonneg 函數(shù)MATLAB：c=0.0372 0.2869;0.6861 0.7071;0.6233 0.6245;0.6344 0.6170;d=0.8587;0.1781;0.0747;0.8405;x=lsqnonneg(c,d)4. 曲線擬合的用戶圖形界面求法Matlab 工具箱提供了命令cftool，該命令給出了一維數(shù)據(jù)擬合的交互式環(huán)境。具體執(zhí)行步驟如下：（1）把數(shù)據(jù)導入到工作空間；（2）運行cfto