创建函数法
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);