matlab学习笔记

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

本科的时候没有接触过MATLAB，现在要用，只能东拼西凑，摸着石头过河，按博主的帖子学习曲线拟合，原博客地址https://blog.csdn.net/u013360881/article/details/44782057

图片摘自百度文库，侵删

不是纯粹复制粘贴，有自己的修改和体会在上面

此博客仅作记录

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

matlab的曲线拟合工具箱cftool，拟合函数及指标：

· Custom Equations：用户自定义的函数类型
· Exponential：指数逼近，有2种类型， a*exp(b*x) 、 a*exp(b*x) + c*exp(d*x)
· Fourier：傅立叶逼近，有7种类型，基础型是 a0 + a1*cos(x*w) + b1*sin(x*w)
· Gaussian：高斯逼近，有8种类型，基础型是 a1*exp(-((x-b1)/c1)^2)
· Interpolant：插值逼近，有4种类型，linear、nearest neighbor、cubic spline、shape-preserving
· Polynomial：多形式逼近，有9种类型，linear ~、quadratic ~、cubic ~、4-9th degree ~
· Power：幂逼近，有2种类型，a*x^b 、a*x^b + c
· Rational：有理数逼近，分子、分母共有的类型是linear ~、quadratic ~、cubic ~、4-5th degree ~；此外，分子还包括constant型
· Smoothing Spline：平滑逼近
· Sum of Sin Functions：正弦曲线逼近，有8种类型，基础型是 a1*sin(b1*x + c1)
· Weibull：只有一种，a*b*x^(b-1)*exp(-a*x^b)

选择好所需的拟合曲线类型及其子类型，并进行相关设置：
——如果是非自定义的类型，根据实际需要点击“Fit options”按钮，设置拟合算法、修改待估计参数的上下限等参数；
——如果选Custom Equations，点击“New”按钮，弹出自定义函数等式窗口，有“Linear Equations线性等式”和“General Equations构造等式”两种标签。

 

一种非线性函数的曲线拟合方法（函数公式： k = A*(T^a)*exp(E/T) ）https://blog.csdn.net/chenyusiyuan/article/details/1942624

感谢大神实例

源码：

 

close all
logT = log(T);
logK = log(K);
daoT = T.^(-1);
C = ones(3);
xMat = [];
lenT = length(T);

for r=1:1:lenT-2
    C(:,2) = [logT(r),logT(r+1),logT(r+2)];
    C(:,3) = [daoT(r),daoT(r+1),daoT(r+2)];
    b = [logK(r),logK(r+1),logK(r+2)];
    x = C\(b');
    xMat = [xMat,x];
end
xMat_1 = xMat';
mean_xMat = mean(xMat_1);
A = 3.8858e-020;
a = mean_xMat(2);
E = mean_xMat(3);

h1=stem(T,K,'bo');  % ‘bo’表示每个点用一个小圆圈表示
    set(h1,'MarkerFaceColor','green');  % 小圆圈内的颜色为绿色
    set(h1,'LineStyle','none');     % 隐藏基线到点的连线
    set(get(h1,'BaseLine'),'LineStyle','none'); % 隐藏基线
    hold on;    % 保持由点集构成的目标曲线，以便和拟合曲线进行对比
    
    t=200:10:1300;
    k = [];
    for i=1:1:111
    k(1,i)=A*(t(i)^a)*exp(E/t(i));
    end
    plot(t,k,'r'); 
    xlabel('T');   ylabel('K');   title('Nonlinear Curve Fitting');