Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
%
4 U: K' i! K! x# e0 A%see also http://www.matlabsky.com
%
: `, p8 m, w( {8 X) F7 s%===得到Bernsein基的次数================
clear;
X=load('data.txt');
4 k' `9 K0 e$ O- P%画出原图形的
* f" B  z, l1 b& _! n- X& ^" P: bplot(X(:,1),X(:,2));
( `* C# J. x, F6 ktitle('原始数据');
5 j% K5 M9 k$ [7 I" g0 f%数据参数化,将X(:,1)变换为[0,1]间数据,
# V8 O8 v* p+ f- r' O+ u# v+ R+ lminX=min(X(:,1));
maxX=max(X(:,1));
paraX=0;
2 r8 y' b* F" |& C- ^7 yfor i=1:length(X)
% z& Z; N$ ^- J: W    paraX(i)=(X(i,1)-minX)/(maxX-minX);
4 W+ Q: h" e! [+ @8 bend
m=input('请输入Bernsein基的次数m(即m+1阶):');
b=OLS(X,paraX,m);
. Y* N" I. v# |7 ~5 sresidu=residual(X,m,paraX,b);
  V3 w& p% o! r7 I8 q) {* Y%循环得读取命令
while 1
    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
9 P" ?6 w( W- x& Q3 V    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
4 X. y% O% a* k    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
    cmd=input('请按照以上提示选择操作:');
( l* ^1 S* z$ p* v# a$ S
8 z$ J6 W( p# M3 |    switch cmd
- w% t  K+ g* z; v; z4 U        case 1

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

4 b. R; \# t- h9 N& z% I5 {6 [

宇仲

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