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 ...