Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
%
2 J9 T5 U/ `  G2 O; ^%see also http://www.matlabsky.com
1 A' J4 b- k3 V8 n( k: X%
5 c* \% ?% a, M" h7 y%===得到Bernsein基的次数================
) o5 X% e& o# t' n/ g( Sclear;
+ X; |  F& K! J0 l- t. T# dX=load('data.txt');
%画出原图形的
; Q2 v. Q( i' E7 xplot(X(:,1),X(:,2));
title('原始数据');
2 W9 I; I  `; Z- T- F; r# e%数据参数化,将X(:,1)变换为[0,1]间数据,
minX=min(X(:,1));
; f, p7 n" }8 smaxX=max(X(:,1));
" n& ^( h3 n9 jparaX=0;
for i=1:length(X)
    paraX(i)=(X(i,1)-minX)/(maxX-minX);
9 n$ F2 ^: J1 G, ], [% Pend
m=input('请输入Bernsein基的次数m(即m+1阶):');
b=OLS(X,paraX,m);
residu=residual(X,m,paraX,b);
! F7 b! `6 B$ }$ [! m%循环得读取命令
while 1
    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
# H# G5 M3 H( E2 [) X+ ?5 ~! H    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
5 d, ?& e8 l% J$ W    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
/ p& r4 g( P1 A1 P& N1 {7 S    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
: D( [+ a' J+ g" @' L    cmd=input('请按照以上提示选择操作:');
! R6 q( k9 ]& A$ A! s" \
; }0 m4 q! G1 E/ B    switch cmd
4 v, s  f& w0 y% R( i- }& U! B        case 1

! O' ~+ |+ a+ ~4 J( Z' K0 t

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

9 i  y3 z  y2 n5 [% j" A

宇仲

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