Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
%
1 f- \) @. {3 S' Z  m%see also http://www.matlabsky.com
3 _: X/ u/ \- v: Z0 h$ F%
%===得到Bernsein基的次数================
* c$ N" b- @; |1 y7 [' oclear;
X=load('data.txt');
%画出原图形的
plot(X(:,1),X(:,2));
+ O' k9 V0 T. V2 utitle('原始数据');
+ G0 y/ }+ ?4 E2 M& t. s%数据参数化,将X(:,1)变换为[0,1]间数据,
9 N- Q2 P9 _$ ?  g- ZminX=min(X(:,1));
maxX=max(X(:,1));
paraX=0;
1 y: N$ W+ S: j$ _for i=1:length(X)
    paraX(i)=(X(i,1)-minX)/(maxX-minX);
end
1 n, x+ x+ P- v8 A- q+ i# vm=input('请输入Bernsein基的次数m(即m+1阶):');
b=OLS(X,paraX,m);
residu=residual(X,m,paraX,b);
& d/ O9 x- [, E* s%循环得读取命令
  r6 C2 K. d3 r+ O( H6 m! swhile 1
' x. m4 X6 s. C9 D& y+ @    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
" m2 a) O7 q; n$ x; k    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
    cmd=input('请按照以上提示选择操作:');

' h7 ~! p* a/ G6 ?. X    switch cmd
        case 1
9 X0 h( J$ D9 S( J
( }  \0 z# L0 H6 {9 `5 E

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

& y% w  T) D8 p7 B& _4 L

宇仲

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