[20040110430]

x=linspace(0,90)
gain=[1/16,1/8,1/4,1/2];
y1=gain(1)-sin(x*gain(1)./180*pi)./sin(x./180*pi); %A端扭矩Ta
y2=gain(2)-sin(x*gain(2)./180*pi)./sin(x./180*pi);
y3=gain(3)-sin(x*gain(3)./180*pi)./sin(x./180*pi);
y4=gain(4)-sin(x*gain(4)./180*pi)./sin(x./180*pi);
plot(x,y1,'-',x,y2,':',x,y3,'o',x,y4,'*');
legend('t=1/16','t=1/8','t=1/4','t=1/2',0)

谢谢楼上的朋友指点

上面是我修改的结果 图像画出来还可以

但是上面的编程好想麻烦点 请问有没有简单的方法赋值给y1,y2,y3,y4呢？呵呵:smile: