Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
/ P* P. y+ u0 ]4 B%
9 T  E' C/ p( V%see also http://www.matlabsky.com
%
%===得到Bernsein基的次数================
" J: T) ]+ E' U" ?* N* A5 q5 Bclear;
X=load('data.txt');
%画出原图形的
plot(X(:,1),X(:,2));
- H' H: H8 _9 a+ atitle('原始数据');
( R+ u$ \, `, `% X8 s%数据参数化,将X(:,1)变换为[0,1]间数据,
minX=min(X(:,1));
) {3 E- N& d' f7 z6 O" M* bmaxX=max(X(:,1));
/ b; s8 E7 f: M1 ^, q" V0 JparaX=0;
1 Q7 Y( b7 o( X- ]& d' A5 C9 D; s: Rfor i=1:length(X)
    paraX(i)=(X(i,1)-minX)/(maxX-minX);
, X. |5 j# r/ n" c5 i6 Mend
- x8 |9 s4 S1 c, ]# Xm=input('请输入Bernsein基的次数m(即m+1阶):');
b=OLS(X,paraX,m);
residu=residual(X,m,paraX,b);
%循环得读取命令
while 1
6 _3 ]' t# t7 l) c4 R6 O- k$ e    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
6 e5 P- e6 D' M, c5 d5 G    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
9 g+ e0 v0 r* F: F! h! }9 D  n    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
    cmd=input('请按照以上提示选择操作:');

    switch cmd
2 P  f' \2 l; ]8 E3 }        case 1

! y8 [0 R' S( }' T6 D! v

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

宇仲

楼主辛苦了，继续加油啊！
/ v( p* l( W6 Z$ f6 R; O: U4 h3 j6 ?# o