怎么用matlab计算螺旋线曲率

Question

目前我想画一条任意的螺旋线，然后求出该螺旋线起始位置点和结束位置点的曲率，然后再推广求该螺旋线上任意一点的曲率。
代码如下

b=100;%螺线间距参数
t = 0*pi:0.01*pi:1.5*pi;%
r=a+b*t;%螺线方程
x = round(r.*cos(t));%获取直角坐标系下x轴的坐标并进行四舍五入
y = round(r.*sin(t));%获取直角坐标系下y轴的坐标
plot(x,y);
k=1+(((b.*sin(t)+r.*cos(t)).^2).^3/2)/((2*b.*sin(t)+2*b.*cos(t)-r.*sin(t)+r.*cos(t))/(b.*cos(t)-r.sin(t)).^3);

目前存在的问题是

此类型的变量不支持使用点进行索引。
出错 (line 8)
k=1+(((b.*sin(t)+r.cos(t)).^2).^3/2)/((2b.sin(t)+2b.*cos(t)-r.*sin(t)+r.*cos(t))/(b.*cos(t)-r.sin(t)).^3)

目前工作因为这个问题停滞，一直查不出原因

Answer 1

最后一句改成这样

k=1+(((b.*sin(t)+r.*cos(t)).^2).^3/2)/((2*b.*sin(t)+2*b.*cos(t)-r.*sin(t)+r.*cos(t))/(b.*cos(t)-r.*sin(t)).^3);

有用麻烦采纳一下，谢谢