數學實驗實驗報告二答案_第1頁
數學實驗實驗報告二答案_第2頁
數學實驗實驗報告二答案_第3頁
數學實驗實驗報告二答案_第4頁
數學實驗實驗報告二答案_第5頁
全文預覽已結束

下載本文檔

版權說明：本文檔由用戶提供并上傳，收益歸屬內容提供方，若內容存在侵權，請進行舉報或認領

文檔簡介

實驗二實驗項目名稱：插值與擬合實驗內容：建立拉格朗日插值的函數M文件lagr.m。Newscript源程序：functiony=lagr(x0,y0,x)n=length(x0);m=length(x);fori=1:mz=x(i);s=0;fork=1:np=1;forj=1:nifj~=kp=p*(z-x0(j))/(x0(k)-x0(j));endends=p*y0(k)+s;endy(i)=s;end保存為larg.m完成機床加工問題中數據增加密度，利用三種插值方法計算加工所需數據，并畫出曲線，比較分析三種插值結果。（參考教材P77）源程序：x0=[0,3,5,7,9,11,12,13,14,15];

y0=[0,1.2,1.7,2.0,2.1,2.0,1.8,1.2,1.0,1.6];

x=0:0.1:15;

y1=lagr(x0,y0,x);

y2=interp1(x0,y0,x);

y3=interp1(x0,y0,x,'spline');

[x',y1',y2',y3'];

subplot(3,1,1)

plot(x0,y0,'k+',x,y1,'r')

grid,title('lagrange')

subplot(3,1,2)

plot(x0,y0,'k+',x,y2,'r')

grid,title('piecewise

linear')

subplot(3,1,3)

plot(x0,y0,'k+',x,y3,'r')

grid,title('spline')實驗結果：完成教材P80習題1中的函數y=sinx的三種插值結果的分析和比較。源程序：n=10;m=100;x=linspace(0,2*pi,100);y=sin(x);z=x*0;x0=linspace(0,2*pi,6);y0=sin(x0);y1=lagr(x0,y0,x);y2=interp1(x0,y0,x);y3=interp1(x0,y0,x,'spline');plot(x,z,x,y,'k',x,y1,x,y2,x,y3,'y')實驗結果：完成教材P81習題2。源程序：x0=[0,1,4,9,16];y0=sqrt(x0);x=5;y=sqrt(5)y1=lagr(x0,y0,x)x0(1)=[];y0=sqrt(x0);y2=lagr(x0,y0,x)x0(4)=[];y0=sqrt(x0);y3=lagr(x0,y0,x)x0(1)=[];y0=sqrt(x0);y4=lagr(x0,y0,x)實驗結果：y=2.2361y1=2.0794y2=2.2540y3=2.2667y4=2.2000理解教材的擬合應用實例-給藥方案，參考該實例完成教材P81習題6。源程序：（τ輸入成t0）t=[0.5,1,2,3,4,5,7,9];v=[6.36,6.48,7.26,8.22,8.66,8.99,9.43,9.63];y=log(10-v);a=polyfit(t,y,1);t0

溫馨提示

  • 1. 本站所有資源如無特殊說明，都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等，如果需要附件，請聯系上傳者。文件的所有權益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙，網頁內容里面會有圖紙預覽，若沒有圖紙預覽就沒有圖紙。
  • 4. 未經權益所有人同意不得將文件中的內容挪作商業或盈利用途。
  • 5. 人人文庫網僅提供信息存儲空間，僅對用戶上傳內容的表現方式做保護處理，對用戶上傳分享的文檔內容本身不做任何修改或編輯，并不能對任何下載內容負責。
  • 6. 下載文件中如有侵權或不適當內容，請與我們聯系，我們立即糾正。
  • 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。

評論

0/150

提交評論