當前位置:首頁 » 編程語言 » 牛頓前差公式C語言編程
擴展閱讀
webinf下怎麼引入js 2023-08-31 21:54:13
堡壘機怎麼打開web 2023-08-31 21:54:11

牛頓前差公式C語言編程

發布時間: 2022-08-26 23:17:22

c語言程序翻譯解釋,這是一段程序中關於牛頓插值的程序,請幫忙把每句解釋一下,printf不用解釋。

{
double x[20],fx[20],f[20][20],xx,Nx,t; //定義雙精度的3個一維數組、一個二維數組、3個變數
int i,j,k,n; //定義4個整型變數
printf(" 請輸入插值點的個數n:");
scanf("%d",&n); //點的個數存入n
printf(" 請輸入插值節點值和插值點函數值:\n");
for(i=0;i<=n-1;i++) //接收各個點的值
{
printf("x%d f(%d)分別為:",i,i);
scanf("%lf%lf",&x[i],&fx[i]); //循環接收各個點的x、y的值並存入
}
for(i=0;i<=n-1;i++)
f[i][0]=fx[i]; //把f[20][20]看成一個表,表的第一列按順序存入因變數的各個值
for(j=1;j<=n-1;j++)
for(i=j;i<=n-1;i++)
f[i][j]=(f[i][j-1]-f[i-1][j-1])/(x[i]-x[i-j]); //公式,結合書上的圖和上面一句話。。
for(k=0;k<20;k++)
{
printf("x=");
scanf("%lf",&xx); //用戶輸入的要求解的自變數值
if(xx==9999)
break; //如果自變數值太大了 就跳出去
Nx=f[0][0];
for(i=1;i<=n-1;i++)
{ t=1.0;
for(j=0;j<=i-1;j++) t=t*(xx-x[j]); //又是公式,就是(x-x0)*(x-x1)*...(x-xn)
Nx=Nx+t*f[i][i]; //公式的盡頭
}
printf("N(%f)=%f\n",xx,Nx); //輸出用戶輸入的自變數對應的因變數的值()
}
}
break;
大概這樣

㈡ 解方程,用牛頓法...(用c語言編寫)

#include<stdio.h>
#include<math.h>
int main()
{
float x1,x,f1,f2;static int count=0;
x1=1.5//定義初始值
do
{
x=x1;
f1=x*(2*x*x-4*x+3)-6;
f2=6*x*x-8*x+3;//對函數f1求導
x1=x-f1/f2; count++;
}while(fabs(x1-x)<=1e-5);
printf("%8.7f\n",x1); printf("%d\n",count);
return 0;
}
//2x3-4x2+3x-6//根據我改了初始值,查看結果,表明:改變初始值得到的結果並不一樣,但是迭代的次數並沒有改變!!

㈢ 牛頓演算法和拉格朗日插值演算法的C語言實現

已經編譯運行確認:

#include<stdio.h>
#include<stdlib.h>
#include<iostream.h>
typedef struct data
{
float x;
float y;
}Data;//變數x和函數值y的結構
Data d[20];//最多二十組數據
float f(int s,int t)//牛頓插值法,用以返回插商
{
if(t==s+1)
return (d[t].y-d[s].y)/(d[t].x-d[s].x);
else
return (f(s+1,t)-f(s,t-1))/(d[t].x-d[s].x);
}
float Newton(float x,int count)
{
int n;
while(1)
{
cout<<"請輸入n值(即n次插值):";//獲得插值次數
cin>>n;
if(n<=count-1)// 插值次數不得大於count-1次
break;
else
system("cls");
}
//初始化t,y,yt。
float t=1.0;
float y=d[0].y;
float yt=0.0;
//計算y值
for(int j=1;j<=n;j++)
{
t=(x-d[j-1].x)*t;
yt=f(0,j)*t;
//cout<<f(0,j)<<endl;
y=y+yt;
}
return y;
}
float lagrange(float x,int count)
{
float y=0.0;
for(int k=0;k<count;k++)//這兒默認為count-1次插值
{
float p=1.0;//初始化p
for(int j=0;j<count;j++)
{//計算p的值
if(k==j)continue;//判斷是否為同一個數
p=p*(x-d[j].x)/(d[k].x-d[j].x);
}
y=y+p*d[k].y;//求和
}
return y;//返回y的值
}
void main()
{
float x,y;
int count;
while(1)
{
cout<<"請輸入x[i],y[i]的組數,不得超過20組:";//要求用戶輸入數據組數
cin>>count;
if(count<=20)
break;//檢查輸入的是否合法
system("cls");
}
//獲得各組數據
for(int i=0;i<count;i++)
{
cout<<"請輸入第"<<i+1<<"組x的值:";
cin>>d[i].x;
cout<<"請輸入第"<<i+1<<"組y的值:";
cin>>d[i].y;
system("cls");
}
cout<<"請輸入x的值:";//獲得變數x的值
cin>>x;
while(1)
{
int choice=3;
cout<<"請您選擇使用哪種插值法計算:"<<endl;
cout<<" (0):退出"<<endl;
cout<<" (1):Lagrange"<<endl;
cout<<" (2):Newton"<<endl;
cout<<"輸入你的選擇:";
cin>>choice;//取得用戶的選擇項
if(choice==2)
{
cout<<"你選擇了牛頓插值計算方法,其結果為:";
y=Newton(x,count);break;//調用相應的處理函數
}
if(choice==1)
{
cout<<"你選擇了拉格朗日插值計算方法,其結果為:";
y=lagrange(x,count);break;//調用相應的處理函數
}
if(choice==0)
break;
system("cls");
cout<<"輸入錯誤!!!!"<<endl;
}
cout<<x<<" , "<<y<<endl;//輸出最終結果

}

