| 5 | 1/1 | 返回列表 |
| 查看: 3872 | 回復(fù): 4 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
wind5fsa金蟲 (小有名氣)
|
[求助]
報錯求教Undefined functionfor input arguments of type 'double'.
|
|
|
大家好,我的matlab程序里面有一部分是我在網(wǎng)上下的一個擬合橢圓的程序,該部分程序如下所示,能將我程序前面部分得到的坐標(biāo)為(X,Y)的一堆數(shù)擬合成橢圓,寫完之后一直順利計算。但是最近忽然報錯: Undefined function 'fitellipse' for input arguments of type 'double'. Error in F0_1Single_4zhang_ZXEDGEDINGDIAN (line 323) PS:323這一行的程序是: fitellipse(nx,ny) 然后再怎么運行也都是這么報錯。求指教,謝謝! if nargin == 0 % Create an ellipse t = linspace(0,2); Rx = 300 Ry = 200 Cx = 250 Cy = 150 Rotation = .4 % Radians x = Rx * cos(t); y = Ry * sin(t); nx = x*cos(Rotation)-y*sin(Rotation) + Cx; ny = x*sin(Rotation)+y*cos(Rotation) + Cy; % Draw it plot(nx,ny,'o'); % Fit it fitellipse(nx,ny) % Note it returns (Rotation - pi/2) and swapped radii, this is fine. return end % normalize data mx = mean(X); my = mean(Y); sx = (max(X)-min(X))/2; sy = (max(Y)-min(Y))/2; x = (X-mx)/sx; y = (Y-my)/sy; % Force to column vectors x = x( ;y = y( ;% Build design matrix D = [ x.*x x.*y y.*y x y ones(size(x)) ]; % Build scatter matrix S = D'*D; % Build 6x6 constraint matrix C(6,6) = 0; C(1,3) = -2; C(2,2) = 1; C(3,1) = -2; % Solve eigensystem [gevec, geval] = eig(S,C); % Find the negative eigenvalue I = find(real(diag(geval)) < 1e-8 & ~isinf(diag(geval))); % Extract eigenvector corresponding to negative eigenvalue A = real(gevec(:,I)); % unnormalize par = [ A(1)*sy*sy, ... A(2)*sx*sy, ... A(3)*sx*sx, ... -2*A(1)*sy*sy*mx - A(2)*sx*sy*my + A(4)*sx*sy*sy, ... -A(2)*sx*sy*mx - 2*A(3)*sx*sx*my + A(5)*sx*sx*sy, ... A(1)*sy*sy*mx*mx + A(2)*sx*sy*mx*my + A(3)*sx*sx*my*my ... - A(4)*sx*sy*sy*mx - A(5)*sx*sx*sy*my ... + A(6)*sx*sx*sy*sy ... ]'; % Convert to geometric radii, and centers thetarad = 0.5*atan2(par(2),par(1) - par(3)); cost = cos(thetarad); sint = sin(thetarad); sin_squared = sint.*sint; cos_squared = cost.*cost; cos_sin = sint .* cost; Ao = par(6); Au = par(4) .* cost + par(5) .* sint; Av = - par(4) .* sint + par(5) .* cost; Auu = par(1) .* cos_squared + par(3) .* sin_squared + par(2) .* cos_sin; Avv = par(1) .* sin_squared + par(3) .* cos_squared - par(2) .* cos_sin; % ROTATED = [Ao Au Av Auu Avv] tuCentre = - Au./(2.*Auu); tvCentre = - Av./(2.*Avv); wCentre = Ao - Auu.*tuCentre.*tuCentre - Avv.*tvCentre.*tvCentre; uCentre = tuCentre .* cost - tvCentre .* sint; vCentre = tuCentre .* sint + tvCentre .* cost; Ru = -wCentre./Auu; Rv = -wCentre./Avv; Ru = sqrt(abs(Ru)).*sign(Ru); Rv = sqrt(abs(Rv)).*sign(Rv); a = [uCentre, vCentre, Ru, Rv, thetarad]; |

金蟲 (小有名氣)

木蟲 (職業(yè)作家)
中國特色

金蟲 (小有名氣)

金蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[基金申請] 面上5B能上會嗎? +4 | redcom 2026-03-29 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿雙一流機械285分求調(diào)劑 +4 | 幸運的三木 2026-03-29 | 5/250 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學(xué)材料與化工方向336分 +9 | 辰沐5211314 2026-03-26 | 9/450 |
|
|
[考研] 085701求調(diào)劑初試286分 +4 | secret0328 2026-03-28 | 4/200 |
|
|
[考研] 一志愿北化085600材料專碩275|有文章專利|求調(diào)劑 +7 | Micky11223 2026-03-25 | 7/350 |
|
|
[考研] 394求調(diào)劑 +3 | 好事多磨靜候佳?/a> 2026-03-26 | 5/250 |
|
|
[考研] 317求調(diào)劑 +6 | 十閑wx 2026-03-24 | 6/300 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +8 | 邱gl 2026-03-27 | 8/400 |
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
[考研] 材料求調(diào)劑 一志愿哈工大總分298分,前三科223分 +5 | dongfang59 2026-03-27 | 5/250 |
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
[考研] 316求調(diào)劑 +5 | Pigcasso 2026-03-24 | 5/250 |
|
|
[考研] 325求調(diào)劑 +5 | 李嘉圖·S·路 2026-03-23 | 5/250 |
|
|
[考研] 290分調(diào)劑求助 +3 | 吉祥止止陳 2026-03-25 | 3/150 |
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
[有機交流]
20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|
[考研] 080500求調(diào)劑 +3 | zzzzfan 2026-03-24 | 3/150 |
|
|
[考研] 277分求調(diào)劑,跨調(diào)材料 +3 | 考研調(diào)劑lxh 2026-03-24 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境,總分308求調(diào)劑 +7 | 墨墨漠 2026-03-23 | 8/400 |
|