Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
%
%see also http://www.matlabsky.com
%
%===得到Bernsein基的次数================
clear;
; C5 s. O. e0 y0 GX=load('data.txt');
* c, U. E% ~- X- K) s6 k%画出原图形的
plot(X(:,1),X(:,2));
7 p9 }7 f, A/ m/ M0 `title('原始数据');
%数据参数化,将X(:,1)变换为[0,1]间数据,
minX=min(X(:,1));
2 w4 U9 O: g6 {0 ^, a4 U; }! CmaxX=max(X(:,1));
paraX=0;
for i=1:length(X)
2 I' d8 J& X) f0 K4 w" X5 u    paraX(i)=(X(i,1)-minX)/(maxX-minX);
0 N8 g* p* @1 ?& @+ nend
m=input('请输入Bernsein基的次数m(即m+1阶):');
b=OLS(X,paraX,m);
+ [0 t! A8 D) Q/ G( Vresidu=residual(X,m,paraX,b);
%循环得读取命令
while 1
    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
7 v2 s- m! ^0 R. ^) O, w% \    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
0 k0 @) Y5 ^- a- Q    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
6 W0 J1 W6 N) k# Q8 U. c    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
    cmd=input('请按照以上提示选择操作:');

    switch cmd
- y. _8 |& P! _: y8 u% [, o( [        case 1
1 z+ ?. u# a3 B

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

" d+ \2 M+ x! `/ n

宇仲

楼主辛苦了，继续加油啊！
: o! |! C8 m. H. P  R: o