Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
%
%see also http://www.matlabsky.com
4 H8 B- a# P5 M  a) Y/ l1 i%
%===得到Bernsein基的次数================
. c0 s7 X7 F" `  |+ G6 w7 mclear;
X=load('data.txt');
%画出原图形的
plot(X(:,1),X(:,2));
title('原始数据');
9 m* n7 i2 L; t" ~" j) W2 n%数据参数化,将X(:,1)变换为[0,1]间数据,
minX=min(X(:,1));
+ |; p9 b# l: ~" Q. t. \maxX=max(X(:,1));
: @  V  p4 w/ w+ W& @2 [8 F/ pparaX=0;
for i=1:length(X)
$ q7 W, c$ E# z9 P' F& U- p    paraX(i)=(X(i,1)-minX)/(maxX-minX);
( @0 \6 l9 ?/ Y; _* z' ?, Uend
$ g) _9 C/ o$ w/ _. u: x; vm=input('请输入Bernsein基的次数m(即m+1阶):');
( g5 W; o2 J% J- W5 ab=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 ');
, {6 N( ^+ @: h; o    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
    cmd=input('请按照以上提示选择操作:');

    switch cmd
        case 1
3 W) n1 J% U' q# I( ?/ }

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

宇仲

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