




下載本文檔
版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、關(guān)于用正交多項(xiàng)式做最小二乘擬合的實(shí)驗(yàn)報(bào)告1. 實(shí)驗(yàn)?zāi)康?用正交多項(xiàng)式做最小二乘擬合及擬合圖形2. 實(shí)驗(yàn)內(nèi)容:編寫用正交多項(xiàng)式做最小二乘擬合的程序,并用于求解一個(gè)任意給定的 數(shù)的3次多項(xiàng)式最小二乘擬合問(wèn)題,在這里給出數(shù)據(jù)如下:X11.31.61.92.22.52.83.13.43.74Y2.7183.6694.956.6869.02512.18216.44522.19829.96440.44754.598對(duì)表格中數(shù)據(jù)用正交多項(xiàng)式做最小二乘擬合在擬合完后作出擬合曲線的圖 形,計(jì)算平方誤差,最后對(duì)它們進(jìn)行分析。程序如下:1) .構(gòu)建的正交多項(xiàng)式最高項(xiàng)次數(shù)為 3時(shí)的程序:>> x= 1:0
2、34;>> y=2.7183.669 4.95 6.686 9.025 12.182 16.445 22.198 29.96440.447 54.598;>> n=3;構(gòu)建的正交多項(xiàng)式最高項(xiàng)次數(shù)為3>> result=inputdlg('請(qǐng)輸入權(quán)向量 w:','charpt-3',1,'1 1 1 1 1 1 1 1 1 1 1');>> w=str2num(char(result);利用str2num函數(shù)將數(shù)值型轉(zhuǎn)化為符號(hào)型>> m=le ngth(x)-1;>> s1=0
3、;>> s2=on es(1,m+1);>> v2=sum(w);>> d(1)=y*w'>> c(1)=d(1)/v2;>> for k=1: nxs=x.*s2.A2*w'a(k)=xs/v2;if(k=1)b(k)=0;elseb(k)=v2/v1;ends3=(x-a(k).*s2-b(k)*s1;v3=s3.A2*w'd(k+1)=y.*s3*w'c(k+1)=d(k+1)/v3;s1=s2;s2=s3;v1=v2;v2=v3;end>> r=y.*y*w'-c*d'
4、;r =0.8918>> alph=zeros(1,n+1)alph =0 0 0 0>> T=zeros(n+1,n+2);>> T(:,2)=ones(n+1,1);>> T(2,3)=-a(1);>>>> if(n>=2)for k=3:n+1for i=3:k+1T(k,i)=T(k-1,i)-a(k-1)*T(k-1,i-1)-b(k-1)*T(k-2,i-2); endendend>> for i=1:n+1for k=i:n+1 alph(n+2-i)=alph(n+2-i)+c(k)*T(
5、k,k+2-i); endend>> xmin=min(x);>> xmax=max(x);>> dx=(xmax-xmin)/(25*m);>> t=(xmin-dx):dx:(xmax+dx);>> s=alph(1);>> for k=2:n+1 s=s.*t+alph(k); end>> plot(x,y,'x',t,s,'-');>> grid on;>> disp(alph);>> disp(r)2) .構(gòu)建的正交多項(xiàng)式最高項(xiàng)次數(shù)為
6、4時(shí),將上述程序中的n改為4即可3) .構(gòu)建的正交多項(xiàng)式最高項(xiàng)次數(shù)為8時(shí),將上述程序中的n改為8即可3. 實(shí)驗(yàn)結(jié)果:1). n=3>> plot(x,y,'x',t,s,'-')Figure匚10" X |E 訂嚼 嗣itnstrt 工阿1卩 Vindov Help>> disp(alph);2.3365 -10.2050 19.2296 -8.9441>> disp(r)0.89182). n=4r =0.0185alph =0 0 0 0 00.5685-3.34859.8348-9.90635.60010.01852). n=89.9721e-007alph =0 0 0 0 0Columns 1 through 6-0.00490.1028 -0.9037 4.4386 -13.1107 24.3727Columns 7 through 9 -26.4044 17.4675 -3.2398 9.9721e-0074. 實(shí)驗(yàn)結(jié)論:從上面的圖像可以看出利用正交多項(xiàng)式做最小二乘擬合的效果很不錯(cuò),誤差很小,且正 交多項(xiàng)式的最高次數(shù)項(xiàng)次數(shù)越高,誤差越
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 印刷機(jī)械行業(yè)智能化市場(chǎng)消費(fèi)者需求變化趨勢(shì)分析考核試卷
- 紡織品售后服務(wù)滿意度分析考核試卷
- 財(cái)務(wù)風(fēng)險(xiǎn)管理與資金安全策略考核試卷
- 光纜在電力系統(tǒng)分布式發(fā)電控制中的應(yīng)用考核試卷
- 航道養(yǎng)護(hù)工程案例分析考核試卷
- 農(nóng)產(chǎn)品市場(chǎng)綠色通道建設(shè)考核試卷
- 部編人教版三年級(jí)語(yǔ)文上冊(cè)閱讀理解專項(xiàng)水平練習(xí)
- 2025年中國(guó)IPL光子嫩膚儀數(shù)據(jù)監(jiān)測(cè)研究報(bào)告
- 2025年中國(guó)LED柔性燈條數(shù)據(jù)監(jiān)測(cè)報(bào)告
- 2025年中國(guó)60度錐磨頭數(shù)據(jù)監(jiān)測(cè)報(bào)告
- 2024-2025北師大版小學(xué)數(shù)學(xué)四年級(jí)上冊(cè)期末考試測(cè)試卷及參考答案(共三套)
- 職業(yè)院校教師人工智能素養(yǎng):內(nèi)涵流變、框架構(gòu)建與生成路徑
- 呼吸閥行業(yè)市場(chǎng)發(fā)展現(xiàn)狀及趨勢(shì)與投資分析研究報(bào)告
- 互聯(lián)網(wǎng)企業(yè)產(chǎn)品經(jīng)理聘用合同
- 課程委托協(xié)議合同范例
- 2023-2024學(xué)年廣東省深圳高級(jí)中學(xué)七年級(jí)(上)期中歷史試卷
- HGT20638-2017化工裝置自控工程設(shè)計(jì)文件深度規(guī)范
- 應(yīng)征公民體格檢查表
- JT-T-1178.2-2019營(yíng)運(yùn)貨車安全技術(shù)條件第2部分:牽引車輛與掛車
- 2024年公務(wù)員考試《言語(yǔ)理解與表達(dá)》題庫(kù)附參考答案【綜合卷】
- 24春國(guó)家開(kāi)放大學(xué)《生活中的數(shù)學(xué)》形考作業(yè)1-4參考答案
評(píng)論
0/150
提交評(píng)論