插值项目

最新推荐文章于 2024-09-12 19:02:20 发布
最新推荐文章于 2024-09-12 19:02:20 发布
阅读量208 收藏 1
点赞数
文章标签： 算法
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qzl19/article/details/107795529
版权

插值项目

一、项目说明

1、面向一维插值和线性拟合，含有牛顿插值、线性插值等四种插值方式；
2、设计了可用鼠标直接操作的GUI界面，比上次贪吃蛇项目的GUI更好；

二、代码说明

1、头文件

#include <graphics.h>              // 引用图形库头文件
#include <conio.h>                 //控制台输入输出
#include <stdio.h>                
#include <windows.h>				//用到了定时函数sleep()
#include <math.h>                  //数学函数库
#include <string.h>                //字符串处理函数

2、插值部分（只解释代码，原理部分的相关资料见文尾）

void nuton()//牛顿插值
{ 
   event=0;
   cleardevice();//清屏（取决于背景色）
   int k,m;
   //算出插商表
   for(int i=0;i<non;i++){
      numc[i][0]=numy[i];
   }
   for(k=1;k<non;k++)
	   for(m=k;m<non;m++)
	   {
		   numc[m][k]=(numc[m-1][k-1]-numc[m][k-1])/(numx[m-k]-numx[m-k+1]);
	   }
   for(int i=0;i<nan;i++)
   {
	   float t=1.0;
	    numry[i]=numc[0][0];
	   for(int j=1;j<=non;j++)
	   {   
		   t=t*(numr[i]-numx[j-1]);
		   numry[i]=numry[i]+numc[j][j]*t;
	   }
   }
   //此后为配合图形界面输出结果
   char msg[40]="所求的插值为";
   for(int i=0;i<nan;i++)
   {  
	   char temp[30];
	   sprintf(temp,"%f",numry[i]);
	   sprintf(msg,"%s%s ",msg,temp);
   }
   MessageBox(NULL,TEXT(msg),TEXT("结果"),MB_OK);
   exit(0);
};

可以直接通过牛顿插值公式写出代码。

void langrage()//拉格朗日插值
{
   cleardevice();//清屏（取决于背景色）
   int k,m;
   float a[30];
   for(int i=0;i<nan;i++){
      for(k=0;k<non;k++)
       {   
	      a[k]=numy[k];
	       for(m=0;m<non;m++)
	      {
		   if(m!=k) a[k]=a[k]*((numr[i]-numx[m])/(numx[k]-numx[m]));
	      }
	      numry[i]+=a[k];
        }
   }
   char msg[40]="所求的插值为";
   for(int i=0;i<nan;i++)
   {   char temp[30];
	   sprintf(temp,"%f",numry[i]);
	   sprintf(msg,"%s%s ",msg,temp);
   }
   MessageBox(NULL,TEXT(msg),TEXT("结果"),MB_OK);
   exit(0);
};

void fenduan()//分段线性
{
	for(int i=0;i<nan;i++)
	{
		for(int j=0;j<non;j++)
		{
			if(numr[i]>numx[j]&&numr[i]<numx[j+1])
				{
					numry[i]=numy[j]+((numy[j+1]-numy[j])/(numx[j+1]-numx[j]))*(numr[i]-numx[j]);
					break;
			}
		}
	}
   char msg[40]="所求的插值为";
   for(int i=0;i<nan;i++)
   {   char temp[30];
	   sprintf(temp,"%f",numry[i]);
	   sprintf(msg,"%s%s ",msg,temp);
   }
   MessageBox(NULL,TEXT(msg),TEXT("结果"),MB_OK);
   exit(0);

};

拉格朗日插值和分段线性插值代码结构与牛顿插值一致。

3、拟合部分

拟合采用了最小二乘法，对一维数据进行拟合

void nihe()//线性拟合
{
   float a,b;
   float up=0,down=0,sumx=0,sumy=0;
   for(int i=0;i<non;i++)
   {
	   sumx+=numx[i];
	   sumy+=numy[i];
	   up+=numx[i]*numy[i];
	   down+=numx[i]*numx[i];
   }
   b=(up-sumx*sumy/non)/(down-sumx*sumx/non);
   a=sumy/non-b*sumx/non;
   for(int i=0;i<nan;i++)
   {
	   numry[i]=numr[i]*b+a;
   }
   char msg[40];
   sprintf(msg,"所求的函数为Y=%fX+%f,所求的点的值为",b,a);
   for(int i=0;i<nan;i++)
   {   char temp[30];
	   sprintf(temp,"%f",numry[i]);
	   sprintf(msg,"%s%s ",msg,temp);
   }
   MessageBox(NULL,TEXT(msg),TEXT("结果"),MB_OK);
   exit(0);

};

