主题:[讨论]help me !!!!
global cta1 L1 L2 L3 L4
omega1=1;L1=20;L2=100;L3=60;L4=120;
cta=[];omega=[];alpha=[];xo=[0.6423,1.2143];
t=0:0.1:2*pi;
for cta1=omega1*t
xo=fsolve('position',xo,optimset('Display','off'));
s1=sin(cta1);s2=sin(xo(1));s3=sin(xo(2));
c1=cos(cta1);c2=cos(xo(1));c3=cos(xo(2));
A=[-L2*s2,L3*s3,0,0
L2*c2 -L3*c3 0 0];
A1=[-L2*c2 L3*c3 0 0
-L2*c2 L3*s3 0 0];
b=[s1;-c1];
b1=[c1;s1];
w=A\(omega1*L1*b);
e=A\(-A1*w.^2+omega1^2*L1*b1);
cta=[cta;xo];
omega=[omega;w'];
alpha=[alpha;e'];
end
subplot(3,1,1);plot(omega1*t*180/pi,cta);axis([0 360 0 3]);grid on;
title('\fontname{宋体}\fontsize{8}各杆角位移曲线');
legend('\fontsize{6}\bf2杆角位移','\fontsize{6}\bf3杆角位移');
xlabel('\fontname{Arial}\fontsize{8}主动件角位移\sl{\theta}_{2}');
ylabel('\fontname{Arial}\fontsize{8}从动件角位移\sl{\theta}');
subplot(3,1,2);plot(omega1*t*180/pi,omega);axis([0 360 -1 1]);grid on;
title('\fontname{宋体}\fontsize{8}各杆角速度曲线');
legend('\fontsize{6}\bf2杆角速度','\fontsize{6}\bf3杆角速度');
xlabel('\fontname{Arial}\fontsize{8}主动件角位移\sl{\theta}_{2}');
ylabel('\fontname{Arial}\fontsize{8}从动件角速度\sl{\omega}');
subplot(3,1,3);plot(omega1*t*180/pi,alpha);axis([0 360 -2 2]);grid on;
title('\fontname{宋体}\fontsize{8}各杆角加速度曲线');
legend('\fontsize{6}\bf2杆角加速度','\fontsize{6}\bf3杆角加速度');
xlabel('\fontname{Arial}\fontsize{8}主动件角位移\sl{\theta}_{2}');
ylabel('\fontname{Arial}\fontsize{8}从动件角加速度\sl{\epsilon}');
function F=position(x)
global cta1 L1 L2 L3 L4
F=[L1*cos(cta1)+L2*cos(x(1))-L3*cos(x(2))-L4;...
L1*sin(cta1)+L2*sin(x(1))-L3*sin(x(2))];
为什么会出现下面的错误提示?谢谢帮忙解答!!!
??? Error using ==> feval
Undefined command/function 'position'.
Error in ==> fsolve at 180
fuser = feval(funfcn{3},x,varargin{:});
请好心的朋友帮忙解答一下
谢谢了[em18][em18][em18][em18][em18][em18][em18]
omega1=1;L1=20;L2=100;L3=60;L4=120;
cta=[];omega=[];alpha=[];xo=[0.6423,1.2143];
t=0:0.1:2*pi;
for cta1=omega1*t
xo=fsolve('position',xo,optimset('Display','off'));
s1=sin(cta1);s2=sin(xo(1));s3=sin(xo(2));
c1=cos(cta1);c2=cos(xo(1));c3=cos(xo(2));
A=[-L2*s2,L3*s3,0,0
L2*c2 -L3*c3 0 0];
A1=[-L2*c2 L3*c3 0 0
-L2*c2 L3*s3 0 0];
b=[s1;-c1];
b1=[c1;s1];
w=A\(omega1*L1*b);
e=A\(-A1*w.^2+omega1^2*L1*b1);
cta=[cta;xo];
omega=[omega;w'];
alpha=[alpha;e'];
end
subplot(3,1,1);plot(omega1*t*180/pi,cta);axis([0 360 0 3]);grid on;
title('\fontname{宋体}\fontsize{8}各杆角位移曲线');
legend('\fontsize{6}\bf2杆角位移','\fontsize{6}\bf3杆角位移');
xlabel('\fontname{Arial}\fontsize{8}主动件角位移\sl{\theta}_{2}');
ylabel('\fontname{Arial}\fontsize{8}从动件角位移\sl{\theta}');
subplot(3,1,2);plot(omega1*t*180/pi,omega);axis([0 360 -1 1]);grid on;
title('\fontname{宋体}\fontsize{8}各杆角速度曲线');
legend('\fontsize{6}\bf2杆角速度','\fontsize{6}\bf3杆角速度');
xlabel('\fontname{Arial}\fontsize{8}主动件角位移\sl{\theta}_{2}');
ylabel('\fontname{Arial}\fontsize{8}从动件角速度\sl{\omega}');
subplot(3,1,3);plot(omega1*t*180/pi,alpha);axis([0 360 -2 2]);grid on;
title('\fontname{宋体}\fontsize{8}各杆角加速度曲线');
legend('\fontsize{6}\bf2杆角加速度','\fontsize{6}\bf3杆角加速度');
xlabel('\fontname{Arial}\fontsize{8}主动件角位移\sl{\theta}_{2}');
ylabel('\fontname{Arial}\fontsize{8}从动件角加速度\sl{\epsilon}');
function F=position(x)
global cta1 L1 L2 L3 L4
F=[L1*cos(cta1)+L2*cos(x(1))-L3*cos(x(2))-L4;...
L1*sin(cta1)+L2*sin(x(1))-L3*sin(x(2))];
为什么会出现下面的错误提示?谢谢帮忙解答!!!
??? Error using ==> feval
Undefined command/function 'position'.
Error in ==> fsolve at 180
fuser = feval(funfcn{3},x,varargin{:});
请好心的朋友帮忙解答一下
谢谢了[em18][em18][em18][em18][em18][em18][em18]