各位板眾大家好 自己摸索學習Matlab 卡在這邊很久了 所以來請教大家 我有一組x.y的數據約600筆 想要對已知"非線性方程"做曲線擬合 方程式為 x=(log((y/b)*((b-c)/(y-c)))-c*(1/y-1/b))/(a*c^2) 想求未知參數a.b.c 並求得R2值 我嘗試用了兩個方法 1. curvefiting 我用Toolboxes裡的Curve Fitting 選擇Custom Equations來輸入自訂函數 可是就是跑不出來 出現 Fit could not be computed due to error: Inf computed by model function, fitting cannot continue. Try using or tightening upper and lower bounds on coefficients. 但是對於參數a.b.c的值 大概有個底 所以對於初始值及上下限都有給值 問題會出現在這邊嗎? 2. lsqcurvefit M.file 輸入 ----------------------------------------------- function F=myfun(x,J); F=(log([(J*(x(2)-x(3)))*diag((x(2)*(J-x(3))).^-1)])-x(3)*(J.^-1-x(2)^-1)) *((x(1)*x(3)^2)^-1) ----------------------------------------------- Command Window輸入 clear;clc xlsFile='data9910031.xls'; x(1,:)=xlsread(xlsFile,'sheet1','A1:A660')'; y(1,:)=xlsread(xlsFile,'sheet1','B1:B660')'; x0=[500000 0.0001 0.00002]; [x,resnorm]=lsqcurvefit(@myfun,x0,x,y) format long; ----------------------------------------------- 擬合的結果很差 也出現根本不可能的負值 還有 請問R2有指令嗎? corrcoef我不會用T^T 我已經困在這裡很久了 請大家有意見的請告訴我~ 謝謝各位!!! --

※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 140.115.150.220