作业帮 > 综合 > 作业

matlab小段代码求解释

来源:学生作业帮 编辑:搜狗做题网作业帮 分类:综合作业 时间:2024/04/29 13:25:51
matlab小段代码求解释
function [R,mx,my,mse]=circle_detect(x,y)
N=length(x);
a1=2*(sum(x)^2-N*sum(x.*x));
a2=2*(sum(x)*sum(y)-N*sum(x.*y));
b1=a2;
b2=2*(sum(y)^2-N*sum(y.*y));
c1=sum(x.*x)*sum(x)-N*sum(x.^3)+sum(x)*sum(y.*y)-N*sum(x.*y.*y);
c2=sum(x.*x)*sum(y)-N*sum(y.^3)+sum(y)*sum(y.*y)-N*sum(x.*x.*y);
mx=(c1*b2-c2*b1)/(a1*b2-a2*b1);
my=(a1*c2-a2*c1)/(a1*b2-a2*b1);
R=sqrt((sum(x.*x)-2*sum(x)*mx+N*mx*mx+sum(y.*y)-2*sum(y)*my+N*my*my)/N);
e=((x-mx).^2+(y-my).^2-R*R)/(R*R);
mse=mean(e.*e);
这是一个检测圆圈的代码,这里面的mx,my,R,
matlab小段代码求解释
实话说,没心思看这样的代码,全是计算.
估计mx和my是某点的坐标(我才是圆心)
R应该就是半径
mse就不知道了.
再问: 。。。就是想知道mse是什么。。。。。。。。。帮忙看看行吗。。给你追加。。。今天早上就要交作业了。。其中一个函数没看太懂。。
再答: e=((x-mx).^2+(y-my).^2-R*R)/(R*R); 跟据这个式子可以判断e的意义和算法 x,y如果是数字的话 mse=mean(e.*e);好像没什么意义 x,y如果是向量的话 mse=mean(e.*e);那么返回值为该矩阵各列向量的均值 剩下的自己看着办吧,我matlab也是自学的、