matlab程式設計問題,已知函式表示式和資料,求表示式中的係數

時間 2021-09-02 08:26:16

1樓:匿名使用者

t=[0.1,1.1,2.3,3.4,4.5,5.7,6.7]';

y=[0,94.0,84.0,77.0,69.0,62.0,56.0]';

ft_ = fittype('a*exp(-beta*t)*cos(w*t+phi)',...

'dependent',,'independent',,...

'coefficients',);

st=[-100 0.01 -pi 5]

[curve, goodness]= fit(t,y,ft_,'startpoint',st)

figure

plot(curve,'predobs',0.95);

hold on,plot(t,y,'b*')

st =

-100.0000 0.0100 -3.1416 5.0000

curve =

general model:

curve(t) = a*exp(-beta*t)*cos(w*t+phi)

coefficients (with 95% confidence bounds):

a = -143.8 (-219.6, -68.05)

beta = 0.135 (0.01584, 0.2541)

phi = -2.131 (-2.379, -1.883)

w = 5.394 (5.259, 5.529)

goodness =

sse: 352.1786

rsquare: 0.9377

dfe: 3

adjrsquare: 0.8754

rmse: 10.8348

2樓:

我知道一個非線性擬合的函式:lsqcurvefit,不知道對你有沒有幫助,可以擬合出表示式的未知引數

用matlab擬合一組資料,要求得出函式表示式

3樓:匿名使用者

^根據給定的資料bai,我們可以假定

du函式表示式為

y=b1+b2*x+b3*x^zhi2+b4*x^3+b5*x^4;

所以上述函dao數版可以用matlab的regress()多權元線性迴歸分析函式來擬合。實現過程如下:

a=[1.75,0.26;2.

25,0.32;2.5,0.

44;2.875,0.57;2.

686,0.50;2.563,0.

46];

x=a(:,1),y=a(:,2)

x=[ones(length(y),1),x,x.^2,x.^3,x.^4];

b=regress(y,x)

x1=linspace(1.5,3,50);

y1=b(1)+b(2)*x1+b(3)*x1.^2+b(4)*x1.^3+b(5)*x1.^4;

plot(x,y,'k*',x1,y1,'r-')

grid on

xlabel('x'),ylabel('y')

執行結果

4樓:昱婷漫步

從資料描bai

點圖來看,它可以用多項式du來擬合zhi

**:clear all;

clc;

close all;

x=[1.75,2.25,2.5,2.875,2.686,2.563];

y=[0.26,0.32,0.44,0.57,0.50,0.46];

plot(x,y,'r*')

[p,s]=polyfit(x,y,5)  %引數改為1就是線性擬dao合

y1=polyval(p,x)

hold on;

plot(x,y1,'b.')

五次版多項式擬合生成的權函式值與原資料基本重合:

五次多項式,從最高次項的係數到常數項,依次為p的值,所以:

y=-3.8276x^5+ 47.6583x^4 -235.8777x^3+579.6790x^2 -706.4245x+341.2326

5樓:匿名使用者

程式如下:

x=[1.75,2.25,2.5,2.875,2.686,2.563];

y=[0.26,0.32,0.44,0.57,0.50,0.46];

plot(x,y,'g.','markersize',25);%%%有這個是先描點,看出大致是什麼圖,這裡看到像二次曲線的回。答

hold on ;

p3=polyfit(x,y,2);

%%%%%%下面這個是擬合後的p3公式,式子為什麼是這個可以由上面擬合程式得出,放在一個圖裡面看效果

x2=0:0.5:4;

y2=0.1701*x2.^2-0.4965*x2+0.6027;

plot(x2,y2,'b');

6樓:匿名使用者

資料bai不多哈

clc;clear all

x=[1.75, 2.25, 2.5, 2.875, 2.686, 2.563];

y=[0.26,0.32,0.44,0.57,0.50,0.46];

p=polyfit(x,y,3);%四次

du擬合

zhiplot(x,y,'r--')

hold on

y1=polyval(p,x);%得出

dao擬合結

回果plot(x,y1,'b-')

xlabel('x')

ylabel('y')

legend('擬合前

答','擬合後')

7樓:科技數碼答疑

擬合函式表示式是咋麼的,是多項式還是指數,是1次函式還是二次函式??

matlab程式設計問題急求,matlab程式設計問題。。。。。。急求!!!!!

從資料庫中取出資料集合,將ip列從字串轉為long型,然後排序即可 從資料庫查詢就不說了吧,取出ip之後,通過以下方式可以轉化,直接比大小即可 public static long parseip string ip catch return l 程式 microsoft sql server 20...

matlab傳遞函式問題,matlab中simulink傳遞函式的建模問題

給你我的一個我曾經編的超前校正的程式,可能有用 對於函式功能不懂的,可以使用help命令 function gc blead g,kc,dpm dpm 50 相位裕度 角度 包含5度的安全裕量 kc 1 設定開環增益 g tf 48e 3 48 1e 9 1e 6 2e 3 1 生成待校正系統傳遞函...

MATLAB程式設計的一點小問題

化學工程 暫時不會,可以參考一下這個 先回答第二個問題 a randn 5,6 i input 請輸入要輸出的行數i if i 1 i 5 b a i,disp b else b a 5,disp b disp error i與行數不符 end執行結果 a 0.6508 1.2205 0.5377 ...