计算机控制实验指导书

1、前 言计算机控制技术是自动化、电气工程及其自动化和测控技术与仪器专业的一门必修课程。本课程的作用是使学生掌握计算机控制技术的基本理论、基本原理、基本方法和计算机控制系统的实现。它要求学生通过课堂教学和实验后，掌握计算机控制系统的基本概念、工作原理和计算机控制系统设计的基本方法，并具有计算机控制系统的初步分析能力。计算机控制技术实验课程是整个课程教学的重要组成部分，既是对理论知识的实践应用，也是对理论课程的补充、延伸和深化。本实验课程的主要目的是通过与课堂教学的密切配合，巩固和扩充课堂讲授的理论知识，加深对课堂教学内容的理解；训练科学实验的基本技能和工程实践的基本方法，养成严谨的科学态度和工作作

2、风，培养应用所学理论知识独立分析、解决实际问题的能力和实际动手能力。本实验课程共有二个实验项目，为设计性实验，实验有难、易不同的要求，为学生留有独立思考和发挥的空间。本实验指导书对每个实验项目的目的、内容、步骤、运行、预习要求、思考题及实验报告都有明确说明。学生应认真按照实验指导书的要求做好每一个实验。本实验指导书适合自动化专业、电气工程及其自动化专业、测控技术与仪器专业和电子信息工程专业的学生使用。目 录实验一：数字PID控制3实验二：最小拍控制系统5实验三：基于Smith预估器的纯滞后控制系统设计5附录 学生实验报告基本内容要求11实验一：数字PID控制实验学时：3实验类型：设计实验要求：

3、必修一实验目的1了解数字PID的一般原理；2掌握数字PID的编程方法及PID参数的调整。二实验器材PC机三实验内容及要求1必做内容针对具体的被控对象，设计PID控制器，其算法采用位置式或增量式算法。设计、编制程序进行仿真分析，给出仿真结果，并讨论比例、积分、微分的作用； 2选做内容采用一到二种改进型算法设计PID控制器，例如：积分分离PID算法、不完全微分PID算法和微分先行PID算法，设计、编制程序进行仿真分析，给出仿真结果；3被控对象设被控对象为二阶惯性环节，其传递函数如下式：，且采样周期T0.1s，输入为单位阶跃信号； 4性能指标要求调整时间；超调量；静态误差。四实验原理图u*(t)e*

4、(t)y(t)Tr(t)e(t)图1 PID闭环控制系统PID控制器TZOHG0(s)G(z)五实验步骤：1将零阶保持器和给定的控制对象构成广义控制对象，写出其z传递函数，并转化成相应的差分方程；2PID控制算法采用位置式算法或者增量式算法 位置式： 增量式：其中：为比例系数，为积分系数，为微分系数；3通过试凑法选择合适的PID参数，使得系统达到给定的性能指标。并讨论比例、积分和微分对系统产生的作用。 记录选择的PID参数； 记录输出的仿真数据； 记录仿真图形曲线。提示: (1) 调整PID参数时，先将I、D参数保持不变，改变P参数，通过曲线图的变化得出P与超调量及调整时间的关系，再固定P、D

5、参数,分析I对超调量及调整时间的影响，以此方法，观察D对超调量及调整时间的影响。PID调整时间超调量不变不变不变不变不变不变不变不变不变不变不变不变分别记录P、I、D三参数对应的仿真曲线(2) 试凑法选择合适的PID参数，使得系统达到给定的性能指标。记录仿真曲线及参数。六预习要求1复习数字PID控制算法的基本原理；2复习MATLAB和VC、C+的编程、调试方法；3用MATLAB或C+设计、编写数字PID控制算法源程序和仿真图形曲线输出的源程序。七实验报告1写出程序清单并进行注释；2分析调试、运行结果。实验二：最小拍控制系统实验学时：3实验类型：设计实验要求：必修一实验目的1了解数字最小拍控制系

