n=900;m=700;a=10*pi/180;holdonyi=0;fork=1:7x=(1:m)+(k-1).*(4*n-2*m).*ones(1,m);y=tan(a)*x+yin=900;m=700;a=10*pi/180;holdonyi=0;fork=1:7x=(1:m)+(k-1).*(4*n-2*m).*ones(1,m);y=tan(a)*x+yi;plot(x,y)x=(m+1:2*n-m)+(k-1).*(4*n-2*m).*ones(1,2*n-m-m)

Question

n=900;m=700;a=10*pi/180; hold on yi=0; for k=1:7 x=(1:m)+(k-1).*(4*n-2*m).*ones(1,m); y=tan(a)*x+yi
n=900;m=700;a=10*pi/180;
hold on
yi=0;
for k=1:7
x=(1:m)+(k-1).*(4*n-2*m).*ones(1,m);
y=tan(a)*x+yi;
plot(x,y)
x=(m+1:2*n-m)+(k-1).*(4*n-2*m).*ones(1,2*n-m-m);
y=sqrt(((n-m)/sin(a)).^2-(x-n).^2)+n*tan(a)-(n-m)/sin(a)*sqrt((tan(a))^2+1)+yi;
plot(x,y)
x=(2*n-m+1:2*n)+(k-1).*(4*n-2*m).*ones(1,m);
y=-tan(a).*(x-2*n)+yi;
plot(x,y)
x=(2*n+1:4*n-2*m)+(k-1).*(4*n-2*m).*ones(1,2*n-2*m);
y=-sqrt(((n-m)/sin(a)).^2-(x-3*n+m).^2)-(n-m)*tan(a)+(n-m)/sin(a)*sqrt((tan(a))^2+1)+yi;
plot(x,y)
yi=y(size(y));
end
Error using ==> plus
Matrix dimensions must agree.
怎么会这样?求教!本人新手
是画一轨迹曲线

Answer 1

额,觉得代码比较奇怪,建议不要将数组和数做加减,这样不好,比如如下几句都是这样
y=tan(a)*x+yi;（yi与x）
x=(m+1:2*n-m)+(k-1).*(4*n-2*m).*ones(1,2*n-m-m);
y=sqrt(((n-m)/sin(a)).^2-(x-n).^2)+n*tan(a)-(n-m)/sin(a)*sqrt((tan(a))^2+1)+yi;
x=(2*n-m+1:2*n)+(k-1).*(4*n-2*m).*ones(1,m);
y=-tan(a).*(x-2*n)+yi;
x=(2*n+1:4*n-2*m)+(k-1).*(4*n-2*m).*ones(1,2*n-2*m);
y=-sqrt(((n-m)/sin(a)).^2-(x-3*n+m).^2)-(n-m)*tan(a)+(n-m)/sin(a)*sqrt((tan(a))^2+1)+yi;
这也就是为什么报错说矩阵加法过程中矩阵维数不一致