当前位置: 首页 > 程序语言 >关于zernike多项式展开的系数问题

关于zernike多项式展开的系数问题

作者 ilovexiaomu
来源: 小木虫 150 3 举报帖子
+关注

各位老师:
最近遇到一个关于zernike多项式展开的系数问题,如下:
function Z = zernike(i, r, theta)

load('zernike_index'); % load the mapping of (n,m) to i
n = zernike_index(i,1);
m = zernike_index(i,2);
if m==0
Z = sqrt(n+1)*zrf(n,0,r);
else
if mod(i,2) == 0 % i is even
Z = sqrt(2*(n+1))*zrf(n,m,r) .* cos(m*theta);
else % i is odd
Z = sqrt(2*(n+1))*zrf(n,m,r) .* sin(m*theta);
end
end
return
% Zernike radial function
function R = zrf(n, m, r)
R = 0;
for s = 0 : ((n-m)/2)
num = (-1)^s * gamma(n-s+1);
denom = gamma(s+1) * gamma((n+m)/2-s+1) ...
* gamma((n-m)/2-s+1);
R = R + num / denom * r.^(n-2*s);
end

上面是代码,下面我简要说一下问题,
我做符号哪里那个load('zernike_index'),需要我load一个zernike_index的文件,这个文件是这样的一个图片,
我怎么把这个图片里面的前三项 (n,m,i)变成能用的zernike_index文件呢?

可以用excel表格么?
望各位大师帮忙。

关于zernike多项式展开的系数问题
QQ图片20170711103100.png 返回小木虫查看更多

今日热帖
  • 精华评论
  • somomo91

    你发帖的时间,已经足够你手动输入 72 个数据很多遍了

  • Mr__Right

    参考这里的公式,自己生成一个矩阵
    替代需要load的文件不就可以了,

猜你喜欢
下载小木虫APP
与700万科研达人随时交流
  • 二维码
  • IOS
  • 安卓