Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
%
%see also http://www.matlabsky.com
%
; o1 C( d# x% P' j0 Y%===得到Bernsein基的次数================
clear;
X=load('data.txt');
: ?, O2 k8 j9 N%画出原图形的
/ F2 o" u$ r0 j0 X" Xplot(X(:,1),X(:,2));
title('原始数据');
%数据参数化,将X(:,1)变换为[0,1]间数据,
minX=min(X(:,1));
maxX=max(X(:,1));
% `. X1 ^+ E- g' N7 s& ^paraX=0;
for i=1:length(X)
    paraX(i)=(X(i,1)-minX)/(maxX-minX);
end
5 e$ y- b. ?' J7 u$ f8 x% Im=input('请输入Bernsein基的次数m(即m+1阶):');
. `6 D" B+ Z# i4 `b=OLS(X,paraX,m);
residu=residual(X,m,paraX,b);
  [8 g& {& m! I7 x3 [: M+ c%循环得读取命令
while 1
* S3 l0 ~' T+ l; ~# t    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
* @4 O. x0 P0 U$ Y$ h    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
+ i, Q0 _9 H) }; c: q) v* m    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
* r# S. ?! E2 i# P    cmd=input('请按照以上提示选择操作:');

' t2 c* Q3 L2 I' f    switch cmd
        case 1

" ?. O; c! c) |* C$ `$ [

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

宇仲

楼主辛苦了，继续加油啊！
& p3 `: E. t; z4 b0 ^) U; w) s. E3 s