Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
%
%see also http://www.matlabsky.com
$ Y) Q9 u- v7 R%
%===得到Bernsein基的次数================
clear;
X=load('data.txt');
%画出原图形的
plot(X(:,1),X(:,2));
title('原始数据');
%数据参数化,将X(:,1)变换为[0,1]间数据,
# v; D; O. F$ f% d1 r5 ~minX=min(X(:,1));
2 j9 n9 Q! l! r( u" I! m; t  PmaxX=max(X(:,1));
3 o7 L+ w, Q1 v0 X) B* eparaX=0;
for i=1:length(X)
0 ?9 B! k  d0 ]  P% z3 b( H% ^" k    paraX(i)=(X(i,1)-minX)/(maxX-minX);
+ y- A" V# q2 I6 L" f  x+ X* K; rend
m=input('请输入Bernsein基的次数m(即m+1阶):');
+ @. g5 O- ]. ~/ B- @; ib=OLS(X,paraX,m);
residu=residual(X,m,paraX,b);
; @3 |4 m0 W. H  i%循环得读取命令
) }1 ^; n% a  f5 u3 x# ~- cwhile 1
9 ?0 s* i. h. h1 [8 ]; `- Y    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
! w" _" t$ z" F( k    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
8 {! m& K/ o+ E) f: R    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
& d- d6 J, r+ ^# f# ^    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
    cmd=input('请按照以上提示选择操作:');

4 H3 C6 e! E; a    switch cmd
        case 1

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

宇仲

楼主辛苦了，继续加油啊！
# q$ w  j/ j) V4 y4 h