 |
ÍÀÓ×ÍÀß ÁÈÁËÈÎÒÅÊÀ - ÐÅÔÅÐÀÒÛ - Âû÷èñëåíèå ïëîùàäè ñëîæíîé ôèãóðû ìåòîäîì èìèòàöèîííîãî ìîäåëèðîâàíèÿ
Âû÷èñëåíèå ïëîùàäè ñëîæíîé ôèãóðû ìåòîäîì èìèòàöèîííîãî ìîäåëèðîâàíèÿ
Ïðèëîæåíèå 1. Òåêñò ïðîãðàììû.
Ôàéë sq.c
/*
Ïpîãpàììà SQ
îñíîâíàÿ
*/
#include <stdio.h>
#include <graphics.h>
#include <conio.h>
#include <math.h>
#include <stdlib.h>
#include "matim.c"
#include "window.c"
#include "f_integr.c"
#include "draft.c"
#include "draft_f.c"
#include "draft_e.c"
int k=20,i=15,l=270,j=140;
void main(void)
{
float
b0,b1,b2,b3,c0,c1,c2,c3,x1,x2,maxb,maxc,minb,minc,min,max,S;
int N;
do{
closegraph();
get_poly(&b3,&b2,&b1,&b0,&c3,&c2,&c1,&c0,&x1,&x2,&N);
f_max(b3,b2,b1,b0,x1,x2,&minb,&maxb);
f_max(c3,c2,c1,c0,x1,x2,&minc,&maxc);
max=(maxb>maxc)?maxb:maxc;
min=(minb<minc)?minb:minc;
S=i_num(b3,b2,b1,b0,c3,c2,c1,c0,x1,x2);
init();
draft_f(b3,b2,b1,b0,c3,c2,c1,c0,x1,x2,min,max,k,i,l,j);
draft_e(b3,b2,b1,b0,c3,c2,c1,c0,x1,x2,min,max,S,k,i+180,l+100,j+160,N);
setcolor(2);
outtextxy(0,340," Press q for exit ");
}
while (( getch()) != 'q');
}
Ôàéë matim.c
/* Ïîäïpîãpàììà ñîäåpæèò ïpîöåäópû ìàòåìàòè÷åñêîé îápàáîòêè
ôóíêöèé*/
#include <stdio.h>
#include <conio.h>
#include <math.h>
#include <stdlib.h>
#include <dos.h>
/* Âû÷èñëåíèå
ìàêñèìàëüíîãî è ìèíèìàëüíîãî
çíà÷åíèÿ
ôóíêöèè íà çàäàíîì èíòåpâàëå */
void f_max(float a3,float a2,float a1,float a0,float
x1,float x2,float *amin,float *amax)
{
float
dx,x,Fx,Fx1,Fmax,Fmin;
dx=(x2-x1)/500;
x=x1;
Fx1=a3*x*x*x+a2*x*x+a1*x+a0;
Fmax=Fx1;
Fmin=Fx1;
do {
x=x+dx;
Fx=a3*x*x*x+a2*x*x+a1*x+a0;
if
(Fx>=Fmax)
Fmax=Fx;
if
(Fx<=Fmin)
Fmin=Fx;
} while ( x<x2 );
*amin=Fmin;
*amax=Fmax;
}
/*Âû÷èñëåíèå êîpíåé êóáè÷íîãî ópàâíåíèÿ */
int f_root(float a0,float a1,float a2,float a3,float
x1,float x2,float e,float *k1,float *k2,float *k3)
{ float ku1,ku2,ku3,x,a,b;
int c=0;
x=x1;
do
{
a=a3*pow(x,3)+a2*pow(x,2)+a1*x+a0;
x+=e;
b=a3*pow(x,3)+a2*pow(x,2)+a1*x+a0;
if (a*b<0)
{ c++;
switch(c) {
case 1: ku1=x;
break;
case 2: ku2=x;
break;
case 3: ku3=x;
break;
default:
printf("\n Âíèìàíèå !!! \n Îøèáêà â matim.c (f_root).");
break;
};
}
} while (x<x2);
*k1=ku1;
*k2=ku2;
*k3=ku3;
return c;
}
float fun(float a3,float a2,float a1,float a0,float x)
{
float s;
s=a3*x*x*x+a2*x*x+a1*x+a0;
return (s);
}
Ôàéë window.c
/* Ïîäïpîãpàììû pàáîòû ñ îêíàìè*/
#include <stdio.h>
#include <conio.h>
#include <math.h>
#include <stdlib.h>
#include <dos.h>
/*ôóíêöèÿ pèñîâàíèÿ pàìêè îêíà */
void border(int sx,int sy,int ex,int ey){
int i;
for
(i=sx+1;i<ex;i++){ //ãîpèçîíòàëüíûå ëèíèè
gotoxy(i,sy);
putch(205);
gotoxy(i,ey);
putch(205);
}
for
(i=sy+1;i<ey;i++){ //âåpòèêàëüíûå ëèíèè
gotoxy(sx,i);
putch(186);
gotoxy(ex,i);
putch(186);
}
gotoxy(sx,sy);putch(201);
gotoxy(sx,ey);putch(200);
gotoxy(ex,sy);putch(187);
gotoxy(ex,ey);putch(188);
}
void talkerror(void)
{
textcolor(15);
textbackground(4);
gotoxy(1,18);
cprintf(" ATTATETION !
DATE ERROR . Press any key to continue... ");
sound(1700); delay(100); nosound();
delay(100);
sound(1400); delay(100); nosound();
getch();
gotoxy(1,18);
textcolor(15);
textbackground(1);
clreol();
}
void get_poly(float *bo3,float *bo2,float *bo1,float
*bo0,float *co3,float *co2,float *co1,float *co0,float *xo1,float *xo2,int *No)
{
float
b3,b2,b1,b0,c3,c2,c1,c0,x1,x2;
int
xb1=5,yb1=4,xb2=76,yb2=22,c,k=3,k1=10,k2=50,N;
char
bc0[5],bc1[5],bc2[5],bc3[5],cc0[5],cc1[5],cc2[5],cc3[5],x1c[5],x2c[5],nc[5];
textbackground(11);
clrscr();
window(xb1,yb1,xb2,yb2);
textcolor(15);
textbackground(1);
clrscr();
do {
textcolor(15);
textbackground(1);
gotoxy(k1,k);
puts("b3= ");
gotoxy(k1,k+1);
puts("b2= ");
gotoxy(k1,k+2);
puts("b1= ");
gotoxy(k1,k+3);
puts("b0= ");
gotoxy(k2,k);
puts("c3= ");
gotoxy(k2,k+1);
puts("c2= ");
gotoxy(k2,k+2);
puts("c1= ");
gotoxy(k2,k+3);
puts("c0= ");
gotoxy(k1,k+6);
puts("x1=");
gotoxy(k2,k+6);
puts("x2=");
gotoxy(k1,k+10); puts("QUANTITY OF ADRESSES TO RNG ");
B3: gotoxy(k1,k); puts("b3= ");
gotoxy(k1+4,k); gets(bc3);
sscanf(bc3,"%f",&b3);
if (fabs(b3)>100) { talkerror(); goto
B3; }
B2:
gotoxy(k1,k+1);
puts("b2=
");
gotoxy(k1+4,k+1); gets(bc2);
sscanf(bc2,"%f",&b2);
if (fabs(b2)>100) { talkerror(); goto
B2; }
B1: gotoxy(k1,k+2); puts("b1= ");
gotoxy(k1+4,k+2); gets(bc1);
sscanf(bc1,"%f",&b1);
if (fabs(b1)>100) { talkerror(); goto
B1; }
B0: gotoxy(k1,k+3); puts("b0= ");
gotoxy(k1+4,k+3); gets(bc0);
sscanf(bc0,"%f",&b0);
if (fabs(b0)>100) { talkerror(); goto
B0; }
C3: gotoxy(k2,k); puts("c3= ");
gotoxy(k2+4,k); gets(cc3);
sscanf(cc3,"%f",&c3);
if (fabs(c3)>100) { talkerror(); goto
C3; }
C2: gotoxy(k2,k+1); puts("c2= ");
gotoxy(k2+4,k+1); gets(cc2);
sscanf(cc2,"%f",&c2);
if (fabs(c2)>100) { talkerror(); goto
C2; }
C1: gotoxy(k2,k+2); puts("c1= ");
gotoxy(k2+4,k+2); gets(cc1); sscanf(cc1,"%f",&c1);
if (fabs(c1)>100) { talkerror(); goto
C1; }
C0:
gotoxy(k2,k+3);
puts("c0=
");
gotoxy(k2+4,k+3); gets(cc0);
sscanf(cc0,"%f",&c0);
if (fabs(c0)>100) { talkerror(); goto
C0; }
X1:
gotoxy(k1,k+6);
puts("x1=
");
gotoxy(k2,k+6); puts("x2= ");
gotoxy(k1+4,k+6); gets(x1c);
sscanf(x1c,"%f",&x1);
if (fabs(x1)>100) { talkerror(); goto
X1; }
X2: gotoxy(k2,k+6); puts("x2= ");
gotoxy(k2+4,k+6); gets(x2c); sscanf(x2c,"%f",&x2);
if (fabs(x2)>100) { talkerror(); goto
X2; }
if (x1>=x2) { talkerror(); goto X1; }
V: R: gotoxy(k1,k+10);
puts("QUANTITY OF ADRESSES TO RNG ");
gotoxy(k1+30,k+10); gets(nc); sscanf(nc,"%d",&N);
if (N>32000) { talkerror(); goto R; }
if (N<1) { talkerror(); goto V; }
textbackground(2);
gotoxy(1,18);
cprintf("
FOR CONFURMATION PRESS 'Y' ");
sound(700); delay(100); nosound(); delay(100);
sound(1400); delay(100); nosound();
delay(100);
sound(700); delay(150); nosound();
gotoxy(1,18);
} while (( getch()) != 'y');
cprintf(" O.K.
WAIT FOR MATIMATITION
");
sound(1000);
delay(200); nosound();
*bo3=b3;
*bo2=b2;
*bo1=b1;
*bo0=b0;
*co3=c3;
*co2=c2;
*co1=c1;
*co0=c0;
*xo1=x1;
*xo2=x2;
*No=N;
}
Ôàéë f_integer.c
#include <stdio.h>
#include <math.h>
#include <time.h>
/* Âû÷èñëåíèå èíòåãpàëà ÷èñëåííûì ìåòîäîì */
float i_num(float a3,float a2,float a1,float a0,float
b3,float b2,float b1,float b0,float x1,float x2)
{
float
xt,sx=0,f1,f2,e=0.01;
xt=x1;
while
(xt<x2)
{
sx=fabs(fun(a3-b3,a2-b2,a1-b1,a0-b0,xt))*e+sx;
xt=xt+e;
};
return (sx);
}
/* Ïpîöåäópà pàñ÷èòûâàþùàÿ ïëîùàäü ñëîæíîé ôèãópû
ñ ïîìîùüþ ìåòîäà
èìèòàöèîííîãî ìîäåëèpîâàíèÿ. Èç-çà ÷åãî âñå íà÷èíàëîñü...*/
float i_rand(float a3,float a2,float a1,float a0,float b3,float
b2,float b1,float b0,float x1,float x2,float fmin,float fmax,int n)
{
float
s,sn=0,f1,f2,min,max,x,y;
int i;
time_t t;
srand((unsigned)
time (&t));
//randomize();
for(i=1;i<n;i++)
{
x=x1+random(x2-x1)+random(100)*0.01;
y=fmin+random(fmax-fmin)+random(100)*0.01;
f1=a3*x*x*x+a2*x*x+a1*x+a0;
f2=b3*x*x*x+b2*x*x+b1*x+b0;
max=(f1>f2)?f1:f2;
min=(f1<f2)?f1:f2;
if (y>=min) {
if (y<=max)
sn++;
//srand((unsigned) time (&t));
}
}
s=(sn*(fmax-fmin)*(x2-x1)/n);
return s;
}
Ôàéë draft.c
/*
Ïîäïpîãpàììà DRAFT
âñå ñâÿçàíîå ñ ãpàôèêîé
*/
#include <stdio.h>
#include <graphics.h>
#include <conio.h>
#include <math.h>
#include <stdlib.h>
extern int k,i,l,j;
/* èíèöèàëèçàöèÿ ãpàôèêè */
void init(void)
{
int
driv,mode,err;
driv=DETECT;
initgraph(&driv,&mode,"");
err=graphresult();
if (err
!=grOk)
{
printf("Îøèáêà
ïpè èíèöèàëèçàöèè ãpàôèêè : %s",grapherrormsg(err));
exit(1);
}
setgraphmode(EGAHI);
return;
}
/*Ââîä ïàpàìåòpîâ ôóíêöèé
F(X)= A3*X^3 + A2*X^2 + A1*X + A0 */
void get_parms(float *a3,float *a2,float *a1,float *a0)
{
printf("Ââåäèòå
êîýôèöèåíòû A3 A2 A1 A0 \n");
scanf("%f
%f %f %f",a3,a2,a1,a0);
}
/*Îáâîäèò íåïpåpûâíûé êîíòóp */
void f_draft(float a0,float a1,float a2,float a3,float
dx,float a,float b,float x1)
{
float xt,y,x;
xt=x1-dx;
y=ceil(a*(a0+a1*x1+a2*x1*x1+a3*x1*x1*x1)+b);
moveto(k,y);
for
(x=k-1;x<l+1;x++)
{
y=a*(a0+a1*xt+a2*xt*xt+a3*xt*xt*xt)+b;
lineto(x,y);
xt+=dx;
delay(0);
}
}
/*Ðèñóåò îñè êîîpäèíàò */
void osi(float x1,float x2,float b)
{
float c;
setcolor(4);
setlinestyle(0,1,1);
settextstyle(2,HORIZ_DIR,4);
setfillstyle(3,13);
line(k-5,b,l+5,b);
c=k-x1*(l-k)/(x2-x1);
line(c,i-5,c,j+5);
/*îñü y */
outtextxy(l+10,b-2,"x");
outtextxy(c+3,i-12,"y");
outtextxy(c-10,b-10,"0");
outtextxy(l,b-3,">");
outtextxy(c-3,i-6,"^");
}
void strout(int f,float a3,float a2,float a1,float a0,int
bx,int by)
{
char s[50];
sprintf(s,"Y%d(X)=(%2.2f)*X^3+(%2.2f)*X^2+(%2.2f)*X+(%2.2f)",f,a3,a2,a1,a0);
outtextxy(bx,by,s);
}
Ôàéë draft_f.c
/*
Ïîäïpîãpàììà DRAFT
âñå ñâÿçàíîå ñ ãpàôèêîé
*/
#include <stdio.h>
#include <graphics.h>
#include <conio.h>
#include <math.h>
#include <stdlib.h>
extern int k,i,l,j;
/* èíèöèàëèçàöèÿ ãpàôèêè */
void init(void)
{
int
driv,mode,err;
driv=DETECT;
initgraph(&driv,&mode,"");
err=graphresult();
if (err
!=grOk)
{
printf("Îøèáêà
ïpè èíèöèàëèçàöèè ãpàôèêè : %s",grapherrormsg(err));
exit(1);
}
setgraphmode(EGAHI);
return;
}
/*Ââîä ïàpàìåòpîâ ôóíêöèé
F(X)= A3*X^3 + A2*X^2 + A1*X + A0 */
void get_parms(float *a3,float *a2,float *a1,float *a0)
{
printf("Ââåäèòå
êîýôèöèåíòû A3 A2 A1 A0 \n");
scanf("%f
%f %f %f",a3,a2,a1,a0);
}
/*Îáâîäèò íåïpåpûâíûé êîíòóp */
void f_draft(float a0,float a1,float a2,float a3,float
dx,float a,float b,float x1)
{
float xt,y,x;
xt=x1-dx;
y=ceil(a*(a0+a1*x1+a2*x1*x1+a3*x1*x1*x1)+b);
moveto(k,y);
for
(x=k-1;x<l+1;x++)
{
y=a*(a0+a1*xt+a2*xt*xt+a3*xt*xt*xt)+b;
lineto(x,y);
xt+=dx;
delay(0);
}
}
/*Ðèñóåò îñè êîîpäèíàò */
void osi(float x1,float x2,float b)
{
float c;
setcolor(4);
setlinestyle(0,1,1);
settextstyle(2,HORIZ_DIR,4);
setfillstyle(3,13);
line(k-5,b,l+5,b);
c=k-x1*(l-k)/(x2-x1);
line(c,i-5,c,j+5);
/*îñü y */
outtextxy(l+10,b-2,"x");
outtextxy(c+3,i-12,"y");
outtextxy(c-10,b-10,"0");
outtextxy(l,b-3,">");
outtextxy(c-3,i-6,"^");
}
void strout(int f,float a3,float a2,float a1,float a0,int
bx,int by)
{
char s[50];
sprintf(s,"Y%d(X)=(%2.2f)*X^3+(%2.2f)*X^2+(%2.2f)*X+(%2.2f)",f,a3,a2,a1,a0);
outtextxy(bx,by,s);
}
Ôàéë draft_e.c
/*
Ïîäïpîãpàììà
DRAFT_N ãpàôèê ïîãpåøíîñòè âû÷èñëåíèÿ èíòåãpàëà pàçëè÷íûìè
ìåòîäàìè
*/
#include <stdio.h>
#include <graphics.h>
#include <conio.h>
#include <math.h>
#include <stdlib.h>
/*Ôóíêöèÿ pèñóåò ãpàôèê ïîëùàäè ñëîæíîé ôèãópû â çàâèñèìîñòè
îò
êîëè÷åñòâà
èñïûòàíèé*/
void draft_e(float b3,float b2,float b1,float b0,float
c3,float c2,float c1,float c0,float x1,float x2,float min,float max,float
Sn,int k,int i,int l,int j,int n)
{
float
dx,x,y,Sr,a,xl,yl,Ss;
int v,nt;
char s[10];
setcolor(4);
setlinestyle(0,1,1);
settextstyle(2,HORIZ_DIR,4);
line(k-5,j,l+5,j);
line(k,i-5,k,j+5);
/*îñü y */
outtextxy(l+10,j-2,"N");
outtextxy(k-8,i,"S");
outtextxy(k-10,j-10,"0");
outtextxy(l,j-3,">");
outtextxy(k-3,i-6,"^");
setbkcolor(15);
setcolor(2);
line(l+50,i+110,l+100,i+110);
outtextxy(l+103,i+107,"Sr-random");
setcolor(1);
line(l+50,i+120,l+100,i+120);
outtextxy(l+103,i+117,"Sn-numeric");
dx=n/10;
a=(i-j)/(2*Sn);
y=a*Sn+j;
line(k+5,y,l-5,y);
settextstyle(2,HORIZ_DIR,4);
setcolor(5);
sprintf(s,"S=%3.2f",Sn);
outtextxy(l+120,i-40,s);
outtextxy(l+50,i-20,"N");
outtextxy(l+120,i-20,"Sr");
outtextxy(l+220,i-20,"Sn-Sr");
xl=k;
yl=j;
for(v=1;v<11;v++){
nt=ceil(v*dx);
Sr=i_rand(b3,b2,b1,b0,c3,c2,c1,c0,x1,x2,min,max,nt);
x=k+v*(l-k)/10;
y=a*Sr+j;
setcolor(2);
line(xl,yl,x,y);
xl=x;
yl=y;
setcolor(4);
settextstyle(2,VERT_DIR,4);
sprintf(s,"%d",nt);
outtextxy(x,j+3,s);
setcolor(8);
settextstyle(2,HORIZ_DIR,4);
outtextxy(l+40,i+(v-1)*10,s);
sprintf(s,"%3.2f",Sr);
outtextxy(l+110,i+(v-1)*10,s);
Ss=100-(Sr*100/Sn);
sprintf(s,"%2.1f%",Ss);
outtextxy(l+205,i+(v-1)*10,s);
}
}
ÌÎÑÊÎÂÑÊÈÉ ÃÎÑÓÄÀÐÑÒÂÅÍÍÛÉ ÈÍÑÒÈÒÓÒ ÝËÅÊÒÐÎÍÍÎÉ ÒÅÕÍÈÊÈ
(ÒÅÕÍÈ×ÅÑÊÈÉ ÓÍÈÂÅÐÑÈÒÅÒ)
Ðàñ÷åò ïëîùàäè
ñëîæíîé ôèãóðû ñ ïîìîùüþ ìåòîäà
èìèòàöåîííîãî
ìîäåëèðîâàíèÿ .
Ëîãâèíåíêî Â.
Ìîñêâà. 1995 ã.
Çàäàíèå: Ðàçðàáîòàòü
ïðîãðàììó, ïîçâîëÿþùóþ ñ ïîìîùüþ ìåòîäà èìèòàöèîííîãî ìîäåëèðîâàíèÿ ðàññ÷èòàòü
ïëîùàäü ñëîæíîé ôèãóðû, îãðàíè÷åííîé ñâåðõó êðèâîé U=Y1(x) , ñíèçó V=Y2(x).
1. Äëÿ ðåøåíèÿ äàííîé çàäà÷è ïðèìåíèì ñëåäóþùèé ìåòîä.
Îãðàíè÷èì çàäàííóþ ôèãóðó
ïðÿìîóãîëüíèêîì, ñòîðîíû êîòîðîãî ïðîõîäÿò:
÷åðåç òî÷êè
ìàêñèìàëüíîãî è ìèíèìàëüíîãî çíà÷åíèÿ ôóíêöèé è ïàðàëëåëüíû îñÿì àáñöèññ;
÷åðåç ëåâóþ è
ïðàâóþ ãðàíè÷íûå òî÷êè îáëàñòè îïðåäåëåíèÿ àðãóìåíòà è ïàðàëëåëüíû îñÿì
îðäèíàò.
Èñïîëüçóÿ äàò÷èê ñëó÷àéíûõ
÷èñåë ðàçûãðûâàþòñÿ êîîðäèíàòû ñëó÷àéíîé òî÷êè èç ýòîãî ïðÿìîóãîëüíèêà . Ïðîâåðÿåì ïîïàäàåòå òî÷êè â çàäàííóþ ôèãóðó.
Çíàÿ ïëîùàäü ïðÿìîóãîëüíèêà è îòíîøåíèå ïîïàâøèõ òî÷åê ê èõ îáùåìó ÷èñëó
ðàçûãðàííûõ, ìîæíî îöåíèòü ïëîùàäü èíòåðåñóþùåé íàñ ôèãóðû.
2. Òåõíè÷åñêèå
õàðàêòåðèñòèêè îáúåêòà èññëåäîâàíèÿ:
2.1. Äèàïàçîí çíà÷åíèé ïàðàìåòðîâ çàäà÷è.
Ìíîæåñòâî êðèâûõ îãðàíè÷èì ïîëèíîìàìè òðåòüåãî ïîðÿäêà, â âèäó òîãî ÷òî
ïîëèíîìû áîëåå âûñîêîãî ïîðÿäêà ñèëüíî óâåëè÷èâàþò âðåìÿ âû÷èñëåíèÿ. Ïðè÷åì äëÿ
íàãëÿäíîñòè ðåøåíèÿ âïîëíå äîñòàòî÷íî ïîðÿäêà "3".
Êîýôôèöèåíòû ïîëèíîìà îãðàíè÷èì äèàïàçîíîì [-100,100] .
Îáëàñòü
îïðåäåëåíèÿ îãðàíè÷èì äèàïàçîíîì [-100,100].
Ýòè
îãðàíè÷åíèÿ ââåäåíû äëÿ áîëåå íàãëÿäíîãî ðåøåíèÿ çàäà÷è, è èçìåíèòü èõ íå ñ
òåõíè÷åñêîé òî÷êè çðåíèÿ íå ñëîæíî.
3. Ðåøåíèå çàäà÷è.
Äàííàÿ çàäà÷à ðåøåíà â ñðåäå Turbo C. Äëÿ ðåøåíèÿ
ïîòðåáîâàëîñü îáùóþ çàäà÷ó ðàçáèòü íà
íåñêîëüêî íåáîëüøèõ çàäà÷ (ïðîöåäóð).
À èìåííî îòäåëüíî( â âèäå ïðîöåäóð) áûëè ðåøåíû
çàäà÷è
-ââîä
ïàðàìåòðîâ; |
ïðîöåäóðà get_poly |
|
-ñîîáùåíèå
îá îøèáêå ïðè ââîäå; | Ôàéë
WINDOW.C
ïðîöåäóðà talkerror
|
|
-ðèñîâàíèå
ðàìêè îêíà; |
ïðîöåäóðà border |
-âû÷èñëåíèå ìèíèìàëüíîãî è |
ìàêñèìàëüíîãî çíà÷åíèè ôóíêöèé ; |
ïðîöåäóðà f_max |
|
-âû÷èñëåíèå
çíà÷åíèÿ ïîëèíîìà â |
çàäàííîé òî÷êå; | Ôàéë
MATIM.C
ïðîöåäóðà fun |
|
-âû÷èñëåíèå êîðíåé êóáè÷íîãî |
óðàâíåíèÿ; |
ïðîöåäóðà
f_root |
-âû÷èñëåíèå
èíòåãðàëà ÷èñëåííûì |
ìåòîäîì; |
ïðîöåäóðà
i_num |
| Ôàéë F_INTEGER.C
-âû÷èñëåíèå èíòåãðàëà ñ ïîìîùüþ |
èìèòàöèîííîãî ìîäåëèðîâàíèÿ; |
ïðîöåäóðà i_rand |
-èíèöèàëèçàöèÿ ãðàôè÷åñêîãî ðåæèìà |
ïðîöåäóðà
init |
|
-îáâîäêà
íåïðåðûâíîãî êîíòóðà | Ôàéë DRAFT.C
ïðîöåäóðà
f_draft |
|
-
âûðèñîâêà îñåé êîîðäèíàò |
ïðîöåäóðà
osi |
-âûðèñîâêè ãðàôèêîâ ôóíêöèé è | Ôàéë DRAFT_F.C
øòðèõîâêà çàäàííîé ïëîùàäè |
ïðîöåäóðà
draft_f |
-âûðèñîâêà
ãðàôèêîâ âû÷èñëåíèÿ |
ïëîùàäè ðàçíûìè ìåòîäàìè è âûâîä | Ôàéë DRAFT_N.C
òàáëèöû ðåçóëüòàòîâ âû÷èñëåíèÿ |
ïðîöåäóðà
draft_n |
Ñõåìà àëãîðèòìà
èìååò âèä:
4. Îïèñàíèå ïðîöåäóð èñïîëüçóåìûé â ïðîãðàììå.
4.1 Ôàéë WINDOW.C.
4.1.1
Ïðîöåäóðà ââîäà ïàðàìåòðîâ.
void get_poly( float *b3,float *b2,float *b1,float
*b0, //-êîýôôèöèåíòû ïîëèíîìà Y1
fliat *c3,float *c2,float *c1,float
*c0, //-êîýôôèöèåíòû ïîëèíîìà Y2
float *x1,float *x2, // îáëàñòü îïðåäåëåíèÿ
[x1,x2]
int *N ) //
êîëè÷åñòâî îáðàùåíèé ê ãåíåðàòîðó //ñëó÷àéíûõ ÷èñåë
4.1.2
Ïðîöåäóðà ðèñîâàíèÿ ðàìêè îêíà.
void border(int sx, int sy, int en, int ey) // ðèñóåò ðàìêó ñ êîîðäèíàòàìè ëåâîãî
âåðõíåãî // óãëà (sx,sy) è êîîðäèíàòàìè ïðàâîãî íèæíåãî // óãëà (ex,ey)
4.1.3 Ïðîöåäóðà
ñîîáùåíèÿ îá îøèáêå ïðè ââîäå.
void talkerror(void) -
Ïðîöåäóðà ïîäàåò çâóêîâîé ñèãíàë è âûâîäèò íà ýêðàí
ñîîáùåíèå îá îøèáêå ïðè ââîäå.
4.2. Ôàéë MATIM.C
4.2.1 Ïðîöåäóðà
âû÷èñëåíèÿ ìàêñèìàëüíîãî è ìèíèìàëüíîãî çíà÷åíèé ôóíêöèé íà çàäàííîì èíòåðâàëå.
void f_max(float b3,float b2,float b1,float b0,
//-êîýôôèöèåíòû ïîëèíîìà Y1
fliat c3,float c2,float c1,float c0, //-êîýôôèöèåíòû ïîëèíîìà Y2
float x1,float x2, // îáëàñòü îïðåäåëåíèÿ [x1,x2]
float *amin, float *amax) // ìèíèìàëüíîå è ìàêñèìàëüíîå çíà÷åíèÿ // ôóíêöèé
4.2.2
Ïðîöåäóðà âû÷èñëåíèÿ çíà÷åíèÿ ïîëèíîìà â äàííîé òî÷êå.
float fun(float b3,float b2,float b1,float b0,
//-êîýôôèöèåíòû ïîëèíîìà
float
x)
Âîçâðàùàåò çíà÷åíèå ïîëèíîìà â òî÷êå õ.
4.2.3 Ïðîöåäóðà
âû÷èñëåíèÿ êîðíåé êóáè÷íîãî óðàâíåíèÿ.
int f_root(float b3,float b2,float b1,float b0,
//-êîýôôèöèåíòû ïîëèíîìà Y1
fliat c3,float c2,float c1,float c0, //-êîýôôèöèåíòû ïîëèíîìà Y2
float x1,float x2, // îáëàñòü îïðåäåëåíèÿ [x1,x2]
float
e, // òî÷íîñòü âû÷èñëåíèÿ êîðíåé
float *k1,float *k2,float *k3) // çíà÷åíèÿ êîðíåé // ôóíêöèé
Âîçâðàùàåò êîëè÷åñòâî äåéñòâèòåëüíûõ êîðíåé íà äàííîì èíòåðâàëå.
4.3. Ôàéë F_INTEGER.C
4.3.1 Ïðîöåäóðà
âû÷èñëåíèÿ ïëîùàäè ñëîæíîé ôèãóðû ÷èñëåííûì ìåòîäîì.
float f_num(float b3,float b2,float b1,float b0,
//-êîýôèöèåíòû ïîëèíîìà Y1
fliat c3,float c2,float c1,float c0, //-êîýôèöèåíòû ïîëèíîìà Y2
float x1,float x2) // îáëàñòü îïðåäåëåíèÿ [x1,x2]
Âû÷èñëÿåò ïëîùàäü ñëîæíîé ôèãóðû.
4.3.2 Ïðîöåäóðà âû÷èñëåíèÿ ïëîùàäè ñëîæíîé ôèãóðû
c ïîìîùüþ ìåòîäà èìèòàöèîííîãî ìîäåëðîâàíèÿ
float f_(float b3,float b2,float b1,float b0,
//-êîýôèöèåíòû ïîëèíîìà Y1
fliat c3,float c2,float c1,float c0, //-êîýôèöèåíòû ïîëèíîìà Y2
float x1,float x2, // îáëàñòü îïðåäåëåíèÿ [x1,x2]
float fmin,float fmax, // ìèíèìàëüíîå è ìàêñèìàëüíîå
çíà÷åíèÿ //ôóíêöèé íà äàííîì
èíòåðâàëå
int n) //
êîëè÷åñòâî îáðàùåíèé ê ãåíåðàòîðó // ñëó÷àéíûé ÷èñåë
Âû÷èñëÿåò ïëîùàäü ñëîæíîé ôèãóðû ñ ïîìîùüþ ìåòîäà
èìèòàöèîííîãî ìîäåëèðîâàíèÿ.
4.4 Ôàéë DRAFT.C
4.4.1
Ïðîöåäóðà èíèöèàëèçàöèè ãðàôè÷åñêîãî ðåæèìà.
void init (void)
4.4.2 Ïðîöåäóðà
îáâîäêè íåïðåðûâíîãî êîíòóðà.
void f_draft (float b0,float b1,float b2,float b3,
//-êîýôèöèåíòû ïîëèíîìà
float x1,float x2) // îáëàñòü îïðåäåëåíèÿ
[x1,x2]
4.4.3 Ïðîöåäóðà
âûðèñîâêè îñåé êîîðäèíàò.
void osi ( float x1, float x2, // îáëàñòü
îïðåäåëåíèÿ ôóíêöèé
float b) //
ìàøòàáíûé êîýôèöèåíò ðàñ÷èòûâàåòñÿ ïî ôîðìóëå
// b= j - Fmin*(i-j) / (Fmax - Fmin)
// ãäå i,j - çàäàþò
ïîëîæåíèå ãðàôèêà íà ýêðàíå
// Fmin,Fmax - ìèíèìàëüíîå è
ìàêñèìàëüíîå çíà÷åíèÿ //ôóíêöèé íà äàííîì èíòåðâàëå
4.5 Ôàéë DRAFT_F.
4.5.1
Ïðîöåäóðà âûðèñîâêè ãðàôèêîâ ôóíêöèé.
void draft_f (float b3,float b2,float b1,float b0,
//-êîýôèöèåíòû ïîëèíîìà Y1
fliat c3,float c2,float c1,float c0, //-êîýôèöèåíòû ïîëèíîìà Y2
float x1,float x2, // îáëàñòü îïðåäåëåíèÿ [x1,x2]
float fmin,float fmax, // ìèíèìàëüíîå è ìàêñèìàëüíîå
çíà÷åíèÿ //ôóíêöèé íà äàííîì
èíòåðâàëå
int k, int i, int l, int j) //
êîîðäèíàòû, çàäàþùèå ïîëîæåíèå //ãðàôèêà
íà ýêðàíå
4.6 Ôàéë DRAFT_N.
4.6.1
Ïðîöåäóðà âûðèñîâêè ãðàôèêîâ çíà÷åíèé ïîëùàäåé ðàñ÷èòàííûõ ÷èñëâûì ìåòîäîì è
ìåòîäîì èìèòàöèîííîãî ìîäåëèðîâàíèÿ â çàâèñèìîñòè îò êîëè÷åñòâà îáðàùåíèé ê
ãåíåðàòîðó ñëó÷àéíûõ ÷èñåë.
void draft_e (float b3,float b2,float b1,float b0,
//-êîýôèöèåíòû ïîëèíîìà Y1
fliat c3,float c2,float c1,float c0, //-êîýôèöèåíòû ïîëèíîìà Y2
float x1,float x2, // îáëàñòü îïðåäåëåíèÿ [x1,x2]
float fmin,float fmax, // ìèíèìàëüíîå è ìàêñèìàëüíîå
çíà÷åíèÿ //ôóíêöèé íà äàííîì
èíòåðâàëå
float Sn, // ïëîùàäü ðàññ÷èòàííàÿ ÷èñëîâûì
ìåòîäîì
int k, int i, int l, int j) //
êîîðäèíàòû, çàäàþùèå ïîëîæåíèå //ãðàôèêà
íà ýêðàíå
4.7 Ôàéë SQ.C
Âñå ôàéëû îáúåäèíåíû â ãëàâíîé ïðîãðàììå SQ.C,
êîòîðàÿ ÿâëÿåòñÿ îñíîâíîé è êîîðäèíèðóåò ðàáîòó ïðîöåäóð.
5 Èñïîëüçîâàíèå ïðîãðàììû.
Äëÿ èñïîëüçîâàíèÿ äàííîé ïðîãðàììû íåîáõîäèìà
îïåðàöèîííàÿ ñðåäà MS DOS,
ôàéë egavega.bgi, è ñîáñòâåííî ñàìà
ñêîìïèëèðîâàííàÿ ïðîãðàììà sq.exe.
6 Èñõîäíûé òåêñò ïðîãðàììû äàí â ïðèëîæåíèè ¹1.
7 Òåñîâûé ïðèìåð ïîêàçàí â ïðèëîæåíèè ¹2.
8 Ñïèñîê èñïîëüçîâàííîé ëèòåðàòóðû.
8.1 ßçûê
ïðîãðàììèðîâàíèÿ Ñè äëÿ ïåðñîíàëüíîãî êîìïüþòåðà .
Ñ.Î. Áî÷êîâ, Ä.Ì. Ñóááîòèí.
8.2 Ñ++ . Îïèñàíèå ÿçûêà ïðîãðàììèðîâàíèÿ.
Áüÿðí Ñòðàóñòðàï.
8.3 TURBO C. User's Guide. Borland International, Inc. 1988.
8.4 TURBO C. Reference Guide. Borland International, Inc. 1988.
9 Çàêëþ÷åíèå.
9.1 Ñîïîñòàâëåíèå
ðåçóëüòàòîâ ðàáîòû ñ òåðáîâàíèåì çàäàíèÿ.
Ñîïîñòàâëÿÿ ðåçóëüòàòû ðàáîòû ñ òðåáîâàíèåì
çàäàíèÿ, ìîæíî ñêàçàòü ÷òî çàäà÷à ðåøåíà â ïîëíîé ìåðå, çà èñêëþ÷åíèåì, áûòü
ìîæåò îáùíîñòè îòíîñèòåëüíî âîçìîæíîñòè ðàñ÷åòà äëÿ ìíîãèå êëàññîâ ôóíêöèé. Íî
ðåøåíèå áîëåå îáùåé çàäà÷è ( ò.å. âîçìîæíîñòü ðàñ÷åòà äëÿ ìíîãèõ êëàññîâ
ôóíêöèé ) ïðåäñòàâëÿåòñÿ çíà÷èòåëüíî áîëåå ãðîìîçäêèì, è âîîáùå ÿâëÿåòñÿ
îòäåëüíîé çàäà÷åé. Ïîýòîìó àâòîð íå ñ÷åë íóæíûì ðàçðàáàòûâàòü àëãîðèòì ââîäà
ìíîãèõ ôóíêöèé è çàîñòðèë âíèìàíèå ñîáñòâåííî íà ñàìîé çàäà÷å - ðàñ÷åòå ïëîùàäè
ñëîæíîé ôèãóðû ñ ïîìîùüþ ìåòîäà èìèòàöèîííîãî ìîäåëèðîâàíèÿ è ñðàâíåíèå ýòîãî
ìåòîäà ñ ÷èñëîâûìè ìåòîäàìè.
9.2 Ðåêîìåíäàöèè
ïî óëó÷øåíèþ ïðîãðàììû.
Ïðè ðàçðàáîòêå ïðîãðàììû àâòîð óïóñòèë âîçìîæíîñòü
ðàáîòû ñ ÷èñëîâûìè ìàññèâàìè. Ïîýòîìó, ìîæíî óëó÷øèòü ïðîãðàììó ïåðåïèñàâ ðÿä
ïðîöåäóð ïîä ìàññèâû , ÷òî ñäåëàåò ïðîãðàììó ìåíåå ìàññèâíîé è áîëåå íàãëÿäíîé.
Øèðîêîå
âîçìîæíîñòè ïî óëó÷øåíèþ
ïðîãðàììû â îáëàñòè ðàçðàáîòêè àëãîðèòìîâ ââîäà ðàçëè÷íûé êëàññîâ
ôóíêöèé.
Äëÿ ðåøåíèÿ çàäà÷è ìåòîäîì èìèòàöèîííîãî ìîäåëèðîâàíèÿ îãðàíè÷èì äàííóþ
/*
Ïîäïpîãpàììà DRAFT âñå ñâÿçàíîå ñ ãpàôèêîé
*/
#include
<stdio.h>
#include
<graphics.h>
#include
<conio.h>
#include
<math.h>
#include
<stdlib.h>
extern
int k,i,l,j;
/*
èíèöèàëèçàöèÿ ãpàôèêè */
void
init(void)
{
int driv,mode,err;
driv=DETECT;
initgraph(&driv,&mode,"");
err=graphresult();
if (err !=grOk)
{
printf("Îøèáêà ïpè
èíèöèàëèçàöèè ãpàôèêè : %s",grapherrormsg(err));
exit(1);
}
setgraphmode(EGAHI);
return;
}
/*Ââîä
ïàpàìåòpîâ ôóíêöèé
F(X)=
A3*X^3 + A2*X^2 + A1*X + A0 */
void
get_parms(float *a3,float *a2,float *a1,float *a0)
{
printf("Ââåäèòå êîýôèöèåíòû A3 A2 A1
A0 \n");
scanf("%f %f %f
%f",a3,a2,a1,a0);
}
/*Îáâîäèò
íåïpåpûâíûé êîíòóp */
void
f_draft(float a0,float a1,float a2,float a3,float dx,float a,float b,float x1)
{
float xt,y,x;
xt=x1-dx;
y=ceil(a*(a0+a1*x1+a2*x1*x1+a3*x1*x1*x1)+b);
moveto(k,y);
for (x=k-1;x<l+1;x++)
{
y=a*(a0+a1*xt+a2*xt*xt+a3*xt*xt*xt)+b;
lineto(x,y);
xt+=dx;
delay(0);
}
}
/*Ðèñóåò
îñè êîîpäèíàò */
void
osi(float x1,float x2,float b)
{
float c;
setcolor(4);
setlinestyle(0,1,1);
settextstyle(2,HORIZ_DIR,4);
setfillstyle(3,13);
line(k-5,b,l+5,b);
c=k-x1*(l-k)/(x2-x1);
line(c,i-5,c,j+5); /*îñü y */
outtextxy(l+10,b-2,"x");
outtextxy(c+3,i-12,"y");
outtextxy(c-10,b-10,"0");
outtextxy(l,b-3,">");
outtextxy(c-3,i-6,"^");
}
void
strout(int f,float a3,float a2,float a1,float a0,int bx,int by)
{
char
s[50];
sprintf(s,"Y%d(X)=(%2.2f)*X^3+(%2.2f)*X^2+(%2.2f)*X+(%2.2f)",f,a3,a2,a1,a0);
outtextxy(bx,by,s);
}
/*
Ïîäïpîãpàììà DRAFT_N ãpàôèê ïîãpåøíîñòè
âû÷èñëåíèÿ èíòåãpàëà pàçëè÷íûìè
ìåòîäàìè
*/
#include
<stdio.h>
#include
<graphics.h>
#include
<conio.h>
#include
<math.h>
#include
<stdlib.h>
/*Ôóíêöèÿ
pèñóåò ãpàôèê ïîëùàäè ñëîæíîé ôèãópû â çàâèñèìîñòè îò
êîëè÷åñòâà
èñïûòàíèé*/
void
draft_e(float b3,float b2,float b1,float b0,float c3,float c2,float c1,float
c0,float x1,float x2,float min,float max,float Sn,int k,int i,int l,int j,int
n)
{
float dx,x,y,Sr,a,xl,yl,Ss;
int v;
long int nt;
char s[10];
setcolor(4);
setlinestyle(0,1,1);
settextstyle(2,HORIZ_DIR,4);
line(k-5,j,l+5,j);
line(k,i-5,k,j+5); /*îñü y */
outtextxy(l+10,j-2,"N");
outtextxy(k-8,i,"S");
outtextxy(k-10,j-10,"0");
outtextxy(l,j-3,">");
outtextxy(k-3,i-6,"^");
setbkcolor(15);
setcolor(2);
line(l+50,i+110,l+100,i+110);
outtextxy(l+103,i+107,"Sr-random");
setcolor(1);
line(l+50,i+120,l+100,i+120);
outtextxy(l+103,i+117,"Sn-numeric");
dx=n/10;
a=(i-j)/(2*Sn);
y=a*Sn+j;
line(k+5,y,l-5,y);
settextstyle(2,HORIZ_DIR,4);
setcolor(5);
sprintf(s,"S=%3.2f",Sn);
outtextxy(l+120,i-40,s);
outtextxy(l+50,i-20,"N");
outtextxy(l+120,i-20,"Sr");
outtextxy(l+220,i-20,"Sn-Sr");
xl=k;
yl=j;
for(v=1;v<11;v++){
nt=ceil(v*dx);
Sr=i_rand(b3,b2,b1,b0,c3,c2,c1,c0,x1,x2,min,max,nt);
x=k+v*(l-k)/10;
y=a*Sr+j;
setcolor(2);
line(xl,yl,x,y);
xl=x;
yl=y;
setcolor(4);
settextstyle(2,VERT_DIR,4);
sprintf(s,"%d",nt);
outtextxy(x,j+3,s);
setcolor(8);
settextstyle(2,HORIZ_DIR,4);
outtextxy(l+40,i+(v-1)*10,s);
sprintf(s,"%3.2f",Sr);
outtextxy(l+110,i+(v-1)*10,s);
Ss=(Sr*100/Sn)-100;
sprintf(s,"%2.1f%",Ss);
outtextxy(l+205,i+(v-1)*10,s);
}
}
/*
Ïîäïpîãpàììà DRAFT_F pèñóåò ãpàôèêè ôóíêöèé
*/
#include
<stdio.h>
#include
<graphics.h>
#include
<conio.h>
#include
<math.h>
#include
<stdlib.h>
void
draft_f(float b3,float b2,float b1,float b0,float c3,float c2,float c1,float
c0,float x1,float x2,float min,float max,int k,int i,int l,int j)
{
float a,b,c,xt,dx,x,y,pby,pcy,k1,k2,k3;
int v,m;
char s[50];
v=f_root(c0-b0,c1-b1,c2-b2,c3-b3,x1,x2,0.1,&k1,&k2,&k3);
a=(i-j)/(max-min);
b=j-a*min;
setbkcolor(15);
osi(x1,x2,b);
dx=(x2-x1)/(l-k);
setcolor(1);
setlinestyle(0,1,1);
pby=a*(b0+b1*x1+b2*x1*x1+b3*x1*x1*x1)+b;
pcy=a*(c0+c1*x1+c2*x1*x1+c3*x1*x1*x1)+b;
x=(pcy>pby)?1:-1;
line(k,pcy+x,k,pby-x);
pby=(a*(b0+b1*x2+b2*x2*x2+b3*x2*x2*x2)+b);
pcy=(a*(c0+c1*x2+c2*x2*x2+c3*x2*x2*x2)+b);
x=(pcy>pby)?1:-1;
line(l,pby-x,l,pcy+x);
f_draft(c0,c1,c2,c3,dx,a,b,x1);
f_draft(b0,b1,b2,b3,dx,a,b,x1);
for(m=0;m<=v;m++)
{
switch(m) {
case 0 : xt=(v>0)?(x1+(k1-x1)/2):(x1+(x2-x1)/2); break;
case 1 : xt=(v>1)?(k1+(k2-k1)/2):(k1+(x2-k1)/2); break;
case 2 : xt=(v>2)?(k2+(k3-k2)/2):(k2+(x2-k2)/2); break;
case 3 : xt=k3+(x2-k3)/2; break;
}
pby=a*(b0+b1*xt+b2*xt*xt+b3*xt*xt*xt)+b;
pcy=a*(c0+c1*xt+c2*xt*xt+c3*xt*xt*xt)+b;
x=(xt-x1)*(l-k)/(x2-x1)+k;
y=(pby>pcy)?(pcy+(pby-pcy)/2):(pby+(pcy-pby)/2);
floodfill(x,y,1);
}
osi(x1,x2,b);
setcolor(1);
sprintf(s,"In
Limits [%2.2f,%2.2f]",x1,x2);
outtextxy(l+20,i+45,s);
line(l+60,i+83,l+110,i+83);
outtextxy(l+113,i+80,"-Y1(x)");
strout(1,b3,b2,b1,b0,l+10,i+5);
strout(2,c3,c2,c1,c0,l+10,i+25);
rectangle(l+60,i+115,l+80,i+125);
floodfill(l+70,i+119,1);
outtextxy(l+83,i+118,"-S(x)");
setcolor(2);
f_draft(c0,c1,c2,c3,dx,a,b,x1);
line(l+60,i+98,l+110,i+98);
outtextxy(l+113,i+95,"-Y2(x)");
}
#include
<stdio.h>
#include
<math.h>
#include
<time.h>
/*
Âû÷èñëåíèå èíòåãpàëà ÷èñëåííûì ìåòîäîì */
float
i_num(float a3,float a2,float a1,float a0,float b3,float b2,float b1,float
b0,float x1,float x2)
{
float xt,sx=0,f1,f2,e=0.1;
xt=x1;
while (xt<x2)
{
sx=fabs(fun(a3-b3,a2-b2,a1-b1,a0-b0,xt))*e+sx;
xt=xt+e;
};
return (sx);
}
/*
Ïpîöåäópà pàñ÷èòûâàþùàÿ ïëîùàäü ñëîæíîé ôèãópû
ñ ïîìîùüþ ìåòîäà èìèòàöèîííîãî
ìîäåëèpîâàíèÿ. Èç-çà ÷åãî âñå íà÷èíàëîñü...*/
float
i_rand(float a3,float a2,float a1,float a0,float b3,float b2,float b1,float
b0,float x1,float x2,float fmin,float fmax,int n)
{
float s,sn=0,f1,f2,min,max,x,y;
int i;
time_t t;
srand((unsigned) time (&t));
//randomize();
for(i=1;i<n;i++)
{
x=x1+random(x2-x1)+random(100)*0.01;
y=fmin+random(fmax-fmin)+random(100)*0.01;
f1=a3*x*x*x+a2*x*x+a1*x+a0;
f2=b3*x*x*x+b2*x*x+b1*x+b0;
max=(f1>f2)?f1:f2;
min=(f1<f2)?f1:f2;
if
(y>=min) {
if (y<=max)
sn++;
//srand((unsigned) time (&t));
}
}
s=(sn*(fmax-fmin)*(x2-x1)/n);
return s;
}
/*
Ïîäïpîãpàììà ñîäåpæèò ïpîöåäópû ìàòåìàòè÷åñêîé îápàáîòêè ôóíêöèé*/
#include
<stdio.h>
#include
<conio.h>
#include
<math.h>
#include
<stdlib.h>
#include
<dos.h>
/* Âû÷èñëåíèå ìàêñèìàëüíîãî è ìèíèìàëüíîãî
çíà÷åíèÿ ôóíêöèè íà çàäàíîì
èíòåpâàëå */
void
f_max(float a3,float a2,float a1,float a0,float x1,float x2,float *amin,float
*amax)
{
float dx,x,Fx,Fx1,Fmax,Fmin;
dx=(x2-x1)/500;
x=x1;
Fx1=a3*x*x*x+a2*x*x+a1*x+a0;
Fmax=Fx1;
Fmin=Fx1;
do {
x=x+dx;
Fx=a3*x*x*x+a2*x*x+a1*x+a0;
if (Fx>=Fmax)
Fmax=Fx;
if (Fx<=Fmin)
Fmin=Fx;
}
while ( x<x2 );
*amin=Fmin;
*amax=Fmax;
}
/*Âû÷èñëåíèå
êîpíåé êóáè÷íîãî ópàâíåíèÿ */
int
f_root(float a0,float a1,float a2,float a3,float x1,float x2,float e,float
*k1,float *k2,float *k3)
{ float
ku1,ku2,ku3,x,a,b;
int c=0;
x=x1;
do
{
a=a3*pow(x,3)+a2*pow(x,2)+a1*x+a0;
x+=e;
b=a3*pow(x,3)+a2*pow(x,2)+a1*x+a0;
if (a*b<0)
{ c++;
switch(c) {
case 1: ku1=x;
break;
case 2: ku2=x;
break;
case 3: ku3=x;
break;
default: printf("\n Âíèìàíèå !!! \n
Îøèáêà â matim.c (f_root).");
break;
};
}
} while (x<x2);
*k1=ku1;
*k2=ku2;
*k3=ku3;
return
c;
}
float
fun(float a3,float a2,float a1,float a0,float x)
{
float
s;
s=a3*x*x*x+a2*x*x+a1*x+a0;
return
(s);
}
/*
Ïpîãpàììà SQ îñíîâíàÿ
*/
#include
<stdio.h>
#include
<graphics.h>
#include
<conio.h>
#include
<math.h>
#include
<stdlib.h>
#include
"matim.c"
#include
"window.c"
#include
"f_integr.c"
#include
"draft.c"
#include
"draft_f.c"
#include
"draft_e.c"
int
k=20,i=15,l=270,j=140;
void
main(void)
{
float
b0,b1,b2,b3,c0,c1,c2,c3,x1,x2,maxb,maxc,minb,minc,min,max,S;
int N;
do{
closegraph();
get_poly(&b3,&b2,&b1,&b0,&c3,&c2,&c1,&c0,&x1,&x2,&N);
f_max(b3,b2,b1,b0,x1,x2,&minb,&maxb);
f_max(c3,c2,c1,c0,x1,x2,&minc,&maxc);
max=(maxb>maxc)?maxb:maxc;
min=(minb<minc)?minb:minc;
S=i_num(b3,b2,b1,b0,c3,c2,c1,c0,x1,x2);
init();
draft_f(b3,b2,b1,b0,c3,c2,c1,c0,x1,x2,min,max,k,i,l,j);
draft_e(b3,b2,b1,b0,c3,c2,c1,c0,x1,x2,min,max,S,k,i+180,l+100,j+160,N);
setcolor(2);
outtextxy(0,340," Press q for
exit ");
} while
(( getch()) != 'q');
}
/*
Ïîäïpîãpàììû pàáîòû ñ îêíàìè*/
#include
<stdio.h>
#include
<conio.h>
#include
<math.h>
#include
<stdlib.h>
#include
<dos.h>
/*ôóíêöèÿ
pèñîâàíèÿ pàìêè îêíà */
void
border(int sx,int sy,int ex,int ey){
int i;
for (i=sx+1;i<ex;i++){ //ãîpèçîíòàëüíûå
ëèíèè
gotoxy(i,sy);
putch(205);
gotoxy(i,ey);
putch(205);
}
for (i=sy+1;i<ey;i++){ //âåpòèêàëüíûå
ëèíèè
gotoxy(sx,i);
putch(186);
gotoxy(ex,i);
putch(186);
}
gotoxy(sx,sy);putch(201);
gotoxy(sx,ey);putch(200);
gotoxy(ex,sy);putch(187);
gotoxy(ex,ey);putch(188);
}
void
talkerror(void)
{
textcolor(15);
textbackground(4);
gotoxy(1,18);
cprintf("
ATTATETION ! DATE ERROR . Press any key to continue... ");
sound(1700); delay(100); nosound(); delay(100);
sound(1400); delay(100); nosound();
getch();
gotoxy(1,18);
textcolor(15);
textbackground(1);
clreol();
}
void
get_poly(float *bo3,float *bo2,float *bo1,float *bo0,float *co3,float
*co2,float *co1,float *co0,float *xo1,float *xo2,int *No)
{
float b3,b2,b1,b0,c3,c2,c1,c0,x1,x2;
int
xb1=5,yb1=4,xb2=76,yb2=22,c,k=3,k1=10,k2=50,N;
char
bc0[5],bc1[5],bc2[5],bc3[5],cc0[5],cc1[5],cc2[5],cc3[5],x1c[5],x2c[5],nc[5];
textbackground(11);
clrscr();
window(xb1,yb1,xb2,yb2);
textcolor(15);
textbackground(1);
clrscr();
do {
textcolor(15);
textbackground(1);
gotoxy(k1,k); puts("b3= ");
gotoxy(k1,k+1); puts("b2= ");
gotoxy(k1,k+2); puts("b1= ");
gotoxy(k1,k+3); puts("b0= ");
gotoxy(k2,k); puts("c3= ");
gotoxy(k2,k+1); puts("c2= ");
gotoxy(k2,k+2); puts("c1= ");
gotoxy(k2,k+3); puts("c0= ");
gotoxy(k1,k+6); puts("x1=");
gotoxy(k2,k+6); puts("x2=");
gotoxy(k1,k+10);
puts("QUANTITY OF ADRESSES TO RNG
");
B3: gotoxy(k1,k); puts("b3=
");
gotoxy(k1+4,k); gets(bc3); sscanf(bc3,"%f",&b3);
if (fabs(b3)>100) { talkerror(); goto B3; }
B2: gotoxy(k1,k+1); puts("b2=
");
gotoxy(k1+4,k+1); gets(bc2);
sscanf(bc2,"%f",&b2);
if (fabs(b2)>100) { talkerror(); goto B2; }
B1: gotoxy(k1,k+2); puts("b1=
");
gotoxy(k1+4,k+2); gets(bc1);
sscanf(bc1,"%f",&b1);
if (fabs(b1)>100) { talkerror(); goto B1; }
B0: gotoxy(k1,k+3); puts("b0=
");
gotoxy(k1+4,k+3); gets(bc0);
sscanf(bc0,"%f",&b0);
if (fabs(b0)>100) { talkerror(); goto B0; }
C3: gotoxy(k2,k); puts("c3=
");
gotoxy(k2+4,k); gets(cc3); sscanf(cc3,"%f",&c3);
if (fabs(c3)>100) { talkerror(); goto C3; }
C2: gotoxy(k2,k+1); puts("c2=
");
gotoxy(k2+4,k+1); gets(cc2);
sscanf(cc2,"%f",&c2);
if (fabs(c2)>100) { talkerror(); goto C2; }
C1: gotoxy(k2,k+2); puts("c1=
");
gotoxy(k2+4,k+2); gets(cc1);
sscanf(cc1,"%f",&c1);
if (fabs(c1)>100) { talkerror(); goto C1; }
C0: gotoxy(k2,k+3); puts("c0=
");
gotoxy(k2+4,k+3); gets(cc0);
sscanf(cc0,"%f",&c0);
if (fabs(c0)>100) { talkerror(); goto C0; }
X1: gotoxy(k1,k+6); puts("x1=
");
gotoxy(k2,k+6); puts("x2=
");
gotoxy(k1+4,k+6); gets(x1c);
sscanf(x1c,"%f",&x1);
if (fabs(x1)>100) { talkerror(); goto X1; }
X2: gotoxy(k2,k+6); puts("x2=
");
gotoxy(k2+4,k+6); gets(x2c);
sscanf(x2c,"%f",&x2);
if (fabs(x2)>100) { talkerror(); goto X2; }
if (x1>=x2) { talkerror(); goto X1; }
V: R: gotoxy(k1,k+10); puts("QUANTITY OF
ADRESSES TO RNG ");
gotoxy(k1+30,k+10); gets(nc); sscanf(nc,"%d",&N);
/*if (N>32000) { talkerror(); goto R; }
if (N<1) { talkerror(); goto V; } */
textbackground(2);
gotoxy(1,18);
cprintf(" FOR CONFURMATION PRESS 'Y' ");
sound(700); delay(100); nosound();
delay(100);
sound(1400); delay(100); nosound(); delay(100);
sound(700); delay(150); nosound();
gotoxy(1,18);
} while
(( getch()) != 'y');
cprintf(" O.K. WAIT FOR MATIMATITION ");
sound(1000); delay(200); nosound();
*bo3=b3;
*bo2=b2;
*bo1=b1;
*bo0=b0;
*co3=c3;
*co2=c2;
*co1=c1;
*co0=c0;
*xo1=x1;
*xo2=x2;
*No=N;
}
Ïðîãðàììà äëÿ âû÷èñëåíèÿ ïëîùàäè íàðèñîâàííîé ôèãóðû. Ðàñ÷¸ò êîîðäèíàò ôèãóðû ïîëó÷åííîé îáúåäèíåíèåì. Ïðîãðàììû äëÿ âû÷èñëåíèÿ ïëîùàäåé ðàçíûõ ôèãóð. Ïðîãðàììà ïî îïðåäåëåíèþ ïëîùàäè ñëîæíûõ ôèãóð. Ïðîãðàììû äëÿ ðàñ÷åòà ïëîùàäè ñëîæíûõ ôèãóð. Ïðîãðàììà âû÷èñëåíèå îáúåìà ñëîæíûõ ôèãóð. Ïðîãðàììà äëÿ âû÷èñëåíèÿ ïëîùàäè ôèãóðû. Ðàññ÷èòàòü ïëîùàäü íåïðàâèëüíîé ôèãóðû. Âû÷èñëåíèå ïëîùàäè ñëîæíîé ôèãóðû. Âû÷èñëåíèå ïëîùàäè ñëîæíîé äåòàëè. Âû÷èñëåíèå îáú¸ìà ñëîæíûõ ôèãóð. Âû÷èñëåíèå îáúåìà ñëîæíûõ ôèãóð. Âû÷èñëåíèå ïëîùàäåé ôèãóð êëàññ. Ðàñ÷åò ïëîùàäè ñëîæíûõ ôèãóð. Ïðîãðàììà âû÷èñëåíèÿ ïëîùàäè.
| |
|
|
