谁能帮忙看个程序

clear all;close all
a=0.005;%定义盲孔半径为a=1；
sigama_1=120;
sigama_2=100;%假设实际的残余应力值；
u=0.28;
E=2.1e11;%定义材料为合金钢，给出材料的泊松比及弹性模量；
r=0.010;%选择参考点为半径为r=2圆上的点；
theta=2*pi*rand(1,100);%选择参考点的角度参数；
A=0.5*(1+u)*a^2/r^2;
B=0.5*(3*(1+u)*a^4/r^4-4*a^2/r^2);%推导公式中的A、B的代数式；
H=zeros(100,2);%定义矩阵H；
Y=zeros(100,1);%定义矩阵Y;
S=zeros(100,1);%定义矩阵S；
for k=1:1:100
H(k,:)=1/E*[A+B*cos(2*theta(k)),A-B*cos(2*theta(k))];
S(k)=1/E*[A+B*cos(2*theta(k)),A-B*cos(2*theta(k))]*[sigama_1;sigama_2];%for循环，实际应变推导公式；
end
coeff=0.01;%定义参考系数；
Y=S+coeff*randn(size(S));%测量得出的应变值；
X=inv(H'*H)*H'*Y