还是那个for循环的问题,我想在图里标坐标表示什么应该在什么地方加y1=60;y4=200;y2=120;y3=150;x1=linspace(0,2*pi,50);A=2*y1*y2*sin(x1);B=2*y2*(y1*cos(x1)-y4);C=y1^2+y2^2+y4^2-2*y1*y4*cos(x1)-y3^2;D=2*y1*y3*sin(x1);E=2*y3*(y1*c

Question

还是那个for循环的问题,我想在图里标坐标表示什么应该在什么地方加
y1=60;y4=200;y2=120;y3=150;
x1=linspace(0,2*pi,50);
A=2*y1*y2*sin(x1);
B=2*y2*(y1*cos(x1)-y4);
C=y1^2+y2^2+y4^2-2*y1*y4*cos(x1)-y3^2;
D=2*y1*y3*sin(x1);
E=2*y3*(y1*cos(x1)-y4);
F=y2^2-y1^2-y3^2-y4^2+2*y1*y4*cos(x1);
for i=1:50
x(i)=2*atan((A(i)+sqrt(A(i)^2+B(i)^2-C(i)^2))/(B(i)-C(i)));
y(i)=2*atan((D(i)-sqrt(D(i)^2+E(i)^2-F(i)^2))/(B(i)-C(i)))
subplot(1,2,1)
plot([x1(i),x(i)])
subplot(1,2,2)
plot([x1(i),y(i)])
endxlabel('凸轮转角')
ylabel('杆2与竖直方向的夹角') 我这样加,图像里什么都没有,

Answer 1

你函数的用法没有错,应该没有问题,我运行了一下是显示的,你再运行一下我给你的这个,如果还不行再联系；y1=60;y4=200;y2=120;y3=150;
x1=linspace(0,2*pi,50);
A=2*y1*y2*sin(x1);
B=2*y2*(y1*cos(x1)-y4);
C=y1^2+y2^2+y4^2-2*y1*y4*cos(x1)-y3^2;
D=2*y1*y3*sin(x1);
E=2*y3*(y1*cos(x1)-y4);
F=y2^2-y1^2-y3^2-y4^2+2*y1*y4*cos(x1);
for i=1:50
x(i)=2*atan((A(i)+sqrt(A(i)^2+B(i)^2-C(i)^2))/(B(i)-C(i)));
y(i)=2*atan((D(i)-sqrt(D(i)^2+E(i)^2-F(i)^2))/(B(i)-C(i)))
subplot(1,2,1)
plot([x1(i),x(i)])
xlabel('凸轮转角')
ylabel('杆2与竖直方向的夹角')
subplot(1,2,2)
plot([x1(i),y(i)])
xlabel('凸轮转角')
ylabel('杆2与竖直方向的夹角')
end