Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
%
( x8 `2 _6 N7 @! s%see also http://www.matlabsky.com
%
3 u* g# v9 U3 c! \, T8 G%===得到Bernsein基的次数================
clear;
! l) ]8 [+ {. ^& C7 D3 t$ W# m7 r, Z" }X=load('data.txt');
) T3 B( _1 K$ v0 z( T( t%画出原图形的
plot(X(:,1),X(:,2));
4 h0 r* M' B+ btitle('原始数据');
%数据参数化,将X(:,1)变换为[0,1]间数据,
minX=min(X(:,1));
) M( \% l# D; f0 I2 v7 u: B2 H% S. K1 [maxX=max(X(:,1));
9 y# `* w3 Z2 {, CparaX=0;
for i=1:length(X)
    paraX(i)=(X(i,1)-minX)/(maxX-minX);
7 g/ e/ d; l* j  I4 @end
4 I5 g/ j" E/ S# l% x1 a& Zm=input('请输入Bernsein基的次数m(即m+1阶):');
1 t$ C( t2 k7 I! _  b* Rb=OLS(X,paraX,m);
4 @; F. _/ E" \* o. k1 ^" |residu=residual(X,m,paraX,b);
%循环得读取命令
! v- n# [" c( B1 {. q+ ~4 ?while 1
    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
$ Q7 ^! S7 `6 _7 K0 F; O& M/ ~4 H5 q    cmd=input('请按照以上提示选择操作:');

    switch cmd
1 l- z, C- I2 R% I! f        case 1

( N. |: m7 E4 ?, G1 Y

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

宇仲

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