1、VB课程设计报告VB课程设计报告Visual Basic程序设计课程设计报告设计题目 直线拟合 题目难度 四星级 姓 名 蒋龙鑫 学 号 *F0309 指导教师 申 屠 德 忠 提交日期 4月15号 一、题目描述与功能要求 1题目描述 本题目重要是用于将实验中得到的具有线性特征的一组数据通过最小二乘法拟合成直线,求出其直线方程和绘制出该直线的图形。也就是寻找一条直线使这些离散的点尽可能平均的位于该直线的两边,它对于寻找数据的规律有很大的帮助,是一个比较实用的程序。2功能要求:(1)运行程序中,首先显示如下图窗体1的“直线拟合”窗口,从中选择要拟合数据的组数。单击“确定”按钮进入“数据输入”窗口
2、。(2)“数据输入”窗口中有10对文本框,用于输入最多10组实验数据。根据在“直线拟合”窗口的选择不同,这10对文本框只有指定数目的文本框可用。输入指定的数据,然后单击“确定”按钮。如果数据未完全输入,则显示错误提示的消息框。(3)如果数据输入正确,在“数据输入”窗口中单击“确定”按钮之后,程序进入“直线参数与图形”窗口,在此窗口中,程序根据输入的数据自动计算拟合直线方程的参数,并将系数和方程显示在窗口顶部。在窗口的绘图区绘制直线图形。二、程序设计(一)计算模型或算法实行假设实验中得到n个点,P1(X1,Y1),P2(X2,Y2),Pn(Xn,Yn),用i表示第i个点.由于各种误差的存在,这些
3、数据点不可能完全在一条直线上。如果我们能找到一条直线,使这n个点到该直线的距离之平方和最小,则该直线就是用最小二乘法拟合出来的直线。设拟合的直线方程为F(X)=Y=bX+a,由于拟合出来的直线斜率是固定的,我们可以用每一个数据点与直线在Y方向上的“残差”来定义,即残差Yi=F(Xi)-Yi,可得: (Yi)2=(F(Xi)-Yi)2=(bXi+a-Yi)2根据最小二乘法原理,(bXi+a-Yi)2取极小值时,即为残差最小,这时确定的直线即为拟合的直线。上式是一个二元函数求极值的问题,可分别对自变量a与b求偏导数,可易得: na+bXi =Yi aXi+b(Xi)2=(XiYi) 从而得到公式:
4、由上公式可知,我们可用n个点的坐标确定直线方程的系数a与b,得到方程并将方程拟合出来.本题目的程序设计由三个窗体组成。在Form1中,有1个框架控件、1个标签控件、2个命令按钮和1个由6个单选按钮组成的控件数组。单选按钮控件数组用来选择要输入的数据组数,当选择好后,单击“确定”按钮,就可以进入下一个窗体。单击“取消”按钮则退出程序Form1的代码见附件。 在Form2中有2个命令按钮、20个文本框(每10个组成一个控件数组)、12个标签控件和许多直线组成。在文本框中输入实验数据后,点击“确定”按钮,若数据没有全部输入,则出现消息框,提醒没有完全输入。若已经完全输入,则显示Form3,隐藏For
5、m2。与此同时,对实验数据计算,得出拟合直线的截距和斜率,并在Form3的标签控件1中体现;根据所得的斜率和截距得到拟合直线方程,并在Form3的标签2中体现;还可以在Form3中绘出拟合直线和每个实验点(用绿色),根据实验数据调整坐标系统,画出拟合直线。Form2的代码见附件。(二)模块设计(窗体界面与代码)本题目的程序设计由3个模块组成,其中窗体模块3个。序号模块名称模块类型模块功能1 Form1窗体模块选择要输入的数据组数,当选择好后,单击“确定”按钮,就可以进入下一个窗体;单击“取消”按钮则退出程序2 Form2窗体模块输入实验数据后,对实验数据计算,得出拟合直线的截距和斜率,根据所得
6、的斜率和截距得到拟合直线方程。 3Form3窗体模块根据实验数据调整坐标系统,画出拟合直线。窗体1的设计 表1 窗体1中控件列表序号控件名称控件类型LeftTopWidthHeightCaption1Label1标签控件14401202055375直线拟合程序2Frame1框架控件48060038551455请选择数据组数3Option1单选框控件3603609752555组4Option2单选框控件3607209752556组5Option3单选框控件36010809752557组6Option4单选框控件25203609752558组7Option5单选框控件25207209752559组
7、8Option6单选框控件2520108097525510组9Command1命令按钮控件72022801575495确定10Command2命令按钮控件252022801575495退出表2 窗体1中事件过程列表序号事件过程功能1Option1-click()选择数据组数为5组2Option2-click()选择数据组数为6组3Option3-click()选择数据组数为7组4Option4-click()选择数据组数为8组5Option5-click()选择数据组数为9组6Option6-click()选择数据组数为10组7Command1-click()进入第二个窗体坐标轴的输入8Com
8、mand2-click()退出拟合值线的数据输入表3 窗体1中通用过程与函数列表序号过程与函数输入参数输出参数功能1MsgBox()以文本向用户显示运行过程中的信息,让用户单击按钮作选择2Chr()10换行符使消息内容在中间某处换行3Chr()13回车符窗体1的代码Form1Private Sub Command1_Click() 显示第二个窗体Form1.HideForm2.ShowForm3.Hidezxnh = MsgBox(你必须完全输入所选的项目! & Chr(10) & Chr(13) & _ 否则程序将会出错!, 48, 直线拟合)End SubPrivate Sub Comma
9、nd2_Click() 结束程序EndEnd Sub使第二个窗体的文本框按照第一个窗体给出的数据组数显示Private Sub Option1_Click()Form2.Label6.Enabled = FalseForm2.Label7.Enabled = FalseForm2.Label8.Enabled = FalseForm2.Label9.Enabled = FalseForm2.Label10.Enabled = FalseFor n = 6 To 10 Form2.Text1(n).Enabled = False Form2.Text1(n).BackColor = &H8000
10、000F Form2.Text2(n).Enabled = False Form2.Text2(n).BackColor = &H8000000FNextEnd SubPrivate Sub Option2_Click()Form2.Label6.Enabled = TrueForm2.Label7.Enabled = FalseForm2.Label8.Enabled = FalseForm2.Label9.Enabled = FalseForm2.Label10.Enabled = FalseForm2.Text1(6).Enabled = TrueForm2.Text1(6).BackC
11、olor = &HFFFFFFForm2.Text2(6).Enabled = TrueForm2.Text2(6).BackColor = &HFFFFFFFor n = 7 To 10 Form2.Text1(n).Enabled = False Form2.Text1(n).BackColor = &H8000000F Form2.Text2(n).Enabled = False Form2.Text2(n).BackColor = &H8000000FNextEnd SubPrivate Sub Option3_Click()Form2.Label6.Enabled = TrueFor
12、m2.Label7.Enabled = TrueForm2.Label8.Enabled = FalseForm2.Label9.Enabled = FalseForm2.Label10.Enabled = FalseFor n = 6 To 7 Form2.Text1(n).Enabled = True Form2.Text1(n).BackColor = &HFFFFFF Form2.Text2(n).Enabled = True Form2.Text2(n).BackColor = &HFFFFFFNext For n = 8 To 10 Form2.Text1(n).Enabled =
13、 False Form2.Text1(n).BackColor = &H8000000F Form2.Text2(n).Enabled = False Form2.Text2(n).BackColor = &H8000000FNextEnd SubPrivate Sub Option4_Click()Form2.Label6.Enabled = TrueForm2.Label7.Enabled = TrueForm2.Label8.Enabled = TrueForm2.Label9.Enabled = FalseForm2.Label10.Enabled = FalseFor n = 6 T
14、o 8 Form2.Text1(n).Enabled = True Form2.Text1(n).BackColor = &HFFFFFF Form2.Text2(n).Enabled = True Form2.Text2(n).BackColor = &HFFFFFFNextFor n = 9 To 10 Form2.Text1(n).Enabled = False Form2.Text1(n).BackColor = &H8000000F Form2.Text2(n).Enabled = False Form2.Text2(n).BackColor = &H8000000FNextEnd
15、SubPrivate Sub Option5_Click()Form2.Label6.Enabled = TrueForm2.Label7.Enabled = TrueForm2.Label8.Enabled = TrueForm2.Label9.Enabled = TrueForm2.Label10.Enabled = FalseForm2.Text1(10).Enabled = FalseForm2.Text1(10).BackColor = &H8000000FForm2.Text2(10).Enabled = FalseForm2.Text2(10).BackColor = &H800
16、0000FFor n = 6 To 9 Form2.Text1(n).Enabled = True Form2.Text1(n).BackColor = &HFFFFFF Form2.Text2(n).Enabled = True Form2.Text2(n).BackColor = &HFFFFFFNextEnd SubPrivate Sub Option6_Click()Form2.Label6.Enabled = TrueForm2.Label7.Enabled = TrueForm2.Label8.Enabled = TrueForm2.Label9.Enabled = TrueFor
17、m2.Label10.Enabled = TrueFor n = 6 To 10 Form2.Text1(n).Enabled = True Form2.Text1(n).BackColor = &HFFFFFF Form2.Text2(n).Enabled = True Form2.Text2(n).BackColor = &HFFFFFFNextEnd Sub窗体2的设计表4 窗体2中控件列标序号控件名称控件类型LeftTopWidthHeightCaption1Label1标签控件480880495255第1点2Label2标签控件4801240495255第2点3Label3标签控件4
18、801600495255第3点4Label4标签控件4801960495255第4点5Label5标签控件4802320495255第5点6Label6标签控件4802680495255第6点7Label7标签控件4803040495255第7点8Label8标签控件4803400495255第8点9Label9标签控件4803760495255第9点10Label10标签控件4804120575255第10点11Label11标签控件1320360975255X坐标值:12Label12标签控件2520360975255Y坐标值:13Text1(1)文本框控件144084073527014T
19、ext1(2)文本框控件1440120073527015Text1(3)文本框控件1440156073527016Text1(4)文本框控件1440192073527017Text1(5)文本框控件1440228073527018Text1(6)文本框控件1440264073527019Text1(7)文本框控件1440300073527020Text1(8)文本框控件1440336073527021Text1(9)文本框控件1440372073527022Text1(10)文本框控件1440408073527023Text2(1)文本框控件264084073527024Text2(2)文本框
20、控件2640120073527025Text2(3)文本框控件2640156073527026Text2(4)文本框控件2640192073527027Text2(5)文本框控件2640228073527028Text2(6)文本框控件2640264073527029Text2(7)文本框控件2640300073527030Text2(8)文本框控件2640336073527031Text2(9)文本框控件2640372073527032Text2(10)文本框控件2640408073527033Command1命令按钮控件48046801335495确定34Command2命令按钮控件204
21、046801335495返回X1Y1X2Y2BorderStyleBorderWidth35Line1直线控件24024036002401-solid136Line2直线控件24024024044401-solid137Line3直线控件3600240360044401-solid138Line4直线控件1200240120044401-solid139Line5直线控件2400240240044401-solid140Line6直线控件24072036007201-solid141Line7直线控件2404440360044401-solid142Line8直线控件2402401200720
22、1-solid1表5 窗体2中事件过程列表序号事件过程功能1Command1-click()进入绘图窗体,拟合成直线2Command2-click()退出拟合直线的数据输入 表6 窗体2中通用过程与函数列表序号过程与函数输入参数输出参数功能1Spc()33个空格组成的字符串返回number个空格组成的字符串 窗体模块代码2的 Form2:Option Base 1Dim m As IntegerDim X(10) As SingleDim Y(10) As SinglePrivate Sub Command1_Click() 进入第三个绘图窗体Form1.HideForm2.HideForm3
23、.ShowDim A As Single 截距Dim B As Single 斜率Dim c As SingleDim d As SingleDim e As SingleDim f As SingleDim W As SingleDim R As SingleDim T As Single讨论点的个数If Form2.Text1(5).Enabled = True Then m = 5If Form2.Text1(6).Enabled = True Then m = 6If Form2.Text1(7).Enabled = True Then m = 7If Form2.Text1(8).E
24、nabled = True Then m = 8If Form2.Text1(9).Enabled = True Then m = 9If Form2.Text1(10).Enabled = True Then m = 10对点的x.y坐标进行赋值If m = 5 Then For n = 1 To 5 X(n) = Text1(n).Text Y(n) = Text2(n).Text NextEnd IfIf m = 6 Then For n = 1 To 6 X(n) = Text1(n).Text Y(n) = Text2(n).Text NextEnd IfIf m = 7 Then
25、For n = 1 To 7 X(n) = Text1(n).Text Y(n) = Text2(n).Text NextEnd IfIf m = 8 Then For n = 1 To 8 X(n) = Text1(n).Text Y(n) = Text2(n).Text NextEnd IfIf m = 9 Then For n = 1 To 9 X(n) = Text1(n).Text Y(n) = Text2(n).Text NextEnd IfIf m = 10 Then For n = 1 To 10 X(n) = Text1(n).Text Y(n) = Text2(n).Tex
26、t NextEnd If For n = 1 To m c = X(n) + c d = X(n) * X(n) + d e = X(n) * Y(n) + e f = Y(n) + fNextA = f/m-(m*e-c*f)*c/m/(m*d-c*c) 计算直线的截距B = (c * f - m * e) / (c * c - m * d) 计算直线的斜率打印数据方程的位置Form3.CurrentX = 200Form3.CurrentY = 350Form3.Print 直线斜率(B)= & BForm3.CurrentX = 200Form3.CurrentY = 580Form3.
27、Print 直线截距(A)= & AForm3.CurrentX = 200Form3.CurrentY = 810Form3.Print 直线方程为:; Spc(3); Y= & B & X+ & A在坐标图上打印出点来For n = 1 To m Form3.PSet (360 + 37 * X(n), 5400 - 37 * Y(n), RGB(0, 0, 0)NextW = 5400 - 3700 * B - 37 * AR = 360 + (3700 - 37 * A) / BT = 360 - 37 * A / B将超出坐标轴的线去掉If A = 0 Then Form3.Line (360, 5400 - 37 * A)-(4060, W) If W 5400 And B 0 Then Form3.Line (T, 5400)-(40
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1