QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2000|回复: 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基构建函数数据(二维)==========================
    : x% |2 ]. C5 O" E%
    # p* |- A3 m" x: f5 m4 |: {%see also http://www.matlabsky.com0 i* f/ D4 V2 N  H
    %( X* i, _. k7 G7 `# Y
    %===得到Bernsein基的次数================1 P( s6 b: ?0 p8 ]
    clear;* k: u* y* v5 _# \
    X=load('data.txt');7 f) \! @* F; ?: _7 {* F: q
    %画出原图形的7 G, Q! a4 d4 r
    plot(X(:,1),X(:,2));6 @& e0 V$ S( H( v
    title('原始数据');/ l+ a9 D8 a5 ~
    %数据参数化,将X(:,1)变换为[0,1]间数据,# e) R) O+ M' j7 X; g8 U- r6 C/ U
    minX=min(X(:,1));: O. k1 _/ S1 \! N$ ^5 _# B
    maxX=max(X(:,1));) p5 W! v% W! T5 T( A8 Y
    paraX=0;' J: f+ I/ H1 e) K9 K  ]  q! o
    for i=1:length(X)
    ' g4 L- Y+ ~' l  ^) O8 [    paraX(i)=(X(i,1)-minX)/(maxX-minX);
    ' }% J4 l4 G% u- A8 gend
    8 D8 n+ F( M! z* [$ [0 Ym=input('请输入Bernsein基的次数m(即m+1阶):');/ U: x+ I" B; ]3 ^. ^
    b=OLS(X,paraX,m);- l! g% L0 n+ A# v, f9 {- T
    residu=residual(X,m,paraX,b);0 t* F9 t# I' G# ^8 b
    %循环得读取命令
    - a) |0 ^4 q/ J: a0 xwhile 19 ^- P  ]# @" M# I! U
        fprintf(' 1.设置基的阶数并拟合\n 2.图形显示原始数据\n 3.图形显示拟合结果\n ');( }0 ?/ ?' F. i) p6 {
        fprintf('4.图形显示残差\n 5.图形显示残差比\n 6.图形显示与平均数残差比\n ');; N2 R& _5 k. L3 X! r2 k0 t
        fprintf('7.图形显示一阶导数\n 8.图形显示二阶导数\n ');$ b9 |5 w! v5 X: D  X" f
        fprintf('9.读取其它数据文件,默认扩展名(.txt),默认文件名(data.txt).\n 输入其它整数退出\n ');5 W" C! a& x. \9 l& T, N0 [
        cmd=input('请按照以上提示选择操作:');
    " v# P$ V* {" ]9 \
    7 D: w) R8 ?  p) W" ^  G    switch cmd
    4 Z  S" \$ F4 _( |+ t        case 1
    - F6 T" e8 Q% I, ]
    0 T; S) W' }  {: B
    游客,如果您要查看本帖隐藏内容请回复

    ! p, {" V6 t; v1 K0 n
    " W+ D# I5 F# ?3 U9 }
    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-21 14:15 , Processed in 0.419120 second(s), 56 queries .

    回顶部