Bernstein基构建函数数据MATLAB代码

建不了的模。

%================Bernstein基构建函数数据(二维)==========================
0 Q$ Y0 `  q8 M2 Y/ w, {+ y%
%see also http://www.matlabsky.com
. e% e: ~1 S% b5 C& o. Z, d%
* V) C8 |# |8 S%===得到Bernsein基的次数================
clear;
X=load('data.txt');
# T" v& C: y' D%画出原图形的
% l. H0 ~0 W5 W+ ?  T7 m7 Xplot(X(:,1),X(:,2));
title('原始数据');
%数据参数化,将X(:,1)变换为[0,1]间数据,
" ~4 _( f4 Z6 t! Z& b# rminX=min(X(:,1));
0 h+ o; ~' M. x# ImaxX=max(X(:,1));
paraX=0;
for i=1:length(X)
* L( M4 G0 ?6 y" E  ^3 p    paraX(i)=(X(i,1)-minX)/(maxX-minX);
end
m=input('请输入Bernsein基的次数m(即m+1阶):');
- o/ V- o) D7 O  ]* D* Fb=OLS(X,paraX,m);
residu=residual(X,m,paraX,b);
%循环得读取命令
while 1
7 Q- U3 N7 O, E( y' P    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
5 M  c# r, A4 F. X    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
" I; W9 {. ]+ E2 W5 y% u9 ?1 ]    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
- N+ I% s% b( ]0 A! N- s( o    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
1 v- B" W5 {0 @! H    cmd=input('请按照以上提示选择操作:');
4 |* x+ B! i. f- M
    switch cmd
5 ?( c7 ]9 @: R+ b4 O" m1 Q& \        case 1
/ h7 i& e& I# s6 k
& B7 ?  s1 u/ Y( F9 O- k% F

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

宇仲

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