Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
%
# c# Q" l( i5 Q9 f" S" ?%see also http://www.matlabsky.com
%
( v, @8 x( J  ?5 e3 s- N- I%===得到Bernsein基的次数================
clear;
2 e4 m( a6 r& F7 j& XX=load('data.txt');
( A% U7 u! i3 S! a: ?" R' W- L. o%画出原图形的
plot(X(:,1),X(:,2));
( ?& s  y/ C0 `/ ]+ k& Dtitle('原始数据');
4 U$ u+ Q, P4 D2 v" |; g5 Y%数据参数化,将X(:,1)变换为[0,1]间数据,
minX=min(X(:,1));
maxX=max(X(:,1));
paraX=0;
* I9 o% d  q4 Q+ t2 u( \for i=1:length(X)
2 z* F8 E; K+ C3 h7 K5 ~* ^# E    paraX(i)=(X(i,1)-minX)/(maxX-minX);
end
8 e% M3 |$ _+ m% J5 _" @6 W  ^& Lm=input('请输入Bernsein基的次数m(即m+1阶):');
b=OLS(X,paraX,m);
residu=residual(X,m,paraX,b);
%循环得读取命令
- ]: s! u. ^, E. Wwhile 1
1 R, c5 ^+ N  _    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
5 Z1 a& C6 T8 L$ a+ b! V( r8 Y* O    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
    cmd=input('请按照以上提示选择操作:');

    switch cmd
# Q  D! P2 B! d, B6 ~8 t) g3 O        case 1

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

- `3 y" g; T& L5 e" M7 N: M

宇仲

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