數(shù)值逼近課程設(shè)計報告(編輯修改稿)
2025-04-19 03:20 本頁面
【文章內(nèi)容簡介】 B yi yi x0 O xi x17 x36 x xi圖1 凸輪模型 圖2 展開曲線 表1：測量數(shù)據(jù)表xix0x1x2x3x4x5x6x7x8yi525492451xiX0x10x11x12x13x14x15x16x17 x36yixi=ih, x0 =0, x36=, h=。是直線段，AB是曲線段，為了數(shù)控加工，需要計算出圓周上任一點處的柱高,試構(gòu)造算法、設(shè)計程序、編程計算。三、問題解決的方法與算法方法：Newton差值和分段線性插值 Newton插值設(shè)y=f(x)在節(jié)點處的函數(shù)值已知，則關(guān)于點的一階差商為：則關(guān)于點的二階差商為：則關(guān)于點的k階差商為：一階差商二階差商三階差商X0f(X0)X1f(X1)f[X0, X1]X2f(X2)f[X1, X2]f[X0,X1, X2]X3f(X3)f[X2, X3]f[X1, X2,X3]f[X0,X1, X2 X3] 插值公式 插值余項 插值函數(shù)設(shè)y=f(x)在節(jié)點處的函數(shù)值為為了提高近似程度，可以考慮用分段線性插值來逼近原函數(shù)，這時的插值函數(shù)為分段函數(shù)：在區(qū)間上的線性函數(shù)為： 插值公式 插值余項第一步：輸入要進行計算的x的值第二步：輸入各組y[i]及x[i]的值，i=0,1,2,…,n第三步：判斷條件：if(x=0amp。amp。x=16*)則計算各階插商，用牛頓插值進行計算if (x=*17amp。amp。x=)則用分段線性插值進行計算，如果兩種情況均不是則輸出輸入錯誤第四步：計算出柱高S第五步：輸出S四、計算程序// : 定義控制臺應(yīng)用程序的入口點。////include includedefine ARRAY_LEN 18double dq(double a){ int i=0。 int j=0。 double mul=1。 double h=。 double x[ARRAY_LEN]={0}。 double y[ARRAY_LEN]={,525,492,451,,,}。 double temp=0。 for(i=0。iARRAY_LEN。i++) { temp=i*h。 x[i]=temp。 } printf(請輸出所有x的值：\n)。 for(i=0。iARRAY_LEN。i++) { printf(%f ,x[i])。 if(i%3==0) printf(\n)。 } printf(\n)。 printf(請輸出所有y的值：\n)。 for(j=0。jARRAY_LEN。j++) { printf(%lf ,y[j])。 if(j%3==0) printf(\n)。 } printf(\n)。 ///////////////////求差商 //printf(各階插商值為：\n)。 for(i=0。iARRAY_LEN。i++) { for(j=ARRAY_LEN1。ji。j) { y[j]=(y[j]y[j1])/(x[j]x[j1i])。 } } ///////////////////////求和 double s=y[0]。 for(j=1。jARRAY_LEN。j++) { mul=1。 for(i=0。i=j1。i++) { mul=mul*(ax[i])。 } s=s+mul*y[j]。 } return s。}double s(double n){ int i=0。 int j=0。 double s=0。 double x[2]={,}。 double y[2]={,}。 printf(請輸出所有x的值：\n)。 for(i=0。i2。i++) { printf
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號-1