EDA365电子论坛网
标题:
fit函数拟合曲线曲面方程的问题请教大神。
[打印本页]
作者:
ccapsemi
时间:
2020-5-15 11:23
标题:
fit函数拟合曲线曲面方程的问题请教大神。
本帖最后由 ccapsemi 于 2020-5-15 13:48 编辑
$ }' w/ H2 i+ j) b8 y) [
" z# ?- ~. ~# h- j
: t( l, p; ]/ M; {
大家好，我最近在使用fit函数拟合曲线曲面方程，如下
: R( D7 j4 ]$ S" q- h
a = [1,2;3,4;5,6;7,8;9,10]
) Z q. X, n0 g4 E% I9 |$ C6 Z
. x$ y$ p* P% S9 H8 J+ W
a =
# i; T7 W/ f: B/ K
0 J/ B* C6 V3 ]: n* N
1 2
! k( ?. L0 ^3 `. I
3 4
4 W; N; B' H/ [% m4 D( v1 H G
5 6
6 Q0 b* U% t: _$ `
7 8
' K6 q- I. Y7 n6 R$ x2 \
9 10
$ L5 Q3 g/ C q& ?" a$ n: g6 u$ Y
`% E4 ]+ w3 N$ S
K>> result = fit(a(:,1),a(:,2),'poly1')
) U" T# `( c' X4 o8 S
, u$ ^$ S+ H& O
result =
/ \( s0 B: C* q) G( S
7 I% O) [" T3 ?) Z) U7 K
Linear model Poly1:
- \+ o) f6 K8 P5 D3 R# Y, G
result(x) = p1*x + p2
9 Z! J- P6 a8 E8 h+ M% O6 f* U
Coefficients (with 95% confidence bounds):
6 _6 j5 L7 M M+ G1 e4 p* z) N
p1 = 1 (1, 1)
u1 ]" j& J7 E' P0 r6 I2 g
p2 = 1 (1, 1)
: ~, C5 U% N! M0 `: a8 t, U
4 R H$ X: a( {* K
虽然fit函数可以拟合出方程，我知道得到的result就相当于f（x，y），但是我不会直接用到方程，我想要的是方程的系数，就是上例中的p1和p2。
4 y# Y+ j0 a' \- j: v" e _/ x
我查过了有关问题，他们是在工具箱显示，我想直接用fit函数，程序不加冒号的时候，他虽然可以显示方程，可以看到系数，但是我需要的是将系数直接输出到一个矩阵或向量里。
: P# g; w1 W' G5 `4 g' O. \
请教大神，非常感谢！
作者:
YNhhU
时间:
2020-5-15 13:45
a = [1,2;3,4;5,6;7,8;9,10];
{9 b/ m, z; v& B; d Y
result = fit(a(:,1),a(:,2),'poly1')
5 n8 a3 U7 S" M- D) d
coefficients = coeffvalues(result)
+ T8 Y6 }: m, s! y6 E- Q
$ N; I ]- j( f1 g& v" S' k/ A
. H! T; ^7 x1 r R' N E
result =
/ B% j; a2 j- x c& \
# Q3 k# l) \9 K; z8 Q
Linear model Poly1:
. t( |8 m: b; `9 V9 N0 r- @. F+ M' T
result(x) = p1*x + p2
2 r( T% F5 e+ K; n7 ^. ~
Coefficients (with 95% confidence bounds):
. K4 ^% u9 j g4 B7 X
p1 = 1 (1, 1)
$ U7 a- ^) Z# E! u5 w
p2 = 1 (1, 1)
7 Y( H5 J( i! f U
" G _: W, g+ b1 i
coefficients =
( ?' R, L) N/ }/ H0 M
9 c5 @. t* h2 s5 U8 s) i2 U3 W
1.000000000000000 1.000000000000000
欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/)
Powered by Discuz! X3.2