用MATLAB進行線性擬合的程序代碼

U=[1.45,1.85,2.02,2.2,2.67,2.98];

V=[200,400,500,600,800,1000];

H=[25.6,43.5,51.4,60.2,79.6,95.2];

[AX,H1,H2]=plotyy(U,V,U,H,'plot');

set(H1,'Marker','o','line','none','color','k')

set(H2,'Marker','*','line','none','color','r')

legend('容積','實際高度')

xlabel('電壓值/V')

yl1=get(AX(1),'Ylabel')

set(yl1,'string','容積/L')

yl2=get(AX(2),'Ylabel')

set(yl2,'string','實際高度/cm')

p1=polyfit(U,V,1);

p2=polyfit(U,H,1);

VV=p1(1)*U+p1(2);

HH=p2(1)*U+p2(2);

hold?on

[AX,HH1,HH2]=plotyy(U,VV,U,HH)

set(HH1,'color','k')

set(HH2,'color','r')

clc

disp(['V=',num2str(p1(1)),'*U',num2str(p1(2))])

disp(['H=',num2str(p2(1)),'*U',num2str(p2(2))])

運行結果

V=513.0775*U-542.8718

H=45.1123*U-39.7716