Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
1 H2 L" X+ e  O%
%see also http://www.matlabsky.com
$ n& }4 q6 U1 k  v9 ]%
1 q* I, i6 ]& j* E7 f%===得到Bernsein基的次数================
3 |2 U3 w' p6 Yclear;
7 l, |/ P: y) Q( U1 M! e: GX=load('data.txt');
%画出原图形的
' _! H: ~* @! j& E3 q1 Z: P' \; Oplot(X(:,1),X(:,2));
3 b  k! A: a/ F8 K4 j- ]title('原始数据');
3 Q2 j) O) x. t( F%数据参数化,将X(:,1)变换为[0,1]间数据,
minX=min(X(:,1));
7 b8 T' L! v% Y  G. @maxX=max(X(:,1));
paraX=0;
for i=1:length(X)
" l' m, U9 H& L. B' j    paraX(i)=(X(i,1)-minX)/(maxX-minX);
end
7 p+ ~6 Z1 O% ?; Mm=input('请输入Bernsein基的次数m(即m+1阶):');
b=OLS(X,paraX,m);
! H6 F" m% F2 ~5 Kresidu=residual(X,m,paraX,b);
%循环得读取命令
while 1
    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
3 i+ {$ x1 a8 c6 V1 w    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
0 K# Y" ^, h- B+ F  }    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
    cmd=input('请按照以上提示选择操作:');
- b. H. o1 q( G+ k1 ]" F0 p
) X9 N5 z; P  h0 [0 l7 x( @1 E. a( L* e$ I    switch cmd
        case 1
5 y3 r% ^9 J7 s! |; B$ Q0 u) y
( k0 F3 \* D4 o! f! i

游客，如果您要查看本帖隐藏内容请回复

宇仲

楼主辛苦了，继续加油啊！