Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
8 }  K8 g2 @# L) b# G! g+ T%
%see also http://www.matlabsky.com
; n, }; v" ]+ b( z% Z, f%
* _* a* u, F. F5 h%===得到Bernsein基的次数================
clear;
X=load('data.txt');
( S4 T" m; n* |1 [2 e! `%画出原图形的
+ S* C0 a" F8 `, d% T8 Lplot(X(:,1),X(:,2));
3 Q, s1 r8 A$ x6 m' V) otitle('原始数据');
%数据参数化,将X(:,1)变换为[0,1]间数据,
* n2 J7 ]( J+ V/ I! I. `' YminX=min(X(:,1));
maxX=max(X(:,1));
1 i; M- l/ s/ O6 y& XparaX=0;
for i=1:length(X)
    paraX(i)=(X(i,1)-minX)/(maxX-minX);
3 S+ O* d8 X( u+ P  @9 Send
1 c( u" b/ n/ x) a' G2 e0 ^m=input('请输入Bernsein基的次数m(即m+1阶):');
8 B, M/ B: N$ j0 Hb=OLS(X,paraX,m);
residu=residual(X,m,paraX,b);
+ B) m- F) p; v% _, \8 d%循环得读取命令
  h: u& r# P" w  y& _6 @2 qwhile 1
+ n! V  l* `: e2 O    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
0 K) V* s3 q6 y5 @% ?/ o    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
7 Z3 \3 b* j8 }5 ~    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
/ W4 f% Y# k2 C1 J2 c    cmd=input('请按照以上提示选择操作:');

' f: J0 c: T3 Z% D, T+ e    switch cmd
        case 1

+ U0 w- R0 V. s

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

# H+ ~% ^- t9 N9 Y4 ]6 U

宇仲

楼主辛苦了，继续加油啊！
1 j( F! G) A3 _- O- ~8 n