Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
6 p: r& @3 t! o0 p7 G5 ~: R%
) L5 [% |& e- M" \%see also http://www.matlabsky.com
%
( D, C/ G* A  e: v9 x%===得到Bernsein基的次数================
clear;
% @5 Y9 R$ T$ Z8 a% y. JX=load('data.txt');
%画出原图形的
plot(X(:,1),X(:,2));
$ y  ~6 W& F" R. \- D* q' P/ ctitle('原始数据');
%数据参数化,将X(:,1)变换为[0,1]间数据,
& ?! a: d! ]5 H* dminX=min(X(:,1));
maxX=max(X(:,1));
paraX=0;
7 u$ v/ R( N$ M1 [2 Wfor i=1:length(X)
    paraX(i)=(X(i,1)-minX)/(maxX-minX);
% ?6 l( _$ d0 e) b0 @! n% ]end
m=input('请输入Bernsein基的次数m(即m+1阶):');
4 K, b1 j7 H9 n# [; tb=OLS(X,paraX,m);
residu=residual(X,m,paraX,b);
8 I7 j8 ]8 P+ }) f+ o4 v%循环得读取命令
* @4 p5 Y$ Q+ {8 i* i) h  z5 |/ A8 ewhile 1
    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
! S; D) x% ?4 V% g) f7 g    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
+ U; x* P. z! {/ @/ T7 f, j    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
    cmd=input('请按照以上提示选择操作:');
$ Z/ U8 U* {2 `! c# `* E5 O1 X5 D" n
    switch cmd
9 o, e7 y3 U' d2 w8 {* Y        case 1

7 U' v0 T. g/ }* k( G

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

$ W8 }$ }7 M( m  W% i( P$ [

宇仲

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