恳请指点关于“函数分布的假设检验”的一段小程序

烦劳各位前辈能否帮晚辈看一看下面的这段程序，是对具体的函数分布进行假设检验的，
在调试的过程中总报错，显示SHM没有定义，但是SHM是一个矩阵，矩阵在这里需要
初始化吗？不好意思，打扰了，恳请帮助！不胜感激，致礼！


fprintf('对该产品的伪寿命进行威布尔分布的假设检验：')%程序输出的文字

phat=weibfit(SHM,alpha);
p2=weibcdf(SHM,phat(1),phat(2));
[H2,s2]=kstest(SHM,[SHM,p2],alpha)
if H2==O
fprintf('该产品的伪寿命服从双参数分别为%f和%f的威布尔分布:',phat(1),phat(2))
fprintf('伪寿命服从威布尔分布的概率为%f; ',s2)
fprintf('n')
[muhat2, sigmahat2]=weibstat(phat(1), phat(2));
else
disp('该产品的伪寿命不服从威布尔分布；');
s2=0;
end


以上程序中，函数weibfit的输入部分：SHM表示伪寿命矩阵；alpha表示指定
的显著性水平。函数weibfit的输出部分：phat表示将伪寿命与weibull分布拟合产生的weibull分布的双参数矢量(注：phat(1)表示该双参数weibull分布的第一个参数——尺度参数；phat(2)表示该双参数weibull分布的第二个参数——形状参数)；
函数weibcdf的输入部分：SHM表示伪寿命矩阵；phat(1)表示尺度参数；phat
(2)表示形状参数。函数weibcdf的输出部分：p2表示WeibuU分布在取尺度参数phat(1)，形状参数phat(2)，取值为SHM时的累积概率分布函数值。
函数kstest的输入部分：SHM表示伪寿命矩阵；[SHM，p2]表示由SHM和p2构
成的矩阵(SHM和p2的含义见上)；alpha表示指定的显著性水平。函数kstest的输出部分：若H2等于0则服从weibull分布；否则不服从该分布；s2为服从weibull分布假设检验的概率。
函数weibstat的输入部分：phat(1)表示尺度参数；phat(2)表示形状参数；函数weibstat的输出部分：muhat2表示双参数分别为phat(1)与phat(2)的weibull的均值估计值；sigmahat2表示双参数分别为phat(1)与phat(2)的weibull的方差估计值。