公式如下：
在这里插入图片描述

4、界面部分（最为繁杂，作为重点讲解）

#define maxwidth 480
#define maxheight 360//规定界面的大小
int event=0;
int r1[4][4]={{40,10,440,50},{160,100,320,140},{160,160,320,200},{160,220,320,260}};
RECT R11={r1[0][0],r1[0][1],r1[0][2],r1[0][3]};
RECT R12={r1[1][0],r1[1][1],r1[1][2],r1[1][3]};
RECT R13={r1[2][0],r1[2][1],r1[2][2],r1[2][3]};
RECT R14={r1[3][0],r1[3][1],r1[3][2],r1[3][3]};//第一个界面的按钮位置参数

int r2[6][4]={{40,10,440,50},{40,100,200,140},{280,100,440,140},{40,200,200,240},{280,200,440,240}};
RECT R21={r2[0][0],r2[0][1],r2[0][2],r2[0][3]};
RECT R22={r2[1][0],r2[1][1],r2[1][2],r2[1][3]};
RECT R23={r2[2][0],r2[2][1],r2[2][2],r2[2][3]};
RECT R24={r2[3][0],r2[3][1],r2[3][2],r2[3][3]};
RECT R25={r2[4][0],r2[4][1],r2[4][2],r2[4][3]};


MOUSEMSG m;//鼠标指针
void win_inter();
void win_enter();
void win_enter1();
void nuton();//牛顿插值
void langrage();//拉格朗日插值
void fenduan();//分段线性插值
void nihe();//线性拟合

//构建了三个设置字体大小的函数，可以直接用，l-large;m-middle;s-small
void textl()
{
    LOGFONT f;//字体样式指针
    gettextstyle(&f);	//获取字体样式
	f.lfHeight=48;
	_tcscpy(f.lfFaceName,_T("宋体"));	//设置字体为宋体
	settextstyle(&f);                     // 设置字体样式
	settextcolor(BLACK);				//BLACK在graphic.h头文件里面被定义为黑色的颜色常量
	
}
void textm()
{
	LOGFONT f;//字体样式指针
    gettextstyle(&f);	//获取字体样式
	f.lfHeight=30;
	_tcscpy(f.lfFaceName,_T("宋体"));	//设置字体为宋体
	settextstyle(&f);                     // 设置字体样式
	settextcolor(BLACK);				//BLACK在graphic.h头文件里面被定义为黑色的颜色

}
void texts()
{
   LOGFONT f;//字体样式指针
    gettextstyle(&f);	//获取字体样式
	f.lfHeight=20;
	_tcscpy(f.lfFaceName,_T("宋体"));	//设置字体为宋体
	settextstyle(&f);                     // 设置字体样式
	settextcolor(BLACK);				//BLACK在graphic.h头文件里面被定义为黑色的颜色
}
//两个判断函数，用于分别在两个界面中判断鼠标是否移动到按钮内
int judge_b1(int x,int y)
{
	if(x>r1[1][0]&&x<r1[1][2]&&y>r1[1][1]&&y<r1[1][3])return 1;
	if(x>r1[2][0]&&x<r1[2][2]&&y>r1[2][1]&&y<r1[2][3])return 2;
	if(x>r1[3][0]&&x<r1[3][2]&&y>r1[3][1]&&y<r1[3][3])return 3;
	else return 0;
}
int judge_b2(int x,int y)
{
	if(x>r2[1][0]&&x<r2[1][2]&&y>r2[1][1]&&y<r2[1][3])return 1;
	if(x>r2[2][0]&&x<r2[2][2]&&y>r2[2][1]&&y<r2[2][3])return 2;
	if(x>r2[3][0]&&x<r2[3][2]&&y>r2[3][1]&&y<r2[3][3])return 3;
	if(x>r2[4][0]&&x<r2[4][2]&&y>r2[4][1]&&y<r2[4][3])return 4;
	else return 0;
}

