Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
* u/ z' [$ ]! h9 @, s; ?%
%see also http://www.matlabsky.com
%
%===得到Bernsein基的次数================
( S* B6 B8 p" }: Q5 d( Rclear;
X=load('data.txt');
%画出原图形的
plot(X(:,1),X(:,2));
) [; c  _" \! _# Q8 Ctitle('原始数据');
%数据参数化,将X(:,1)变换为[0,1]间数据,
minX=min(X(:,1));
maxX=max(X(:,1));
# t4 H9 g* K: cparaX=0;
$ U  o3 i  H7 J2 {for i=1:length(X)
    paraX(i)=(X(i,1)-minX)/(maxX-minX);
( z) I, u: J: x3 q% ?end
, A; |/ c* |* Hm=input('请输入Bernsein基的次数m(即m+1阶):');
b=OLS(X,paraX,m);
residu=residual(X,m,paraX,b);
%循环得读取命令
/ B( O; V% A. H- Zwhile 1
    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
8 k/ S; J/ q. K, v1 l  ]    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
' ]6 e+ G" ]2 z/ m  T    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
0 K+ r9 e7 O" P1 L7 \5 L    cmd=input('请按照以上提示选择操作:');

    switch cmd
& u1 Y# e8 W- Y7 k        case 1
- C7 U& N- z: n, E+ U/ r

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

% s, T' ^" b7 V, E
' @0 x2 R1 ]- M

宇仲

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