Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
%
! i0 U* U; r& ^" y: W2 y* P%see also http://www.matlabsky.com
%
7 m' f+ R' M% U/ g! Y. m%===得到Bernsein基的次数================
clear;
3 c* L* F6 {) }: G1 U- QX=load('data.txt');
# |5 ~% h$ I% B3 |%画出原图形的
- G, l: S4 b& W9 |plot(X(:,1),X(:,2));
title('原始数据');
$ x* @- g6 }. V1 S/ g$ p, U%数据参数化,将X(:,1)变换为[0,1]间数据,
minX=min(X(:,1));
maxX=max(X(:,1));
" Q6 [, N/ x+ c1 D0 |paraX=0;
6 o1 h8 L% U4 A1 i! H$ c3 cfor i=1:length(X)
    paraX(i)=(X(i,1)-minX)/(maxX-minX);
6 e2 Z- ?. s/ Q0 Tend
m=input('请输入Bernsein基的次数m(即m+1阶):');
b=OLS(X,paraX,m);
residu=residual(X,m,paraX,b);
$ j6 s) I! v* h' A( a%循环得读取命令
while 1
/ E; m" d; n8 x, m' J- L4 O    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
- ]( g: e- N( o1 V; t    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
) Q' U. G$ b' d7 D7 G9 x    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
    cmd=input('请按照以上提示选择操作:');
( Z8 ~) A7 e' s
* r# K% j; ^; @0 K    switch cmd
        case 1

3 c6 d5 g) O8 ~

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

+ g% D$ T" ~8 d& W! i

宇仲

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