关于#matlab#的问题：出错 Untitled4 (第 21 行)polar(ct,rou' , 'b-'),hold on

救命啊，大 哥们，看看到底哪里出错了

右标感觉没有用啊

ct=s./s(10000)*2*pi;
rou=e+E.*sin(atan(L./(R+m)));
rou1=e+E.*sin(atan(L./(R+m)))-2;
polar(ct,rou' , 'b-'),hold on;

for i=1:10000
    x(i)=rou(i)*cos(ct(i));
    y(i)=rou(i)*sin(ct(i));
    z(i)=0;
end
tulun=[x;y;z];
writematrix(tulun1','tulun1.txt');

THETA 和 RHO 的大小必须相同。

出错 Untitled4 (第 21 行)
polar(ct,rou' , 'b-'),hold on;
照着别人写的，为什么报错啊

我想要达到的结果，如果你需要快速回答，请尝试 “付费悬赏”

polar(ct,rou' , 'b-'); 这里是分号
hold on;

完整代码如下：

clc;
clear;
M=80;
m=10;
E=20;
L=120;
e=50;
XRY=textread('guiji2.0.txt');
X=XRY(:,2);
Y=XRY(:,3);
R=XRY(:,5);
plot(X,Y,'b-'),hold on,axis equal;
s(1)=0;
for  i=2:10000
    s(i)=s(i-1)+sqrt((X(i)-X(i-1))^2+(Y(i)-Y(i-1))^2)*(R(i)-M/2)/R(i);
end
ct=s./s(10000)*2*pi;
rou=e+E.*sin(atan(L./(R+m)));
rou1=e+E.*sin(atan(L./(R+m)))-2;
polar(ct,rou' ,'b-');
hold on;
for i=1:10000
    x(i)=rou(i)*cos(ct(i));
    y(i)=rou(i)*sin(ct(i));
    z(i)=0;
end
tulun=[x;y;z];
writematrix(tulun','tulun1.txt');