6、统的一般原理；2掌握数字最小拍控制系统的编程方法。二实验器材PC机三实验内容1必做内容针对具体的被控对象，设计系统存在波纹的最小拍数字控制器，编制程序进行仿真分析，给出仿真结果；2选做内容采用阻尼因子法设计最小拍数字控制器，编制程序进行仿真分析，给出仿真结果；3被控对象设被控对象为二阶惯性环节，其传递函数如下式：，且采样周期T1s，输入分别为单位阶跃信号、单位速度信号和单位加速度信号。四实验原理图u*(t)e*(t)y(t)Tr(t)e(t)图1 最小拍闭环控制系统TZOHG0(s)G(z)五实验步骤1将零阶保持器和给定的控制对象构成广义控制对象，写出其传递函数，并化成相应的差分方程；2确定希

7、望的数字控制器： 单位阶跃信号： 单位速度信号： 单位加速度信号： 3通过MATLAB或C+编程进行仿真，验证系统是否在规定拍数达到稳态。分别对单位阶跃信号、单位速度信号和单位加速度信号： 记录调整的最小拍控制参数； 记录输出的仿真数据； 记录仿真图形曲线，标明拍数六预习要求1复习数字最小拍控制系统的基本原理； 2复习MATLAB和C+的编程、调试方法；3用MATLAB或C+设计、编写数字最小拍控制算法源程序和仿真图形曲线输出的源程序。七实验报告1写出程序清单并进行注释；2分析调试、运行结果。实验三：基于Smith预估器的纯滞后控制系统设计实验学时：3实验类型：设计实验要求：必修一实验目的1了

8、解smith预估控制器算法的一般原理；2掌握smith预估控制器算法的编程方法；3. 与PID控制算法进行比较，加深对该控制算法的掌握和理解。二实验器材PC机三实验内容1根据给定对象特性，设计smith预估控制器算法，利用Matlab（Simulink）软件或C+软件编制程序进行仿真分析，给出仿真结果；2设广义被控对象为： 控制系统框图为：取T=1、=2、T1=2.88，经采样（T=1s）保持后，其广义对象z传递函数为， 转换为2个单位迟延控制器参数：Kp=0.5，Ki=0.2，Kd=0。四实验步骤1设计smith预估控制算法，作给定值扰动和外部扰动响应实验，并绘制控制器输出P和系统输出y响应

9、曲线；2被控对象不变，采用理想PID进行给定值扰动和外部扰动响应实验，并绘制控制器输出P和系统输出y响应曲线。五预习要求1复习smith预估控制算法和PID算法的基本原理； 2复习MATLAB和C+的编程、调试方法；3用MATLAB或C+设计、编写smith预估控制算法源程序和仿真图形曲线输出的源程序。六思考题1分析两类控制算法对带迟延对象的控制效果；2根据实验分析Smith预估控制算法的优点是什么，若采用PID算法解决同类问题效果如何？ 七实验报告1写出程序清单并进行注释；2分析调试、运行结果。Matlab辅助设计软件:具体操作步骤：1、 启动Matlab；2、单击工具栏中的Simulink

10、仿真图标 ，进入Simulink仿真环境新建模块库3、新建仿真结构图，寻找模块，拖动到新建仿真结构图中所涉及模块的位置：加法器Sum：在Simulink/Math Operations子库中。离散PID控制器：在SimPowerSystems/Extra Library/DiscreteControl Blocks子库中。离散传递函数Discrete Transfer Fcn：在Simulink/Discrete子库中。示波器Scope：在Simulink/Sinks模型库中。阶跃信号Step：在Simulink/Sources模型库中。4、 修改模块参数。双击模块，在出现的窗口中设置参数。5

