Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
( U9 C' w& _. o! ^%
6 H8 A# I" }: u& h. T%see also http://www.matlabsky.com
%
: H$ k! M% O) K: o1 B%===得到Bernsein基的次数================
* }  H0 O6 }$ k' Q3 N- {clear;
X=load('data.txt');
: V7 [5 \0 N5 s/ o2 h%画出原图形的
+ C9 K" P* z* d& F2 Aplot(X(:,1),X(:,2));
- f4 t2 V* z$ K/ o: c4 gtitle('原始数据');
%数据参数化,将X(:,1)变换为[0,1]间数据,
# d' A( X& B- Y1 f( jminX=min(X(:,1));
7 p; p3 z; G+ {7 KmaxX=max(X(:,1));
paraX=0;
- Y+ n$ `8 k. _2 j4 E. N$ tfor i=1:length(X)
    paraX(i)=(X(i,1)-minX)/(maxX-minX);
end
9 Z2 w+ n+ G4 }5 c+ Vm=input('请输入Bernsein基的次数m(即m+1阶):');
b=OLS(X,paraX,m);
residu=residual(X,m,paraX,b);
$ b  Z! z! v% w: j7 J4 U! p%循环得读取命令
. a! o6 _$ x1 h& mwhile 1
5 d/ d2 }; d6 l    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
1 _5 I  J7 v$ O2 g7 c    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
    cmd=input('请按照以上提示选择操作:');
; [5 ?: a1 r9 A1 l# X8 F, \8 P+ S
    switch cmd
        case 1
" B+ Q4 O* b8 `1 a' p

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

* {$ Y" H/ V' ?) \! h4 G& {

宇仲

楼主辛苦了，继续加油啊！
& c# G* b2 k3 I* n8 i: z( b( |7 x