| 5 | 1/1 | 返回列表 |
| 查看: 3860 | 回復(fù): 4 | ||
wind5fsa金蟲 (小有名氣)
|
[求助]
報(bào)錯(cuò)求教Undefined functionfor input arguments of type 'double'.
|
|
大家好,我的matlab程序里面有一部分是我在網(wǎng)上下的一個(gè)擬合橢圓的程序,該部分程序如下所示,能將我程序前面部分得到的坐標(biāo)為(X,Y)的一堆數(shù)擬合成橢圓,寫完之后一直順利計(jì)算。但是最近忽然報(bào)錯(cuò): Undefined function 'fitellipse' for input arguments of type 'double'. Error in F0_1Single_4zhang_ZXEDGEDINGDIAN (line 323) PS:323這一行的程序是: fitellipse(nx,ny) 然后再怎么運(yùn)行也都是這么報(bào)錯(cuò)。求指教,謝謝! 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è)作家)
中國(guó)特色

金蟲 (小有名氣)

金蟲 (小有名氣)

金蟲 (小有名氣)
| 5 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿天津大學(xué)化學(xué)工藝專業(yè)(081702)315分求調(diào)劑 +11 | yangfz 2026-03-17 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +3 | Mqqqqqq 2026-03-19 | 3/150 |
|
|
[考研] 0817調(diào)劑 +3 | 沒有答案_ 2026-03-14 | 3/150 |
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院,一志愿華科070200學(xué)碩,346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
[考研] 材料專碩英一數(shù)二306 +5 | z1z2z3879 2026-03-18 | 5/250 |
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
[考研] 298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑 +7 | 手機(jī)用戶 2026-03-17 | 7/350 |
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
[考研]
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 一志愿南京大學(xué),080500材料科學(xué)與工程,調(diào)劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語(yǔ) 2026-03-15 | 5/250 |
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
[考研] 招收0805(材料)調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|