課程設計名稱： 設計二：數(shù)據(jù)擬合 指導教師： 張莉 課程設計時數(shù)： 6 課程設計設備：安裝了Matlab 、C ++軟件的計算機 課程設計日期： 實驗地點： 第五教學樓北902 課程設計目的：

1. 了解最小二乘擬合的原理，掌握用MA TLAB 作最小二乘擬合的方法； 2. 學會利用曲線擬合的方法建立數(shù)學模型。

課程設計準備：

1. 在開始本實驗之前，請回顧相關(guān)內(nèi)容；

2. 需要一臺準備安裝Windows XP Professional操作系統(tǒng)和裝有數(shù)學軟件的計算機。

課程設計內(nèi)容及要求

要求：設計過程必須包括問題的簡要敘述、問題分析、實驗程序及注釋、實驗數(shù)據(jù)及結(jié)果分析和實驗結(jié)

論幾個主要部分。

1. 用切削機床進行金屬品加工時，為了適當?shù)卣{(diào)整機床，需要測定刀具的磨損速度，在一定的時間測量刀具的厚度，得數(shù)據(jù)如表所示，請選用合適的函數(shù)來描述切削時間與刀具厚度的關(guān)系。

首先對數(shù)據(jù)進行分析，畫出離散的點，觀察點近似的曲線：

t=0:1:15;

y=[30.0 29.1 29.8 28.1 28.0 27.7 27.5 27.2 27.0 26.8 26.5 26.3 26.1 25.7 25.3 24.8]; plot(t,y,'r*')

判斷出曲線是近似直線函數(shù)，所以對數(shù)據(jù)進行測試可以做三次函數(shù)擬合：

t=0:1:15;

y=[30.0 29.1 29.8 28.1 28.0 27.7 27.5 27.2 27.0 26.8 26.5 26.3 26.1 25.7 25.3 24.8]; %plot(t,y,'r*') A=polyfit(t,y,3) z=polyval(A,t); plot(t,y,'r*',t,z, 'b' )

051015

擬合結(jié)果： A =

-0.3099 29.5676

擬合函數(shù)為：y=-0.3099t+29.5676

2. Malthus 人口指數(shù)增長模型

用以上數(shù)據(jù)檢驗馬爾薩斯人口指數(shù)增長模型，根據(jù)檢驗結(jié)果進一步討論馬爾薩斯人口模型的改進。 首先對數(shù)據(jù)進行分析，畫出離散的點，觀察點近似的曲線：

x=1790:10:2000;

y=[3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.4 38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 251.4 281.4]; plot(x,y,'r*')

[**************]

判斷出曲線是近似拋物線函數(shù)，所以對數(shù)據(jù)進行二次函數(shù)擬合：

x=1:10:220;

y=[3.9 5.3 7.2 9.6 12.9 17.1 23.2 31.4 38.6 50.2 62.9 76.0 92.0 106.5 123.2 131.7 150.7 179.3 204.0 226.5 251.4 281.4]; %plot(x,y,'r*')

A=polyfit(x,y,2)； z=polyval(A,x); plot(x,y,'r*',x,z, 'b' ) y=poly2str(A,'x' )

[**************]

擬合函數(shù)結(jié)果： y =

0.0067365 x^2 - 0.13636 x + 6.4212

3．價格預測

美國舊車價格的調(diào)查數(shù)據(jù)

分析用什么形式的曲線來擬合數(shù)據(jù)，并預測使用4、5年后的舊車平均價格大致為多少。 首先對數(shù)據(jù)進行分析，畫出離散的點，觀察點近似的曲線：

x=1:1:10;

y=[2615 1943 1494 1087 765 538 484 290 226 204];

plot(x,y,'r*')

判斷出曲線是近似多次函數(shù)，所以對數(shù)據(jù)進行三次函數(shù)擬合：

x=1:1:10;

y=[2615 1943 1494 1087 765 538 484 290 226 204]; %plot(x,y,'r*') A=polyfit(x,y,3); z=polyval(A,x); plot(x,y,'r*',x,z, 'b' ) y=poly2str(A,'x' )

擬合結(jié)果：

y =

-2.6548 x^3 + 79.9522 x^2 - 852.8324 x + 3380.1 預測4、5年后的舊車價格為：

>>y= -2.6548 *4^3 + 79.9522* 4^2 - 852.8324* 4 + 3380.1 y =

1.0781e+003

>> y= -2.6548 *5^3 + 79.9522* 5^2 - 852.8324* 5 + 3380.1 y = 782.8930

4．用最小二乘法求一個形如y =a +bx 2的經(jīng)驗公式，數(shù)據(jù)如下：

x=[19,25,31,38,44];

y=[19.0,32.3,49.0,73.3,97.8];

fun1=inline('c(1)+c(2)*x.^2', 'c' , 'x' ); %擬合函數(shù) k=[0 0];

c=lsqcurvefit(fun1,k,x,y) %求擬合系數(shù)

['y=',num2str(c(1)),'+',num2str(c(2)),'x^2']%擬合函數(shù)表達式 y1=c(1)+c(2)*x.^2; %擬合函數(shù)表達式值

plot (x,y,'b*', x, y1,'r-' ); %繪制x ，y 點與擬合函數(shù)，比較擬合程度

10090

[***********]

[1**********]5

實驗結(jié)果如下：

Optimization terminated: first-order optimality less than OPTIONS.TolFun, and no negative/zero curvature detected in trust region model. c =

0.9726 0.0500

ans =

y=0.97258+0.050035x^2

課程設計總結(jié)： 這次實驗主要講的是數(shù)據(jù)的擬合，對數(shù)據(jù)的處理首先是觀察圖像在近似擬合，有的函數(shù)掌握的還不夠徹底，需要查閱資料才能完成，以后遇到這種問題會更加仔

細，有些錯誤的地方還希望老師指出，再加以改進，謝謝！

總結(jié)

以上是生活随笔為你收集整理的matlab数据拟合语句,Matlab数据拟合程序 - 范文中心的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯，歡迎將生活随笔推薦給好友。