[討論] 擬合分式係數的問題
各位大大好,我最近在嘗試擬合分式的係數,但一直不知道問題在哪裡
例如:
%產生數據
function y=mydata(x)
y=1./(x-2);
%誤差函數
function E=myerr(a,x,y)
Y=a(1)./a(2)*x+a(3);
E=sum((y-Y).^2);
%最小誤差
clc
clear all
u=100;
for s=1:u
x(s)=(s-1)*0.1;
end
y=mydata(x)
x=x'
y=y'
a0=[0.1,0.1,0.1]
for i=1:10
a=fminsearch('myerr',a0,[],x,y)
a0=a
end
但一直顯示這個問題,擬合不出係數
Exiting: Maximum number of function evaluations has been exceeded
- increase MaxFunEvals option.
Current function value: Inf
在擬合一般多項式的時候都沒什麼問題,請問各位大大能否直接擬合出分式的係數,
像能否直接擬合出y=(a(1)*x^2+a(2)*x+a(3))/(a(4)*x^2+a(5)*x+a(6)) 係數a
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 114.35.135.39