函数文件由函数声明行、H1行、在线帮助文本区、编写和修改记录、函数主体等几个部分组成。
注意:
1.在线帮助文档区必须在H1行下一行,中间不能空行。在线帮助文档区中每行前面都要有注释号%
,空行也要有
2.H1简要说明该函数的功能
3.在线帮助文档区写明输入代表的是什么、输出代表的是什么。如果函数有不同的调用形式则要说明每个调用的方式的功能
示例:
function [ p ] = arch( obj_filename )
% arch 求牙齿模型的三次拟合牙弓
% 输入参数:
% obj_filename: 【切掉底座】的牙齿模型的文件名(obj格式)
% 输出参数:
% p:三次拟合曲线的四个参数,整条曲线可以用y=polyval(p,x)还原出来,x是牙齿模型的宽度
% 调用示例:
% p=arch('data-toothAndGum.obj');得到三次拟合的四个多项式系数
% arch('data-toothAndGum.obj');将画图输出部分取消注释,得到二维牙齿与牙弓的图
%读取obj得到点集和面集
[v,f]=read_obj(obj_filename);
%得到曲率
tri.faces=f';
tri.vertices=v';
[Cmean,Cgaussian,Dir1,Dir2,Lambda1,Lambda2,vertex_curv]=patchcurvature(tri,true);
% 提取曲率>0,5的点
C=Cmean';
x=v(1,:);
y=v(3,:);
x1=x(:,C(1,:)>0.5);
y1=y(:,C(1,:)>0.5);
p = polyfit(x1,y1,3);%p是三次拟合曲线的四个参数,后面可以用ployval还原出来
end
然后就可以用help查看自己的函数啦!