请教matlab变化圆心画圆具体函数如下:R0=100;r=27.5;L1=180;L2=105;du=180/pi;wt3=0:1:75;Y0=(acos((L1^2+L2^2-(R0+r)^2)/(2*L1*L2)))*du;Y1=23.9123*(pi*wt3/75-0.5*sin(2*pi*wt3/75))/pi;Y=Y0+Y1;R=sqrt(L1^2+L2^2-2*L1*L2*cos(Y/du));%以上为公
来源:学生作业帮助网 编辑:作业帮 时间:2024/09/11 04:13:06
![请教matlab变化圆心画圆具体函数如下:R0=100;r=27.5;L1=180;L2=105;du=180/pi;wt3=0:1:75;Y0=(acos((L1^2+L2^2-(R0+r)^2)/(2*L1*L2)))*du;Y1=23.9123*(pi*wt3/75-0.5*sin(2*pi*wt3/75))/pi;Y=Y0+Y1;R=sqrt(L1^2+L2^2-2*L1*L2*cos(Y/du));%以上为公](/uploads/image/z/4098310-70-0.jpg?t=%E8%AF%B7%E6%95%99matlab%E5%8F%98%E5%8C%96%E5%9C%86%E5%BF%83%E7%94%BB%E5%9C%86%E5%85%B7%E4%BD%93%E5%87%BD%E6%95%B0%E5%A6%82%E4%B8%8B%EF%BC%9AR0%3D100%3Br%3D27.5%3BL1%3D180%3BL2%3D105%3Bdu%3D180%2Fpi%3Bwt3%3D0%3A1%3A75%3BY0%3D%28acos%28%28L1%5E2%2BL2%5E2-%28R0%2Br%29%5E2%29%2F%282%2AL1%2AL2%29%29%29%2Adu%3BY1%3D23.9123%2A%28pi%2Awt3%2F75-0.5%2Asin%282%2Api%2Awt3%2F75%29%29%2Fpi%3BY%3DY0%2BY1%3BR%3Dsqrt%28L1%5E2%2BL2%5E2-2%2AL1%2AL2%2Acos%28Y%2Fdu%29%29%3B%25%E4%BB%A5%E4%B8%8A%E4%B8%BA%E5%85%AC)
请教matlab变化圆心画圆具体函数如下:R0=100;r=27.5;L1=180;L2=105;du=180/pi;wt3=0:1:75;Y0=(acos((L1^2+L2^2-(R0+r)^2)/(2*L1*L2)))*du;Y1=23.9123*(pi*wt3/75-0.5*sin(2*pi*wt3/75))/pi;Y=Y0+Y1;R=sqrt(L1^2+L2^2-2*L1*L2*cos(Y/du));%以上为公
请教matlab变化圆心画圆
具体函数如下:
R0=100;
r=27.5;
L1=180;
L2=105;
du=180/pi;
wt3=0:1:75;
Y0=(acos((L1^2+L2^2-(R0+r)^2)/(2*L1*L2)))*du;
Y1=23.9123*(pi*wt3/75-0.5*sin(2*pi*wt3/75))/pi;
Y=Y0+Y1;
R=sqrt(L1^2+L2^2-2*L1*L2*cos(Y/du));
%以上为公式
x=zeros(1,360);
y=zeros(1,360);
m=zeros(1,360);
n=zeros(1,360);
%以上为缩小范围
for i=0:1:75 %画曲线
x(i+1)=R(i+1)*cos(i*pi/180.0);
y(i+1)=R(i+1)*sin(i*pi/180.0);
for t=0:1:2*pi %画圆
m=x(i+1)+r*cos(t);
n=y(i+1)+r*sin(t);
plot(m,n,'r')
end
end
plot(x,y,'r')
axis equal
hold on
程序运行后是个曲线,我想显示出曲线上每隔i=1,以(x,y)为圆心,r为半径画圆,这样会有76个圆,为什么显示不出来呢,我哪里出错了?
请教matlab变化圆心画圆具体函数如下:R0=100;r=27.5;L1=180;L2=105;du=180/pi;wt3=0:1:75;Y0=(acos((L1^2+L2^2-(R0+r)^2)/(2*L1*L2)))*du;Y1=23.9123*(pi*wt3/75-0.5*sin(2*pi*wt3/75))/pi;Y=Y0+Y1;R=sqrt(L1^2+L2^2-2*L1*L2*cos(Y/du));%以上为公
第一种和第二种是动态的,第三种是静态的
1、3D版
clear all;
clc;
R0=100;
r=27.5;
L1=180;
L2=105;
du=180/pi;
wt3=0:1:75;
Y0=(acos((L1^2+L2^2-(R0+r)^2)/(2*L1*L2)))*du;
Y1=23.9123*(pi*wt3/75-0.5*sin(2*pi*wt3/75))/pi;
Y=Y0+Y1;
R=sqrt(L1^2+L2^2-2*L1*L2*cos(Y/du));
%以上为公式
%hold on;% 加上hold on使得所有图形都在一张图上,按照matlab顺序执行的特点,hold on不能放在多重绘图语句的最后面,一般放在第一个绘图语句的前面或者第一第二个语句的中间都行,如:a注释所示
%axis equal% 是一种对当前图形句柄的修正命令,放在任意位置都行,注意只对当前图形句柄有效,对于动态图,其用法跟hold on类似
x=zeros(1,76);% 删除m,n的预开辟内存空间的语句,因为他们的长度是由0:0.1:2*pi的长度决定
y=zeros(1,76);% x,y的长度由0:75决定,所以为76
%以上为缩小范围
for i=0:1:75 %画曲线
x(i+1)=R(i+1)*cos(i*pi/180.0);
y(i+1)=R(i+1)*sin(i*pi/180.0);
k=1;
for t=0:0.1:2*pi %画圆
m(k)=x(i+1)+r*cos(t);
n(k)=y(i+1)+r*sin(t);
k=k+1;
end
z=repmat(i,1,k-1);
plot3(m,n,z,'r');
axis([-20 220 -40 200 0 78]);
hold on
plot3(x(1:i+1),y(1:i+1),0:i,'b');
pause(0.1);
% 将绘图语句放在循环外
% hold on % a 如:把hold on语句放在此处也行
end
2、2D版
clear all;
clc;
R0=100;
r=27.5;
L1=180;
L2=105;
du=180/pi;
wt3=0:1:75;
Y0=(acos((L1^2+L2^2-(R0+r)^2)/(2*L1*L2)))*du;
Y1=23.9123*(pi*wt3/75-0.5*sin(2*pi*wt3/75))/pi;
Y=Y0+Y1;
R=sqrt(L1^2+L2^2-2*L1*L2*cos(Y/du));
%以上为公式
hold on;% 加上hold on使得所有图形都在一张图上,按照matlab顺序执行的特点,hold on不能放在多重绘图语句的最后面,一般放在第一个绘图语句的前面或者第一第二个语句的中间都行,如:a注释所示
axis equal% 是一种对当前图形句柄的修正命令,放在任意位置都行,注意只对当前图形句柄有效,对于动态图,其用法跟hold on类似
x=zeros(1,76);% 删除m,n的预开辟内存空间的语句,因为他们的长度是由0:0.1:2*pi的长度决定
y=zeros(1,76);% x,y的长度由0:75决定,所以为76
%以上为缩小范围
for i=0:1:75 %画曲线
x(i+1)=R(i+1)*cos(i*pi/180.0);
y(i+1)=R(i+1)*sin(i*pi/180.0);
k=1;
for t=0:0.1:2*pi %画圆
m(k)=x(i+1)+r*cos(t);
n(k)=y(i+1)+r*sin(t);
k=k+1;
end
plot(m,n,'r');
axis([-20 220 -40 200]);
hold on
plot(x(1:i+1),y(1:i+1),'b');
pause(0.1);
% 将绘图语句放在循环外
% hold on % a 如:把hold on语句放在此处也行
end
3、静态2D版
clear all;
clc;
R0=100;
r=27.5;
L1=180;
L2=105;
du=180/pi;
wt3=0:1:75;
Y0=(acos((L1^2+L2^2-(R0+r)^2)/(2*L1*L2)))*du;
Y1=23.9123*(pi*wt3/75-0.5*sin(2*pi*wt3/75))/pi;
Y=Y0+Y1;
R=sqrt(L1^2+L2^2-2*L1*L2*cos(Y/du));
%以上为公式
hold on;% 加上hold on使得所有图形都在一张图上,按照matlab顺序执行的特点,hold on不能放在多重绘图语句的最后面,一般放在第一个绘图语句的前面或者第一第二个语句的中间都行,如:a注释所示
axis equal% 是一种对当前图形句柄的修正命令,放在任意位置都行,注意只对当前图形句柄有效,对于动态图,其用法跟hold on类似
x=zeros(1,76);% 删除m,n的预开辟内存空间的语句,因为他们的长度是由0:0.1:2*pi的长度决定
y=zeros(1,76);% x,y的长度由0:75决定,所以为76
%以上为缩小范围
for i=0:1:75 %画曲线
x(i+1)=R(i+1)*cos(i*pi/180.0);
y(i+1)=R(i+1)*sin(i*pi/180.0);
k=1;
for t=0:0.1:2*pi %画圆
m(k)=x(i+1)+r*cos(t);
n(k)=y(i+1)+r*sin(t);
k=k+1;
end
plot(m,n,'r');
axis([-20 220 -40 200]);
hold on
plot(x(1:i+1),y(1:i+1),'b');
% 将绘图语句放在循环外
% hold on % a 如:把hold on语句放在此处也行
end
%试一下
R0=100;
r=27.5;
L1=180;
L2=105;
du=180/pi;
wt3=0:1:75;
Y0=(acos((L1^2+L2^2-(R0+r)^2)/(2*L1*L2)))*du;
Y1=23.9123*(pi*wt3/75-0.5*sin(2*pi*wt3/75))/pi;
Y=Y0+Y1;
全部展开
%试一下
R0=100;
r=27.5;
L1=180;
L2=105;
du=180/pi;
wt3=0:1:75;
Y0=(acos((L1^2+L2^2-(R0+r)^2)/(2*L1*L2)))*du;
Y1=23.9123*(pi*wt3/75-0.5*sin(2*pi*wt3/75))/pi;
Y=Y0+Y1;
R=sqrt(L1^2+L2^2-2*L1*L2*cos(Y/du));
%以上为公式
x=zeros(1,360);
y=zeros(1,360);
m=zeros(1,360);
n=zeros(1,360);
%以上为缩小范围
for i=0:1:75 %画曲线
x(i+1)=R(i+1)*cos(i*pi/180.0);
y(i+1)=R(i+1)*sin(i*pi/180.0);
%%%%%%%%%%%
%改变t间隔
for t=0:0.1:2*pi %画圆
m=x(i+1)+r*cos(t);
n=y(i+1)+r*sin(t);
%%%%%%%%%%%
%改变画图指令及hold on
plot(m,n,'-r')
hold on
end
end
%%%%%%%%%%%
%最后的零点不出现
plot(x(1:76),y(1:76),'r')
axis equal
收起
for i=0:1:75 %画曲线
x(i+1)=R(i+1)*cos(i*pi/180.0);
y(i+1)=R(i+1)*sin(i*pi/180.0);
for t=0:1:2*pi %画圆
m=x(i+1)+r*cos(t);
n=y(i+1)+r*sin(t);
plot(m,n,'r'...
全部展开
for i=0:1:75 %画曲线
x(i+1)=R(i+1)*cos(i*pi/180.0);
y(i+1)=R(i+1)*sin(i*pi/180.0);
for t=0:1:2*pi %画圆
m=x(i+1)+r*cos(t);
n=y(i+1)+r*sin(t);
plot(m,n,'r')
hold on
end
end
plot(x,y,'r')
axis equal
hold on
不知道是不是你hold on的问题
收起