㈣ 求用c語言編寫牛頓插值法

牛頓插值法:

#include<stdio.h>
#include<alloc.h>
float Language(float *x,float *y,float xx,int n)
{
int i,j;
float *a,yy=0.0;
a=(float *)malloc(n*sizeof(float));
for(i=0;i<=n-1;i++)
{
a[i]=y[i];
for(j=0;j<=n-1;j++)
if(j!=i)a[i]*=(xx-x[j])/(x[i]-x[j]);
yy+=a[i];
}
free(a);
return yy;
}
void main()
{
float x[4]={0.56160,0.5628,0.56401,0.56521};
float y[4]={0.82741,0.82659,0.82577,0.82495};
float xx=0.5635,yy;
float Language(float *,float *,float,int);
yy=Language(x,y,xx,4);
printf("x=%f,y=%f\n",xx,yy);
getchar();
}
‍2.牛頓插值法#include<stdio.h>
#include<math.h>
#define N 4
void Difference(float *x,float *y,int n)
{
float *f;
int k,i;
f=(float *)malloc(n*sizeof(float));
for(k=1;k<=n;k++)
{
f[0]=y[k];
for(i=0;i<k;i++)
f[i+1]=(f[i]-y[i])/(x[k]-x[i]);
y[k]=f[k];
}
return;
}
main()
{
int i;
float varx=0.895,b;
float x[N+1]={0.4,0.55,0.65,0.8,0.9};
float y[N+1]={0.41075,0.57815,0.69675,0.88811,1.02652};
Difference(x,(float *)y,N);
b=y[N];
for(i=N-1;i>=0;i--)b=b*(varx-x[i])+y[i];
printf("Nn(%f)=%f",varx,b);
getchar();
}
留下個郵箱,我發給你:牛頓插值法的程序設計與應用

㈤ C語言編程中,牛頓迭代法是什麼

牛頓迭代法是一種常用的計算方法,這個大學大三應該學過。
具體為:設r是f(x) = 0的根,選取x0作為r初始近似值,過點(x0,f(x0))做曲線y = f(x)的切線L,L的方程為y = f(x0)+f'(x0)(x-x0),求出L與x軸交點的橫坐標 x1 = x0-f(x0)/f'(x0),稱x1為r的一次近似值。過點(x1,f(x1))做曲線y = f(x)的切線,並求該切線與x軸交點的橫坐標 x2 = x1-f(x1)/f'(x1),稱x2為r的二次近似值。重復以上過程,得r的近似值序列,其中x(n+1)=x(n)-f(x(n))/f'(x(n)),稱為r的n+1次近似值,上式稱為牛頓迭代公式。
你把這段文字認真仔細慢慢讀一遍,把給的方程式寫出來,然後照這個在紙上畫出圖形,就會明白牛頓迭代法的概要了。
你講的xopint?root?float?這些都是自己定義的函數。float是c語言中定義浮點型變數的寫法。
#include <iostream>
#include <math.h>
void main()
{
float f(float);
float xpoint(float,float);
float root(float,float);
float x,x1,x2,f1,f2;
do
{
printf("輸入x1,x2\n\n");
scanf("%f%f",&x1,&x2);
f1=f(x1);
f2=f(x2);
}while(f1*f2>0);
x=root(x1,x2);
printf("方程在1.5附近的根為:%f\n\n",x);
}
float f(float x)//定義一個f函數,返回值y
{
float y;
y=2*x*x*x-4*x*x+3*x-6;
return(y);
}
float xpoint(float x1,float x2)//定義一個帶返回值的函數即y,也就是求y的函數,main()中調用
{
float y;
y=(x1*f(x2)-x2*f(x1))/(f(x2)-f(x1));
return(y);
}
float root(float x1,float x2)//這也是定義一個函數,是求根的函數,利用了上面自己定義的函數
{
float x,y,y1;
y1=f(x1);
do
{
x=xpoint(x1,x2);
y=f(x);
if(y*y1>0)
{
y1=y;
x1=x;
}
else
x2=x;
}while(fabs(y)>1e-4);
return(x);
}
建議你看看c 語言教程,上面講得很詳細噢。