void win_first()//第一个界面
{
	for(int i=0;i<256;i++)
	{
	  setbkcolor(RGB(i,i,i));
	  cleardevice();
	  Sleep(5);
	}
	textl();
	//现在规定的按钮区域中写出文字
	drawtext("Math Tool",&R11,DT_CENTER | DT_VCENTER | DT_SINGLELINE);//在矩形区域R1内输入文字，水平居中，垂直居中，单行显示
	textm();
	drawtext("插值",&R12,DT_CENTER | DT_VCENTER | DT_SINGLELINE);//在矩形区域R2内输入文字，水平居中，垂直居中，单行显示
	drawtext("拟合",&R13,DT_CENTER | DT_VCENTER | DT_SINGLELINE);//在矩形区域R3内输入文字，水平居中，垂直居中，单行显示
	drawtext("退出",&R14,DT_CENTER | DT_VCENTER | DT_SINGLELINE);//在矩形区域R3内输入文字，水平居中，垂直居中，单行显示
	setlinecolor(BLACK);
	//然后画出矩形，形成按钮
	rectangle(r1[1][0],r1[1][1],r1[1][2],r1[1][3]);
	rectangle(r1[2][0],r1[2][1],r1[2][2],r1[2][3]);
	rectangle(r1[3][0],r1[3][1],r1[3][2],r1[3][3]);
	while(true)
	{
		m=GetMouseMsg();//读取鼠标操作
		switch(m.uMsg)
		{
	    case WM_MOUSEMOVE:
			    setrop2(R2_XORPEN);
				setlinecolor(LIGHTCYAN);//线条颜色为亮青色
				setlinestyle(PS_SOLID, 3);//设置画线样式为实现，10磅
				setfillcolor(WHITE);//填充颜色为白色
				if(judge_b1(m.x,m.y)!=0)
				{
					if(event != judge_b1(m.x,m.y))
					{
						event =judge_b1(m.x,m.y);//记录这一次触发的按钮
						fillrectangle(r1[event][0],r1[event][1],r1[event][2],r1[event][3]);//有框填充矩形（X1,Y1,X2,Y2）
					}
				}
				else
				{
					if(event != 0)//上次触发的按钮未被修正为原来的颜色
					{
						fillrectangle(r1[event][0],r1[event][1],r1[event][2],r1[event][3]);//两次同或为原来颜色
						event = 0;
					}
				}
				break;
		case WM_LBUTTONDOWN:
				setrop2(R2_NOTXORPEN);//二元光栅——NOT(屏幕颜色 XOR 当前颜色)
				for(int i=0;i<=10;i++)//当鼠标点击时会产生水波特效
				{
					setlinecolor(RGB(25*i,25*i,25*i));//设置圆颜色
					circle(m.x,m.y,2*i);
					Sleep(30);//停顿30ms
					circle(m.x,m.y,2*i);//抹去刚刚画的圆
				}
				switch(judge_b1(m.x,m.y))//switch函数，规定每个按钮的去向
				{
				case 1:  win_inter();break;
				case 2:  win_enter1();
				         nihe();
					     break;
				case 3:  closegraph();exit(0);
				default:
				FlushMouseMsgBuffer();//清空鼠标消息缓存区
	            break;
				}
	    break;
	 }
	}
}

void win_inter()//第二个界面，与第一个界面基本一致
{  
	event=0;
	setbkcolor(RGB(255,255,255));
	cleardevice();//清屏（取决于背景色）
	textl();
	drawtext("插值",&R21,DT_CENTER|DT_VCENTER|DT_SINGLELINE);
	textm();
	drawtext("牛顿",&R22,DT_CENTER | DT_VCENTER | DT_SINGLELINE);//在矩形区域R2内输入文字，水平居中，垂直居中，单行显示
	drawtext("拉格朗日",&R23,DT_CENTER | DT_VCENTER | DT_SINGLELINE);//在矩形区域R3内输入文字，水平居中，垂直居中，单行显示
	drawtext("线性插值",&R24,DT_CENTER | DT_VCENTER | DT_SINGLELINE);//在矩形区域R3内输入文字，水平居中，垂直居中，单行显示
	drawtext("返回",&R25,DT_CENTER | DT_VCENTER | DT_SINGLELINE);//在矩形区域R2内输入文字，水平居中，垂直居中，单行显示

	setlinecolor(BLACK);
	rectangle(r2[1][0],r2[1][1],r2[1][2],r2[1][3]);
	rectangle(r2[2][0],r2[2][1],r2[2][2],r2[2][3]);
	rectangle(r2[3][0],r2[3][1],r2[3][2],r2[3][3]);
	rectangle(r2[4][0],r2[4][1],r2[4][2],r2[4][3]);
	while(true)
	{
		m=GetMouseMsg();//读取鼠标操作
		switch(m.uMsg)
		{
	    case WM_MOUSEMOVE:
			    setrop2(R2_XORPEN);
				setlinecolor(LIGHTCYAN);//线条颜色为亮青色
				setlinestyle(PS_SOLID, 3);//设置画线样式为实现，10磅
				setfillcolor(WHITE);//填充颜色为白色
				if(judge_b2(m.x,m.y)!=0)
				{
					if(event != judge_b2(m.x,m.y))
					{
						event =judge_b2(m.x,m.y);//记录这一次触发的按钮
						fillrectangle(r2[event][0],r2[event][1],r2[event][2],r2[event][3]);//有框填充矩形（X1,Y1,X2,Y2）
					}
				}
				else
				{
					if(event != 0)//上次触发的按钮未被修正为原来的颜色
					{
						fillrectangle(r2[event][0],r2[event][1],r2[event][2],r2[event][3]);//两次同或为原来颜色
						event = 0;
					}
				}
				break;
		case WM_LBUTTONDOWN:
				setrop2(R2_NOTXORPEN);//二元光栅——NOT(屏幕颜色 XOR 当前颜色)
				for(int i=0;i<=10;i++)
				{
					setlinecolor(RGB(25*i,25*i,25*i));//设置圆颜色
					circle(m.x,m.y,2*i);
					Sleep(30);//停顿30ms
					circle(m.x,m.y,2*i);//抹去刚刚画的圆
				}
				switch(judge_b2(m.x,m.y))
				{
				case 1:  win_enter();
					     nuton();
						 break;
				case 2:  win_enter();
					     langrage();
					     break;
				case 3:  win_enter();
					     fenduan();
					     break;
				case 4:  cleardevice();//清屏（取决于背景色）
					     event=0;
					     win_first();
				default:
				FlushMouseMsgBuffer();//清空鼠标消息缓存区
	            break;
				}
	    break;
	 }
	}
}

