本文目錄一覽:
- 1、用C語言實現拉格朗日插值、牛頓插值、等距結點插值算法
- 2、用C語言編寫一個線性插值程序
- 3、c語言中插值排序法怎麼打程序
- 4、用拉格朗日插值法(二次插值)來計算某個f(x)的值,我用了C語言編,但是編譯的時候很多錯誤,求修改啊!
- 5、兩元lagrange插值如何用C語言表示
- 6、求雙線性插值法的C語言程序!幫幫忙!拜託各位了!
用C語言實現拉格朗日插值、牛頓插值、等距結點插值算法
#includestdio.h
#includestdlib.h
#includeiostream.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次插值):”;//獲得插值次數
cinn;
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;
//coutf(0,j)endl;
y=y+yt;
}
return y;
}
float lagrange(float x,int count)
{
float y=0.0;
for(int k=0;kcount;k++)//這兒默認為count-1次插值
{
float p=1.0;//初始化p
for(int j=0;jcount;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組:”;//要求用戶輸入數據組數
cincount;
if(count=20)
break;//檢查輸入的是否合法
system(“cls”);
}
//獲得各組數據
for(int i=0;icount;i++)
{
cout”請輸入第”i+1″組x的值:”;
cind[i].x;
cout”請輸入第”i+1″組y的值:”;
cind[i].y;
system(“cls”);
}
cout”請輸入x的值:”;//獲得變量x的值
cinx;
while(1)
{
int choice=3;
cout”請您選擇使用哪種插值法計算:”endl;
cout” (0):退出”endl;
cout” (1):Lagrange”endl;
cout” (2):Newton”endl;
cout”輸入你的選擇:”;
cinchoice;//取得用戶的選擇項
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;
}
coutx” , “yendl;//輸出最終結果
}
用C語言編寫一個線性插值程序
#include stdio.h
double Lerp(double x0,double y0,double x1,double y1,double x)
{
double dy = y1 – y0;
if(dy == 0){
printf(“除0錯誤!\n”);
return 0;
}
return x * (x1 – x0) / dy;
}
int main()
{
double x0,x1,y1,y0,x,y;
printf(“Inptu x0 y0 x1 y1 x:”);
scanf(“%lf %lf %lf %lf %lf”,x0,y0,x1,y1,x);
y = Lerp(x0,y0,x1,y1,x);
printf(“y = %lf\n”,y);
return 0;
}
c語言中插值排序法怎麼打程序
(1)“冒泡法”
冒泡法大家都較熟悉。其原理為從a[0]開始,依次將其和後面的元素比較,若a[0]a[i],則交換它們,一直比較到a[n]。同理對a[1],a[2],…a[n-1]處理,即完成排序。下面列出其代碼:
void bubble(int *a,int n) /*定義兩個參數:數組首地址與數組大小*/
{
int i,j,temp;
for(i=0;in-1;i++)
for(j=i+1;jn;j++) /*注意循環的上下限*/
if(a[i]a[j]) {
temp=a[i];
a[i]=a[j];
a[j]=temp;
}
}
冒泡法原理簡單,但其缺點是交換次數多,效率低。
下面介紹一種源自冒泡法但更有效率的方法“選擇法”。
(2)“選擇法”
選擇法循環過程與冒泡法一致,它還定義了記號k=i,然後依次把a[k]同後面元素比較,若a[k]a[j],則使k=j.最後看看k=i是否還成立,不成立則交換a[k],a[i],這樣就比冒泡法省下許多無用的交換,提高了效率。
void choise(int *a,int n)
{
int i,j,k,temp;
for(i=0;in-1;i++) {
k=i; /*給記號賦值*/
for(j=i+1;jn;j++)
if(a[k]a[j]) k=j; /*是k總是指向最小元素*/
if(i!=k) { /*當k!=i是才交換,否則a[i]即為最小*/
temp=a[i];
a[i]=a[k];
a[k]=temp;
}
}
}
選擇法比冒泡法效率更高,但說到高效率,非“快速法”莫屬,現在就讓我們來了解它。
(3)“快速法”
快速法定義了三個參數,(數組首地址*a,要排序數組起始元素下標i,要排序數組結束元素下標j). 它首先選一個數組元素(一般為a[(i+j)/2],即中間元素)作為參照,把比它小的元素放到它的左邊,比它大的放在右邊。然後運用遞歸,在將它左,右兩個子數組排序,最後完成整個數組的排序。下面分析其代碼:
void quick(int *a,int i,int j)
{
int m,n,temp;
int k;
m=i;
n=j;
k=a[(i+j)/2]; /*選取的參照*/
do {
while(a[m]kmj) m++; /* 從左到右找比k大的元素*/
while(a[n]kni) n–; /* 從右到左找比k小的元素*/
if(m=n) { /*若找到且滿足條件,則交換*/
temp=a[m];
a[m]=a[n];
a[n]=temp;
m++;
n–;
}
}while(m=n);
if(mj) quick(a,m,j); /*運用遞歸*/
if(ni) quick(a,i,n);
}
(4)“插入法”
插入法是一種比較直觀的排序方法。它首先把數組頭兩個元素排好序,再依次把後面的元素插入適當的位置。把數組元素插完也就完成了排序。
void insert(int *a,int n)
{
int i,j,temp;
for(i=1;in;i++) {
temp=a[i]; /*temp為要插入的元素*/
j=i-1;
while(j=0tempa[j]) { /*從a[i-1]開始找比a[i]小的數,同時把數組元素向後移*/
a[j+1]=a[j];
j–;
}
a[j+1]=temp; /*插入*/
}
}
(5)“shell法”
shell法是一個叫 shell 的美國人與1969年發明的。它首先把相距k(k=1)的那幾個元素排好序,再縮小k值(一般取其一半),再排序,直到k=1時完成排序。下面讓我們來分析其代碼:
void shell(int *a,int n)
{
int i,j,k,x;
k=n/2; /*間距值*/
while(k=1) {
for(i=k;in;i++) {
x=a[i];
j=i-k;
while(j=0xa[j]) {
a[j+k]=a[j];
j-=k;
}
a[j+k]=x;
}
k/=2; /*縮小間距值*/
}
}
上面我們已經對幾種排序法作了介紹,現在讓我們寫個主函數檢驗一下。
#includestdio.h
/*別偷懶,下面的”…”代表函數體,自己加上去哦!*/
void bubble(int *a,int n)
{
…
}
void choise(int *a,int n)
{
…
}
void quick(int *a,int i,int j)
{
…
}
void insert(int *a,int n)
{
…
}
void shell(int *a,int n)
{
…
}
/*為了打印方便,我們寫一個print吧。*/[code]
void print(int *a,int n)
{
int i;
for(i=0;in;i++)
printf(“%5d”,a[i]);
printf(“\n”);
}
main()
{ /*為了公平,我們給每個函數定義一個相同數組*/
int a1[]={13,0,5,8,1,7,21,50,9,2};
int a2[]={13,0,5,8,1,7,21,50,9,2};
int a3[]={13,0,5,8,1,7,21,50,9,2};
int a4[]={13,0,5,8,1,7,21,50,9,2};
int a5[]={13,0,5,8,1,7,21,50,9,2};
printf(“the original list:”);
print(a1,10);
printf(“according to bubble:”);
bubble(a1,10);
print(a1,10);
printf(“according to choise:”);
choise(a2,10);
print(a2,10);
printf(“according to quick:”);
quick(a3,0,9);
print(a3,10);
printf(“according to insert:”);
insert(a4,10);
print(a4,10);
printf(“according to shell:”);
shell(a5,10);
print(a5,10);
}
用拉格朗日插值法(二次插值)來計算某個f(x)的值,我用了C語言編,但是編譯的時候很多錯誤,求修改啊!
#include stdio.h
float lagrange(float x[],float y[],float xx) // 不寫 i,不要 int n
{
int i,j;
float a[3],yy=0.0; // 定態數組 a[3]
for(i=0;i=2;i++) { a[i]=y[i];
for(j=0;j=2;j++)
if(j!=i) a[i]*=(xx-x[j])/(x[i]-x[j]);
yy+=a[i]; }
return yy; }
void main() {
int i;
// 去掉函數原型聲明
float x[3],y[3],xx,yy; // 定態數組 x[3],y[3]
for(i=0;i=2;i++) {
printf(“x[%d]:”,i); scanf(“%f”,x[i]); }
printf(“\n”);
for(i=0;i=2;i++)
{
printf(“y[%d]:”,i);scanf(“%f”,y[i]);
}
printf(“\n”); printf(“Input xx:”);
scanf(“%f”,xx);
yy=lagrange(x,y,xx);
printf(“x=%f,y=%f\n”,xx,yy);
}
兩元lagrange插值如何用C語言表示
//Lagrange插值多項式
//算法描述:
// 1、輸入:插值節點控制數n,插值點序列(x_i,y_i),i=0,1,…n,要計算的函數點x。
// 2、for(i=0,1,2,…,n)
// {
// temp=1;
// for(j=0,1,…i,i+1,…n)
// !x要事先給定
// temp=temp*(x-x_j)/(x_i-x_j);
// }
// fx=fx+temp*y_i;
// }
#includestdio.h
#includestring.h
#define MAX_n 20
typedef struct tagPOINT
{
double x;
double y;
}POINT;
double Lagrange()
{
int n,i,j;
double x,temp,fx=0;
POINT points[MAX_n];
printf(“Now,please input the n value: \n”);
scanf(“%d”,n);
if(n=1||nMAX_n)
{
printf(“The value of n should be between 2 and %d\n”,MAX_n);
return 1;
}
printf(“Now,please input the (x_i,y_i),i=0,…,%d\n”,n-1);
for(i=0;in;i++)
scanf(“%lf%lf”,points[i].x,points[i].y);
printf(“Now,please input the x value:\n”);
scanf(“%lf”,x);
for(i=0;in;i++)
{
temp=1;
for(j=0;jn;j++)
if(i==j)continue;
else temp=temp*(x-points[j].x)/(points[i].x-points[j].x);
fx=fx+temp*points[i].y;
}
printf(“So,when x=%lf,the Lagrange(%lf)=%lf\n”,x,x,fx);
}
int main()
{
char s[10];
Lagrange();
gets(s);
while(strcmp(s,”exit”))
{
if(strcmp(s,”con”)==0)
{
Lagrange();
getchar();
}
printf(“繼續(輸入con),退出(輸入exit)!\n”);
gets(s);
}
return 0;
}
求雙線性插值法的C語言程序!幫幫忙!拜託各位了!
a b
t
c d
就是兩次線性插值,先在x方向插出t上下方的_t1、_t2,然後再用它們插出t來
float test(float x,float y)
{
float _t1,_t2,t;
_t1 = a+(b-a)*(x-ax)/(bx-ax);
_t2 = c+(d-c)*(x-cx)/(dx-cx);
t = _t1 +(_t2-_t1)*(y – ay);
return t;
}
原創文章,作者:小藍,如若轉載,請註明出處:https://www.506064.com/zh-hant/n/153376.html