找回密码
 注册
关于网站域名变更的通知
查看: 374|回复: 1
打印 上一主题 下一主题

fit函数拟合曲线曲面方程的问题请教大神。

[复制链接]
  • TA的每日心情
    开心
    2019-11-29 15:38
  • 签到天数: 4 天

    [LV.2]偶尔看看I

    跳转到指定楼层
    1#
    发表于 2020-5-15 11:23 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式

    EDA365欢迎您登录!

    您需要 登录 才可以下载或查看,没有帐号?注册

    x
    本帖最后由 ccapsemi 于 2020-5-15 13:48 编辑 # _; @% i# n4 D0 F) V0 @7 r7 o
    4 }7 @8 {/ Y% F) {$ t

    ( X" i- A0 I" o, z4 F; ^大家好,我最近在使用fit函数拟合曲线曲面方程,如下) R! }7 ~" `. J4 c8 j
    a = [1,2;3,4;5,6;7,8;9,10]0 s  @6 Y; [, [1 d- h8 `/ b
    5 d" y* u1 R8 z+ Y& P% u' z/ X5 |' W
    a =
    $ J* H8 X: o* F( L8 c0 `; X0 o8 ?$ J, H0 X: e' C! s/ w5 H
         1     2& @* r4 O! [2 c2 ^" r5 E  `
         3     4
    # @2 f/ s0 Y" x2 V* b0 x     5     6
    * R; D1 V7 J9 k3 R5 b) ]     7     8. T* w) c% V; m4 Q, M
         9    10  Z* H+ r3 M) B" _- K( m7 D

    + Z' g5 Z/ v8 {: s. LK>> result = fit(a(:,1),a(:,2),'poly1')( b) v. N' \1 m) C1 o! t- e9 H
    ; i6 k5 @0 X3 s* h7 ?3 s
    result =
    " V' b- M0 N1 M' X" a( h7 H) `! G: J+ h7 F' _( T6 c- _
         Linear model Poly1:
    7 n5 b) e( s* r% B: ^6 R; J: J7 h8 e     result(x) = p1*x + p2
    3 X9 _+ \' c+ S" f# J! Z7 @     Coefficients (with 95% confidence bounds):
    " G. E) h7 o+ p8 Y4 i: m       p1 =           1  (1, 1)
    4 _- ?  A* C3 c4 z, e       p2 =           1  (1, 1)) u/ d! L* l/ J' M1 }. a4 b  h0 `
    , x# K/ C- g: e- v1 y4 @
    虽然fit函数可以拟合出方程,我知道得到的result就相当于f(x,y),但是我不会直接用到方程,我想要的是方程的系数,就是上例中的p1和p2。4 [+ V) p3 o3 b  g  z0 I
    我查过了有关问题,他们是在工具箱显示,我想直接用fit函数,程序不加冒号的时候,他虽然可以显示方程,可以看到系数,但是我需要的是将系数直接输出到一个矩阵或向量里。
    / ]- e7 P& j& K- b7 ^7 F4 ?请教大神,非常感谢!

    该用户从未签到

    2#
    发表于 2020-5-15 13:45 | 只看该作者
    a = [1,2;3,4;5,6;7,8;9,10];! T' ~) c# B0 _, Z" t7 K6 w9 L
    result = fit(a(:,1),a(:,2),'poly1')" I# P; U1 {# y/ g/ Y! |9 |8 g
    coefficients = coeffvalues(result)
    " |, Q8 s7 L: u9 W- Q; v2 ~4 q/ H- w' p: `& s- L: G

    " m( ~( ^+ i5 i% @% Presult =/ \) o+ M% y* F$ \$ U! r5 h* {% u

    . S  ^: B9 O5 [# _- S; `) `/ u     Linear model Poly1:
    ) W8 O: y, s5 T. }% T: q% m# A     result(x) = p1*x + p2
    : f! W$ S: N# |0 e7 `8 [     Coefficients (with 95% confidence bounds):/ s+ W( m3 a/ a- i/ E% x2 d  n
           p1 =           1  (1, 1)
      ?& r- j; ?6 X9 k3 h; n1 T1 ]       p2 =           1  (1, 1)0 ^0 d( B- d6 H) \
    - C/ }8 w, Q% s( z8 t0 v
    coefficients =
    7 I' [4 S6 n7 N" _# [
    ' X6 ~( y  e- m# {: K   1.000000000000000   1.000000000000000
    您需要登录后才可以回帖 登录 | 注册

    本版积分规则

    关闭

    推荐内容上一条 /1 下一条

    EDA365公众号

    关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

    GMT+8, 2025-6-25 20:33 , Processed in 0.062500 second(s), 23 queries , Gzip On.

    深圳市墨知创新科技有限公司

    地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

    快速回复 返回顶部 返回列表