Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
+ m/ V% g9 ^7 ?* A/ x& A& b8 i+ X3 c' e%
  C! Y% Z/ E  v%see also http://www.matlabsky.com
: B6 n* _; J# c. {6 e# @7 ~%
%===得到Bernsein基的次数================
3 [' f- W8 T, j  E8 m7 G" Pclear;
! ^- @9 T3 Y7 {; |, dX=load('data.txt');
%画出原图形的
+ T5 W6 i" W8 g+ x; n7 ^9 P3 Fplot(X(:,1),X(:,2));
, W# d& d8 J" t% Atitle('原始数据');
# b" k+ B+ F5 ~) e# [%数据参数化,将X(:,1)变换为[0,1]间数据,
  y/ @" e$ r  ZminX=min(X(:,1));
maxX=max(X(:,1));
paraX=0;
for i=1:length(X)
8 I* q/ E! f$ `- o1 K5 J' V    paraX(i)=(X(i,1)-minX)/(maxX-minX);
0 G; Z# }, `% e# Zend
m=input('请输入Bernsein基的次数m(即m+1阶):');
b=OLS(X,paraX,m);
6 h& {" L* A  xresidu=residual(X,m,paraX,b);
$ p3 ?$ N' H, w4 f' o. C%循环得读取命令
while 1
    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
* S% c% @9 y8 h" I7 v    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
    cmd=input('请按照以上提示选择操作:');

$ A5 k9 B! D+ _* d" @    switch cmd
        case 1

2 E" `0 F3 |- V; s9 u

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

! X. ^# F' k3 Y  Z

宇仲

楼主辛苦了，继续加油啊！
3 |" \1 J2 \: k5 d- W  g/ i9 }