论坛 › K01：MatLab › 请高手指教一下本人编的一个四杆机构仿真程序

darkbluexi 发表于 2007-6-5 15:31:15

请高手指教一下本人编的一个四杆机构仿真程序

程序如下：
figure(1)
axis([-100,200,-300,100]);
hold on
xd=89.21;

yd=-11.91;
theta0=0;
theta1=pi/6;
theta5=5*pi/36;
theta0=0; l1=34; xb=l1*cos(theta0);
yb=l1*sin(theta0);
l2=66;l4=80; l0=((l1).^2+(l4).^2-2*l1*l4*cos(theta1+theta0)).^(1/2);
l3=90;theta3=acos(((l0).^2+(l3).^2-(l2).^2)/(2*(l0)*(l3)));
theta2=acos(((l4).^2+(l0).^2-(l1).^2)/(2*(l4)*(l0)));
theta3=acos(((l0).^2+(l2).^2-(l3).^2)/(2*(l0)*(l2)));
xc=89.21+(l3)*cos(31*pi/36+theta2+theta3);
yc=-11.91+(l3)*sin(31*pi/36+theta2+theta3);
l5=154;
xm=l1*cos(theta0)+l5*cos(theta3-theta2+theta1-theta5);
ym=l1*sin(theta0)+l5*sin(theta3-theta2+theta1-theta5);
jijia=line(,[-11.91,0],′color′,′r′,′linewidth′,3);
handle_a=line(0,0,′color′,′b′,′marker′,′.′,...
′erasemode′,′xor′,′markersize′,20);
handle_b=line(xb,yb,′color′,′b′,′marker′,′.′,...
′erasemode′,′xor′,′markersize′,20)
gan1=line(,,′color′,′r′,′linewidth′,3.5,...
   ′erasemode′,′xor′);
handle_c = line (xc, yc,′color′,′k′,′marker′,′.′,...
   ′erasemode′,′xor′,′markersize′,5);
gan2=line(,,′color′,′r′,′linewidth′,3.5,...
′erasemode′,′xor′);
handle_d=line(89.21,-11.91,′color′,′b′, ′marker′,′.′,...
′erasemode′,′xor′,′markersize′,20);
gan3=line(,,′color′,′r′,′linewidth′,3.5,...
   ′erasemode′,′xor′);
handle_m = line (xm, ym,′color′,′k′,′marker′,′.′,...
   ′erasemode′,′xor′,′markersize′,5);
gan4=line(,,′color′,′r′,′linewidth′,3.5,...
′erasemode′,′xor′);

运行后出现：
??? jijia=line(,[-11.91,0],′color′,′r′,′linewidth′,3);
                                    |
Error: Missing variable or function.

怎么调都不行，很郁闷啊，哪位高手指点一下

shunfly 发表于 2007-6-5 15:52:58

jijia=line(,[-11.91,0],'color','r','linewidth',3);
不知道你是在哪里抄得程序，依此类推

yizi0000 发表于 2007-6-5 16:22:39

:) 楼主,好像你的输入法不对吧.单引号哪有那样的啊?

alexqxp 发表于 2007-6-5 19:11:18

估计是作本科毕业设计。四杆仿真的程序机械原理里面有详细的公式说明。先把杆组的点坐标，速度，加速度公式看明白了，再找个VB的程序看看，最后参考VB程序转成MATLAB的。另外，MATLAB里面有现成的SimMechanics例子，但是恐怕你看不了。

lishuai8 发表于 2007-6-14 09:10:19

为什么不用 matlab 里的simechnic呢？demo里就有四连杆机构。非常好学，建议一用。

qhm2lhy 发表于 2007-6-22 16:27:36

呵呵对阿
   simechnic   我用它做过 一个RRR-RRP的6杆机构的运动仿真

lin2009 发表于 2010-10-28 16:32:40

可以运行，但结果不对。
看了plot中的参数   ...'erasemode', 'xor');应该是做动画的。
但是动不起来，可能只Copy的部分的程序。
figure(1)
axis([ -100, 200, -300, 100 ]);
hold on
xd = 89.21;

yd = -11.91;
theta0 = 0;
theta1 = pi/6;
theta5 = 5*pi/36;

