Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
%
( A* Z! v2 l' p1 ^9 Y& l* y%see also http://www.matlabsky.com
; [* T( f) q+ u' O4 ]8 d%
%===得到Bernsein基的次数================
clear;
5 A( J$ [/ {0 F3 ^X=load('data.txt');
' j, D8 d3 g. H9 y% _6 s2 T3 F%画出原图形的
plot(X(:,1),X(:,2));
title('原始数据');
7 s( `1 j0 I& A2 w  G$ o%数据参数化,将X(:,1)变换为[0,1]间数据,
1 `: T( ^  H/ [9 r6 y, F. dminX=min(X(:,1));
8 @. o& ]) Q- {! _+ a! ]maxX=max(X(:,1));
. Z7 v: m$ u  x; Y3 Q3 a; tparaX=0;
for i=1:length(X)
$ n* P) Q1 e2 B; {2 H7 l/ r    paraX(i)=(X(i,1)-minX)/(maxX-minX);
end
3 s* G3 X, i- r& sm=input('请输入Bernsein基的次数m(即m+1阶):');
$ ]( e" k" E, P2 F" |; }* ~b=OLS(X,paraX,m);
residu=residual(X,m,paraX,b);
%循环得读取命令
while 1
    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 ');
' g$ R# p& s* i" ?3 X7 r    cmd=input('请按照以上提示选择操作:');

    switch cmd
        case 1
& ?! F' P  t3 }/ s9 [! g

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

" d0 f0 Z4 p( E
# P% j9 w" k: @4 {

宇仲

楼主辛苦了，继续加油啊！
4 K( }+ ~" C' i' H$ [. S4 n0 o  E