11、、 连接模块。将光标移到一个模块的输出端(>)按下鼠标左键拖动鼠标到另一个模块的输入端(>)，松开鼠标左键就可以完成两个模块的连接。6、 设置仿真参数，进行仿真。 在Simulation下拉菜单中设置仿真参数，单击工具栏中的Start Simulation 图标 ），计算机开始仿真，示波器可显示出仿真曲线 。7、通过曲线，分析系统性能。附录参考程序1#include "stdafx.h"#include "1111.h"#include "1111Doc.h"#include "1111View.h"#

12、ifdef _DEBUG#define new DEBUG_NEW#undef THIS_FILEstatic char THIS_FILE = _FILE_;#endifIMPLEMENT_DYNCREATE(CMy1111View, CView)BEGIN_MESSAGE_MAP(CMy1111View, CView)ON_COMMAND(ID_FILE_PRINT, CView:OnFilePrint)ON_COMMAND(ID_FILE_PRINT_DIRECT, CView:OnFilePrint)ON_COMMAND(ID_FILE_PRINT_PREVIEW, CView:OnF

13、ilePrintPreview)END_MESSAGE_MAP()CMy1111View:CMy1111View()CMy1111View:CMy1111View()BOOL CMy1111View:PreCreateWindow(CREATESTRUCT& cs)return CView:PreCreateWindow(cs);/ CMy1111View drawingvoid CMy1111View:OnDraw(CDC* pDC) /求PID参数int i=0; double Kp=0.5,Ki=0.045,Kd=1;double r100;double y100,u100,e1

