Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
9 v& \5 W# `! S4 W! `%
. s8 {0 p: l6 E! R; S) h%see also http://www.matlabsky.com
%
%===得到Bernsein基的次数================
clear;
X=load('data.txt');
, J; L+ b* K( s3 l) @- O) D, j4 z" S%画出原图形的
; {. I& r  \* Mplot(X(:,1),X(:,2));
/ P! Y$ e0 O0 J9 Q4 otitle('原始数据');
  Y$ M" p0 ^) i& O%数据参数化,将X(:,1)变换为[0,1]间数据,
minX=min(X(:,1));
maxX=max(X(:,1));
5 C" z4 q9 K  r, ?0 N8 N/ c6 wparaX=0;
for i=1:length(X)
    paraX(i)=(X(i,1)-minX)/(maxX-minX);
end
; t9 E" y# \0 cm=input('请输入Bernsein基的次数m(即m+1阶):');
b=OLS(X,paraX,m);
residu=residual(X,m,paraX,b);
( i! s, p) d' x1 I( u. g%循环得读取命令
8 n+ L4 I1 d: B- L. vwhile 1
    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
7 u/ E( N% n% ~6 L- y7 j( `( W    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
  t$ q8 ?' D! R' c. q# ?    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
4 I% r, e1 }8 F# l% H: d) a* Y    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
    cmd=input('请按照以上提示选择操作:');
! P" R% j- Z0 P6 o% f  j1 k( X
    switch cmd
        case 1
' t! y6 u+ g+ Q6 O2 D( a
% @0 N4 M# r4 @) T0 x# S, T

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

: }' N8 i8 V6 s$ u$ F1 q/ G% {
; [6 L8 O7 _0 ]& }( Z$ P# t

宇仲

楼主辛苦了，继续加油啊！
( s: j/ [0 c0 v, @5 B