void win_enter()//输入参数时的界面
{
	InputBox(s,30,"请输入点的个数");//弹出对话框
	strcpy(c,s);
	non=c[0]-'0';
	char str[50];
	for(int i=0;i<non;i++)
	{ 
		sprintf(str,"请输入第%d个点的x坐标和y坐标,用空格隔开,支持浮点数",i+1);
		InputBox(s,30,str);
		sscanf(s,"%f%f",&numx[i],&numy[i]);
	}
	InputBox(s,30,"请输入想要计算的插值点的个数");
	strcpy(c,s);
	nan=c[0]-'0';
	for(int i=0;i<nan;i++)
	{ 
		sprintf(str,"请输入所要求的第%d个插值点的x坐标",i+1);
		InputBox(s,30,str);
		sscanf(s,"%f",&numr[i]);
	}

}

void win_enter1()
{
	InputBox(s,30,"请输入点的个数");
	strcpy(c,s);
	non=c[0]-'0';
	char str[50];
	for(int i=0;i<non;i++)
	{ 
		sprintf(str,"请输入第%d个点的x坐标和y坐标,用空格隔开,支持浮点数",i+1);
		InputBox(s,30,str);
		sscanf(s,"%f%f",&numx[i],&numy[i]);
	}
	InputBox(s,30,"请输入想要计算的点的个数");
	strcpy(c,s);
	nan=c[0]-'0';
	for(int i=0;i<nan;i++)
	{ 
		sprintf(str,"请输入所要求的第%d个点的x坐标",i+1);
		InputBox(s,30,str);
		sscanf(s,"%f",&numr[i]);
	}

}

三、所遇困难和现存不足

1、GUI仍然不够美观
2、只能处理一维数据
3、离群点未判断（难搞，放弃了）
4、最后显示数据的对话框不能弹出到最前，会被数据录入窗口遮挡，需要手动挪开第一个窗口才能看到数据输出。

四、参考资料

插值原理:https://www.cnblogs.com/duye/p/8671820.html
界面设计：https://blog.csdn.net/weixin_44044411/article/details/104276757
感谢所有为我提供帮助的人。

五、源代码

https://download.csdn.net/download/qzl19/12708237

