c語言二元插值,二維線性插值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

(0)
打賞 微信掃一掃 微信掃一掃 支付寶掃一掃 支付寶掃一掃
小藍的頭像小藍
上一篇 2024-11-14 03:05
下一篇 2024-11-14 03:05

相關推薦

  • AES加密解密算法的C語言實現

    AES(Advanced Encryption Standard)是一種對稱加密算法,可用於對數據進行加密和解密。在本篇文章中,我們將介紹C語言中如何實現AES算法,並對實現過程進…

    編程 2025-04-29
  • 學習Python對學習C語言有幫助嗎?

    Python和C語言是兩種非常受歡迎的編程語言,在程序開發中都扮演着非常重要的角色。那麼,學習Python對學習C語言有幫助嗎?答案是肯定的。在本文中,我們將從多個角度探討Pyth…

    編程 2025-04-29
  • Python被稱為膠水語言

    Python作為一種跨平台的解釋性高級語言,最大的特點是被稱為”膠水語言”。 一、簡單易學 Python的語法簡單易學,更加人性化,這使得它成為了初學者的入…

    編程 2025-04-29
  • OpenJudge答案1.6的C語言實現

    本文將從多個方面詳細闡述OpenJudge答案1.6在C語言中的實現方法,幫助初學者更好地學習和理解。 一、需求概述 OpenJudge答案1.6的要求是,輸入兩個整數a和b,輸出…

    編程 2025-04-29
  • Python按位運算符和C語言

    本文將從多個方面詳細闡述Python按位運算符和C語言的相關內容,並給出相應的代碼示例。 一、概述 Python是一種動態的、面向對象的編程語言,其按位運算符是用於按位操作的運算符…

    編程 2025-04-29
  • Python實現一元線性回歸模型

    本文將從多個方面詳細闡述Python實現一元線性回歸模型的代碼。如果你對線性回歸模型有一些了解,對Python語言也有所掌握,那麼本文將對你有所幫助。在開始介紹具體代碼前,讓我們先…

    編程 2025-04-29
  • Python語言由荷蘭人為中心的全能編程開發工程師

    Python語言是一種高級語言,很多編程開發工程師都喜歡使用Python語言進行開發。Python語言的創始人是荷蘭人Guido van Rossum,他在1989年聖誕節期間開始…

    編程 2025-04-28
  • Python語言設計基礎第2版PDF

    Python語言設計基礎第2版PDF是一本介紹Python編程語言的經典教材。本篇文章將從多個方面對該教材進行詳細的闡述和介紹。 一、基礎知識 本教材中介紹了Python編程語言的…

    編程 2025-04-28
  • Python語言實現人名最多數統計

    本文將從幾個方面詳細介紹Python語言實現人名最多數統計的方法和應用。 一、Python實現人名最多數統計的基礎 1、首先,我們需要了解Python語言的一些基礎知識,如列表、字…

    編程 2025-04-28
  • Python作為中心語言,在編程中取代C語言的優勢和挑戰

    Python一直以其簡單易懂的語法和高效的編碼環境而著名。然而,它最近的發展趨勢表明Python的使用範圍已經從腳本語言擴展到了從Web應用到機器學習等廣泛的開發領域。與此同時,C…

    編程 2025-04-28

發表回復

登錄後才能評論