[未答复] 这个错误是什么呀,像求解f点关于sita的函数,

[复制链接]
IRHH 发表于 2022-5-14 16:51:56
sym x_f
sym y_f
l_ef=16;
x_e=-30;
y_e=-36;
l_bc=43;
l_cd=20;
y_d=0;x_d=0;

F=[x_f,y_f];
E=[x_e,y_e];
eny1=(y_f-y_e)*(y_f-y_c)/(x_f-x_e)/(x_f-x_c)+1;
eny2=(F-E)'*(F-E)-l_ef^2;
F_0=[-10,-15];
for sita=1:0.01:2*pi
    x_c=x_d+l_cd*cos(sita);
    y_c=y_d+l_cd*sin(sita);
    ENY=[eny1;eny2];
    n=1;
    while norm(ENY)>=0.01
        J(1,1)=diff(eny1,x_f);J(1,2)=diff(eny1,y_f);diff(eny2,x_f);diff(eny2,y_f);
        F=F_0-inv(J)*ENY;
        ENY=[eny1,eny2];
        n=n+1;
        F_0=F;
        if n>=1000
            n;
            break
        end
    end
end

ans =

x_f


ans =

y_f

错误使用 sym/cat>checkDimensions (第 68 行)
CAT arguments dimensions not consistent.

出错 sym/cat>catMany (第 33 行)
[resz, ranges] = checkDimensions(sz,dim);

出错 sym/cat (第 25 行)
    ySym = catMany(dim, args);

出错 sym/vertcat (第 19 行)
    ySym = cat(1,args{:});

>>

您需要登录后才可以回帖 登录 | 注册

本版积分规则

相关帖子
相关文章
热门教程
站长推荐
快速回复 返回顶部 返回列表