Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
, [3 b& y% C( a%
%see also http://www.matlabsky.com
0 L! A  [7 J  |7 Z% D, r( L%
%===得到Bernsein基的次数================
$ [/ {# p2 D2 F7 M; D* sclear;
X=load('data.txt');
%画出原图形的
  E2 P; N% X: [" f3 a" Eplot(X(:,1),X(:,2));
title('原始数据');
/ c& D( x$ N0 P' F%数据参数化,将X(:,1)变换为[0,1]间数据,
. t& ~8 s, E3 q! U. O5 tminX=min(X(:,1));
8 R! `) e: q. N" F! nmaxX=max(X(:,1));
+ m+ V2 B8 v1 f; S  o; d' uparaX=0;
& z; i, |" ]- O  {8 a6 d* `; k0 efor i=1:length(X)
" N8 m* q4 ^& k# q; R% E" I% ]) M    paraX(i)=(X(i,1)-minX)/(maxX-minX);
$ H3 i) r, i/ Q. `2 z, \: Eend
7 j3 W3 O" |* e& xm=input('请输入Bernsein基的次数m(即m+1阶):');
3 s1 q1 s# U' D3 c/ k7 e4 s, M: `# C: x9 ]b=OLS(X,paraX,m);
7 `7 s! J! M7 F1 S" c- Vresidu=residual(X,m,paraX,b);
%循环得读取命令
' L/ R: u4 ]" p6 e6 c" h. }while 1
    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
( z# t& e8 z* b; F; e    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
& k; M5 S+ D! d% n    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
$ p0 |# q" p& ]! T! B" e    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
1 Y6 P7 M0 P( I+ Y5 C" u" w    cmd=input('请按照以上提示选择操作:');

% Y% R+ f0 o. M* H$ O; l    switch cmd
; ~! w7 P4 ?* c) s. y- _: Z/ F( Z        case 1
  x# |$ @0 l+ u7 I* F- q. _

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

宇仲

楼主辛苦了，继续加油啊！
4 s, q' X" Q* k