EDA365电子论坛网

标题: MATLAB线性拟合问题 [打印本页]
作者: elephant_    时间: 2022-8-17 13:46
标题: MATLAB线性拟合问题
怎么把两个线拟合到一个图里，然后显示y轴标题为“电阻值”，x轴标题为“个数”# ?* H0 F7 d2 n! m8 f$ M
1 `. U" B1 R! k0 _  k) [' Z; M! M
x=[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16];. x2 W+ A$ B- K1 b' N- l
y1=[2.23,2.53,4.52,3.01,2.75,3.86,2.25,2.39,0.19,2.62,-4.68,-4.76,-10.75,-7.84,-5.25,-9.96];
  W/ C$ @9 W! _& }! uy2=[2.37,2.53,3.91,3.01,2.75,3.86,2.25,2.39,0.19,2.62,1.07,-4.76,-10.75,-7.84,-5.25,-9.96];$ X6 m2 y) C& K* A

* D, G0 Q3 d2 X: r6 afigure(1);
4 ~' x3 N# B9 T/ {3 w%plot(x,y1,'r',x,y2,'b);
# ~- p. m/ K# n- _plot(x,y,'r');
1 R) J1 L, i0 ]2 f, E8 wtitle('线性拟合')# H5 C; ]8 `0 z; x3 D4 Q1 A0 U
hold on;
+ z. R) _" C( D[a,b]=Linear_fitting(x,y);6 m4 O. u: }' D1 e

7 E" t" n- ?1 kt=1:1:16;
5 l3 J1 v  j2 V. b% T2 O6 M* q9 p' K% lz=a*t+b;3 f1 V2 i7 F! n& D7 E3 a
plot(t,z);
# W* H0 Y  K: u5 b7 K. E2 c
5 Q: x- C1 T( X' T! r- e( Mfunction [a,b]=Linear_fitting(x,y)
/ M: E0 m: S; w( Z* t. t4 c" Yn=size(x,2);
1 R6 D$ p9 R$ dsum_xy=sum(x.*y);
( M1 V4 `8 L& isum_x=sum(x);' k! N  I4 c* _/ S0 z7 C4 \" v
sum_y=sum(y);9 p  l  e& `2 J2 ~' v; F: V9 a- a& |  X
sum_xx=sum(x.*x);! q1 c/ M9 A9 K0 w0 n
# K+ u2 f, M) C" [: D" B7 m: g1 a
a=(n*sum_xy-sum_x*sum_y)/(n*sum_xx-sum_x^2);4 e& v/ j. H; s
b=(sum_xy-sum_xx*sum_y/sum_x)/(sum_x-n*sum_xx/sum_x);' w% O) s- `. V4 f
; o6 o* B. ]/ U2 P7 v
end
- Q( I" }7 z8 ?3 ^' _
作者: purpose_857    时间: 2022-8-17 14:31
能否先给出一个你认为已达成目的的图呢，它长什么样子。
作者: elephant_    时间: 2022-8-17 14:38

" d5 H; A" {3 Z' P3 V2 W
作者: whatever_    时间: 2022-8-17 14:49
你是要 这种效果吧？
拟合4次

9 Z1 Z( U" ?& C7 f. K5 X% z' ^. S4 I! S$ G' o
拟合2次
8 u. i0 |8 Q' b% H: k' Z6 F
" e5 C2 u, ^7 _. {3 o
测试例程
2 b) V! B* c6 \7 w, N

4 t* N* n- o& B* L3 x# `$ H0 ^8 U" x5 ~/ @6 ~




欢迎光临 EDA365电子论坛网 (https://bbs.eda365.com/) Powered by Discuz! X3.2