㈥ C語言編程牛頓法解方程(只限C語言)!~急~!!~!3Q!

你可以看一下譚浩強變得c教程,裡面有

㈦ 牛頓的插值法用C語言怎麼編寫怎麼編啊

程序代碼如下。
希望能幫助到你!
牛頓插值法
#include<stdio.h>
#include<math.h>
#define
n
4
void
difference(float
*x,float
*y,int
n)
{
float
*f;
int
k,i;
f=(float
*)malloc(n*sizeof(float));
for(k=1;k<=n;k
)
{
f[0]=y[k];
for(i=0;i<k;i
)
f[i
1]=(f[i]-y[i])/(x[k]-x[i]);
y[k]=f[k];
}
return;
}
main()
{
int
i;
float
varx=0.895,b;
float
x[n
1]={0.4,0.55,0.65,0.8,0.9};
float
y[n
1]={0.41075,0.57815,0.69675,0.88811,1.02652};
difference(x,(float
*

㈧ 求 C語言 牛頓迭代法 程序~~~~~~~~~~~~

給你一點提示。
牛頓迭代法要計算
(1)
y1=f(x)

x
的函數值
(2)
d1=f(x)
的一階導數

x
的值
你可以寫兩個函數,分別計算y1,d1
如果一階導數有解析解,則可用賦值語句,否則要寫數值解子程序。
步驟:
設解的精度,例
float
eps=0.000001;
設x初值,x1;
算y1=f(x1);
迭代循環開始
算一階導數

x1
的值
d1
用牛頓公式
算出
x2;
[x2
=
x1
-
y1
/
d1]
如果
fabs(x2-x1)
>
eps
則從新迭代
--
用新的函數值和一階導數值推下一個
新x.

㈨ 急,用C語言表示牛頓柯特斯公式

將Simpson一元函數積分**************
//***************拓展到任意一個牛頓柯特斯積分****

#include<iostream>
#include<math.h>
using namespace std;

//設置全局數組——牛頓 科特斯公式系數表

double C[6][7]={{1.0/2,1.0/2},{1.0/6,4.0/6,1.0/6},{1.0/8,3.0/8,3.0/8,1.0/8},

{7.0/90,16.0/45,2.0/15,16.0/45,7.0/90},

{19.0/288,25.0/96,25.0/144,25.0/144,25.0/96,19.0/188},

{41.0/840,9.0/35,9.0/280,34.0/105,9.0/280,9.0/35,41.0/840}};

int main()

{

double a=0.0,b=0.0,Cotes=0.0;

int n=0;

cout<<"請分別輸入積分段的下屆和上屆:"<<endl;

cin>>a>>b;

cout<<"請輸入您想設置的分段數(節點數-1):"<<endl;

cin>>n;

//檢測輸入

while(!(n>=1&&n<=6))
{

cout<<"分段數最多為6,請重新輸入"<<endl;

cin>>n;

}

//計算科特斯公式的值

for(int j=0;j<=n;j++)

Cotes=Cotes+C[n-1][j]*log((j*(b-a)/n)+a); //函數f(x)為f(x)=ln x 這里你可以自己改為你想要的函數 當然也可以像你的代碼一樣通過函數指針傳遞f(x)表達式~~

Cotes=(b-a)*Cotes;

cout<<"牛頓—柯特斯公式計算積分的結果是"<<Cotes<<endl;

return 0;