Ⅰ 高分悬赏 求三维数据点c语言插值计算程序
问题补充，因字数限制，挪到这
1.拉格朗日插值简介：
对给定的n个插值节点x1,x2,…,xn,及其对应的函数值y1=f(x1), y2=f(x2),…, yn=f(xn);使用拉格朗日插值公式，计算在x点处的对应的函数值f(x);
2.一维拉格朗日插值c语言程序：
Int lagrange(x0, y0, n, x, y)
Float xo[], yo[], x;
Int n;
Float *y
{
Int i, j;
Float p;
*y=0;
If (n>1)
{
For(i=0;i<n;i++)
{
P=1;
For(j=1;j<n;j++)
{
If(i!=J)
P=p*(x-x0[j]/x0[i]-x0[j]);
}
*y=*y+p*y0[i];
Return(0);
}
Else
Return(-1);
}
3.例题。已知函数如下表所示，求x=0.472处的函数值：
X 0.46 0.47 0.48 0.49
Y 0.484655 0.4903745 0.502750 0.511668
计算这个问题的c语言程序如下：
#minclude stdio
#includeM<nath.h>
Main()
{
Float x0[4]={ 0.46, 0.47,0.48,0.49};
Float y0[4]={ 0.484655 ,0.4903745 ,0.502750 ,0.511668};
Float x, y;
Int n, rtn;
N=4;
X=0.472;
Rth=lagrange(x0,y0,n,x,&y);
If(rtn=0)
{
Prinf(“Y(0.472)=:%f\n”,y);
}
Else
{
Prinf(“n must be larger than 1.\n”);
}
}
计算结果：Y(0.472)=:0.495553
4.问题补充
我的问题与上面的例子类似，计算三维空间一点(x,y,z)对应的函数值（Vx,Vy,Vz）.不同的是自变量(point_coordinate.txt)为三维空间散乱点（不是正方体的顶点），因变量(point_data.txt)为矢量（向量 ）。插值算法比较多，常数法，拉格朗日插值，埃特金插值，三阶样条插值等。最简单的就是常数法，查找离目标点（x,y,z）距离最近的已知自变量(Xi,Yi,Zi)，把该点的函数值赋给目标点做函数值，求高手帮忙写写。
Ⅱ c语言 函数名: bar3d 怎么用
bar3d是Turbo C工具中绘制长方体（也可理解为绘制三维条形图）的一个图形函数。其基本语法如下：
bar3d(int left, int top, int right, int bottom, int depth, int topflag)；
其中：int left, int top构成第一个坐标点，即屏幕中长方体前面（即离我们最近的长方形）左上角的一个点；
int right, int bottom构成第二个坐标点，即屏幕中长方体前面（即离我们最近的长方形）右下角的一个点；
由以上两个坐标点就可以确定一个长方形，然后结合int depth（用于确定长方体的厚度，即长方体前后两个面的间距），就可以完全确定长方体的形状。
int topflag 用于对长方体的顶部显示方式。
举例：例如在Turbo C++工具中编写的C程序：
# include <graphics.h>
# include <stdio.h>
void main()
{
int graphdriver=DETECT, graphmode;
initgraph(&graphdriver,&graphmode,"C:\TURBOC3\BGI"); //*图形系统初始化*/
bar3d(0,0,250,250,10,1);
getch() ;
closegraph(); // 关闭图形模式//
}
其运行结果为：
运行结果
Ⅲ 求一个C语言图形变换程序，先进行图形绘制，后能进行图形复制、移动、缩放、旋转，跪求啦。。。
从这一部分开始，进入了图形编程的比较烦琐的部分，要真正对图形编程有所了解，这一部分的内容是必须要掌握的。
在计算机绘图过程中，经常需要进行绘图变换，主要包括二维图形变换和三维图形变换。这一部分讨论二维图形变换，其内容有用户坐标到屏幕坐标的变换、图形的比例变换、对称变换、错切变换、旋转变换、平移变换和复合变换等。后面讲到了二维剪裁，即线段裁剪与多边形裁剪。
第一节 用户坐标到屏幕坐标变换
假设纸上有一个图形，要用计算机把它在屏幕上画出来。那么首先遇到的问题是，纸上的图形采用的坐标是实数域域中的直角坐标系或是极坐标系，统称为用户坐标系。而屏幕上采用的坐标系是整数域中直角坐标系，这类坐标系统称为设备坐标系。因此用户坐标系中图形需要经过变换才能绘制在屏幕上，显然这个变换的内容包括： 1)将用户坐标系中任意范围区域转换到屏幕某个范围区域，从而用户坐标系此范围区域内的图形也转换到屏幕上该范围区域内。 2)用户坐标系此区域内图形上的坐标值转换到屏幕上该范围区域内后不一定是整数，取整后才成为该范围区域内的屏幕坐标值。 3)用户坐标右手系到屏幕坐标左手系的坐标轴方向变换。 4)当屏幕坐标系水平方向与垂直方向刻度不等(即像素间距不等)时，为保持图形不走样，还要进行比例变换。下面介绍这些内容的具体计算问题。
1.窗口到视口的变换
更确切地说，是实际图形到屏幕图形的转换。有时也称为数据规格化。
在用户坐标系中，指定一矩形域以确定要显示(或绘制)的图形部分，这个矩形区域称为窗口。在屏幕上可任选一矩形域以显示(或绘制)窗口内的图形，该域称为视口。如图2-1所示。
一般视窗口的四条边界分别为：
左边界 x=x1、右边界 x=x2.下边界 y=y1,上边界y=y2。
视口的四条边界分别为：
左边界sx=sx1,右边界sx=sx2,上边界sy=sy1,下边界sy=sy2。
经变换后应有，窗口的上边界线段(或下边界线段)长x2-x1变换成视口上边界线段(或下边界线段)长sx2-sx1。设其比例变换因子为k1，则可得
k1*(x2-x1)=sx2-sx1
k1=(sx2-sx1)/(x2-x1)
对窗口内任一x坐标(x1<=x<=x2)变换后为视口内水平方向sx坐标(sx1<=sx<=sx2)。由上述有：
k1*(x-x1)=sx-sx1
sx=sx1+k1*(x-x1)
=sx1+(x-x1)*(sx2-sx1)/(x2-x1)
同样，经变换后窗口的左边界线段(或右边界线段)长y2-y1变换成视口左边界线段(或右边界线段)长sy2-sy1。设其比例变换因子为k2，则可得
k2*(y2-y1)=sy2-sy1
k2=(sy2-sy1)/(y2-y1)
对窗口内任一y坐标(y1<=y<=y2)变换后为视口内垂直sy坐标(sy1<=sy<=sy2)，应有
k2*(y-y1)=sy-sy1
sy=sy1+k2*(y-y1)
=sy1+(y-y1)*(sy2-sy1)/(y2-y1)
于是对窗口内图形上任一点坐标(x,y)变换到屏幕上视口内成为(sx,sy)，则
sx=sx1+(x-x1)*(sx2-sx1)/(x2-x1)
sy=sy1+(y-y1)*(sy2-sy1)/(y2-y1)
写成简式
sx=k1*x+a
sy=k2*y+b
这里
a=sx1-k1*x1
b-sy1-k2*y1
k1=(sx2-sx1)/(x2-x1)
k2=(sy2-sy1)/(y2-y1)
2. 实型值到整型值的变换
上面对窗口内图形上任一点坐标(x,y)变换到屏幕上视口内成为(sx,sy)，
sx=k1*x+a
sy=k2*y+b k1,k2,a,b同上
这样计算出来的sx,sy一般是实型值，而屏幕上视口内屏幕坐标是整型值，因此要将sx,sy实型值转换成屏幕坐标系的整型值。这可以通过四舍五入的方法将实型值的绝对值圆整化。由于C语言中已经替我们想到了这点，它提供的函数可以自动取整，因此用户在调用标准函数在屏幕上绘图时一般不需要考虑这个问题。当然也可以用赋值的类型转换规则来实现实型值到整型值的变换。
3. y坐标值方向变换
一般屏幕坐标系是直角左手系，y轴方向向下为正，原点在屏幕的左上角，如图2-2所示。
窗口内图形上任一点(x,y)变换到视口内成为(sx,xy)，而(x,y)是相对用户坐标系(直角右手系)的。(sx,sy)是相对屏幕坐标系(直角左手系)的，因此y轴方向相反。为使窗口内图形变换到视口上图形其形状一致，需将视口上图形y轴方向变换成窗口内图形y轴方向。这只要将求得的视口内各点的sy整型坐标均用sy2去减，即sy2-sy(整型)代替sy(整型)即可，经这样的坐标轴方向变换后得到的视口内图形与窗口内图形一致。
4.长宽比例变换
屏幕坐标系x方向与y方向上的刻度可能不一样，这取决于水平方向像素间距与垂直方向偈素间距大小是否一致。如果两个方向的刻度不相等，那么用户坐标系下一个正方形将显示(或绘制)成为一个长方形有，一个圆将成为一个椭圆。
为保持原图形的长宽比。使图形显示(或绘制)后不走样，需求出屏幕上两侍标轴刻度的比值(即纵横比)。可以用函数getaspectratio()(见前文所述)返回x方向和y方向的比例数，从而求得这个比值。再瘵原图形y方向坐标乘以该比值，这样显示(或绘制)出来的图形应不走样。若不考虑图形的走样，就不必作这个变换。
第二节 二维几何变换
图形的几何变换一般是指对图形的几何信息经过变换后产生新的图形，图形几何变换既可以看作是坐标系不动而图形变动，变动后的图形在坐标系中的坐标值发生变化；出可以看作图形不动而坐标系变动，变动后的图形在新坐标系下具有新的坐标值。这两种情况本质上都是一样的，都是图形由新的坐标值表示，因此是新产生的图形。图形几何变换包括比例变换、对称变换、错切变换、旋转变换、平移变换及其复合变换。图形上所有的点在几何变换前后的坐标关系一般用解析几何方法可以求得，但这些几何关系用矩阵方法表示，运算更为方便。
一、基本变换
图形基本几何变换是指比例变换、对称变换、错切变换、旋转变换和平移变换等。除平移变换外，这里其它四种几何变换都可以用组成图形的点向量(或称1×2阶矩阵)和2×2阶变换矩阵相乘表示，而平移变换需引入新方法来实现。
1、比例变换
设图形上一点P(x,y)，经比例变换后成为新的菜上一点P'(x',y')，即有
x'=a*x
y'=d*y
式中a,d为比例因子
将此比例变换式写成矩阵式得
a 0
[x' y']=[x y] = [x y] * T
0 d
a 0
这里 T= 叫做比例变换矩阵。若a=d，则x,y坐标按同一比例变换。
0 d
当a=d>1时，图形放大；当0<a=d<1时，图形缩小。
若a≠d，则x,y坐标按各自不同比例变换。
3 0
例 1： 设有比例变换矩阵 T= , 三角形abc经过比例变换成为三角形a'b'c'。
0 1
如图2-3所示。
3 0
a [1 2] = [3 2] a'
0 1
3 0
b [2 2] = [6 2] b'
0 1
3 0
c [2 3] = [6 3] c'
0 1
2. 对称变换
图形上一点P(x,y)经关于原点对称变换后成为新图形上一点P'(x',y')，则
x' = -x
y' = -y
写成矩阵形式成为
-1 0
[x' y'] = [x y] = [x y] * T
0 -1
-1 0
这里 T = 为关于原点对称变换矩阵。
0 -1
若关于x轴对称，则对称变换的矩阵表示为
1 0
[x' y'] = [x y] = [x y] * T
0 -1
1 0
于是关于x轴对称变换矩阵 T =
0 -1
若关于y轴对称，则对称变换的矩阵表示为
-1 0
[x' y'] = [x y] = [x y] * T
0 1
-1 0
于是关于y轴对称变换矩阵 T =
0 1
若关于直线y = -x对称，则对称变换矩阵表示为
0 -1
[x' y'] = [x y] = [x y] * T
-1 0
0 1
于是关于直线 y = x对称变换矩阵 T =
1 0
各种对称变换的图形均可由实例程序绘出，参见实例程序图形。
3. 错切变换
对图形的任一点P(x,y)，作线性变换如下
x' = x + by
y' = y + dx
式中b,d为不全为零的常 数，点P'(x',y')为新图形上相应的点，这个变换称为图形的错切变换。
错切变换的矩阵表示为
1 d
[x' y'] = [x y] = [x y] * T
b 1
1 d
T = 叫做错切变换矩阵(b,d不全为零)。
b 1
① 当d=0时，x'=x+by,y'=y，这时图形的y坐标不变，x坐标值随(x,y)及系数b作线性变化。若b>0时，图形沿x轴作错切位移；若b<0，图形沿x轴负向作错切位移。
② 当b=0时，x'=x,y'=dx+y，此时图形的x坐标不变y坐标随(x,y)及系数d作线性变化。如d>0，图形沿y轴正向作错切位移；如d<0，图形沿y轴负向作错切位移。
③ 当b≠0且d≠0时，x'=x+by,y'=y+dx，图形沿x,y两个方向作错切位移。
1 2
例 2： 设有错切变换 矩阵 T = ，正方形abcd经此错切变换成为四边形a'b'c'd'，
0 1
如图2-4所示。
1 2
a [0 0] = [0 0] a'
0 1
1 2
b [1 0] = [1 2] b'
0 1
1 2
c [1 1] = [1 3] c'
0 1
1 2
d [0 1] = [0 1] d'
0 1
4. 旋转变换
设图形上一点P(x,y)绕原点逆时针旋转θ角后成为新的图形上一点P'(x',y')，则由解析几何方法可得
x' = xcosθ + ysinθ
y' = -xsinθ + ycosθ
用矩阵表示为
cosθ -sinθ
[x' y'] = [x y] = [x y] * T
sinθ cosθ
cosθ -sinθ
这里 T = 为绕原点逆时针变换矩阵。若顺时针旋转时，θ角为负值。
sinθ cosθ
5. 平移变换
若图形上一点P(x,y)沿x轴平移 l距离，沿y轴平移m距离后成为新的图形上一点P'(x',y')，则有
x' = x + l
y' = y + m
式中l,m不全为零，这称为平移变换。但此变换无法用组成图形的点向量和2×2阶变换矩阵相乘来实现。
用二维点向量和2×2阶矩阵相乘不能表示图形的平移变换，那么自然会想到用三维点向量和3×3阶矩阵相乘来实现图形的平移变换。因此对图形上二个坐标的点向量需要添加一个坐标，使之成为三维点向量以便与三阶矩阵相乘，进而实现用矩阵表示平移变换。实际上就是对上面的二个坐标变换式添加第三个坐标变换式，即成为
x' = x + l
y' = y + m
k = k
这第三个坐标变换式(即k=k)必须是恒等式，因为不需作变换，本质上是为了进行矩阵运算而引入的。
将此三个变换式(仍然是图形的平移变换，不妨将k = k取成1=1)写成矩阵得
1 0 0
[x' y' l] = [x y l] 0 1 0 = [x y 1] * T
l m 1
1 0 0
显然 T = 0 1 0 为图形的平移变换矩阵。
l m 1
这里通过对原图形上二维点向量引进第三个坐标成为三维点向量，从而使原图形的平移变换 能用矩阵表示。同样其它基本变换也可以如此用矩阵表示。因此图形的基本变换都可以在这样的三维点向量下统一、整齐用矩阵表示。这样的三维点向量称为齐次点向量，也叫三维齐次坐标点，简称三维齐次坐标。只有在三维齐次坐标下，二维几何变换才都可以用矩阵表示。下面再进一步讨论一下齐次坐标的优点。
引用齐次坐标后，可将上面各种基本变换矩阵统一在一个三阶矩阵中。即
a b 0
T = c d 0
l m 1
式中左上角二阶矩阵实现比例、对称、错切、旋转等变换，左下角1×2阶矩阵实现平移变换，其中a,b,c,d,l,m只要赋以相应的值，并建立图形上点的齐次坐标(即在图形上点的坐标后引入第三个坐标1)，这样就可以用图形上点的三维齐次坐标与此三阶矩阵相乘来表示三维图形的基本几何变换了。而变换后，不用考虑第三个坐标1，前面两个坐标就反映了图形的整个变换情况。
用齐次坐标表示一个图形上的点，可以有多种表示，如(6,8，1)、(12,16,2)、(30,40,5)等均表示图形上同一个点(6,8)。这样，齐次坐标可以表示计算机无法容纳的数。例如当计算机的字长为16位时，它能表示的最大整数为216-1=32767。若点坐标为(80 000,40 000)，则计算机无法表示。但用齐次坐标可表示为(20 000,10 000,1/4)，经过处理后再用第三个坐标支除前面两个坐标，从而得到原来通常的坐标。
齐次坐标优点很多，在计算机绘图中都采用这种表示来处理图形。下面介绍的图形复合几何变换就是如此。
二、复合变换
图形的复合几何变换是指图形作一次以上的基本几何变换，变换结果是每次基本变换矩阵的乘积。图殂的复合几何变换简称复合变换。
1. 复合平移
若对图形首先作平移变换 T1，然后再作平移变换T2，相应的平移变换矩阵分别为
1 0 0
T1 = 0 1 0
l1 m1 1
1 0 0
T2 = 0 1 0
l2 m2 1
则变换结果为复合平移变换T，其复合平移变换矩阵为
T = T1 * T2
1 0 0 1 0 0
= 0 1 0 * 0 1 0
l1 m1 1 l2 m2 1
1 0 0
= 0 1 0
l1+l2 m1|m2 1
Ⅳ 如何用C语言实现三维
最简单的做法是在结构体里存函数指针，然后初始化结构体的时候将函数的地址赋值给它。
复杂一点的，就得按照c++的多态原理，用虚函数表（其实就是能记录函数名对应的函数地址就好）记录每种结构体的成员函数地址，然后每个结构体里多存一个虚函数表的地址。。。然后调用函数的时候查到函数地址强转成函数类型然后调用。。。这个可以自己研究一下c++的虚函数表机制。
Ⅳ 用c语言来画1个正方体，一个简单的程序。。
void far bar3d(int x1, int y1, int x2, int y2,int depth,int topflag);当
topflag为非0时, 画出一个三维的长方体。当topflag为0时，三维图形不封顶，
实际上很少这样使用。
void far setfillstyle(int pattern, int color); color的值是当前屏幕图形
模式时颜色的有效值,SOLID_FILL 1 以实填充
void far floodfill(int x, int y, int border);
其中：x, y为封闭图形内的任意一border为边界的颜色，也就是封闭图形轮廓的
颜色。调用了该函数后，将用规定的颜色和图模填满整个封闭图形。
#include<stdlib.h>
#include<graphics.h>
main()
{
int gdriver, gmode;
struct fillsettingstype save;
gdriver=DETECT;
initgraph(&gdriver, &gmode, "");
setbkcolor(BLUE);
cleardevice();
setcolor(LIGHTRED);
setlinestyle(0,0,3);
setfillstyle(1,14); /*设置填充方式*/
bar3d(100,200,400,350,200,1); /*画长方体并填充*/
floodfill(450,300,LIGHTRED);
/*填充长方体另外两个面*/
floodfill(250,150, LIGHTRED);
getch();
closegraph();
}