Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
%
8 d0 _5 o( g9 q  i4 q: y/ D! ^0 H%see also http://www.matlabsky.com
%
5 B1 V, [- B! o5 |%===得到Bernsein基的次数================
clear;
! r8 [8 e' f, t& g7 u- e0 ZX=load('data.txt');
%画出原图形的
plot(X(:,1),X(:,2));
title('原始数据');
%数据参数化,将X(:,1)变换为[0,1]间数据,
3 y& l* s5 B9 g* e8 TminX=min(X(:,1));
' q% `9 v2 \- d- @maxX=max(X(:,1));
) l" Y! _# ]8 C, yparaX=0;
for i=1:length(X)
    paraX(i)=(X(i,1)-minX)/(maxX-minX);
" a- w+ Z* b7 _( a2 g/ g5 Mend
m=input('请输入Bernsein基的次数m(即m+1阶):');
b=OLS(X,paraX,m);
3 ?( q4 q7 [" Yresidu=residual(X,m,paraX,b);
%循环得读取命令
7 _# o8 B: b. gwhile 1
    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
: g, z. ?+ ~& L  Y$ y    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
    cmd=input('请按照以上提示选择操作:');
# y0 ]2 T, j8 L5 h6 i3 a) {
    switch cmd
4 o+ X8 Z2 M/ Z4 E, ^; f4 |/ {( s, n        case 1

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

6 ^+ o, D+ P/ I  [

宇仲

楼主辛苦了，继续加油啊！
. X9 F$ A) P6 l1 W8 W