14、00,datau100;y0=0;e0=0;u0=0;e1=0;y1=0;u1=0;for(i=0;i<100;i+)ri=4.01;for(i=2;i<100;i+)ei=ri-yi-1;dataui=(Kp+Ki+Kd)*ei-(Kp+2*Kd)*ei-1+ei-2; ui=ui-1+dataui; yi=1.724*yi-1-0.741*yi-2+0.0453*ui-1+0.04095*ui-2;CPen hPen1(PS_SOLID,1,RGB(0,0,0);CPen hPen2(PS_SOLID,1,RGB(255,0,0); CPen hPen3(PS_SOLID,1,

15、RGB(0,255,0);CPen * pen0;CRect rect0;char s50; CFont font; font.CreateFont(10,6, 0, 0, 30, FALSE,FALSE, FALSE,DEFAULT_CHARSET, OUT_DEFAULT_PRECIS,CLIP_DEFAULT_PRECIS, PROOF_QUALITY, DEFAULT_PITCH|FF_ROMAN, "Arial");CFont * font0;pen0=pDC->SelectObject(&hPen1); / 画X-Y坐标抽pDC->MoveT

16、o(50,610);pDC->LineTo(950,610);pDC->MoveTo(50,610);pDC->LineTo(50,10);pDC->SelectObject(pen0); pen0=pDC->SelectObject(&hPen1);pDC->MoveTo(950,610);pDC->LineTo(945,605);pDC->MoveTo(950,610);pDC->LineTo(945,615);pDC->MoveTo(50,10);pDC->LineTo(45,15);pDC->MoveTo(

17、50,10);pDC->LineTo(55,15);pDC->SelectObject(pen0);for(i=0;i<15;i+)/画X轴刻度pen0=pDC->SelectObject(&hPen1);pDC->MoveTo(100+i*50,610); pDC->LineTo(100+i*50,605); pDC->SelectObject(pen0);rect0.left=90+i*50;rect0.right=rect0.left+20;rect0.top=615;rect0.bottom=630;sprintf(s,"%d

18、",(i+1)*10);font0=pDC->SelectObject(&font); pDC->SetTextColor(RGB(0,0,0); pDC->SetBkMode(TRANSPARENT);pDC->DrawText(s,-1,&rect0,DT_CENTER);pDC->SelectObject(font0);for(i=0;i<5;i+)/画Y轴刻度pen0=pDC->SelectObject(&hPen1);pDC->MoveTo(50,510-i*100); pDC->LineTo(5

19、5,510-i*100);pDC->SelectObject(pen0);rect0.left=30;rect0.right=rect0.left+20;rect0.top=510-i*100-5;rect0.bottom=rect0.top+20;sprintf(s,"%d",i+1);font0=pDC->SelectObject(&font); pDC->SetTextColor(RGB(0,0,0); pDC->SetBkMode(TRANSPARENT);pDC->DrawText(s,-1,&rect0,DT_CEN

20、TER);pDC->SelectObject(font0); pen0=pDC->SelectObject(&hPen2); /画曲线rtfor(i=0;i<99;i+)pDC->MoveTo(50+i*5),(610-ri*100); pDC->LineTo(50+i*5+5),(610-ri+1*100); pDC->SelectObject(pen0);pen0=pDC->SelectObject(&hPen3); /画曲线ytfor(i=0;i<99;i+)pDC->MoveTo(50+i*5),(610-yi*10

21、0); pDC->LineTo(50+i*5+5),(610-yi+1*100); pDC->SelectObject(pen0); hPen1.DeleteObject();hPen2.DeleteObject(); hPen3.DeleteObject();font.DeleteObject();/ CMy1111View printingBOOL CMy1111View:OnPreparePrinting(CPrintInfo* pInfo)return DoPreparePrinting(pInfo);void CMy1111View:OnBeginPrinting(CDC

22、*, CPrintInfo*)/ TODO: add extra initialization before printingvoid CMy1111View:OnEndPrinting(CDC*, CPrintInfo*)/ TODO: add cleanup after printing/ CMy1111View diagnostics#ifdef _DEBUGvoid CMy1111View:AssertValid() constCView:AssertValid();void CMy1111View:Dump(CDumpContext& dc) constCView:Dump(

23、dc);CMy1111Doc* CMy1111View:GetDocument()/ non-debug version is /inlineASSERT(m_pDocument->IsKindOf(RUNTIME_CLASS(CMy1111Doc);return (CMy1111Doc*)m_pDocument;#endif /_DEBUG参考程序2#include "stdafx.h"#include "1111.h"#include "1111Doc.h"#include "1111View.h"#if

24、def _DEBUG#define new DEBUG_NEW#undef THIS_FILEstatic char THIS_FILE = _FILE_;#endifIMPLEMENT_DYNCREATE(CMy1111View, CView)BEGIN_MESSAGE_MAP(CMy1111View, CView)ON_COMMAND(ID_FILE_PRINT, CView:OnFilePrint)ON_COMMAND(ID_FILE_PRINT_DIRECT, CView:OnFilePrint)ON_COMMAND(ID_FILE_PRINT_PREVIEW, CView:OnFil

25、ePrintPreview)END_MESSAGE_MAP()CMy1111View:CMy1111View()CMy1111View:CMy1111View()BOOL CMy1111View:PreCreateWindow(CREATESTRUCT& cs)return CView:PreCreateWindow(cs);void CMy1111View:OnDraw(CDC* pDC)int i;double r1000;double y1000;double e1000;double u1000;for(i=0;i<1000;i+)ri=1;e0=0;u0=0;y0=0;

26、e1=1;u1=1/3.68;y1=0;for( i=2;i<1000;i+)yi=1.368*yi-1-0.368*yi-2+3.68*ui-1+2.64224*ui-2;ei=ri-yi;ui=-0.718*ui-1+(ei-0.368*ei-1)/3.68; CPen hPen1(PS_SOLID,1,RGB(0,0,0);CPen hPen2(PS_SOLID,1,RGB(255,0,0); CPen hPen3(PS_SOLID,1,RGB(0,255,0);CPen * pen0; CRect rect0;char s50; CFont font; font.CreateFo

27、nt(10,6, 0, 0, 30, FALSE,FALSE, FALSE,DEFAULT_CHARSET, OUT_DEFAULT_PRECIS, CLIP_DEFAULT_PRECIS, PROOF_QUALITY, DEFAULT_PITCH|FF_ROMAN, "Arial");CFont * font0;/ 画X-Y坐标抽pen0=pDC->SelectObject(&hPen1);pDC->MoveTo(50,610);pDC->LineTo(950,610);pDC->MoveTo(50,610);pDC->LineTo(

28、50,10);pDC->SelectObject(pen0); pen0=pDC->SelectObject(&hPen1);pDC->MoveTo(950,610);pDC->LineTo(945,605);pDC->MoveTo(950,610);pDC->LineTo(945,615);pDC->MoveTo(50,10);pDC->LineTo(45,15);pDC->MoveTo(50,10);pDC->LineTo(55,15);pDC->SelectObject(pen0);/画X轴刻度for(i=0;i&

29、lt;9;i+)pen0=pDC->SelectObject(&hPen1);/画的刻度线pDC->MoveTo(50+i*100,610); pDC->LineTo(50+i*100,605); pDC->SelectObject(pen0);rect0.left=40+i*100;rect0.right=rect0.left+20;rect0.top=615;rect0.bottom=630;sprintf(s,"%d",i);font0=pDC->SelectObject(&font); pDC->SetTextCo

30、lor(RGB(0,0,0); pDC->SetBkMode(TRANSPARENT);pDC->DrawText(s,-1,&rect0,DT_CENTER);pDC->SelectObject(font0);/画Y轴刻度for(i=0;i<2;i+)pen0=pDC->SelectObject(&hPen1);/画刻度线pDC->MoveTo(50,410-i*200); pDC->LineTo(55,410-i*200);pDC->SelectObject(pen0);/rect0.left=30;/画刻度rect0.rig

31、ht=rect0.left+20;rect0.top=410-i*200-5;rect0.bottom=rect0.top+20;sprintf(s,"%d",i+1);font0=pDC->SelectObject(&font); pDC->SetTextColor(RGB(0,0,0); pDC->SetBkMode(TRANSPARENT);pDC->DrawText(s,-1,&rect0,DT_CENTER);pDC->SelectObject(font0); /画曲线rt pen0=pDC->SelectObj

32、ect(&hPen2);for(i=0;i<8;i+)pDC->MoveTo(50+i*100),(610-ri*200); pDC->LineTo(50+i*100+100),(610-ri+1*200); pDC->SelectObject(pen0);/画曲线ytpen0=pDC->SelectObject(&hPen3);for(i=0;i<8;i+)pDC->MoveTo(50+i*100),(610-yi*200); pDC->LineTo(50+i*100+100),(610-yi+1*200); pDC->S

33、electObject(pen0); hPen1.DeleteObject();hPen2.DeleteObject(); hPen3.DeleteObject();font.DeleteObject();BOOL CMy1111View:OnPreparePrinting(CPrintInfo* pInfo)return DoPreparePrinting(pInfo);void CMy1111View:OnBeginPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/)void CMy1111View:OnEndPrinting(CDC* /*pDC*/

34、, CPrintInfo* /*pInfo*/)#ifdef _DEBUGvoid CMy1111View:AssertValid() constCView:AssertValid();void CMy1111View:Dump(CDumpContext& dc) constCView:Dump(dc);CMy1111Doc* CMy1111View:GetDocument() / non-debug version is inlineASSERT(m_pDocument->IsKindOf(RUNTIME_CLASS(CMy1111Doc);return (CMy1111Doc

35、*)m_pDocument;#endif /_DEBUG参考程序3#include "stdafx.h"#include "1111.h"#include "1111Doc.h"#include "1111View.h"#ifdef _DEBUG#define new DEBUG_NEW#undef THIS_FILEstatic char THIS_FILE = _FILE_;#endifIMPLEMENT_DYNCREATE(CMy1111View, CView)BEGIN_MESSAGE_MAP(CMy111

36、1View, CView)ON_COMMAND(ID_FILE_PRINT, CView:OnFilePrint)ON_COMMAND(ID_FILE_PRINT_DIRECT, CView:OnFilePrint)ON_COMMAND(ID_FILE_PRINT_PREVIEW, CView:OnFilePrintPreview)END_MESSAGE_MAP()CMy1111View:CMy1111View()CMy1111View:CMy1111View()BOOL CMy1111View:PreCreateWindow(CREATESTRUCT& cs)return CView

37、:PreCreateWindow(cs);void CMy1111View:OnDraw(CDC* pDC)int i;double r1000;double y1000;double e1000;double u1000;for(i=0;i<1000;i+)ri=i;e0=0;u0=0;y0=0;e1=1;u1=0.543;y1=0;for( i=2;i<1000;i+)yi=1.368*yi-1-0.368*yi-2+3.68*ui-1+2.64224*ui-2;ei=ri-yi;ui=0.282*ui-1+0.718*ui-2+0.543*ei-0.471324*ei-1+0

38、.099912*ui-2; CPen hPen1(PS_SOLID,1,RGB(0,0,0);CPen hPen2(PS_SOLID,1,RGB(255,0,0); CPen hPen3(PS_SOLID,1,RGB(0,255,0);CPen * pen0; CRect rect0;char s50; CFont font; font.CreateFont(10,6, 0, 0, 30, FALSE,FALSE, FALSE,DEFAULT_CHARSET, OUT_DEFAULT_PRECIS,CLIP_DEFAULT_PRECIS, PROOF_QUALITY, DEFAULT_PITC

39、H|FF_ROMAN, "Arial");CFont * font0;/ 画X-Y坐标抽pen0=pDC->SelectObject(&hPen1);pDC->MoveTo(50,610);pDC->LineTo(950,610);pDC->MoveTo(50,610);pDC->LineTo(50,10);pDC->SelectObject(pen0); pen0=pDC->SelectObject(&hPen1);pDC->MoveTo(950,610);pDC->LineTo(945,605);pD

40、C->MoveTo(950,610);pDC->LineTo(945,615);pDC->MoveTo(50,10);pDC->LineTo(45,15);pDC->MoveTo(50,10);pDC->LineTo(55,15);pDC->SelectObject(pen0);/画X轴刻度for(i=0;i<9;i+)pen0=pDC->SelectObject(&hPen1);/画的刻度线pDC->MoveTo(50+i*100,610); pDC->LineTo(50+i*100,605); pDC->Sel

41、ectObject(pen0);rect0.left=40+i*100; /画刻度rect0.right=rect0.left+20;rect0.top=615;rect0.bottom=630;sprintf(s,"%d",i);font0=pDC->SelectObject(&font); pDC->SetTextColor(RGB(0,0,0); pDC->SetBkMode(TRANSPARENT);pDC->DrawText(s,-1,&rect0,DT_CENTER);pDC->SelectObject(font0)

42、;for(i=0;i<10;i+)/画Y轴刻度pen0=pDC->SelectObject(&hPen1);/画刻度线pDC->MoveTo(50,560-i*50); pDC->LineTo(55,560-i*50);pDC->SelectObject(pen0);rect0.left=30;/画刻度rect0.right=rect0.left+20;rect0.top=560-i*50-5;rect0.bottom=rect0.top+20;sprintf(s,"%d",i+1);font0=pDC->SelectObject

43、(&font); pDC->SetTextColor(RGB(0,0,0); pDC->SetBkMode(TRANSPARENT);pDC->DrawText(s,-1,&rect0,DT_CENTER);pDC->SelectObject(font0); / pen0=pDC->SelectObject(&hPen2); /画曲线rtfor(i=0;i<8;i+)pDC->MoveTo(50+i*100),(610-ri*50); pDC->LineTo(50+i*100+100),(610-ri+1*50); pDC

44、->SelectObject(pen0);pen0=pDC->SelectObject(&hPen3); /画曲线ytfor(i=0;i<8;i+)pDC->MoveTo(50+i*100),(610-yi*50); pDC->LineTo(50+i*100+100),(610-yi+1*50); pDC->SelectObject(pen0); hPen1.DeleteObject();hPen2.DeleteObject(); hPen3.DeleteObject();font.DeleteObject();BOOL CMy1111View:On

45、PreparePrinting(CPrintInfo* pInfo)return DoPreparePrinting(pInfo);void CMy1111View:OnBeginPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/)void CMy1111View:OnEndPrinting(CDC* /*pDC*/, CPrintInfo* /*pInfo*/)#ifdef _DEBUGvoid CMy1111View:AssertValid() constCView:AssertValid();void CMy1111View:Dump(CDumpCo

46、ntext& dc) constCView:Dump(dc);CMy1111Doc* CMy1111View:GetDocument() / non-debug version is inlineASSERT(m_pDocument->IsKindOf(RUNTIME_CLASS(CMy1111Doc);return (CMy1111Doc*)m_pDocument;#endif /_DEBUG参考程序4#include "stdafx.h"#include "1111.h"#include "1111Doc.h"#in

47、clude "1111View.h"#ifdef _DEBUG#define new DEBUG_NEW#undef THIS_FILEstatic char THIS_FILE = _FILE_;#endifIMPLEMENT_DYNCREATE(CMy1111View, CView)BEGIN_MESSAGE_MAP(CMy1111View, CView)ON_COMMAND(ID_FILE_PRINT, CView:OnFilePrint)ON_COMMAND(ID_FILE_PRINT_DIRECT, CView:OnFilePrint)ON_COMMAND(ID_

48、FILE_PRINT_PREVIEW, CView:OnFilePrintPreview)END_MESSAGE_MAP()CMy1111View:CMy1111View()CMy1111View:CMy1111View()BOOL CMy1111View:PreCreateWindow(CREATESTRUCT& cs)return CView:PreCreateWindow(cs);void CMy1111View:OnDraw(CDC* pDC)int i;float r1000;double y1000;double e1000;double u1000;for(i=0;i&l

49、t;1000;i+)ri=i*i/2;y0=0;e0=0;u0=0;y1=0;e1=0.5;u1=1.5/3.68;y2=1.5;e2=0.5;u2=1.282*1.5/3.68+(3*0.5-4.104*0.5)/3.68;for( i=3;i<1000;i+)yi=1.368*yi-1-0.368*yi-2+3.68*ui-1+2.64224*ui-2;ei=ri-yi;ui=1.282*ui-1+0.436*ui-2-0.718*ui-3+(3*ei-4.104*ei-1+2.104*ui-2-0.368*ui-3)/3.68; CPen hPen1(PS_SOLID,1,RGB(

50、0,0,0);CPen hPen2(PS_SOLID,1,RGB(255,0,0); CPen hPen3(PS_SOLID,1,RGB(0,255,0);CPen * pen0; CRect rect0;char s50; CFont font; font.CreateFont(10,6, 0, 0, 30, FALSE,FALSE, FALSE,DEFAULT_CHARSET, OUT_DEFAULT_PRECIS,CLIP_DEFAULT_PRECIS, PROOF_QUALITY, DEFAULT_PITCH|FF_ROMAN, "Arial");CFont * f

51、ont0;pen0=pDC->SelectObject(&hPen1); / 画X-Y坐标抽pDC->MoveTo(50,610);pDC->LineTo(950,610);pDC->MoveTo(50,610);pDC->LineTo(50,10);pDC->SelectObject(pen0); pen0=pDC->SelectObject(&hPen1);pDC->MoveTo(950,610);pDC->LineTo(945,605);pDC->MoveTo(950,610);pDC->LineTo(945,615);pDC->MoveTo(50,10);pDC