Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
0 \$ @$ E5 e* M5 \4 r" I%
9 V" a8 U4 Y( d8 d  d8 E%see also http://www.matlabsky.com
; P+ Y, x( E7 e; _  m%
7 i' }, h, {# B# D+ D( [# ^%===得到Bernsein基的次数================
clear;
X=load('data.txt');
%画出原图形的
$ s) V. d" Q3 M1 O1 K7 M) kplot(X(:,1),X(:,2));
. k' N- Z& w! {2 |7 }; Ttitle('原始数据');
1 _8 }4 P0 _  x2 v2 F%数据参数化,将X(:,1)变换为[0,1]间数据,
minX=min(X(:,1));
; J9 B- D/ c# w: k: ?+ s. d! q" NmaxX=max(X(:,1));
paraX=0;
for i=1:length(X)
    paraX(i)=(X(i,1)-minX)/(maxX-minX);
7 x: y' v2 E" m  Fend
4 f  `! t# O; a# H9 lm=input('请输入Bernsein基的次数m(即m+1阶):');
b=OLS(X,paraX,m);
5 J8 Q0 Q. _1 Uresidu=residual(X,m,paraX,b);
% p8 i; j$ Y; L0 _& Y3 L%循环得读取命令
while 1
0 Y) ~% u6 Q7 P$ X    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
" V) T4 d/ @* s- D$ v    cmd=input('请按照以上提示选择操作:');

$ ^6 E2 P  x' N/ [# E    switch cmd
. `- D6 \6 ^! b2 A5 h5 D- d        case 1

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

宇仲

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