版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

返回列表

楊小胖

金蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 1007
散金: 100
紅花: 6
帖子: 936
在線: 379小時(shí)
蟲號(hào): 454210
注冊(cè): 2007-11-08
性別: MM
專業(yè): 測(cè)量與地圖學(xué)

[求助] matlab曲線擬合 basic fitting 與polyfit

想請(qǐng)高手幫忙改個(gè)代碼。
其實(shí)基本實(shí)現(xiàn)了，但怎么越改越亂。
問(wèn)題
<1>我不知道怎么在圖中把擬合的公式顯示出來(lái)。
<2>當(dāng)我只用一組數(shù)據(jù)時(shí)直接在命令窗口輸入，那如果我有很多組這樣的數(shù)據(jù)，那要怎么改呢？
x=[ 10.5 13 13.8 16 18 19 20.8];
y=[1190 1040 1000 800 600 430 200];
<3>為什么這樣擬合出來(lái)的結(jié)果是曲線，但如果通過(guò)普通的polyfit擬合，出來(lái)是折線？

想要實(shí)現(xiàn)的是設(shè)置不一樣的擬合次數(shù)，將擬合的系數(shù)都顯示出來(lái)，并在每一個(gè)擬合系數(shù)求出來(lái)以后把表達(dá)式顯示出來(lái)。

這是我從matlab寶典一書中找到的代碼。

CODE:

function [th,err,yi]=polyfits(x,y,N,xi,r)

%x,y  :數(shù)據(jù)點(diǎn)系列

%N    :多項(xiàng)式擬合的系統(tǒng)

%r    :加權(quán)系數(shù)的逆矩陣

M=length(x);

x=x(:);

y=y(:);

%判斷調(diào)用函數(shù)的格式，考慮權(quán)陣

if nargin==4

%當(dāng)用戶調(diào)用函數(shù)格式為(x,y,N,r)

    if length(xi)==M

        r=xi;

        xi=x;

%當(dāng)用戶調(diào)用函數(shù)格式為(x,y,N, xi)，權(quán)陣為1

    else r=1;

    end

%當(dāng)用戶調(diào)用函數(shù)格式為(x,y,N)，不考慮權(quán)陣

elseif nargin==3

    xi=x;

    r=1;

end

%求解系數(shù)矩陣

A(:,N+1)=ones(M,1);

for n=N:-1:1

    A(:,n)=A(:,n+1).*x;

end

if length(r)==M

    for m=1:M

        A(m,:)=A(m,:)/r(m);

        y(m)=y(m)/r(m);

    end

end

%計(jì)算擬合系數(shù)

th=(A\y)';

ye=polyval(th,x);

err=norm(y-ye)/norm(y);

yi=polyval(th,xi);

CODE:

clear;

clf;

x=input('請(qǐng)輸入x值:');

y=input('請(qǐng)輸入y值:');

% x=[ 10.5  13  13.8  16  18  19 20.8];

% y=[1190  1040 1000 800 600 430 200];

% eb=[0.2*ones(5,1);ones(5,1)];

[x,i]=sort(x);

y=y(i);

xi=min(x)+[0:100]/100*(max(x)-min(x));

 for i=1:4

 N=2*i-1;

[th,err1,yi]=polyfits(x,y,N,xi);

  subplot(220+i);

plot(x,y,'k*')

hold on 

plot(xi,yi,'g:','LineWidth',1.5)

title(['The ',num2str(N),'th Polynomial Curve Fitting'])

grid on

 end

用的是這組數(shù)據(jù)
x=[ 10.5 13 13.8 16 18 19 20.8];
y=[1190 1040 1000 800 600 430 200];

運(yùn)行結(jié)果

[ Last edited by 楊小胖 on 2011-6-15 at 16:26 ]

回復(fù)此樓

» 猜你喜歡

求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過(guò) 已經(jīng)有9人回復(fù)
289材料與化工（085600）B區(qū)求調(diào)劑已經(jīng)有4人回復(fù)
一志愿武理材料工程348求調(diào)劑已經(jīng)有6人回復(fù)
08工科 320總分求調(diào)劑已經(jīng)有11人回復(fù)
一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑已經(jīng)有6人回復(fù)
323求調(diào)劑已經(jīng)有5人回復(fù)
284求調(diào)劑已經(jīng)有5人回復(fù)
317求調(diào)劑已經(jīng)有16人回復(fù)
287求調(diào)劑已經(jīng)有9人回復(fù)
求調(diào)劑已經(jīng)有6人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