l1 = 34;
xb = l1*cos(theta0);
yb = l1*sin(theta0);
l2 = 66;
l4 = 80;
l0 = ((l1).^2 + (l4).^2 - 2*l1*l4*cos(theta1+theta0)).^(1/2);
l3 = 90;
theta2 = acos(((l4).^2 + (l0).^2 - (l1).^2)/(2*(l4) * (l0)));
theta3 = acos(((l0).^2 + (l2).^2 - (l3).^2)/(2*(l0) * (l2)));
xc = 89.21 + (l3)*cos(31*pi/36 + theta2 + theta3);
yc = -11.91 + (l3)*sin(31*pi/36 + theta2 + theta3);
l5 = 154;
xm = l1*cos(theta0) + l5*cos(theta3 - theta2 + theta1 - theta5);
ym = l1*sin(theta0) + l5*sin(theta3 - theta2 + theta1 - theta5);
jijia = line([ 0, 89.21 ], [ -11.91, 0 ], 'color', 'r', 'linewidth', 3);
handle_a = line(0, 0, 'color', 'b', 'marker', '.', 'erasemode', 'xor', 'markersize', 20);
handle_b = line(xb, yb, 'color', 'b', 'marker', '.', 'erasemode', 'xor', 'markersize', 20);
gan1 = line([ 0, xb ], [ 0, yb ], 'color', 'r', 'linewidth', 3.5, 'erasemode', 'xor');
handle_c = line (xc, yc, 'color', 'k', 'marker', '.', 'erasemode', 'xor', 'markersize', 5);
gan2 = line([ xb, xc ], [ yb, yc ], 'color', 'r', 'linewidth', 3.5, 'erasemode', 'xor');
handle_d = line(89.21, -11.91, 'color', 'b', 'marker', '.', 'erasemode', 'xor', 'markersize', 20);
gan3 = line([ xc, 89.21 ], [ yc, -11.91 ], 'color', 'r', 'linewidth', 3.5, 'erasemode', 'xor');
handle_m = line (xm, ym, 'color', 'k', 'marker', '.', 'erasemode', 'xor', 'markersize', 5);
gan4 = line([ xc, xm ], [ yc, ym ], 'color', 'r', 'linewidth', 3.5, 'erasemode', 'xor');

lin2009 发表于 2010-10-28 16:34:50

可以运行，但结果不对。
看了plot中的参数   ...'erasemode', 'xor');应该是做动画的。
但是动不起来，可能只Copy的部分的程序。
figure(1)
axis([ -100, 200, -300, 100 ]);
hold on
xd = 89.21;

yd = -11.91;
theta0 = 0;
theta1 = pi/6;
theta5 = 5*pi/36;

l1 = 34;
xb = l1*cos(theta0);
yb = l1*sin(theta0);
l2 = 66;
l4 = 80;
l0 = ((l1).^2 + (l4).^2 - 2*l1*l4*cos(theta1+theta0)).^(1/2);
l3 = 90;
theta2 = acos(((l4).^2 + (l0).^2 - (l1).^2)/(2*(l4) * (l0)));
theta3 = acos(((l0).^2 + (l2).^2 - (l3).^2)/(2*(l0) * (l2)));
xc = 89.21 + (l3)*cos(31*pi/36 + theta2 + theta3);
yc = -11.91 + (l3)*sin(31*pi/36 + theta2 + theta3);
l5 = 154;
xm = l1*cos(theta0) + l5*cos(theta3 - theta2 + theta1 - theta5);
ym = l1*sin(theta0) + l5*sin(theta3 - theta2 + theta1 - theta5);
jijia = line([ 0, 89.21 ], [ -11.91, 0 ], 'color', 'r', 'linewidth', 3);
handle_a = line(0, 0, 'color', 'b', 'marker', '.', 'erasemode', 'xor', 'markersize', 20);
handle_b = line(xb, yb, 'color', 'b', 'marker', '.', 'erasemode', 'xor', 'markersize', 20);
gan1 = line([ 0, xb ], [ 0, yb ], 'color', 'r', 'linewidth', 3.5, 'erasemode', 'xor');
handle_c = line (xc, yc, 'color', 'k', 'marker', '.', 'erasemode', 'xor', 'markersize', 5);
gan2 = line([ xb, xc ], [ yb, yc ], 'color', 'r', 'linewidth', 3.5, 'erasemode', 'xor');
handle_d = line(89.21, -11.91, 'color', 'b', 'marker', '.', 'erasemode', 'xor', 'markersize', 20);
gan3 = line([ xc, 89.21 ], [ yc, -11.91 ], 'color', 'r', 'linewidth', 3.5, 'erasemode', 'xor');
handle_m = line (xm, ym, 'color', 'k', 'marker', '.', 'erasemode', 'xor', 'markersize', 5);
gan4 = line([ xc, xm ], [ yc, ym ], 'color', 'r', 'linewidth', 3.5, 'erasemode', 'xor');

查看完整版本: 请高手指教一下本人编的一个四杆机构仿真程序