MATLAB-2: 拟合实例

最新推荐文章于 2024-03-16 22:17:55 发布

SemiconductorPhysics

最新推荐文章于 2024-03-16 22:17:55 发布

阅读量599

点赞数

分类专栏：编程文章标签： matlab 经验分享

本文链接：https://blog.csdn.net/Jelly_Zhou/article/details/122973867

版权

编程专栏收录该内容

51 篇文章 1 订阅

订阅专栏

函数

function [f, g, o, y_fit]=jfit(x,y)

if(size(x,2)>1)
    x=x';
end

if(size(y,2)>1)
    y=y';
end

k=diff(y)./diff(x);

jL=[-max(abs(k)) -max(abs(y))];
jU=[max(abs(k)) max(abs(y))];
jSt=[mean(k) mean(y)];

Method='NonlinearLeastSquares';
Robust='off';

MaxFunEvals=100000; 
MaxIter=100000;
Algorithm='Trust-Region';
TolX=1e-9;
TolFun=1e-9; 

jType={'a*x+b'};

s = fitoptions(  'Method',Method,...  
                'Robust',Robust,...
                'Lower',jL,...
                'Upper',jU,...
                'MaxFunEvals',MaxFunEvals,...
                'MaxIter',MaxIter,...
                'TolFun',TolFun,...
                'Algorithm',Algorithm,...
                'TolX',TolX,...
                'Startpoint',jSt);        

g1=fittype(jType{1},'options',s);

[f, g, o]=fit(x,y,g1);
y_fit=f(x);
plot(x,y_fit,'-');
hold on;


end

调用

clear
close 
clc

x = 0:20;
y = -3*x+10+2*randn(size(x));
figure
plot(x,y,'o');
hold on

[f,g,o] = fit(x',y','poly1');
y_fit = f(x);
plot(x,y_fit,'-');
hold on;

f
r2=g.rsquare
coeff=coeffvalues(f)
confid=confint(f)

输出


f = 

     Linear model Poly1:
     f(x) = p1*x + p2
     Coefficients (with 95% confidence bounds):
       p1 =      -2.994  (-3.127, -2.862)
       p2 =       11.35  (9.805, 12.9)

r2 =

    0.9916


coeff =

   -2.9943   11.3520


confid =

   -3.1267    9.8049
   -2.8620   12.8990