Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
! }: ~. Y/ f+ t4 _$ o3 ~" m%
/ o/ s8 ]$ N$ B& `- p%see also http://www.matlabsky.com
%
3 J- t# M! S9 M% j3 b. p$ P+ W%===得到Bernsein基的次数================
clear;
X=load('data.txt');
$ E0 b* q- f" f3 Y5 f$ B%画出原图形的
plot(X(:,1),X(:,2));
' ^- K0 \2 W1 v% [  a' ititle('原始数据');
- J+ O5 i4 }, m0 B5 `- u%数据参数化,将X(:,1)变换为[0,1]间数据,
$ z# E6 c+ d# n5 Y7 I" B5 g- }minX=min(X(:,1));
maxX=max(X(:,1));
" t1 h  C8 j! C# c, t% JparaX=0;
. n/ L7 W+ q8 E7 h' Cfor i=1:length(X)
( s, q% j- d% j  v    paraX(i)=(X(i,1)-minX)/(maxX-minX);
end
m=input('请输入Bernsein基的次数m(即m+1阶):');
$ _1 s; Z6 c0 q! i" m" Ub=OLS(X,paraX,m);
residu=residual(X,m,paraX,b);
4 N1 s+ }2 I; J%循环得读取命令
; _! g$ u! \! U: x" Ywhile 1
    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
! e  K4 J8 v5 e% u5 U% E    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
# |7 n# c9 c5 F$ J    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
1 n1 ?4 c6 w3 u# P$ s    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
( d# @* f  K5 E2 X. T    cmd=input('请按照以上提示选择操作:');

9 b/ @# |! T  Z1 o6 H    switch cmd
6 O  ?0 n  t" p. f+ p        case 1

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

宇仲

楼主辛苦了，继续加油啊！
6 C2 n( D1 i5 X0 P4 ~+ q7 V