作业帮关于MATLAB,哪位大神帮看看这个函数哪里错了function a=fit_ellipse(x,y)D=[x.*x x.*y y.*y x y ones(size(x))];S=D'*D;C(6,6)=0;C(1,3)=2;C(2,2)=-1;C(3,1)=2;[gevec,geval]=eig(inv(S)*C);[PosR,PosC]=find(geval>0 & ~isinf(geval));a=gevec(:,PosC
来源：学生作业帮助网 编辑：作业帮 时间：2024/05/04 18:44:50
![关于MATLAB,哪位大神帮看看这个函数哪里错了function a=fit_ellipse(x,y)D=[x.*x x.*y y.*y x y ones(size(x))];S=D'*D;C(6,6)=0;C(1,3)=2;C(2,2)=-1;C(3,1)=2;[gevec,geval]=eig(inv(S)*C);[PosR,PosC]=find(geval>0 & ~isinf(geval));a=gevec(:,PosC](/uploads/image/z/2693910-30-0.jpg?t=%E5%85%B3%E4%BA%8EMATLAB%2C%E5%93%AA%E4%BD%8D%E5%A4%A7%E7%A5%9E%E5%B8%AE%E7%9C%8B%E7%9C%8B%E8%BF%99%E4%B8%AA%E5%87%BD%E6%95%B0%E5%93%AA%E9%87%8C%E9%94%99%E4%BA%86function+a%3Dfit_ellipse%28x%2Cy%29D%3D%5Bx.%2Ax+x.%2Ay+y.%2Ay+x+y+ones%28size%28x%29%29%5D%3BS%3DD%27%2AD%3BC%286%2C6%29%3D0%3BC%281%2C3%29%3D2%3BC%282%2C2%29%3D-1%3BC%283%2C1%29%3D2%3B%5Bgevec%2Cgeval%5D%3Deig%28inv%28S%29%2AC%29%3B%5BPosR%2CPosC%5D%3Dfind%28geval%3E0+%26+%7Eisinf%28geval%29%29%3Ba%3Dgevec%28%3A%2CPosC)
关于MATLAB,哪位大神帮看看这个函数哪里错了function a=fit_ellipse(x,y)D=[x.*x x.*y y.*y x y ones(size(x))];S=D'*D;C(6,6)=0;C(1,3)=2;C(2,2)=-1;C(3,1)=2;[gevec,geval]=eig(inv(S)*C);[PosR,PosC]=find(geval>0 & ~isinf(geval));a=gevec(:,PosC
关于MATLAB,哪位大神帮看看这个函数哪里错了
function a=fit_ellipse(x,y)
D=[x.*x x.*y y.*y x y ones(size(x))];
S=D'*D;
C(6,6)=0;C(1,3)=2;C(2,2)=-1;C(3,1)=2;
[gevec,geval]=eig(inv(S)*C);
[PosR,PosC]=find(geval>0 & ~isinf(geval));
a=gevec(:,PosC);
该函数是要根据已知点集求最佳椭圆一般方程的系数
每次输入点集后调用函数都说第5行[gevec,geval]=eig(inv(S)*C);有错误
到底是怎么回事呢?
关于MATLAB,哪位大神帮看看这个函数哪里错了function a=fit_ellipse(x,y)D=[x.*x x.*y y.*y x y ones(size(x))];S=D'*D;C(6,6)=0;C(1,3)=2;C(2,2)=-1;C(3,1)=2;[gevec,geval]=eig(inv(S)*C);[PosR,PosC]=find(geval>0 & ~isinf(geval));a=gevec(:,PosC
如果你的输入是行向量,那么你的D的维数是6乘以你的向量个数,比如你的行向量个数是9的话,那D就是54×54 的矩阵,而C是6×6的矩阵维数不一样,矩阵相乘出错,如果你的输入是列向量,那么D的维数是6×6的矩阵,维数匹配,不会有问题