Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
%
%see also http://www.matlabsky.com
%
%===得到Bernsein基的次数================
clear;
- P- J, u, v/ G6 b% T+ f3 CX=load('data.txt');
%画出原图形的
, D7 r# T& C- v8 aplot(X(:,1),X(:,2));
title('原始数据');
%数据参数化,将X(:,1)变换为[0,1]间数据,
minX=min(X(:,1));
maxX=max(X(:,1));
paraX=0;
+ d  z" v  `2 T$ Z3 v) Sfor i=1:length(X)
4 Z- w7 R+ d2 n' T" g9 t7 o0 D    paraX(i)=(X(i,1)-minX)/(maxX-minX);
  q  U/ e* T+ y. w! h; I' j& ?end
0 T- j# n) \" ~m=input('请输入Bernsein基的次数m(即m+1阶):');
4 U) ?( Y6 ^; [, S1 Ub=OLS(X,paraX,m);
residu=residual(X,m,paraX,b);
%循环得读取命令
while 1
    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
( ~5 P/ q0 A1 y, K2 x' c    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
    cmd=input('请按照以上提示选择操作:');

! z* x* u3 o$ A( g0 b, F    switch cmd
) R2 Y: ^/ ^0 h- f6 Z( Y  k! v        case 1

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

/ [  q$ t' r& I2 t. j7 s! f) L! d

宇仲

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