matlab 已經(jīng)有5人回復(fù)
-COOH與-NH2的縮合反應(yīng)可以只用EDC，不加NHS嗎？求這類文獻(xiàn)，謝謝！已經(jīng)有3人回復(fù)
魚腥草提取物與魚腥草素鈉已經(jīng)有6人回復(fù)
關(guān)于Matlab作圖的問(wèn)題已經(jīng)有14人回復(fù)
DL_POLY報(bào)錯(cuò)，自由度太低已經(jīng)有10人回復(fù)
matlab矩陣數(shù)據(jù)分塊已經(jīng)有3人回復(fù)
【【吐血求助】】Matlab微分方程作圖-報(bào)錯(cuò)報(bào)錯(cuò) 已經(jīng)有4人回復(fù)
請(qǐng)幫忙用matlab求矩陣特征值已經(jīng)有8人回復(fù)
matlab曲線擬合的一個(gè)小問(wèn)題已經(jīng)有7人回復(fù)
急切求助用matlab曲線擬合色散方程已經(jīng)有10人回復(fù)
matlab工具箱中curve fitting中的圖如何添加橫縱坐標(biāo)軸的名稱已經(jīng)有5人回復(fù)
【討論】matlab量化程序已經(jīng)有4人回復(fù)
【求助】polyfit的擬合效果不好，應(yīng)該怎么修改？【已解決】已經(jīng)有9人回復(fù)
【求助】最小二乘法編寫指數(shù)擬合方程已經(jīng)有19人回復(fù)

人生中最輝煌的不是功成名就的時(shí)候，而是在失敗和挫折中看到希望并為之奮斗的日子

1樓 2011-06-15 16:15:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

楊小胖

金蟲 (正式寫手)

應(yīng)助: 2 (幼兒園)
金幣: 1007
散金: 100
紅花: 6
帖子: 936
在線: 379小時(shí)
蟲號(hào): 454210
注冊(cè): 2007-11-08
性別: MM
專業(yè): 測(cè)量與地圖學(xué)

自己頂貼

回復(fù)此樓

人生中最輝煌的不是功成名就的時(shí)候，而是在失敗和挫折中看到希望并為之奮斗的日子

2樓2011-06-15 19:06:31

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主楊小胖的主題更新

返回列表

普通表情龍兔虎貓

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 一志愿武理材料工程348求調(diào)劑 +4	￣^￣゜汗 2026-03-19	6/300	2026-03-22 17:44 by luoyongfeng
[考研] 08工科 320總分求調(diào)劑 +11	梨花珞晚風(fēng) 2026-03-17	11/550	2026-03-22 17:42 by luoyongfeng
[考研] 環(huán)境學(xué)碩288求調(diào)劑 +6	皮皮皮123456 2026-03-22	6/300	2026-03-22 16:52 by i_cooler
[考研] 275求調(diào)劑 +6	shansx 2026-03-22	8/400	2026-03-22 15:27 by barlinike
[考研] 招08考數(shù)學(xué) +3	laoshidan 2026-03-20	7/350	2026-03-22 12:54 by 在風(fēng)落中
[考研] 286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！ +4	lemonzzn 2026-03-17	8/400	2026-03-21 22:49 by lemonzzn
[考研] 一志愿東華大學(xué)控制學(xué)碩320求調(diào)劑 +3	Grand777 2026-03-21	3/150	2026-03-21 19:23 by 簡(jiǎn)之-
[考研] 0805 316求調(diào)劑 +3	大雪深藏 2026-03-18	3/150	2026-03-21 18:55 by 學(xué)員8dgXkO
[考研] 311求調(diào)劑 +3	勇敢的小吳 2026-03-20	3/150	2026-03-21 17:40 by ColorlessPI
[考研] 279分求調(diào)劑一志愿211 +14	chaojifeixia 2026-03-19	15/750	2026-03-21 13:24 by zhukairuo
[考研] 一志愿山大07化學(xué) 332分四六級(jí)已過(guò) 本科山東雙非求調(diào)劑！ +3	不想理你 2026-03-16	3/150	2026-03-21 03:59 by JourneyLucky
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3	kkcoco25 2026-03-16	4/200	2026-03-21 03:52 by JourneyLucky
[考研] 265求調(diào)劑 +3	Jack?k?y 2026-03-17	3/150	2026-03-21 03:17 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] 一志愿西南交通專碩材料355 本科雙非求調(diào)劑 +5	西南交通專材355 2026-03-19	5/250	2026-03-20 21:10 by JourneyLucky
[考研] 0817 化學(xué)工程 299分求調(diào)劑有科研經(jīng)歷有二區(qū)文章 +22	rare12345 2026-03-18	22/1100	2026-03-20 20:39 by zhukairuo
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考博] 26申博 +4	八6八68 2026-03-16	4/200	2026-03-17 13:00 by 輕松不少隨
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产 高清 中文字幕,99re热久久亚洲综合精品成人,熟妇 一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

24小時(shí)熱門版塊排行榜

楊小胖

[求助] matlab曲線擬合 basic fitting 與polyfit

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

楊小胖

亭亭五月天在线观看,亭亭五月天在线观看,国产最新av一区二区,国产高清中文字幕,99re热久久亚洲综合精品成人,熟妇一区二区三区,一级做a爰片性色毛片武则天,美女的骚穴视频播放,国产美女午夜免费视频

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: