C语言机械原理编程,连杆运动分析图线

1、编程大作业基于 C 语言的机械原理分析图：只有库函数包含头文件：graphics.hconio.h才能观看该程序效果；若已经安装 VC+，可以打开文件夹里面的 EasyX压缩包， 并安装 EasyX文件， 也可以观看程序效果；一、角度与l3 与角 1 的关系；#include #include #include#include#define SZ 2#define PI 3.1415926float cSZSZ=0;void danwei()char s5;int i=0;int k=-5;while(k= -400)sprintf(s,%0.2f,i);outtextxy(610,k,s);

2、k-=40;i+=0.05;outtextxy(610, -420, l3/m/s);void danwei1()char s10;int i= -60;int k=-10;while(k= -400)sprintf(s,%d,i);outtextxy( -30,k,s);k-=20;i+=10;outtextxy( -85,-420, 2、 3/ 度);int main()/ 创建大小为 800 * 600 的绘图窗口initgraph(800, 600);/ 设置原点 (0, 0) 为屏幕中央（ Y 轴默认向下为正）setorigin(100, 500);/ 使用白色填充背景setbkco

3、lor(WHITE);cleardevice();setlinecolor(RGB(7, 190, 234);rectangle(0, -400,600,0);int i,j;/ 网格线for(i=0;i600;i=i+20)line(i,0,i, -400);for(j=0;j400;j=j+20)line(0, -j,600, -j);/ 输出屏幕提示settextcolor(BLACK);setfont(18,0, 宋体 );danwei();danwei1();settextcolor(RED);danwei2();settextcolor(RGB(201, 62, 0);outtex

4、txy(240, -450, 角度关系图 );float o3,o1,o2,w1=10;float temp1;float l1=0.100,l2=0.040,l4=0.200;float l3,m,m1,A,B,C,D;float a22=0,b21=0,a122=0;float k,t,m10=0,t10= -16.2602*40/20 -120,m11=0,t11=(-(16.2602+90)*40/20 -120,m12=0,t12= -0.22*40/(0.05) -80;float n=0;while(n=360)o1=n*PI/180;A=l1*sin(o1) -l4;B=l1*

5、cos(o1);C=-l2;D=(A+sqrt(pow(A,2)+pow(B,2) -pow(C,2)/(B -C);o2=2*atan(D);m=o2*180/PI;l3=(l1*cos(o1) -l2*cos(o2)/sin(o2);k=m;m1=n*600/360;t=-k*40/20 -120;setcolor(BLACK);line(m10,t10,m1,t);m10=m1,t10=t;k=m+90;m1=n*600/360;t=-k*40/20 -120;setcolor(GREEN);line(m11,t11,m1,t);m11=m1,t11=t;k=l3;m1=n*600/36

6、0;t=-k*40/(0.05) -80;setcolor(MAGENTA);line(m12,t12,m1,t);m12=m1,t12=t;n+=0.1;settextcolor(GREEN);outtextxy(160, -80, 2);outtextxy(160, -180, L3);outtextxy(160, -260, 3);outtextxy(500,60, 莫智斌 /2014124066);getch();closegraph();return 0;二、角速度、 vl3 与角一的关系；#include #include #include#include#define SZ 2#

7、define PI 3.1415926float cSZSZ=0;void matrixmultiplication(float a_222,float b_221,float y21)int i,j,k,a_column=2,a_row=2,b_row=1;float temp;for(i=0;ia_column;i+)for(j=0;jb_row;j+) temp=0;for(k=0;ka_row;k+)temp+=a_2ik*b_2kj;yij=temp;void danwei()/ 横坐标char s5;int i=0;int k=-5;while(k= -400)sprintf(s,

8、%0.2f,i);outtextxy(610,k,s);k-=40;i+=0.5;outtextxy(610, -420, Vl3/m/s);void danwei1() / 左边边的纵坐标char s10;int i= -14;int k=-10;while(k= -400)sprintf(s,%d,i);outtextxy( -30,k,s);k-=40;i+=2;outtextxy( -85,-420, w2&w3/rad/s);void qiuni(float a_122,float ni22)int N=2;int i,j,k;float cSZSZ=0;float aSZSZ=0;

9、for(i=0;i2;i+)for(j=0;j2;j+)aij=a_1ij;float bSZ2*SZ=0;float t;for(i=0;iN;i+)for(j=0;jN;j+)bij=aij;for(i=0;iN;i+)for(j=N;j2*N;j+)bij=0;for(i=0;iN;i+)biN+i=1;for(k=0;kk)bik=t;for(j=0;jk;j+)t=bkj;bkj=bij;bij=t;for(j=k+1;j2*N;j+)t=bkj;bkj=bij;bij=t;for(i=k+1;i=k;j -)bij -=(bkj/bkk)*bik;for(j=2*N -1;j=k;

10、j -)bkj/=bkk;/ 第二步转化 k=N-1;while(k0)for(i=0;i=k;j -)bij -=bik*bkj;k-;for(i=0;iN;i+)for(j=0;jN;j+)cij=biN+j;for(i=0;i2;i+)for(j=0;j2;j+)niij=cij;int main()/ 创建大小为 800 * 600 的绘图窗口initgraph(800, 600);/ 设置原点 (0, 0) 为屏幕中央（ Y 轴默认向下为正）setorigin(100, 500);/ 使用白色填充背景setbkcolor(WHITE);cleardevice();setlinecol

11、or(RGB(7, 190, 234);rectangle(0, -400,600,0);int i,j;/ 网格线for(i=0;i600;i=i+20)line(i,0,i, -400);for(j=0;j400;j=j+20)line(0, -j,600, -j);/ 输出屏幕提示settextcolor(BLACK);setfont(18,0, 宋体 );danwei();danwei1();settextcolor(RED);danwei2();settextcolor(RGB(201, 62, 0);outtextxy(240, -450, 角速度、速度关系图);float o3,

12、o1,o2,w1=10;float temp1;float l1=0.100,l2=0.040,l4=0.200;float l3,m,m1,A,B,C,D;float a22=0,b_121=0,a122=0;float ni22=0;float x21=0;float k,t,m10=0,t10= -1.2727*40/2 -280,m12=0,t12=0.909*40/(0.5) -160;float n=0;while(n=360)o1=n*PI/180;A=l1*sin(o1) -l4;B=l1*cos(o1);C=-l2;D=(A+sqrt(pow(A,2)+pow(B,2) -p

13、ow(C,2)/(B -C);o2=2*atan(D);m=o2*180/PI;l3=(l1*cos(o1) -l2*cos(o2)/sin(o2);a00=sin(o2);a01=l3*cos(o2) -l2*sin(o2);/ 赋值矩阵a10= -cos(o2);a11=l3*sin(o2)+l2*cos(o2);b_100= -l1*sin(o1); b_110=l1*cos(o1);qiuni(a,ni);matrixmultiplication(ni,b_1,x);for(i=0;i2;i+)xi0=10*xi0;k=x10;m1=n*600/360;t=-k*40/2 -280;s

14、etcolor(BLACK);line(m10,t10,m1,t);m10=m1,t10=t;k=x00;m1=n*600/360;t=-k*40/(0.5) -160;setcolor(MAGENTA);line(m12,t12,m1,t);m12=m1,t12=t;n+=0.1;settextcolor(GREEN);outtextxy(130, -110, V3);outtextxy(120, -200, w3 、 w4);outtextxy(500,60, 莫智斌 /2014124066);getch();closegraph();return 0;三、角加速度、al3 与角一的关系；

15、#include #include #include#include#define SZ 2#define PI 3.1415926float cSZSZ=0;void matrixmultiplication(float a_222,float b_221,float y21)int i,j,k,a_column=2,a_row=2,b_row=1;float temp;for(i=0;ia_column;i+)for(j=0;jb_row;j+) temp=0;for(k=0;ka_row;k+)temp+=a_2ik*b_2kj;yij=temp;void danwei()/ 横坐标ch

16、ar s5;int i=0;int k=-5;while(k= -400)sprintf(s,%d,i);outtextxy(610,k,s);k-=40;i+=3;outtextxy(610, -420, al3/m/s);void danwei1() / 左边边的纵坐标char s10;int i= -150;int k=-10;while(k= -400)sprintf(s,%d,i);outtextxy( -38,k,s);k-=40;i+=33;outtextxy( -85,-420, a2&a3/rad/s);void qiuni(float a_122,float ni22)in

17、t N=2;int i,j,k;float cSZSZ=0;float aSZSZ=0;for(i=0;i2;i+)for(j=0;j2;j+)aij=a_1ij;float bSZ2*SZ=0;float t;for(i=0;iN;i+)for(j=0;jN;j+)bij=aij;for(i=0;iN;i+)for(j=N;j2*N;j+)bij=0;for(i=0;iN;i+)biN+i=1;for(k=0;kk)bik=t;for(j=0;jk;j+)t=bkj;bkj=bij;bij=t;for(j=k+1;j2*N;j+)t=bkj;bkj=bij;bij=t;for(i=k+1;i

18、=k;j -)bij -=(bkj/bkk)*bik;for(j=2*N -1;j=k;j -)bkj/=bkk;/ 第二步转化 k=N-1;while(k0)for(i=0;i=k;j -)bij -=bik*bkj;k-;for(i=0;iN;i+)for(j=0;jN;j+)cij=biN+j;for(i=0;i2;i+)for(j=0;j2;j+)niij=cij;int main()/ 创建大小为 800 * 600 的绘图窗口initgraph(800, 600);/ 设置原点 (0, 0) 为屏幕中央（ Y 轴默认向下为正）setorigin(100, 500);/ 使用白色填充

19、背景setbkcolor(WHITE);cleardevice();setlinecolor(RGB(7, 190, 234);rectangle(0, -400,600,0);int i,j;/ 网格线for(i=0;i600;i=i+20)line(i,0,i, -400);for(j=0;j400;j=j+20)line(0, -j,600, -j);/ 输出屏幕提示settextcolor(BLACK);setfont(18,0, 宋体 );danwei();danwei1();settextcolor(RED);danwei2();settextcolor(RGB(201, 62,

20、0);outtextxy(240, -450, 角加速度、加速度关系图);float a22=0;float b_121=0;float ni22=0;float o3,o1,o2,w1=10;float temp;float l1=0.100,l2=0.040,l4=0.200;float l3,m,m1,A,B,C,D;float x21=0;float w2,l3_1;floata_A22=0,a_C21=0,b_A21=0,b_B21=0,b_C21=0,jieguo210=0;float k,t,m10=0,t10=32.8234*40/33-181.818,m12=0,t12=3.

21、75657*40/(3) -106.667;float n=0;while(n=360)o1=n*PI/180;A=l1*sin(o1) -l4;B=l1*cos(o1);C=-l2;D=(A+sqrt(pow(A,2)+pow(B,2) -pow(C,2)/(B -C);o2=2*atan(D);m=o2*180/PI;l3=(l1*cos(o1) -l2*cos(o2)/sin(o2);a00=sin(o2);a01=l3*cos(o2) -l2*sin(o2);a10= -cos(o2);a11=l3*sin(o2)+l2*cos(o2);b_100= -l1*sin(o1); b_110=l1*cos(o1);qiuni(a,ni);matrixmultiplication(ni,b_1,x);for( i=0;i2;i+)xi0=w1*xi0;w2=x10;l3_1=x00;a