matlab中对于两个变量的求和函数代码该怎么写？

h=2.5;d=6;q=0.8;c=1.1;
c1=0.1;n=16;vv=[];
for theta = (10:20).*pi/180
 for b=2:.1:3
        v=DWS(b,theta);
        vv=[vv,v];%这里生成了一个1*121的一个矩阵
 end
end
b=10:20;
theta = 2:.1:3;
[bb,theta2]=meshgrid(b,theta);%把b和θ矩阵形式由1*11变成11*11
ww=reshape(vv,11,11);%把1*121的矩阵变成了一个11*11的矩阵
mesh(bb,theta2,ww)
ylabel('b')
xlabel('theta')
%用函数代码求最高点
%区间变小 调大精度
%代码简化 一个矩阵变成多列矩阵，能不能直接生成多列矩阵，或直接画图
%没有把tan（α)和tan（β）分开
%用python来写

function v=DWS(b,theta)%定义一个函数，输入b，θ就输出v
h=2.5;
d=6;
q=0.8;
c=1.1;
n=16;
k=1:n;
alpha=atand((h.*(d+(k-1).*q))./ ...%度数用atand，弧度用atan
    ((d+(k-1).*q).^2+(b-c+h-(k-1) ...
    .*q.*tan(theta)).*(b-c-(k-1).*q.*tan(theta))));
m=sum(alpha)/n;%度数用atand，弧度用atan
s=std(alpha(:),0);%std求无偏估计的方差
v=m/s;
end

这是数学建模中影院的视角与仰角的我写的代码，有点粗糙，但结果运行和课本是一样的，但是这里求和函数的处理我用的是分别输入B，θ，求得一共121个数据，再生成表格再画图，我觉得这一步有些复杂，运行起来需要很长时间，有没有更好的办法来处理这个求和函数