TA的每日心情 | 开心 2022-1-29 15:07 |
---|
签到天数: 2 天 [LV.1]初来乍到
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
怎么把两个线拟合到一个图里,然后显示y轴标题为“电阻值”,x轴标题为“个数”
; _8 W& Q2 \/ b7 E
$ y( y; u- E" ?: Nx=[1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16];7 d# R: @' i2 h% ~0 ?8 u
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];& C6 w$ |$ M" o& n% ~0 n
y2=[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];
) u: H1 P S9 {/ X: L& ^# _% w$ o W" r
figure(1);. B7 Z" k4 [1 k! v
%plot(x,y1,'r',x,y2,'b);
6 q. N0 `+ c8 z% [plot(x,y,'r');
6 w4 [5 P3 B* H/ W) n' @ Xtitle('线性拟合')# {- b' W6 n+ C
hold on;
4 ~0 W" T) U8 J[a,b]=Linear_fitting(x,y);
9 @( Y. V$ A+ y
7 f1 F% x5 g. }8 ^" w! ]) F# O9 {9 ct=1:1:16;" A" c1 r; C7 b9 ]0 _3 A
z=a*t+b;* j$ h; W$ T P2 k/ i: q) C _
plot(t,z);
+ }5 I! v- Q ?2 D. b8 k$ L: |, m% h1 W
function [a,b]=Linear_fitting(x,y)! O7 d) q2 m s( ] h0 a9 D
n=size(x,2);
8 j9 t: z) P& B- dsum_xy=sum(x.*y);1 H, B7 W5 m) i/ s/ V* w
sum_x=sum(x);
) ?4 \- ]# V& |3 Y7 Z. f2 {sum_y=sum(y);
: i' k, _7 U% A/ Esum_xx=sum(x.*x);
6 A% }* i* y( S2 H9 E. {+ E+ w. ?6 Z2 P
a=(n*sum_xy-sum_x*sum_y)/(n*sum_xx-sum_x^2);
: C/ H3 d7 s: k3 S6 Pb=(sum_xy-sum_xx*sum_y/sum_x)/(sum_x-n*sum_xx/sum_x);
5 l+ r- _# |1 M2 @. H* y
3 D* D. D1 K: s( }* {end
( R2 |4 \! B( W |
|