計(jì)算物理學(xué)——線性與非線性擬合

計(jì)算物理學(xué)——線性與非線性擬合

ID:34604611

大?。?01.54 KB

頁數(shù):3頁

時(shí)間:2019-03-08

計(jì)算物理學(xué)——線性與非線性擬合_第1頁
計(jì)算物理學(xué)——線性與非線性擬合_第2頁
計(jì)算物理學(xué)——線性與非線性擬合_第3頁
資源描述:

《計(jì)算物理學(xué)——線性與非線性擬合》由會員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在教育資源-天天文庫。

1、計(jì)算物理學(xué)實(shí)驗(yàn)三——多元線性擬合&非線性擬合一、實(shí)驗(yàn)題目試編程完成下題:1.(多元線性擬合)已知一組實(shí)驗(yàn)數(shù)據(jù)如下所示,試用最小二乘法求它的多元線性擬合曲線.數(shù)據(jù)產(chǎn)生A0=2.2;A1=4;A2=3;A3=5;A4=4.3;X=rand(4,9);fori=1:9y(i)=A0+A1*X(1,i)+A2*X(2,i)-A3*X(3,i)+A4*X(4,i);endX=fix(X.*1000)/1000;y=fix(y.*1000)/1000;2.(非線性擬合)已知一組實(shí)驗(yàn)數(shù)據(jù)如下表所示,試用最小二乘法求它的二次多項(xiàng)式擬合

2、曲線.13456789101054211234二、流程圖三、代碼及結(jié)果11)線性擬合函數(shù)lfit.mfunctionA=lfit(x,y)%線性擬合x=x';[~,n]=size(x);X=zeros(n+1);Y=zeros(n+1,1);X(1,:)=[1,mean(x)];X(2:n+1,1)=(mean(x))';Y(1)=mean(y);forl=1:nX(l+1,2:n+1)=mean(x.*(x(:,l)*ones(1,n)));%求XY(l+1)=mean(y'.*x(:,l));%求YendA=XY

3、;%求系數(shù)主程序linefit.mclear;A0=2.2;A1=4;A2=3;A3=5;A4=4.3;x=rand(4,9);fori=1:9y(i)=A0+A1*x(1,i)+A2*x(2,i)-A3*x(3,i)+A4*x(4,i);endx=fix(x.*1000)/1000;y=fix(y.*1000)/1000;B=lfit(x,y);%多元線性擬合fprintf('多元線性擬合函數(shù):y=(%4.4f)+(%4.4f*)*x1+(%4.4f)*x2+(%4.4f)*x3+(%4.4f)*x4',B);[

4、m,n]=size(x);X=[ones(n,1)x'];B=regress(y',X)'%驗(yàn)證運(yùn)行結(jié)果:多元線性擬合函數(shù):y=(2.2016)+(4.0032*)*x1+(2.9990)*x2+(-4.9977)*x3+(4.2946)*x4B=2.20164.00322.9990-4.99774.29462)非線性擬合函數(shù)nlfit.mfunctionA=nlfit(x,y,n)%非線性擬合x=x(:);y=y(:);X(:,1)=ones(length(x),1);forl=1:nX(:,l+1)=x.*X(:,

5、l);%求XendA=Xy;%求超定方程主程序nlinefit.mclear;x=[1345678910];y=[1054211234];plot(x,y,'-');axis([110010]);holdon;A=nlfit(x,y,2);%二次非線性擬合x=1:10;y=(A(1))+(A(2))*x+(A(3))*x.^2;plot(x,y,'--');title('二次擬合函數(shù)曲線');%二次擬合曲線fprintf('二次擬合函數(shù):y=(%4.2f)+(%4.2f)*x+(%4.2f)*x^2',A);B=

6、polyfit(x,y,2)%驗(yàn)證運(yùn)行結(jié)果:二次擬合函數(shù):y=(13.46)+(-3.61)*x+(0.27)*x^22B=0.2676-3.605313.45973

當(dāng)前文檔最多預(yù)覽五頁,下載文檔查看全文

此文檔下載收益歸作者所有

當(dāng)前文檔最多預(yù)覽五頁,下載文檔查看全文
溫馨提示:
1. 部分包含數(shù)學(xué)公式或PPT動畫的文件,查看預(yù)覽時(shí)可能會顯示錯(cuò)亂或異常,文件下載后無此問題,請放心下載。
2. 本文檔由用戶上傳,版權(quán)歸屬用戶,天天文庫負(fù)責(zé)整理代發(fā)布。如果您對本文檔版權(quán)有爭議請及時(shí)聯(lián)系客服。
3. 下載前請仔細(xì)閱讀文檔內(nèi)容,確認(rèn)文檔內(nèi)容符合您的需求后進(jìn)行下載,若出現(xiàn)內(nèi)容與標(biāo)題不符可向本站投訴處理。
4. 下載文檔時(shí)可能由于網(wǎng)絡(luò)波動等原因無法下載或下載錯(cuò)誤,付費(fèi)完成后未能成功下載的用戶請聯(lián)系客服處理。