QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 1991|回复: 1
打印 上一主题 下一主题

[代码资源] Bernstein基构建函数数据MATLAB代码

[复制链接]
字体大小: 正常 放大

2620

主题

162

听众

1万

积分

升级  0%

  • TA的每日心情
    开心
    2015-3-12 15:35
  • 签到天数: 207 天

    [LV.7]常住居民III

    社区QQ达人 发帖功臣 新人进步奖 优秀斑竹奖 金点子奖 原创写作奖 最具活力勋章 助人为乐奖 风雨历程奖

    群组第六届国赛赛前冲刺培

    群组国赛讨论

    群组2014美赛讨论

    群组2014研究生数学建模竞

    群组数学中国试看培训视频

    跳转到指定楼层
    1#
    发表于 2015-1-20 13:50 |只看该作者 |倒序浏览
    |招呼Ta 关注Ta
    %================Bernstein基构建函数数据(二维)==========================% D9 g) o, ~, r9 A1 C, Y0 m
    %4 o5 j! F1 J0 R5 O
    %see also http://www.matlabsky.com
    1 V' S( D2 y# w- |/ R& K%( G0 }/ O0 D2 U8 t/ {8 [
    %===得到Bernsein基的次数================
    " L8 }" F4 ~4 F/ N7 C1 [# c* Mclear;1 }1 O$ C/ L8 d# K+ V* x
    X=load('data.txt');
    4 w6 I6 W- g+ R* t/ F, V2 `, M%画出原图形的7 s" E6 U- O# }
    plot(X(:,1),X(:,2));" S7 g. Z/ Y9 {9 `& `5 M$ D% c
    title('原始数据');5 A( D; t0 B  Q/ @
    %数据参数化,将X(:,1)变换为[0,1]间数据,/ E& U8 q: E' f' t( {8 p
    minX=min(X(:,1));0 m2 |1 f0 B1 b( w! |
    maxX=max(X(:,1));3 v2 F8 q' r: Z& N
    paraX=0;1 H* e: C: I! r6 ~
    for i=1:length(X)
    $ J* q" r1 E, N+ N. A    paraX(i)=(X(i,1)-minX)/(maxX-minX);& V$ {& C  R4 I9 K% u. F8 i0 P
    end9 c: t4 k8 r3 x+ m
    m=input('请输入Bernsein基的次数m(即m+1阶):');2 o; k0 V& L' [1 M4 U
    b=OLS(X,paraX,m);
    $ a# \+ K: k# b) m" _residu=residual(X,m,paraX,b);/ A9 n) a$ Q+ T6 X: j8 y. y+ z
    %循环得读取命令5 V7 B& e/ x. M* t4 }; H
    while 1
    " Y, m; n+ M6 k( {" V" z    fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');
    % G1 _8 Q; ]: b* F4 D0 Y6 Q% N5 L    fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');
    # U% {# z5 B4 w2 V7 {' O    fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');( P6 }! M) D/ F5 ?: o) [5 x" Q
        fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');8 w7 d# z8 D- I5 g
        cmd=input('请按照以上提示选择操作:');7 `) E$ G: }3 Y8 a$ ^; d9 Z
    & f  i! b% y0 T1 @
        switch cmd' u6 O3 N* c; ~$ o
            case 16 p; e. |9 V+ ]6 N' i2 I
    ; a1 `9 \3 r1 B
    游客,如果您要查看本帖隐藏内容请回复
    1 O. @! Y! X, I. \  }0 s' n4 K% F1 y

    ) w3 ?$ q* q4 j. z% V
    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏0 支持支持0 反对反对0 微信微信
    宇仲        

    2

    主题

    13

    听众

    1263

    积分

    升级  26.3%

  • TA的每日心情
    开心
    2015-4-10 07:33
  • 签到天数: 36 天

    [LV.5]常住居民I

    社区QQ达人

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-13 17:35 , Processed in 0.358489 second(s), 56 queries .

    回顶部