用matlab擬合函式,關係式y a b xc

2021-05-02 11:24:26 字數 4153 閱讀 2639

1樓:匿名使用者

clc;clear all

xy=[3.1,25;6.06,100;7.

99,225;9.52,400;10.81,625;11.

93,900;12.93,1225;13.84,1600;14.

65,2025;15.41,2500;16.14,3025;16.

82,3600;17.47,4225;18.08,4900;18.

65,5625;19.22,6400;19.75,7225;20.

25,8100,;20.73,9025;21.18,10000];

x=xy(:,1);

y=xy(:,2);

plot(x,y)

f=@(a,x) a(1)./(a(2)+x.^(-a(3))); %擬合形式,陣列a即代表[a,b,c]

a0=[1,0.02,-5];  % 假設初始值r=lsqcurvefit(f,a0,x,y);

hold on

plot(x,r(1)./(r(2)+x.^(-r(3))),'r--')

legend('擬合前','擬合後')

disp(r)

計算出:

a、b、c=

0.2082   -0.0000    3.3837

2樓:匿名使用者

[x,i]=sort(x);y=y(i);%由於你的數字是雜亂的,我排列一下

f=@(a,x)a(1).*x.^(a(2))+a(3); %擬合形式,陣列a即代表[a,b,c]

a0=[1,0.02,-5]; % 假設初始值r=lsqcurvefit(f,a0,x,y);

%非線性擬合

%擬合結果儲存在r=[a,b,c]

plot(x,y)

hold on

plot(x,r(1)*x.^r(2)+r(3),'r--')legend('擬合前','擬合後')

disp(r)

求用y=a*exp(b*x)+c的方程用matlab來擬合,求引數值。 25

3樓:我行我素

從圖上可以看出,選用y=a*exp(b*x)+c進行擬合是不對的,應改用polynomial中的

cubic polynomial,我試著擬合得:

linear model poly3:

f(x) = p1*x^3 + p2*x^2 + p3*x + p4

where x is normalized by mean 2007 and std 3.028

coefficients (with 95% confidence bounds):

p1 =  1.903e+006  (-3.354e+007, 3.734e+007)

p2 =  2.888e+007  (5.686e+005, 5.719e+007)

p3 =  7.472e+007  (1.333e+007, 1.361e+008)

p4 =  1.204e+008  (8.642e+007, 1.543e+008)

goodness of fit:

sse: 5.048e+015

r-square: 0.922

adjusted r-square: 0.883

rmse: 2.901e+007

畫出圖如下,

求助:函式模型:y=a*e^-x/m+b*e^-x/n+c,還有x和y的兩組資料,求用matlab擬合出引數a、b、c、m、n的程式 20

4樓:匿名使用者

我這裡用抄x向量代表a,b,c,m,n,樓主可

襲對號入座

%%close all;clear;clc;

f=@(x,xdata)(x(1)*exp(-xdata/x(2))+x(3)*exp(-xdata/x(4))+x(5));

xdata=[1 2 3 4 5];

ydata=[5 4 3 2 1];

x0=[1 1 1 1 1];

x=lsqcurvefit(f,x0,xdata,ydata)

%%close all;clear;clc;

xdata=[1 2 3 4 5];

ydata=[5 4 3 2 1];

f=@(x)(x(1)*exp(-xdata/x(2))+x(3)*exp(-xdata/x(4))+x(5)-ydata);

x0=[1 1 1 1 1];

x=lsqnonlin(f,x0)

兩個**結果都是 5.5394 10.1713 8.3698 10.8647 -7.5688

需要提醒的是,如果擬合的函式選取不當,或者資料不夠,或者資料雜亂無章的話,起點不同結果差別比較大

用matlab擬合函式,關係式y=a/(b+x^(-c)),求引數a b c的值 5

5樓:匿名使用者

xy=[25,3.1;100,6.06;225,7.

99;400,9.52;625,10.81;900,11.

93;1225,12.93;1600,13.84;2025,14.

65;2500,15.41;3025,16.14;3600,16.

82;4225,17.47;4900,18.08;5625,18.

65;6400,19.22;7225,19.75;8100,20.

25;9025,20.73;10000,21.18];

y=xy(:,1);

x=xy(:,2);

plot(x,y)

f=@(a,x) a(1)./(a(2)+x.^(-a(3))); %擬合形式,陣列a即代表[a,b,c]

a0=[0,8,1];  % 假設初始值

r=lsqcurvefit(f,a0,x,y);

hold on

plot(x,r(1)./(r(2)+x.^(-r(3))),'r--')

legend('擬合前','擬合後')

disp(r)

我有2013a的訪問密碼 c27f

6樓:匿名使用者

123456

78910

11[x,i]=sort(x);y=y(i);%由於你的數字是雜亂的,我排列一下

f=@(a,x)a(1).*x.^(a(2))+a(3); %擬合形式,陣列a即代表[a,b,c]

a0=[1,0.02,-5]; % 假設初始值r=lsqcurvefit(f,a0,x,y);

%非線性擬合

%擬合結果儲存在r=[a,b,c]

plot(x,y)

hold on

plot(x,r(1)*x.^r(2)+r(3),'r--')legend('擬合前','擬合後')

disp(r)

7樓:傲慢的

>> xy=[25,3.1;100,6.06;225,7.

99;400,9.52;625,10.81;900,11.

93;1225,12.93;1600,13.84;2025,14.

65;2500,15.41;3025,16.14;3600,16.

82;4225,17.47;4900,18.08;5625,18.

65;6400,19.22;7225,19.75;8100,20.

25;9025,20.73;10000,21.18];

>> x=xy(:,1);

>> y=xy(:,1);

>> syms t;

f=fittype('a/(b+t^(-c))','independent','t','coefficients',);

cfun=fit(x,y,f)

結果:cfun =

general model:

cfun(t) = a/(b+t^(-c))

coefficients (with 95% confidence bounds):

a = 1 (1, 1)

b = -2.647e-014 (-9.633e-014, 4.339e-014)

c = 1 (1, 1)

這裡括號裡的值是置信區間

怎樣從一兩個座標求出函式關係式

通常,要copy先知道什麼函式bai,寫出表示式,如 y ax 2 bx c 若給du出zhi x1,y1 x2,y2 x3,y3 三組數分別代入上式dao y1 ax1 2 bx1 c y2 ax2 2 bx2 c y3 ax3 2 bx3 c 解出a,b,c的具體數值後代入原式就得到所要求的函式...

在二次函式中,知道點的座標,怎樣求函式關係式

令y ax 2 bx c 代入座標 0 4a 2b c 1 0 a b c 2 8 4a 2b c 3 三個未知數,三個方程,這樣可以求出a,b,c 3 1 4b 8 b 2往 1 2 中代入b 2 4a 4 c 0 4 a 2 c 0 5 4 5 3a 6 0 a 2代入 5 c 4 y 2x 2...

三角函式的關係式例題求三角函式大題30道及答案,要簡單點的

同角三角函式的基本關係 式 倒數關係 商的關係 平方關係 tan cot 1 sin csc 1 cos sec 1 sin cos tan sec csc cos sin cot csc sec sin2 cos2 1 1 tan2 sec2 1 cot2 csc2 誘導公式 sin sin co...