qzl19
关注
为什么要知道.rar加#include stdafx_h_百度
10-11
为什么要加#include stdfx——h
MFC 改变控件字体大小
bury_EE的专栏
06-25 726
LOGFONT: 一个结构体，用来逻辑上表达一个字体   ,在Windows内部，字体是以一个名为LOGFONT的结构来表示 typedef struct tagLOGFONT { LONG lfHeight; LONG lfWidth; LONG lfEscapement; LONG lfOrientation; LONG lfWeight; BYTE lfItalic; BYTE
数值计算方法（一）——插值
朱坚强的博客
09-27 3618
数学原理 此处的拉格朗日插值均为多项式插值，固定下节点，多项式立刻确定下来。 n+1个互异节点满足插值条件的n次拉格朗日插值多项式为： 代码实现 算法实现过程： 1.获取节点个数 2.将节点的x和y存入两个数组中 3.在进行拉格朗日插值中，利用两层循环，外循环累加每一项lk(x)的的值，内循环计算各项lk(x)的值（i!=j时进行累乘） #include<iostream> #i...
C++ jpg内存中绘制图片，字体并以char*的形式输出
Billy Jue的博客
11-12 934
1.Gdi的调用 链接附加依赖项：gdiplus.lib #include <comdef.h>//（不加这个可能会报错） #include <GdiPlus.h> //开始调用 ULONG_PTR g_pGdiToken = NULL; Gdiplus::GdiplusStartupInput g_gdiplus...
Interpolation:MAP4384 的多元插值项目
07-08
标题"Interpolation:MAP4384的多元插值项目"暗示我们将探讨MAP4384课程中的一个具体插值问题，可能是学生作业或研究项目。 首先，我们要理解多元插值的基本原理。在二维情况下，我们有线性插值、双线性插值等方法；...
VideoInterpolation_Capstone：云游戏流的初始视频插值项目
02-16
"VideoInterpolation_Capstone：云游戏流的初始视频插值项目" 这个标题表明我们正在探讨一个与视频处理相关的项目，具体是关于视频插值技术在云游戏流中的应用。"Capstone"通常指的是毕业设计或重大项目，暗示这可能...
MATLAB sinc插值源码
03-07
在信号处理领域，插值是一种常见的技术，用于在已有数据点之间估计未知点的值。MATLAB 是一种广泛使用的编程环境，特别适用于...通过阅读和运行提供的源码，你可以更深入地了解这一概念，并将其应用到自己的项目中。
插值查数_插值labview_
10-01
通过研究“插值查数.vi”，初学者可以了解到如何在LabVIEW中构建插值函数，理解数据处理流程，并将其应用到自己的项目中。实践是学习LabVIEW和其他编程语言的关键，因此，动手操作这个范例程序并进行调试将是提升...
Labview 插值预测_插值预测_
09-30
在这个“Labview 插值预测”项目中，重点是利用插值方法来预测阀门开度与流量之间的关系。插值是一种数学技术，它允许我们通过已知的一系列数据点构建一个函数，这个函数可以精确地通过所有这些点，并用于估计在这些...
解析LOGFONT
独爱糕
11-27 3064
首先我们先来了解LOGFONT这个结构体：LOGFONT是Windows内部字体的逻辑结构，主要用于设置字体格式，其定义如下：typedef struct tagLOGFONTA{LONG lfHeight;LONG lfWidth;LONG lfEscapement;LONG lfOrientation;LONG lfWeight;BYTE lfItalic;BYTE lfUnderline
Windows的字体LOGFONT
weixin_34101784的博客
06-04 859
LOGFONT是Windows内部字体的逻辑结构，主要用于设置字体格式，其定义如下：typedef struct tagLOGFONTA{LONGlfHeight;LONGlfWidth;LONGlfEscapement;LONGlfOrientation;LONGlfWeight;BYTE...
插值方法
weixin_34006468的博客
01-08 207
插值方法有很多，对应到MATLAB里面也有很多函数，和参数选择。 对应的原理分析，嘻嘻，我毕竟不是数学专业的学生，有时间能弄懂就弄懂吧，^_^ 分段线性插值，拉格朗如插值，样条插值(B样条插值好熟悉的名字啊，图形学里面有哦)。 MATLAB实现： 插值函数 语法：y = interp1(x0,y0,x,'method'); 参数method方法： 　　'neares...
算法刷题：300. 最长递增子序列、674. 最长连续递增序列、718. 最长重复子数组、1143. 最长公共子序列
最新发布
qq_61936886的博客
09-12 607
如果text1[i - 1] 与 text2[j - 1]不相同，那就看看text1[0, i - 2]与text2[0, j - 1]的最长公共子序列 和 text1[0, i - 1]与text2[0, j - 2]的最长公共子序列，取最大的。但dp[i][0] 和dp[0][j]要初始值，因为 为了方便递归公式dp[i][j] = dp[i - 1][j - 1] + 1;2.递推公式：当A[i - 1] 和B[j - 1]相等的时候，dp[i][j] = dp[i - 1][j - 1] + 1;
Python实现坐标线性插值计算
"Python坐标线性插值应用实现" 在IT领域，线性插值是一种常见的数值分析技术，常...通过理解这段代码，开发者可以学习到如何处理实际场景中的数据插值问题，这对于数据分析、GIS应用或其他相关领域的项目非常有价值。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值