TA的每日心情 | 开心 2022-1-29 15:07 |
---|
签到天数: 2 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
怎么把两个线拟合到一个图里,然后显示y轴标题为“电阻值”,x轴标题为“个数”
$ P$ i+ f/ f$ [. y
* _ b9 a" G" w/ M, Vx=[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16];
+ j) d+ \# n! D: E7 Z& ]8 wy1=[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];
# M# K$ P+ M# Z( ]- Py2=[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];7 N8 B+ }, o2 z g! X# y* d- |
, g( f! g' \3 l3 G4 W
figure(1);4 t) D4 p3 u) H2 h9 P$ a
%plot(x,y1,'r',x,y2,'b);! @/ F! u+ H6 J! L
plot(x,y,'r');
0 Z8 G- d2 s gtitle('线性拟合')" j0 X- H+ ?2 Z2 m( A6 ?9 ^
hold on;
; A! y3 W7 k. L# ? s( A$ f[a,b]=Linear_fitting(x,y);8 m( `# x" A3 v/ F- ^
) A4 ]) ]( n9 B% ?t=1:1:16;7 w. w; W2 S* L# o$ U- S: A
z=a*t+b;: E5 v- Q% T5 S
plot(t,z);' L$ p: ~2 U v, a. [! S
3 z' B: K' m }% |function [a,b]=Linear_fitting(x,y)3 }4 n/ ^4 _( g M( @ k
n=size(x,2);- u) o7 y8 Z" y2 w4 _
sum_xy=sum(x.*y);# W2 E h- G5 X2 m! k
sum_x=sum(x);0 s" {2 R6 G7 h+ \% M9 `
sum_y=sum(y);4 |- }: N9 V6 _! T3 D0 o/ k
sum_xx=sum(x.*x);
. T' a6 ]0 w/ j9 b% k& R
7 M' e1 i" S: x7 d5 ]& ga=(n*sum_xy-sum_x*sum_y)/(n*sum_xx-sum_x^2);0 d/ Y8 [; f$ @
b=(sum_xy-sum_xx*sum_y/sum_x)/(sum_x-n*sum_xx/sum_x);
9 Q( i8 Y& g9 }$ I
; }5 P8 @* P: ]' y B; V/ _: [* Eend: D( F0 N& L2 Z2 V) C
|
|