数学建模社区-数学中国

标题: Bernstein基构建函数数据MATLAB代码 [打印本页]

---

作者: 建不了的模。    时间: 2015-1-20 13:50
标题: Bernstein基构建函数数据MATLAB代码
%================Bernstein基构建函数数据(二维)==========================
1 t) u3 l; X3 d  Q2 ]  \%
" |: g- {2 O& d% w4 h5 q# m& n%see also http://www.matlabsky.com
%
5 l) ~1 t% g/ V- ?6 Y8 ^4 M! \+ q) c%===得到Bernsein基的次数================
clear;
X=load('data.txt');
- g  _* j8 Q$ t7 Z7 s- x. r%画出原图形的
plot(X(:,1),X(:,2));
title('原始数据');
5 \& a5 `. P0 ]  H; O8 T8 E3 _& j%数据参数化,将X(:,1)变换为[0,1]间数据,
minX=min(X(:,1));
maxX=max(X(:,1));
* I& E: V, l, p1 FparaX=0;
for i=1:length(X)
    paraX(i)=(X(i,1)-minX)/(maxX-minX);
end
$ ^" V2 t2 Y) ^  v( |+ wm=input('请输入Bernsein基的次数m(即m+1阶):');
# z! d$ \0 ~: @4 Pb=OLS(X,paraX,m);
; E2 x5 S' Q# N3 n0 K' X1 W8 Iresidu=residual(X,m,paraX,b);
& ^0 ^0 J7 J/ ?4 u7 l6 V%循环得读取命令
while 1
    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
* k0 S, f6 J" @* K8 v7 p    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');
. }. d! o% Q' z7 `    fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');
    cmd=input('请按照以上提示选择操作:');

* W: T+ X1 m( U0 u1 j3 [  i    switch cmd
        case 1

' I% P+ O) u# z
+ m# a% H9 _1 }8 b! g* h, }* k: C: [

---

作者: 宇仲    时间: 2015-1-20 14:28
楼主辛苦了，继续加油啊！
3 F. T2 ~9 ^6 }/ `$ L0 E: C

---

欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)

Powered by Discuz! X2.5