机械机构计算机辅助分析与设计试题(新)

机械机构计算机辅助 分析与设计 院院 系：系： 机电工程学院机电工程学院 专专 业：业：农业生物环境与能源功能农业生物环境与能源功能 姓姓 名：名： 刘海洋刘海洋 学学 号：号： 20102050201020502424 指导指导教教师：师： 王春光王春光 时时 间：间： 20102010 年年 7 7 月月 2020 日日 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 1 机械机构计算机辅助分析作业机械机构计算机辅助分析作业 第一题第一题 二分法二分法 文件名：文件名： erfenfaerfenfa 程序：程序： clear; clear; clcclc; ; n=0; n=0; a=input(a=input(请输入图形下限：请输入图形下限：x1=x1=); ); b=input(b=input(请输入图形上限：请输入图形上限：x2=x2=);); x=linspace(a,b,200);x=linspace(a,b,200); y=x.*xy=x.*x- -x x- -2;2; plot(x,y,plot(x,y,- -b); b); xlabel(Xxlabel(X 值值);); ylabel(Yylabel(Y 值值); ); title(title(f(x)=f(x)=xx2 2- -x x- -2 2 的图形的图形); ); grid ongrid on; ; c=input(c=input(请输入区间下限：请输入区间下限：a =a =);); d=input(d=input(请输入区间上限：请输入区间上限：b =b =); ); jingdu= input(jingdu= input(请输入精度：请输入精度：jingdu =jingdu =); ); f1=c2f1=c2- -c c- -2;2; f2=d2f2=d2- -d d- -2; 2; cha=abs(f2cha=abs(f2- -f1);f1); while chajingdu while chajingdu fen=(c+d)/2; fen=(c+d)/2; f3=fen*fenf3=fen*fen- -fenfen- -2; 2; chen=f1*f3; chen=f1*f3; if if chenjingdu y(n)=exp(1)y(n)=exp(1)- -exp(x(n); exp(x(n); Y(n)=Y(n)=- -exp(x(n); exp(x(n); n=n+1;n=n+1; x(n)=x(nx(n)=x(n- -1)1)- -y(ny(n- -1)/Y(n1)/Y(n- -1); 1); endend X X = = x(n) x(n) 运行结果：运行结果： 请输入图形下限：请输入图形下限：x1 = 0x1 = 0 请输入图形上限：请输入图形上限：x2 = 2x2 = 2 请输入起始位置：请输入起始位置：x0 = 0.7x0 = 0.7 请输入精度：请输入精度：jingdu = 0.001jingdu = 0.001 X =X = 1.00001.0000 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 4 第三题第三题 梯形法梯形法 文件名：文件名： tixingfatixingfa 程序：程序： c clearlear; ; c clclc; ; x1=input(x1=input(请输入积分下限：请输入积分下限：x1 =x1 =);); x2=input(x2=input(请输入积分上限：请输入积分上限：x2 =x2 =);); x=linspace(x1,x2,100); x=linspace(x1,x2,100); y=sin(x);y=sin(x); plot(x,y,plot(x,y,- -m); m); grid on; grid on; legend(sin(x);legend(sin(x); title(title( f(x)=sin(f(x)=sin(x x) )的图形的图形);); xlabel(Xxlabel(X 值值);); ylabel(Yylabel(Y 值值); ); h=(x2h=(x2- -x1)/100;x1)/100; Sum=2*sum(y)Sum=2*sum(y)- -y(1)y(1)- -y(100);y(100); area=area= Sum*h/2Sum*h/2 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 5 运行结果：运行结果： 请输入积分下限：请输入积分下限：x1 = 0x1 = 0 请输入积分上限：请输入积分上限：x2 = 1x2 = 1 area =area = 0.45970.4597 第四题第四题：欧拉法欧拉法 文件名：文件名： oulafaoulafa 程序：程序： clear; clear; clc;clc; x1=input(x1=input(请输入区间下限：请输入区间下限：x1 =x1 =);); x2=input(x2=input(请输入区间上限：请输入区间上限：x2 =x2 =);); h=input(h=input(请输入步长：请输入步长：h =h =);); x0=input(x0=input(请输入初始请输入初始 X X 值：值：x0 =x0 =);); y0=input(y0=input(请输入初始请输入初始 Y Y 值：值：y0 =y0 =);); N=(x2N=(x2- -x1)/h;x1)/h; x(1)=x0;x(1)=x0; y(1)=y0;y(1)=y0; for n=1:Nfor n=1:N x(n+1)=x(n)+h;x(n+1)=x(n)+h; 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 6 y(n+1)=y(n)+h*(y(n)y(n+1)=y(n)+h*(y(n)- -2*x(n)/y(n);2*x(n)/y(n); endend n=input(n=input(请输入拟合多项式的最高次幂：请输入拟合多项式的最高次幂：n =n =); ); p=polyfit(x,y,n)p=polyfit(x,y,n); ; yi=polyval(p,x);yi=polyval(p,x); plot(x,y,*r,x,yi,plot(x,y,*r,x,yi,- -b); b); legend(legend(采样点采样点, ,拟合曲线拟合曲线); ); xlabel(Xxlabel(X 值值);); ylabel(Yylabel(Y 值值); ); grid on;grid on; f=poly2sym(p,X); f=poly2sym(p,X); F=vpa(f,4) F=vpa(f,4) 运行结果：运行结果： 请输入区间下限：请输入区间下限：x1 = 0x1 = 0 请输入区间上限：请输入区间上限：x2 = 2x2 = 2 请输入步长：请输入步长：h = 0.1h = 0.1 请输入初始请输入初始 X X 值：值：x0 = 0x0 = 0 请输入初始请输入初始 Y Y 值：值：y0 = 1y0 = 1 请输入拟合多项式的最高次幂：请输入拟合多项式的最高次幂：n = 4n = 4 p =p = - -0.0063 0.1143 0.0063 0.1143 - -0.3309 1.0059 1.00160.3309 1.0059 1.0016 F =F = - -.6276e.6276e- -2*X4+.1143*X32*X4+.1143*X3- -.3309*X2+1.006*X+1.002.3309*X2+1.006*X+1.002 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 7 第第五五题题 改进欧拉法改进欧拉法 文件名：文件名： gaijinoulafagaijinoulafa 程序：程序： clear; clear; clc; clc; x1=input(x1=input(请输入区间下限：请输入区间下限：x1 =x1 =);); x2=input(x2=input(请输入区间上限：请输入区间上限：x2 =x2 =); ); h=input(h=input(请输入步长：请输入步长：h =h =);); x0=input(x0=input(请输入初始请输入初始 X X 值：值：x0 =x0 =);); y0=input(y0=input(请输入初始请输入初始 Y Y 值：值：y0 =y0 =);); p(1)=0;p(1)=0; N=(x2N=(x2- -x1)/h;x1)/h; x(1)=x0;x(1)=x0; y(1)=y0;y(1)=y0; for n=1:Nfor n=1:N x(n+1)=x(n)+h;x(n+1)=x(n)+h; y(n+1)=y(n)+h*(y(n)y(n+1)=y(n)+h*(y(n)- -2*x(n)/y(n);2*x(n)/y(n); p(n+1)=y(n)+h*(y(n+1)p(n+1)=y(n)+h*(y(n+1)- -2*x(n)/y(n+1);2*x(n)/y(n+1); y(n+1)=(y(n+1)+p(n+1)/2;y(n+1)=(y(n+1)+p(n+1)/2; endend n=input(n=input(请输入拟合多项式的最高次幂：请输入拟合多项式的最高次幂：n =n =);); p=polyfit(x,y,n)p=polyfit(x,y,n); ; yi=polyval(p,x);yi=polyval(p,x); plot(x,y,*r,x,yi,plot(x,y,*r,x,yi,- -b); b); legend(legend(采样点采样点, ,拟合曲线拟合曲线); ); xlabel(Xxlabel(X 值值);); ylabel(Yylabel(Y 值值); ); grid on;grid on; f=poly2sym(p,X); f=poly2sym(p,X); F=vpa(f,4) F=vpa(f,4) 运行结果：运行结果： 请输入区间下限：请输入区间下限：x1 = 0x1 = 0 请输入区间上限：请输入区间上限：x2 = 2x2 = 2 请输入步长：请输入步长：h = 0.1h = 0.1 请输入初始请输入初始 X X 值：值：x0 = 0x0 = 0 请输入初始请输入初始 Y Y 值：值：y0 = 1y0 = 1 请输入拟合多项式的最高次幂：请输入拟合多项式的最高次幂：n = 4n = 4 p = p = 0.0272 0.0843 0.0272 0.0843 - -0.2612 1.0439 1.00190.2612 1.0439 1.0019 F =F = .2721e.2721e- -1*X4+.8428e1*X4+.8428e- -1*X31*X3- -.2612*X2+1.044*X+1.002.2612*X2+1.044*X+1.002 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 8 第第六六题题 黄金分割法黄金分割法 文件名：文件名： huangjinfengefahuangjinfengefa 程序：程序： clear; clear; clc;clc; x1=input(x1=input(请输入图形下限：请输入图形下限：x1 =x1 =);); x2=input(x2=input(请输入图形上限：请输入图形上限：x2 =x2 =);); x=linspace(x1,x2,1000); x=linspace(x1,x2,1000); p=3,p=3,- -16,30,16,30,- -24,8; 24,8; v=polyval(p,x); v=polyval(p,x); plot(x,v,plot(x,v,- -r); r); xlabel(Xxlabel(X 值值);); ylabel(Yylabel(Y 值值); ); grid on;grid on; a=input(a=input(请输入区间下限：请输入区间下限：a =a =);); b=input(b=input(请输入区间上限：请输入区间上限：b =b =); ); while while abs(babs(b- -a)0.001a)0.001 x1=a+0.618*(bx1=a+0.618*(b- -a);a); x2=a+0.382*(bx2=a+0.382*(b- -a);a); f1= 3*x14f1= 3*x14- -16*x13+30*x1216*x13+30*x12- -24*x1+8;24*x1+8; 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 9 f2= 3*x24f2= 3*x24- -16*x23+30*x2216*x23+30*x22- -24*x2+8;24*x2+8; if f1f2if f1f2 b=x1;b=x1; x1=x2;x1=x2; x2=a+0.382*(bx2=a+0.382*(b- -a);a); elseelse a=x2;a=x2; x2=x1;x2=x1; x1=a+0.618*(bx1=a+0.618*(b- -a); a); endend endend X X=(x1+x2)/2=(x1+x2)/2 F F= 3*= 3*X X 44- -16*16*X X3+30*3+30*X X22- -24*24*X X+8+8 运行结果：运行结果： 请输入图形下限：请输入图形下限：x1 = 0x1 = 0 请输入图形上限：请输入图形上限：x2 = 2.6x2 = 2.6 请输入区间下限：请输入区间下限：a = 1.8a = 1.8 请输入区间上限：请输入区间上限：b = 2.2b = 2.2 x =x = 2.00022.0002 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 10 第第七七题题 斐波那契法斐波那契法 文件名：文件名： feibonaqifafeibonaqifa 程序：程序： clear;clear; clc; clc; a=input(a=input(请输入区间下限：请输入区间下限：a =a =);); b=input(b=input(请输入区间上限：请输入区间上限：b =b =);); x=linspace(a,b,150); x=linspace(a,b,150); y=x+20./x;y=x+20./x; plot(x,y,plot(x,y,- -c); c); xlabel(Xxlabel(X 值值);); ylabel(Yylabel(Y 值值); ); grid on;grid on; jingdu=input(jingdu=input(请输入精度：请输入精度：jingdujingdu = =); ); n=4;n=4; F(1)=1;F(1)=1; F(2)=1;F(2)=1; F(3)=2;F(3)=2; F(4)=3;F(4)=3; x1=a+F(3)/F(4)*(bx1=a+F(3)/F(4)*(b- -a);a); x2=a+F(2)/F(4)*(bx2=a+F(2)/F(4)*(b- -a); a); while abs(x2while abs(x2- -x1)jingdux1)jingdu n=n+1;n=n+1; F(n)=F(nF(n)=F(n- -1)+F(n1)+F(n- -2);2); x1=a+F(nx1=a+F(n- -1)/F(n)*(b1)/F(n)*(b- -a);a); x2=a+F(nx2=a+F(n- -2)/F(n)*(b2)/F(n)*(b- -a); a); f1=x1+20/x1;f1=x1+20/x1; f2=x2+20/x2; f2=x2+20/x2; if f1if f1f2f2 b=x1;b=x1; elseelse a=x2;a=x2; endend endend X=(x1+x2)/2X=(x1+x2)/2 F=X+20/XF=X+20/X 运行结果：运行结果： 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 11 请输入区间下限：请输入区间下限：a = 3.7a = 3.7 请输入区间上限：请输入区间上限：b = 5b = 5 请输入精度：请输入精度：jingdu = 0.001jingdu = 0.001 X =X = 4.47254.4725 F =F = 8.94438.9443 第第八八题题 进退法进退法 文件名：文件名： jintuifajintuifa 程序：程序： clear; clear; clc;clc; x1=input(x1=input(请输入图形下限：请输入图形下限：x1 =x1 =); ); x2=input(x2=input(请输入图形上限：请输入图形上限：x2 =x2 =);); x=linspace(x1,x2,1000); x=linspace(x1,x2,1000); p=3,p=3,- -16,30,16,30,- -2,8; 2,8; v=polyval(p,x); v=polyval(p,x); plot(x,v,plot(x,v,- -r); r); 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 12 xlabel(Xxlabel(X 值值);); ylabel(Yylabel(Y 值值); ); grid on;grid on; h=input(h=input(请输入步长：请输入步长：h =h =);); x0=input(x0=input(请输入起始位置：请输入起始位置：x0 =x0 =);); jingdu=input(jingdu=input(请输入精度：请输入精度：jingdu =jingdu =);); x1=x0;x1=x0; x2=x1+h; x2=x1+h; f1= 3*x14f1= 3*x14- -16*x13+30*x1216*x13+30*x12- -2*x1+8;2*x1+8; f2= 3*x24f2= 3*x24- -16*x23+30*x2216*x23+30*x22- -2*x2+8; 2*x2+8; while while abs(f2abs(f2- -f1)f1)jingdujingdu if f1f2if f1f2 x1=x2;x1=x2; x2=x1+h; x2=x1+h; elseelse x1=x2;x1=x2; x2=x1x2=x1- -h/4; h/4; endend f1= 3*x14f1= 3*x14- -16*x13+16*x13+30*x1230*x12- -2*x1+82*x1+8 f2= 3*x24f2= 3*x24- -16*x23+30*x2216*x23+30*x22- -2*x2+82*x2+8 endend X X=(x1+x2)/2=(x1+x2)/2 F F=3*=3*X X44- -16*16*X X3+30*3+30*X X22- -2*2*X X+8+8 运行结果：运行结果： 请输入图形下限：请输入图形下限：x1 = x1 = - -1 1 请输入图形上限：请输入图形上限：x2 = 3x2 = 3 请输入步长：请输入步长：h = 0.01h = 0.01 请输入起始位置：请输入起始位置：x0 = x0 = - -0.50.5 请输入精度：请输入精度：jingdu = 0.002jingdu = 0.002 X =X = 0.03500.0350 F =F = 7.96617.9661 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 13 第第九九题题 牛顿法牛顿法 文件名：文件名：niudunfaniudunfa 程序：程序： clear; clear; clc; clc; a=input(a=input(请输入图形下限：请输入图形下限：x1 =x1 =);); b=input(b=input(请输入图形上限：请输入图形上限：x2 =x2 =);); m=linspace(a,b,200); m=linspace(a,b,200); n=m.2n=m.2- -16*m+35;16*m+35; plot(m,n,plot(m,n,- -c)c); ; xlabel(Xxlabel(X 值值);); ylabel(Yylabel(Y 值值); ); grid ongrid on; ; x0=input(x0=input(请输入起始位置：请输入起始位置：x0 =x0 =);); jingdu=input(jingdu=input(请输入精度：请输入精度：jingdu =jingdu =);); x(1)=x0;x(1)=x0; y(1)=2*x(1)y(1)=2*x(1)- -16; 16; x(2)=x(1)x(2)=x(1)- -y(1)/2; y(1)/2; n=2;n=2; while abs(x(n)while abs(x(n)- -x(nx(n- -1)jingdu1)jingdu 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 14 y(n)=2*x(n)y(n)=2*x(n)- -16;16; n=n+1;n=n+1; x(n)=x(nx(n)=x(n- -1)1)- -y(ny(n- -1)/2; 1)/2; endend X X=x(n) =x(n) Y=Y= X X.2.2- -16*16*X X+35+35 运行结果：运行结果： 请输入图形下限：请输入图形下限：x1 = 5x1 = 5 请输入图形上限：请输入图形上限：x2 = 10x2 = 10 请输入起始位置：请输入起始位置：x0 = 7x0 = 7 请输入精度：请输入精度：jingdu = 0.001jingdu = 0.001 X =8X =8 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 15 第第十十题题 筛分筛分机构机构 文件名：文件名： shaifenshaifen 程序：程序： clear;clear; clcclc; ; while while 1 1 angle=input(angle=input(请输入角度：请输入角度：angleangle =);=); if abs(angle)30if abs(angle)30 break;break; endend figure;figure; alpha=angle*pi/180;alpha=angle*pi/180; theta=linspace(0.0001,30*pi/180,100);theta=linspace(0.0001,30*pi/180,100); beta=30*(pi/180);beta=30*(pi/180); k1=sin(betak1=sin(beta- -alpha)./cos(thetaalpha)./cos(theta- -alpha+beta);alpha+beta); k2=sin(beta+alpha)./cos(thetak2=sin(beta+alpha)./cos(theta- -alphaalpha- -beta);beta); k3=cos(alpha)./abs(sin(thetak3=cos(alpha)./abs(sin(theta- -alpha);alpha); theta1=theta*180/pi;theta1=theta*180/pi; plot(theta1,k1,plot(theta1,k1,- -b,theta1,k2,:g,theta1,k3,b,theta1,k2,:g,theta1,k3,- -r)r); ; xlabel(thetaxlabel(theta 值值) ); ; ylabel(kylabel(k 值值) ); ; title(k1,ktitle(k1,k2,k32,k3 曲线图曲线图) ); ; legend(k1legend(k1 左移左移,k2,k2 右移右移,k3,k3 跃起跃起) ); ; axis(0 30 0 8)axis(0 30 0 8); ; grid on;grid on; endend 运行结果：运行结果： 请输入角度：请输入角度：angle = angle = - -2525 请输入角度：请输入角度：angle = 25angle = 25 请输入角度：请输入角度：angle = angle = - -1010 请输入角度：请输入角度：angle = 10angle = 10 请输入角度：请输入角度：angle = 0angle = 0 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 16 当当 = =- -2525 时：时： 当当 =25=25 时：时： 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 17 当当 = =- -1010 时：时： 当当 = =1010 时：时： 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 18 当当 = =0 0 时：时： 第十第十一一题题 清选清选机构机构 文件名：文件名： qingxuanqingxuan 程序：程序： clear; clear; clc; clc; w=input(w=input(请输入风速请输入风速: : windspeed =windspeed =);); H=input(H=input(请输入清选仓的高度请输入清选仓的高度: : H =H =);); L=input(L=input(请输入清选仓的宽度请输入清选仓的宽度: : L =L =);); v0=input(v0=input(请输入物料的初速度请输入物料的初速度: : v0 =v0 =);); a=10*pi/180; a=10*pi/180; b=a;b=a; g=9.8; g=9.8; t=0.001;t=0.001; h=w/2000; h=w/2000; for j=1:1:4for j=1:1:4 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 19 u=0;u=0; if if j j =1 1 xL=24;xL=24; elseif j=2elseif j=2 xL=10xL=10; ; elseif j=3elseif j=3 xL=6.5; xL=6.5; elseelse xL=3; xL=3; end end for i=1:1:2000for i=1:1:2000 vx=v0*cos(a); vx=v0*cos(a); vy=vy=- -v0*sin(a);v0*sin(a); u=u+h; u=u+h; x1=0; x1=0; y1=0;y1=0; while x10 a(i)=atan(sqrt(1a(i)=atan(sqrt(1- -N2)/N);N2)/N); elseelse a(i)=atan(sqrt(1a(i)=atan(sqrt(1- -N2)/N)+pi;N2)/N)+pi; endend a1(i)=a(i)+aC(i);a1(i)=a(i)+aC(i); xA(i)=L1*cos(a1(i);xA(i)=L1*cos(a1(i); yA(i)=L1*sin(a1(i);yA(i)=L1*sin(a1(i); if (xC(i)if (xC(i)- -xA(i)=0xA(i)=0 if (yC(i)if (yC(i)- -yA(i)0yA(i)0 a2(i)=pi/2;a2(i)=pi/2; elseelse a2(i)=3*pi/2;a2(i)=3*pi/2; endend elseelse a2(i)=atan(yC(i)a2(i)=atan(yC(i)- -yA(i)/(xC(i)yA(i)/(xC(i)- -xA(i);xA(i); if (xC(i)if (xC(i)- -xA(i)0xA(i)0 if (yC(i)if (yC(i)- -yA(i)yA(i)0)0 a2(i)=a2(i);a2(i)=a2(i); elseelse a2(i)=a2(i)+2*pi;a2(i)=a2(i)+2*pi; endend elseelse a2(i)=a2(i)+pi;a2(i)=a2(i)+pi; endend endend a3(i)=a2(i)a3(i)=a2(i)- -a21; a21; xB(i)=xA(i)+L3*cos(a3(i);xB(i)=xA(i)+L3*cos(a3(i); yB(i)=yA(i)+L3*sin(a3(i);yB(i)=yA(i)+L3*sin(a3(i); endend plot(real(xB),real(yB); plot(real(xB),real(yB); title(Btitle(B 点运动轨迹点运动轨迹);); ylabel(Yylabel(Y 值值);); xlabel(Xxlabel(X 值值);); grid on;grid on; 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 2626 运行结果：运行结果： 请输入圆心的横坐标请输入圆心的横坐标: x1 = : x1 = - -400400 请输入圆心的纵坐标请输入圆心的纵坐标: y1 = : y1 = - -400400 请输入圆心的半径请输入圆心的半径: R = 30: R = 30 请输入请输入 CC 点转速点转速: w = 1: w = 1 请输入请输入 OAOA 杆的长度杆的长度: L1 = 300: L1 = 300 请输入请输入 ACAC 杆的长度杆的长度: L2 = 50: L2 = 50 请输入请输入 ABAB 杆的长度杆的长度: L3 = 60: L3 = 60 请输入请输入 ACAC 杆和杆和 ABAB 杆的夹角杆的夹角: a21= 45: a21= 45 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 2727 第十第十四四题题 曲柄摇杆式分插机构曲柄摇杆式分插机构 文件名：文件名： fenchafencha 程序：程序： clear;clear; clc;clc; xC=input(xC=input(请输入请输入 CC 点横坐标点横坐标(cm):xC =);(cm):xC =); yC=input(yC=input(请输入请输入 CC 点纵坐标点纵坐标(cm):yC =);(cm):yC =); R=input(R=input(请输入主动杆请输入主动杆 OAOA 长度长度(cm):R =);(cm):R =); m3=input(m3=input(请输入主动杆请输入主动杆 OAOA 的转速的转速( (转转/ /分分):w =);):w =); L2=input(L2=input(请输入从动杆请输入从动杆 ABAB 长度长度(cm):LAB =);(cm):LAB =); L3=input(L3=input(请输入从动杆请输入从动杆 BCBC 长度长度(cm):LBC =);(cm):LBC =); LAE=input(LAE=input(请输入从动杆请输入从动杆 AEAE 长度长度(cm):LAE =);(cm):LAE =); LDE=input(LDE=input(请输入从动杆请输入从动杆 DEDE 长度长度(cm):LDE =);(cm):LDE =); m1=input(m1=input(请输入从动杆请输入从动杆 ABAB 和和 AEAE 的夹角的夹角( (度度):a0 =);):a0 =); m2=input(m2=input(请输入从动杆请输入从动杆 DEDE 和和 AEAE 的夹角的夹角( (度度):a5 =);):a5 =); n=input(n=input(请输入输出图形形式请输入输出图形形式(0:(0:多图多图; 1:; 1:单图单图) )：n =);n =); a0=m1*pi/180;a0=m1*pi/180; a5=m2*pi/180;a5=m2*pi/180; w=m3*pi/30;w=m3*pi/30; for i=1:300for i=1:300 t=linspace(0,2*pi/t=linspace(0,2*pi/w,300);w,300); 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 2828 xA(i)=R*cos(w*t(i);xA(i)=R*cos(w*t(i); yA(i)=R*sin(w*t(i);yA(i)=R*sin(w*t(i); b(i)=atan(yA(i)b(i)=atan(yA(i)- -yC)/(xA(i)yC)/(xA(i)- -xC)+pi;xC)+pi; N=(L32+(xA(i)N=(L32+(xA(i)- -xC)2+(yA(i)xC)2+(yA(i)- -yC)2yC)2- -L22)/(2*L3*sqrt(xA(i)L22)/(2*L3*sqrt(xA(i)- -xC)2+(yA(i)xC)2+(yA(i)- -yC)2);yC)2); if N0if N0 a3(i)=b(i)+atan(sqrt(1a3(i)=b(i)+atan(sqrt(1- -N2)/N);N2)/N); elseelse a3(i)=pi+b(i)+atan(sqrt(1a3(i)=pi+b(i)+atan(sqrt(1- -N2)/N);N2)/N); endend xB(i)=xC+L3*cos(a3(i);xB(i)=xC+L3*cos(a3(i); yB(i)=yC+L3*sin(a3(i);yB(i)=yC+L3*sin(a3(i); u=atan(yB(i)u=atan(yB(i)- -yA(i)/(xB(i)yA(i)/(xB(i)- -xA(i);xA(i); if xB(i)if xB(i)- -xA(i)=0xA(i)=0 if yB(i)if yB(i)- -yA(i)0yA(i)0 a2(i)=pi/2;a2(i)=pi/2; else else a2(i)=3*pi/2;a2(i)=3*pi/2; endend elseelse if xB(i)if xB(i)- -xA(i)0xA(i)0 if yB(i)if yB(i)- -yA(i)0yA(i)0 a2(i)=u;a2(i)=u; else else a2(i)=ua2(i)=u- -2*pi;2*pi; endend else else a2(i)=u+pi;a2(i)=u+pi; endend endend c(i)=atan(yA(i)c(i)=atan(yA(i)- -yB(i)/(xA(i)yB(i)/(xA(i)- -xB(i);xB(i); LAD=sqrt(LDE2+LAE2LAD=sqrt(LDE2+LAE2- -2*LAE*LDE*cos(a0);2*LAE*LDE*cos(a0); a4=acos(LAD2+LAE2a4=acos(LAD2+LAE2- -LDE2)/(2*LAE*LDE);LDE2)/(2*LAE*LDE); a6(i)=a5+a4+c(i);a6(i)=a5+a4+c(i); xD(i)=xA(i)+LAD*cos(a6(i);xD(i)=xA(i)+LAD*cos(a6(i); yD(i)=yA(i)+LAD*sin(a6(i);yD(i)=yA(i)+LAD*sin(a6(i); dxA(i)=dxA(i)=- -R*w*sin(w*t(i);R*w*sin(w*t(i); dyA(i)=R*w*cos(w*t(i);dyA(i)=R*w*cos(w*t(i); da2(i)=(dxA(i)*cos(a3(i)+dyA(i)*sin(a3(i)/(L2*sin(a2(i)da2(i)=(dxA(i)*cos(a3(i)+dyA(i)*sin(a3(i)/(L2*sin(a2(i)- -a3(i);a3(i); da3(i)=(dxA(i)*cos(a2(i)+dyA(i)*sin(a2(i)/(L3*sin(a2(i)da3(i)=(dxA(i)*cos(a2(i)+dyA(i)*sin(a2(i)/(L3*sin(a2(i)- -a3(i);a3(i); 机械机构计算机辅助分析作业（机械机构计算机辅助分析作业（2010 级）级） 刘海洋刘海洋 2929 endend if n=0if n=0 figure(1);figure(1); axis=(axis=(- -400,400,- -150,150,- -400,400);400,400); plot(real(xD),plot(real(xD),real(yD),g);real(yD),g); grid on;grid on; xlabel(Xxlabel(X 值值);); ylabel(Yylabel(Y 值值) ;) ; title(Dtitle(D 点运动轨迹曲线点运动轨迹曲线);); figure(2);figure(2); plot(dxA,dyA,k);plot(dxA,dyA,k); title(Atitle(A 点速度曲线点速度曲线);); xlabel(Xxlabel(X 值值);); ylabel(Yylabel(Y 值值);); grid on;grid on; figure(3);figure(3); plot(w*t,real(da2),b);plot(w*t,real(d