Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
%
%see also http://www.matlabsky.com
%
' x" V0 s+ y& A1 P) \%===得到Bernsein基的次数================
clear;
X=load('data.txt');
%画出原图形的
' F0 a) o) z: ^0 R( ~5 qplot(X(:,1),X(:,2));
  ]% }# j' P9 F: {5 ^& n7 ^title('原始数据');
' S5 [' y: w  i) d; i%数据参数化,将X(:,1)变换为[0,1]间数据,
; }" G2 C+ c# N2 {0 VminX=min(X(:,1));
maxX=max(X(:,1));
* n) p! s5 o# \4 s0 h7 B/ {paraX=0;
for i=1:length(X)
    paraX(i)=(X(i,1)-minX)/(maxX-minX);
+ m6 Z0 @5 `. ?end
m=input('请输入Bernsein基的次数m(即m+1阶):');
b=OLS(X,paraX,m);
7 N" v  K3 u, t4 Xresidu=residual(X,m,paraX,b);
%循环得读取命令
4 O/ b$ ?3 b% k+ E1 V, iwhile 1
2 W& {( N+ r7 S$ K& M1 x2 T    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
6 H' C- E- A$ L& U$ Q    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
# O. k4 O7 K+ j6 i& d* _    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
9 q$ U. w2 ^+ \7 v7 [4 y# n    cmd=input('请按照以上提示选择操作:');

& T! O1 x# ?+ {7 u7 v, [! n- M( s    switch cmd
: Y( e! [: F( O        case 1

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

1 ?7 |: u0 X6 ^, N8 y' o
% y+ [$ {4 t: e5 w: R

宇仲

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