插值法c语言程序例子

学号2131388 姓名 范宇超

程序：

#include<stdio.h>

#include<math.h>

#define N 6

float sub(float a[],float b[],float x,float e); void main(void)

{

float u[N]={100,121,144,169,196,225}; float v[N]={10,11,12,13,14,15}; float x,y,e,*p1,*p2;

printf("Input number x E=:"); scanf("%f%e",&x,&e);

p1=u;

p2=v;

y=sub(p1,p2,x,e);

printf("y=%f\n",y);

}

float sub(float *pp1,float *pp2,float x,float e) {

float a[N],b[N],t[N],y,y1,c; int i,k;

for(i=0;i<N;i++,pp1++) {

a[i]=*pp1;

printf("%12.6f",a[i]); }

#include #include

#include #define N 10

using namespace std; int n; //定义n为全局变量

double lagrange(double a[],double b[],double x ); int main() {

double x; double y; double a[N]; double b[N];

cout<<\输入x0,y0的个数\ cin>>n;

cout<<\输入x0的值\ for(int i=0;i<=n-1;i++) cin>>a[i];

cout<<\输入y0的值\ for(int i=0;i<=n-1;i++) cin>>b[i];

cout<<\输入x的值\ cin>>x;

y=lagrange(a,b,x);

cout<

getch();}

double lagrange(double a[],double b[],double x) {

do

插值法的程序实现

一 实验目的

1.熟悉Matlab编程；

2.学习插值方法及程序设计算法

二 实验题目

分别用拉格朗日插值、牛顿插值、自然样条函数对1910、1965、2002的

人口进行估算。

三 实验原理与理论基础

1.拉格朗日插值算法设计

①利用已知条件得到xi，yi，i=0，1，2，… ②由Lk(x)=((x-x0)*…*(x-x(k-1))*(x-x(k+1))…(x-xn))/ ((xk-x0)*…*(xk-x(k-1))*(xk-x(k+1))…(xk-xn))得出Li(x);

③由Y=y1* L1(x)+…+yn*Ln(x)得出Y关于x的表达式。 ④带值计算即可。 2. 牛顿插值算法设计

①利用已知条件得到xi，yi，i=0，1，2，… ②利用差商公式

f[x0,…xk]=(f[x0,…,x(k-2),xk]-f[x0,…,x(k-1)])/(xk-x(k-1))各阶差商。

③利用牛顿插值公式

f(x)=f(x0)-f[x0,x1]*(x-x0)+…f[x0,x1,…xn]*(x-x0)*…(x-x(n-1)).

数值分析,拉格朗日插值法C语言的实现

实验 一 .拉格朗日插值法C语言的实现

1.实验目的：

进一步熟悉拉格朗日插值法。

掌握编程语言字符处理程序的设计和调试技术。

2.实验要求：

已知：某些点的坐标以及点数。

输入：条件点数以及这些点的坐标 。

输出：根据给定的点求出其对应的拉格朗日插值多项式的值 。

3.程序流程：

（1）输入已知点的个数；

（2）分别输入已知点的X坐标；

（3）分别输入已知点的Y坐标；

（4）通过调用函数lagrange函数，来求某点所对应的函数值。

拉格朗日插值多项式如下：

Ln(xj) yklk(xj) yjj 0,1,……n k 0n

其中lk(x)

(x x0)……(x-xk-1)(x-xk+1) …(x-xn)(xk x0)……(xk-xk-1)(xk-xk+1) …(xk-xn)k 0,1,……,n

程序流程图：

数值分析,拉格朗日插值法C语言的实现

↓

程序如下：

#include <iostream>

#include <conio.h>

#include <malloc.h>

float lagrange(float *x,float *y,float xx,int n) /*拉格朗日插值算法*/

C

语言程序设计 —基础篇

第 7 章 简单程序设计之二

——基础篇

本章内容提要

多分支选择结构—switch语句 循环嵌套 break语句和continue语句 程序举例

2014年4月11日

C语言程序设计

——基础篇

教学目标在三种基本结构的基础上，结合本章介绍流 程 控 制 的 其 他 几 个 语 句 ： switch 语 句 、 break语句和continue语句及循环的嵌套， 实现复杂的算法，设计出结构清晰、功能强大 的程序，对C语言程序开发设计有比较深刻的 认识。

2014年4月11日

C语言程序设计

7.1 多分支选择结构程序设计1、语句形式：

——基础篇

switch (表达式) { case 常量1: 语句1; 数值型 或 case 常量2: 语句2; 字符型 case 常量3: 语句3; ....... 常量表达式的值必须 case 常量n: 语句n; 互不相同，否则执行 default : 语句n+1;时将出现矛盾，即同 一个开关值，将对应 } 多种执行方案。 2014年4月11日 C语言程序设计

根据表达式的不同值 ，选择不同的程序分 千万不 支，又称开关语句。 能接(;) 号

——基础篇

2、语句执行流程：

计算表达式值

=简单语句或 复合

数值逼近课程的内容,用C语言实现的拉格朗日插值方法,可以随时增加一个点,也可以随时算出某个点的值,还可以算出插值多项式,功能很全

#include<stdio.h>

#include <stdlib.h>

#define M 20;

int n=0;

int p=1;

int num=0;

double *x;

double *y;

double Calculate(double tt) ;

void Insert(int m);

void Print( );

void NewTon(int m)

{

double tt;

Insert(m);

Print( );

printf("是否继续进行插值、计算还是结束？继续插值请输入1，结束请输入0，求值计算请输入2；p=");

scanf("%d",&p);

printf("\n");

while(p!=0)

{

if(p==1)

{

printf("请输入再次插值点个数num=");

scanf("%d",&num);

NewTon(num);

}

else if(p==2)

{

printf("请输入x=")

MATLAB拉格郎日插值法与牛顿插值法构造插值多项式

姓名：樊元君 学号：2012200902 日期：2012.10.25

1.实验目的：

掌握拉格郎日插值法与牛顿插值法构造插值多项式。

2.实验内容：

分别写出拉格郎日插值法与牛顿插值法的算法，编写程序上机调试出结果，要求所编程序适用于任何一组插值节点，即能解决这一类问题，而不是某一个问题。实验中以下列数据验证程序的正确性。 已知下列函数表

求x=0.5635时的函数值。

MATLAB拉格郎日插值法与牛顿插值法构造插值多项式

3.程序流程图：

● 拉格朗日插值法流程图：

MATLAB拉格郎日插值法与牛顿插值法构造插值多项式

●牛顿插值法流程图：

MATLAB拉格郎日插值法与牛顿插值法构造插值多项式

4.源程序：

● 拉格朗日插值法：

function [] = LGLR(x,y,v)

x=input('X数组=：');

y=input('Y数组=');

v=input('插值点数值=：');

n=length(x);

u=0;

for k=1:n

t=1;

for j=1:n

if j~=k

t=t*(v-x(j))/(x(k)-x(j));

end

end

u=u+t*y(k);

end

disp('插值结果=');

经典C语言程序设计实例讲解例子 【程序1】

题目:有1、2、3、4个数字,能组成多少个互不相同且无重复数字的三位数?都是多少?

1.程序分析:可填在百位、十位、个位的数字都是1、2、3、4。组成所有的排列后再去

掉不满足条件的排列。 2.程序源代码: main( { int i,j,k; printf(\

for(i=1;i<5;i++/*以下为三重循环*/ for(j=1;j<5;j++ for (k=1;k<5;k++ {

if (i!=k&&i!=j&&j!=k /*确保i、j、k三位互不相同*/ printf(\}

}

【程序2】

题目:企业发放的奖金根据利润提成。利润(I低于或等于10万元时,奖金可提10%;利润高

于10万元,低于20万元时,低于10万元的部分按10%提成,高于10万元的部分,可可提

成7.5%;20万到40万之间时,高于20万元的部分,可提成5%;40万到60万之间时高于

40万元的部分,可提成3%;60万到100万之间时,高于60万元的部分,可提成1.5%,高于

100万元时,超过100万元的部分按1%提成,从键盘输入当月利润I,求应发放奖金总数?

1.程序分析:请利用数轴来分界,定位。注意定义时需把奖金定义成长整型。 2.程序源代码: main( { long int i;

int bonus1,bonus2,bo

顺序结构程序设计

1、已知三角形的底边为20，高为10，求其面积。 main() {float a,h,s; a=20;

h=10;

s=1.0/2*a*h; /*1.0/2不能写作1/2*/ printf(\}

2、从键盘输入三角形的底边及高的长度，求其面积。 main() {float a,h,s;

scanf(\

s=1.0/2*a*h; /*1.0/2不能写作1/2*/ printf(\}

3、已知圆的半径为10，求圆的面积。 解法1：

#define PI 3.14159 main() {float r,s;

r=10; s=PI*r*r;

printf(\}

解法2：

#include #define PI 3.14159 main() {float r,s; r=10;

s=PI*pow(r,2); printf(\}

4、从键盘输入圆的半径值，求圆的面积。 #define PI 3.14159 main() {float r,s;

scanf(\ s=PI*r*r;

printf(\

}

5、从键盘输入球体的半径，求其体积和表面积。 解法1：

#define PI 3.1

插值法计算实际利率

插 值 法 计 算 实 际 利 率

“插值法”计算实际利率。在08年考题中涉及到了实际利率的计算，其原理是根据比例关系建立一个方程，然后，解方程计算得出所要求的数据。

例如：假设与A1对应的数据是B1，与A2对应的数据是B2，现在已知与A对应的数据是B，A介于A1和A2之间，即下对应关系：

A1 B1

A(?) B

A2 B2

则可以按照（A1-A)/(A1-A2)=(B1-B)/(B1-B2)计算 得出A的数值，其中A1、A2、B1、B2、B都是已知数据。根本不必记忆教材中的公式，也没有任何规定必须B1＞B2 验证如下： 根据：（A1-A)/(A1-A2)=(B1-B)/(B1-B2)可知：

（A1-A)＝(B1-B)/(B1-B2)×（A1- A2）

A＝A1－(B1-B)/(B1-B2)×（A1-A2）＝A1＋(B1-B)/(B1- B2)×（A2-A1） 考生需理解和掌握相应的计算。

例如：

某人向银行存入5000元，在利率为多少时才能保证在未来10年中每年末收到750元？ 5000/750=6.667 或 750*m=5000

查年金现值表，期数为10，利率i=8%时，系数为6.710；利率