#include <dos.h>
#include <time.h>
#include <math.h>
#include <stdio.h>
#include <conio.h>
#include <string.h>
#include <stdlib.h>
#include <graphics.h>
#include "main.h"
Rect win;
char pr_fm[80],str[80],KBinput[80];
int dir[2]={0,1};
int MaxX,MaxY;
char *par_name[parS]={"F","omiga","k1","k3","damp","ICx1","ICx2","upper value of CP","step of CP"};
char *set_name[SettingS]={"sin0cos1","transient cycles","steady cycles","wave cycles","Poincare cycles","Bifur Ampl","segments per cycle","control parameter"};
/* GRAPH.C */
/* NO.5 function */
clear_rec(int l,int t,int r,int b)
{
struct viewporttype v;/*定义的结构变量viewport*/
getviewsettings(&v);/*获得关于现行窗口的信息,并将其存于viewporttype定义的结构变量viewport*/
setviewport(l,t,r,b,0);/*设置图形视口函数
设定一个以(l,t)象元点为左上角, (r,b)象元为右下角的图形窗口, 其中l,t,r,b是相对于整个屏幕的坐标。图形窗口以外可以接触。不做裁剪 窗口调用 */
clearviewport();/* 清除现行图形窗口的内容。*/
setviewport(v.left,v.top,v.right,v.bottom,v.clip);/*viewporttype的结构*/
}
/* NO.6 function 在屏幕上指定位置以指定格式输出一个浮点数*/
print_value(int x,int y,DATA v,char *format)
{
int w,h;
char str[80];/*定义一个字符数组*/
struct textsettingstype ty;
gettextsettings(&ty);/*将当前字体、显示方向、字符大小和水平垂直方向文本排齐方式等设置值填入textsettingstype型结构里，从而可从中获取文本信息。*/
sprintf(str,format,v); /*将数字转化为字符串 format是调用Pr_Fm(*s) sprintf(pr_fm,"%s%s",s,PR_FM); */
w=textwidth(str);/*计算出字符串宽度的像素数,1个字符占8个像素数*/
h=textheight(str);/*计算出字符串高度的像素数,占8个像素数*/
y+=h;
if (ty.direction==HORIZ_DIR) /*判断字符串输出是否从左到右*/
clear_rec(x,y,x+w,y+h);
else clear_rec(x-h,y,x,y+w);
outtextxy(x,y,str);/*在(x,y)输出str字符串*/
}
/* NO.7 function 在屏幕上指定位置输出一个字符串 */
print_str(int x,int y,char *str)
{ int w,h;
struct textsettingstype ty;
gettextsettings(&ty);/*将当前字体、显示方向、字符大小和水平垂直方向文本排齐方式等设置值填入textsettingstype型结构里，从而可从中获取文本信息。*/
w=textwidth(str);
h=textheight(str);
y+=h;
if (ty.direction==HORIZ_DIR)
clear_rec(x,y,x+w,y+h);
else clear_rec(x-h,y,x,y+w);
outtextxy(x,y,str);/*在(x,y)输出str字符串*/
}
/* NO.8 function */
axis_name(Rect *rec,char *name,int axis,int dir,int w1_c0)/*在指定窗口边缘显示或擦除坐标轴名称,w1_c0为1显示,0擦除*/
{ int x,y,w;
if (!w1_c0) /*w1_c0=0或1 判断真假*/
setcolor(0); else setcolor(15);
if (!name) sprintf(pr_fm,"x%d",axis+1); /*name字符串,值为t(s),str,NULL,par_name[par]*/
else strcpy(pr_fm,name);/*字符串复制*/
w=textwidth(pr_fm);
settextstyle(0,dir,1);/*设置图形文本当前字体8*8点阵字型、文本显示方向（水平显示或垂直显示）以及字符大小。dir取0,1 0从做到右 1从底向上*/
if (dir) { x=rec->l-10; y=rec->cy-w/2; clear_rec(x-10,y,x,y+w);}
else if (rec->r-rec->l<100)
{ y=rec->b+4;
x=rec->cx-w/2;
clear_rec(x-2,y-2,x+w,y+8);}
else {y=rec->b+10;
x=rec->cx-w/2;
clear_rec(x-2,y-2,x+w,y+8);}
outtextxy(x,y,pr_fm);
if (dir) settextstyle(0,0,1);
if (!w1_c0) setcolor(15);
}
/* NO.9 function */
Init_rect(Rect *win) /*计算矩形窗口win的中心点图形坐标*/
{
win->cx=(win->l+win->r)/2;
win->cy=(win->t+win->b)/2;
}
/* NO.10 function */
char *Pr_Fm(char *s)
{
sprintf(pr_fm,"%s%s",s,PR_FM); /*输出字符串*/
return pr_fm;
}
/* ******NO.11 function */
int lever(DATA f)/*判断浮点数的浮点级*/
{
int i=0;
DATA a=1.0;
if (ABS(f)<1.0)
{
while(1.-ABS(f/a)>1.e-1)
{
i--;
a/=10.0;
}
}
else
{
while(ABS(f/a)-1.0>1.e-14)
{
i++;
a*=10.0;
}
}
return (i);
}
/* NO.12 function */
massage(char *str,int dx,int flag,int w1_c0)/*字符串信息str显示或擦除.flag=0,1分别在屏幕底部倒数第一,第二行显示.w1_c0赋1为显示 0为擦除 */
{
if (w1_c0)/*判断真值*/
{
if (flag<0) print_str(2+dx,win.b+8,str);
else if (!flag)
print_str(2+dx,win.b+18,str);
else print_str(2+dx,win.b+28,str);/*在(2+dx,win.b+8)输出str*/
}
else if (str)/*字符串str显示*/
{
setcolor(0);
if (flag<0)
print_str(2+dx,win.b+8,str);
else if (!flag)
print_str(2+dx,win.b+18,str);
else print_str(2+dx,win.b+28,str);/*在(2+dx,win.b+8)输出str*/
setcolor(7);
}
else
{
if (flag<0)
clear_rec(2+dx,win.b+16,MaxX-1,win.b+24);
else if (!flag)
clear_rec(2+dx,win.b+26,MaxX-1,win.b+34);
else clear_rec(2+dx,win.b+36,MaxX-1,win.b+44);
}
}
/* NO.13 function */
massage_value(DATA d,char *fm,int dx,int flag,int w1_c0)/*浮点信息d显示或擦除.flag=0,1分别在屏幕底部倒数第一,第二显示.we_c0赋1为显示,0为擦除*/
{
if (w1_c0)
{
if (!flag)
print_value(2+dx,win.b+18,d,Pr_Fm(fm));
else
print_value(2+dx,win.b+28,d,Pr_Fm(fm));
}
else if (fm)
{
setcolor(0);
if (!flag)
print_value(2+dx,win.b+18,d,Pr_Fm(fm));
else
print_value(2+dx,win.b+28,d,Pr_Fm(fm));
setcolor(7);
}
else
{
if (!flag)
clear_rec(2+dx,win.b+17,MaxX-1,win.b+27);
else clear_rec(2+dx,win.b+27,MaxX-1,win.b+37);
}
}
/* ******NO.14 function */
getstr(int K,int lim)/*从键盘获取字符并顺次存入KBinput中.K为字符ASCAII码,lim为最多输入的字符数*/
{ static int i=0; int j;
if (K>=14&&K<=126)
{ KBinput[i]=K;
if (i>lim)
{ massage("The buffer is full!",win.l+50,1,1);i--;}
i++; KBinput[i]=0;
massage( KBinput,win.l+50,1,1);
}
switch (K)
{ case 8:massage( KBinput,win.l+50,1,0);
if (--i<0)
{i=0;break;}
KBinput[i]=0; massage( KBinput,win.l+50,1,1); break;
case 13:case F1:case F2:
if(i&&i<lim)
massage(NULL,0,1,0);
else if (!i)
KBinput[0]=0;
i=0;
break;
}
}
/* ******NO.15 function */
Getstring(char *str)/*调用函数gestr()从键盘获取字符串于KBinput中 str为提示信息*/
{
int K;
putch(7);
massage(NULL,0,0,0);
massage(str,0,0,1);
do
{
if (kbhit())
{
K=getch();
if (!K)
K=getch()+256;
getstr(K,20);
}
if (K==ENT||K==F1||K==F2)
break;/*F1 F2 ENT停止*/
}
while(1);
}
/* NO.16 function */
Init_fig(DATA xL,DATA xU,DATA yL,DATA yU,Rect *wv,int flag)
/*初始化图形显示窗口wv xL，xU，yL，yU，分别为水平，铅锤轴所代表*/
/*物理量的上，下限.flag选择窗口:flag=0为波形窗口1;flag=1为波形*/
/*窗口2;flag>1为分岔图窗口.*/
{
*wv=win;
if (!flag)
wv->b=wv->t+(wv->b-wv->t)/2-15;
else if (flag==1)
wv->t=wv->t+(wv->b-wv->t)/2+15;
else
{
wv->l=70;
wv->t=flag;
}
Init_rect(wv);
wv->region[0][0]=xL;
wv->region[0][1]=xU;
wv->region[1][0]=yL;
wv->region[1][1]=yU;
wv->factor[0]=(wv->r-wv->l)/(wv->region[0][1]-wv->region[0][0]);
wv->factor[1]=(wv->b-wv->t)/(wv->region[1][1]-wv->region[1][0]);
}
/* NO.17 function */
int plot_point(DATA *x,Rect *win,int p)/*在win窗口所代表的相空间(dir[0],dir[1])平面上以颜色p画由状态向量x定义的相点.p=0为擦除点,p=15画小十字,p<0画小圆圈*/
{ int xg,yg; p=p%16;
xg=win->l+(x[dir[0]]-win->region[dir[0]][0])*win->factor[dir[0]];
yg=win->b-(x[dir[1]]-win->region[dir[1]][0])*win->factor[dir[1]];
if (xg<win->l)xg=win->l+1; if (xg>win->r)xg=win->r-1;
if (yg<win->t)yg=win->t+1; if (yg>win->b)yg=win->b-1;
if (p>=0)
{ putpixel(xg,yg,p);
if (p==15) /*画小十字*/
{ putpixel(xg-1,yg,10); putpixel(xg+1,yg,10);
putpixel(xg,yg-1,10); putpixel(xg,yg+1,10); }
}
else if (p<0)/*画小圆圈*/
{ int d=0.1*MIN(win->factor[dir[0]],win->factor[dir[1]])/2;
setcolor(9+p);
circle(xg,y
评论0