實(shí)驗(yàn)二-數(shù)據(jù)建模

1、年級、專業(yè) 數(shù)學(xué)與應(yīng)用數(shù)學(xué) 2012 級 姓名 張旭 學(xué)號 名單序號 16實(shí)驗(yàn)時(shí)間 2014 年 4 月 13 日 使用設(shè)備、軟件PC, MA TLAB注： （ 1）實(shí)驗(yàn)報(bào)告的最后一部分是實(shí)驗(yàn)小結(jié)與收獲（ 2）本實(shí)驗(yàn)的數(shù)據(jù)均在 Excel工作簿： 實(shí)驗(yàn)二 數(shù)據(jù)建模實(shí)驗(yàn) 實(shí)驗(yàn)報(bào)告數(shù)據(jù) .xls中的 某個(gè)工作表中，工作表名稱在各題中給出實(shí)驗(yàn)二 數(shù)據(jù)建模插值、擬合與回歸1. 編寫一個(gè)函數(shù)實(shí)現(xiàn) Lagrangian 插值。 function f=lagrangian(x,y,x0) syms m;if length(x)=length(y) n=length(x);else disp('x 與

2、 y維數(shù)不等 '); return;end f=0.0;for k=1:nl=y(k);for j=1:k-1l=l*(m-x(j)/(x(k)-x(j); end for j=k+1:nl=l*(m-x(j)/(x(k)-x(j); end f=f+l;simplify(f);if k=nif nargin=3 f=subs(f,'m',x0); else f=collect(f); f=vpa(f,6);endendend x0=0:0.01:6;x=1,2,3,4,5;y=4,2,5,3,1;y0= lagrangian(x,y,x0);plot(x,y,'

3、;o',x0,y0)2. 在下面的函數(shù)中選擇兩個(gè)，在n 個(gè)結(jié)點(diǎn)上（ n 不要太大，如 5 11）分別用 Lagrangian 、1/ 72014 春 數(shù)學(xué)實(shí)驗(yàn) 實(shí)驗(yàn)二 數(shù)據(jù)建模插值、擬合與回歸年級、專業(yè) 數(shù)學(xué)與應(yīng)用數(shù)學(xué) 2012 級 姓名 張旭 學(xué)號 名單序號 16實(shí)驗(yàn)時(shí)間 2014 年 4 月 13 日 使用設(shè)備、軟件PC, MA TLAB注： （ 1）實(shí)驗(yàn)報(bào)告的最后一部分是實(shí)驗(yàn)小結(jié)與收獲（ 2）本實(shí)驗(yàn)的數(shù)據(jù)均在 Excel工作簿： 實(shí)驗(yàn)二 數(shù)據(jù)建模實(shí)驗(yàn) 實(shí)驗(yàn)報(bào)告數(shù)據(jù) .xls中的 某個(gè)工作表中，工作表名稱在各題中給出實(shí)驗(yàn)二 數(shù)據(jù)建模插值、擬合與回歸分段線性、 三次樣條插值三種插值方

4、法， 計(jì)算 m個(gè)插值點(diǎn)的函數(shù)值 （m要適中，如 50 100）。通過數(shù)值和圖形輸出， 將三種插值結(jié)果與精確值進(jìn)行比較。 適當(dāng)增加 n, 再做比較， 由此做初步分析。（1） y sin（ x）,0 x 2 ;（2） y1 x2 , 1 x 1;10(3) y cos x, 2 x 2;( 4) y e(1) x=1,2,3,4,5,6,7;y=sin(x); x0=0:0.01:7;y0= lagrangian(x,y,x0); plot(x,y,'o',x0,y0) hold on y0=interp1(x,y,x0);plot(x0,y0) y0=interp1(x,y,x0

5、,'spline');plot(x0,y0)legend(' 拉格朗日 ','分段線性插值 ','三次樣條插值 ') ； xlabel('x0 軸 ');ylabel('y0 軸 ');當(dāng) n=14 時(shí) x=1:0.5:7; y=sin(x);x0=0:0.01:8;y0= lagrangian (x,y,x0); plot(x,y,'o',x0,y0) hold on y0=interp1(x,y,x0);plot(x0,y0) y0=interp1(x,y,x0,'spl

6、ine');plot(x0,y0)legend(' 拉格朗日 ',' 分段線性插值 ',' 三次樣條差值 ') xlabel('x0 軸 ');ylabel('y0 軸 ');(2) x=-1:0.2:1;2/ 72014 春 數(shù)學(xué)實(shí)驗(yàn) 實(shí)驗(yàn)二 數(shù)據(jù)建模插值、擬合與回歸3/ 7年級、專業(yè) 數(shù)學(xué)與應(yīng)用數(shù)學(xué) 2012 級 姓名 張旭 學(xué)號 名單序號 16實(shí)驗(yàn)時(shí)間 2014 年 4 月 13 日 使用設(shè)備、軟件PC, MA TLAB注： （ 1）實(shí)驗(yàn)報(bào)告的最后一部分是實(shí)驗(yàn)小結(jié)與收獲（ 2）本實(shí)驗(yàn)的數(shù)據(jù)均在 Exc

7、el工作簿： 實(shí)驗(yàn)二 數(shù)據(jù)建模實(shí)驗(yàn) 實(shí)驗(yàn)報(bào)告數(shù)據(jù) .xls中的 某個(gè)工作表中，工作表名稱在各題中給出實(shí)驗(yàn)二 數(shù)據(jù)建模插值、擬合與回歸y=sqrt(1-x.2);x0=-1:0.01:1; y0= lagrangian (x,y,x0); plot(x,y,'o',x0,y0) hold on y0=interp1(x,y,x0);plot(x0,y0)y0=interp1(x,y,x0,'spline');plot(x0,y0)legend(' 拉格朗日 ',' 分段線性插值 ',' 三次樣條插值 ') xlabe

8、l('x0 軸 ');ylabel('y0 軸 ');(3)x=-2:0.2:2;y=(cos(x).10;x0=-1:0.01:1;y0= lagrangian (x,y,x0); plot(x,y,'o',x0,y0) hold on y0=interp1(x,y,x0); plot(x0,y0) y0=interp1(x,y,x0,'spline');plot(x0,y0)legend(' 拉格朗日 ',' 分段線性插值 ',' 三次樣條插值 ') xlabel('x0

9、 軸 ');ylabel('y0 軸 ');x=-2:0.2:2;y=exp(-(x).2);x0=-1:0.01:1;y0= lagrangian (x,y,x0); plot(x,y,'o',x0,y0) hold on y0=interp1(x,y,x0);plot(x0,y0)y0=interp1(x,y,x0,'spline');plot(x0,y0)2014 春 數(shù)學(xué)實(shí)驗(yàn) 實(shí)驗(yàn)二 數(shù)據(jù)建模插值、擬合與回歸年級、專業(yè) 數(shù)學(xué)與應(yīng)用數(shù)學(xué) 2012 級 姓名 張旭 學(xué)號 名單序號 16實(shí)驗(yàn)時(shí)間 2014 年 4 月 13 日 使用設(shè)備

10、、軟件PC, MA TLAB注： ( 1)實(shí)驗(yàn)報(bào)告的最后一部分是實(shí)驗(yàn)小結(jié)與收獲( 2)本實(shí)驗(yàn)的數(shù)據(jù)均在 Excel工作簿： 實(shí)驗(yàn)二 數(shù)據(jù)建模實(shí)驗(yàn) 實(shí)驗(yàn)報(bào)告數(shù)據(jù) .xls中的 某個(gè)工作表中，工作表名稱在各題中給出實(shí)驗(yàn)二 數(shù)據(jù)建模插值、擬合與回歸legend(' 拉格朗日 ',' 分段線性插值 ',' 三次樣條插值 ') xlabel('x0 軸 ');ylabel('y0 軸 ');3. 某天的氣溫變化如第 3 題數(shù)據(jù)表，試用最小二乘法找出這一天的氣溫變化規(guī)律。第 3 題數(shù)據(jù)表(工作表： 氣溫變化 )時(shí)刻 t(h)0

11、123456789101112溫度 T( )15141414141516182022232528時(shí)刻 t(h)131415161718192021222324溫度 T( )313231292725242220181716解：程序： x=0:24;y=15 14 14 14 14 15 16 18 20 22 23 25 28 31 32 31 29 27 25 24 22 20 18 17 16 c=lsqcurvefit(inline('c(1)*x.2+c(2)*x+c(3)','c','x'),0 0 0,x,y)結(jié)果： c = -0.09

12、36 2.5943 8.4157故溫度和時(shí)間的函數(shù)關(guān)系為 y=-0.0936*x2+2.5943*x+8.41574. 用電壓 V=10v 的電池給電容器充電，電容器上 t 時(shí)刻的電壓為 u(t) V (V V0)et/其中 V0是電容器的初始電壓，是充電常數(shù)。 試由第 4 題數(shù)據(jù)表中的一組 t,V 數(shù)據(jù)確定第 4 題數(shù)據(jù)表(工作表： 電容器 )t(s)0.51234579V(v)6.46.57.38.28.799.49.6解： = 原函數(shù)變型為 v-u(t)=(v-v0)*exp(-t/c), 兩邊取對數(shù)得： log(v-u(t)=-t/c+log(v-v0)令 x= v-u(t),a=-1

13、/c,y=log(v-v0), 這樣函數(shù)就變?yōu)橐淮尉€性函數(shù) x=a*t+y, 運(yùn)用 polyfit 求 解.t=0.5 1 2 3 4 5 7 9;v=6.4 6.5 7.3 8.2 8.7 9 9.4 9.6;x=10-v;m=log(x);t1=linspace(0,9,100);p=polyfit(t,y,1);pp =-0.27351.4458所以 y=1.4458 a=-0.27354/ 72014 春 數(shù)學(xué)實(shí)驗(yàn) 實(shí)驗(yàn)二 數(shù)據(jù)建模插值、擬合與回歸年級、專業(yè) 數(shù)學(xué)與應(yīng)用數(shù)學(xué) 2012 級 姓名 張旭 學(xué)號 名單序號 16實(shí)驗(yàn)時(shí)間 2014 年 4 月 13 日 使用設(shè)備、軟件PC, M

14、A TLAB注： ( 1)實(shí)驗(yàn)報(bào)告的最后一部分是實(shí)驗(yàn)小結(jié)與收獲( 2)本實(shí)驗(yàn)的數(shù)據(jù)均在 Excel工作簿： 實(shí)驗(yàn)二 數(shù)據(jù)建模實(shí)驗(yàn) 實(shí)驗(yàn)報(bào)告數(shù)據(jù) .xls中的 某個(gè)工作表中，工作表名稱在各題中給出實(shí)驗(yàn)二 數(shù)據(jù)建模插值、擬合與回歸y=log(10-v0) 得到 v0=5.7548 c=1/0.2735=3.6563;u=10-(10-5.7548)*exp(-0.2735*t)5. 給定數(shù)據(jù)見第 5 題數(shù)據(jù)表：第 5 題數(shù)據(jù)表(工作表： 三次樣條 )x0.250.300.390.450.53y0.50000.54770.62450.67080.7280分別就下列邊界條件求三次樣條函數(shù)S(x) 并作

15、圖。(1) S (0.25) 1,S (0.53) 0.6868;(2) S (0.25) S (0.53) 0.解：(1)x=0.25 0.30 0.39 0.45 0.53;y=0.5000 0.5477 0.6245 0.6708 0.7280;m=csape(x,y,'complete',1,0.6868);fnplt(s,'r');hold onplot(x,y,'o')m.coefs1.8863-1.01431.00000.50000.7952-0.73140.91270.54770.6320-0.51670.80040.62450.

16、3151-0.40290.74520.6708ans =(2)x=0.25 0.30 0.39 0.45 0.53; y=0.5000 0.5477 0.6245 0.6708 0.7280;m=csape(x,y,'variational');fnplt(m,'r'); hold on plot(x,y,'o') m.coefs ans =-6.26520.00000.96970.50001.8813-0.93980.92270.54772014 春 數(shù)學(xué)實(shí)驗(yàn) 實(shí)驗(yàn)二 數(shù)據(jù)建模插值、擬合與回歸5/ 7年級、專業(yè) 數(shù)學(xué)與應(yīng)用數(shù)學(xué) 2012 級 姓

17、名 張旭 學(xué)號 名單序號 16實(shí)驗(yàn)時(shí)間 2014 年 4 月 13 日 使用設(shè)備、軟件PC, MA TLAB注： ( 1)實(shí)驗(yàn)報(bào)告的最后一部分是實(shí)驗(yàn)小結(jié)與收獲( 2)本實(shí)驗(yàn)的數(shù)據(jù)均在 Excel工作簿： 實(shí)驗(yàn)二 數(shù)據(jù)建模實(shí)驗(yàn) 實(shí)驗(yàn)報(bào)告數(shù)據(jù) .xls中的 某個(gè)工作表中，工作表名稱在各題中給出實(shí)驗(yàn)二 數(shù)據(jù)建模插值、擬合與回歸-0.4600 -0.43180.79920.62452.1442-0.51460.74240.67086 題數(shù)據(jù)表：6. 某商品的需求量與消費(fèi)者的平均收入及商品價(jià)格的統(tǒng)計(jì)數(shù)據(jù)見第第 6 題數(shù)據(jù)表(工作表：需求量收入價(jià)格)需求量10075807050659010011060收入

18、10006001200500300400130011001300300價(jià)格5766875439(1) 建立回歸模型并進(jìn)行檢驗(yàn)；(2) 預(yù)測平均收入為 1000，價(jià)格為 6 時(shí)的商品需求量。解： box=xlsread(' 實(shí)驗(yàn)二 數(shù)據(jù)建模實(shí)驗(yàn) 實(shí)驗(yàn)報(bào)告數(shù)據(jù) .xls','Prob12 需求量輸入價(jià)格 ') x1=box(2,:)'x2=box(3,:)'x=x1,x2;x=ones(10,1),x;y=box(1,:)' b,bint,r,rint,stats=regress(y,x,0.05)9.95235.0477-5.7188-5.7109-8.4750-2.0929-4.33681.33442014 春 數(shù)學(xué)實(shí)驗(yàn) 實(shí)驗(yàn)二b =111.69180.0143-7.1882bint =56.0503167.3334-0.01200