【MATLAB】变形高斯函数积分求解与多项式拟合（善用矩阵）

创建函数法

inline(expr)
inline(expr,arg1,arg2,…)
.

f=inline('1/(sqrt(2*pi)).*exp(-(x.^2)/2','x');
x=0:100;
y=f(x);

.

概率密度函数积分

高斯分布函数的变形使得函数无法通过分析积分求解（该函数不可积）；采用矩阵计算的形式轻松解决了变量H1的问题，得到该概率密度函数的积分。

ni=100;
hsss=linspace(0,3,ni);
hi=ni;
H1=hsss/sigma;
f=@(x)(1/(sqrt(2*pi)).*(x-H1).^(3/2).*exp(-(x.^2)/2));
%求概率密度函数积分
x0=H1;
xu=15*ones(1,hi);
N=1000;
htrap=((xu-x0)/N);
area=zeros(1,hi);
while(x0<xu)
    area=area+(htrap/2).*(f(x0)+f(x0+htrap));
    x0=x0+htrap;
end

.

多项式拟合

采用多项式拟合函数polyfit、polyval

%曲线拟合
[p1,S]=polyfit(H1,area,4);%S 中normr为残差的范数（评价拟合的参数）
y1=polyval(p1,H1,S);
deltam=max((abs(y1-area))/area);