版塊導(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í)熱門(mén)版塊排行榜

返回列表

下一頁(yè)

本帖產(chǎn)生 1 個(gè) 程序強(qiáng)帖，點(diǎn)擊這里進(jìn)行查看

yangger

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1.1
紅花: 1
帖子: 20
在線(xiàn): 17.2小時(shí)
蟲(chóng)號(hào): 1226591
注冊(cè): 2011-03-09
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[求助] 急切求助用matlab曲線(xiàn)擬合色散方程

現(xiàn)在有如下色散方程：

y^2 = A + B/(x^2-C) - D x^2 ----------------Eq. 1

其中A，B, C, D為待定系數(shù)，現(xiàn)在有(x,y)的一組實(shí)驗(yàn)值；如何得到最佳的A、B、C、D

我的思路如下：

step1 通過(guò)四對(duì)(x_i,y_i) 得到 A，B，C，D 值，
step2 代入Eq.1 輸入x實(shí)驗(yàn)值，得到y(tǒng)擬合值，
step3 判斷y擬合值和y實(shí)驗(yàn)值差值，
step4 如果差值大于10^-4則調(diào)節(jié)A,B,C,D;
step5 用新的ABCD，返回到step2

現(xiàn)在遇到的問(wèn)題是 1.不收斂
2. 調(diào)節(jié)過(guò)的A，B，C，D有時(shí)候使得Eq.1方程右邊小于零

求助：要是用matlab如何處理，

[ Last edited by yangger on 2011-6-3 at 12:53 ]

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

liuweihong1980

wrve

» 猜你喜歡

298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑已經(jīng)有5人回復(fù)
297求調(diào)劑已經(jīng)有5人回復(fù)
302求調(diào)劑已經(jīng)有10人回復(fù)
266求調(diào)劑已經(jīng)有8人回復(fù)
化學(xué)工程321分求調(diào)劑已經(jīng)有17人回復(fù)
344求調(diào)劑已經(jīng)有5人回復(fù)
312求調(diào)劑已經(jīng)有9人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專(zhuān)業(yè)（081702）315分求調(diào)劑已經(jīng)有9人回復(fù)
288求調(diào)劑，一志愿華南理工大學(xué)071005 已經(jīng)有4人回復(fù)
304求調(diào)劑已經(jīng)有13人回復(fù)

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

高金求助matlab解微分方程組已經(jīng)有12人回復(fù)
【求助】matlab 求解微分方程中的未知參數(shù) 已經(jīng)有20人回復(fù)
【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問(wèn)題已經(jīng)有21人回復(fù)
【求助】matlab 解矩陣方程重金懸賞解決方案已經(jīng)有10人回復(fù)
【求助】積分微分方程matlab求解已經(jīng)有6人回復(fù)
【原創(chuàng)】多項(xiàng)式和非多項(xiàng)式曲線(xiàn)擬合的MATLAB實(shí)現(xiàn)視頻教程已經(jīng)有10人回復(fù)
求助曲線(xiàn)擬合已經(jīng)有6人回復(fù)
【求助】求教matlab解非線(xiàn)性方程組已經(jīng)有9人回復(fù)
求助：用Origin做微分方程的曲線(xiàn)擬合已經(jīng)有7人回復(fù)

1樓 2011-06-03 12:43:12

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

回帖支持 ( 顯示支持度最高的前 50 名 )

dubo

金蟲(chóng) (著名寫(xiě)手)

程序強(qiáng)帖: 4
應(yīng)助: 23 (小學(xué)生)
貴賓: 0.779
金幣: 569.2
散金: 3220
紅花: 31
沙發(fā): 1
帖子: 1821
在線(xiàn): 349.5小時(shí)
蟲(chóng)號(hào): 559371
注冊(cè): 2008-05-17
性別: GG
專(zhuān)業(yè): 高分子組裝與超分子結(jié)構(gòu)

【答案】應(yīng)助回帖

★ ★
余澤成(金幣+2, 程序強(qiáng)帖+1): 謝謝參與應(yīng)助！ 2011-06-04 19:30:05
yangger(金幣+3): 非常感謝~~@ 2011-06-17 19:02:33

或者你仿照這個(gè)代碼改改
function y = fitme_diheral( c,phi )
% A1=c(1);
% k1=c(2);
% x1=c(3);
% kBT=443*8.3145/1000;
aa=c(1);
bb=c(2);
cc=c(3);
dd=c(4);
ee=c(5);
ff=c(6);
psi=(phi-pi);
% y=0.5*(f1*(1+cosx)+f2*(1-cos(2*x))+f3*(1+cos(3*x))+f4*(1-cos(4*x)));
y=aa-bb.*cos(psi)+cc.*(cos(psi)).^2-dd.*(cos(psi)).^3+ee.*(cos(psi)).^4-ff.*(cos(psi)).^5;
end

clear;
clc;
load diherals.txt;
data=diherals;
%%avge the theta<0 and theta>0
a=data(1:180,2);
b=data(361:-1:182,2);
avg=0.5*(a+b);
data(1:180,2)=avg;
data(361:-1:182,2)=avg;

p=[24.5156799887984,-1.90007004205358,-1.66217137679272,3.06265999012061,2.33333792446086,-2.48454162655902;];
% data(:,2)=acf(data(:,2));
data3=data;
data3(:,1)=data3(:,1)*pi/180;

for i=1:50
% [p,res]=lsqcurvefit('dihe',p,data3(:,1),data3(:,2));
[p,res]=lsqcurvefit('fitme_diheral',p,data3(:,1),data3(:,2));
%注意修改函數(shù)名dihe
plot(data3(:,1),data3(:,2),'ro');
hold on
plot(data3(:,1),fitme_diheral(p,data3(:,1)));
%注意修改函數(shù)名dihe(p,data3(:,1)
end

贊一下(2人)

回復(fù)此樓

5樓2011-06-04 14:39:52

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

普通回帖

smithjohn

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 4 (幼兒園)
金幣: 1334.9
散金: 372
紅花: 4
帖子: 673
在線(xiàn): 148.8小時(shí)
蟲(chóng)號(hào): 1030430
注冊(cè): 2010-05-28

【答案】應(yīng)助回帖

★
dubo(金幣+1): 謝謝參與 2011-06-04 14:34:19

加上幾十上百組數(shù)不就over了？
為什么要一個(gè)一個(gè)的加啊，就算驗(yàn)證的話(huà)也不用這樣吧？

贊一下(1人)

回復(fù)此樓

2樓2011-06-03 23:28:10

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

楊小胖

金蟲(chóng) (正式寫(xiě)手)

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

★
dubo(金幣+1): 謝謝參與 2011-06-04 14:34:31

你說(shuō)的這step1到step4其實(shí)在matlab中直接用polyfit和polyval這個(gè)函數(shù)不就可以實(shí)現(xiàn)了嗎？
而且step4其實(shí)就是誤差棒。
你再找找help里面的例子吧。我現(xiàn)在電腦上沒(méi)有相關(guān)的代碼。
回頭貼出來(lái)。

贊一下(1人)

回復(fù)此樓

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

3樓2011-06-04 13:36:21

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

dubo

金蟲(chóng) (著名寫(xiě)手)

程序強(qiáng)帖: 4
應(yīng)助: 23 (小學(xué)生)
貴賓: 0.779
金幣: 569.2
散金: 3220
紅花: 31
沙發(fā): 1
帖子: 1821
在線(xiàn): 349.5小時(shí)
蟲(chóng)號(hào): 559371
注冊(cè): 2008-05-17
性別: GG
專(zhuān)業(yè): 高分子組裝與超分子結(jié)構(gòu)

【答案】應(yīng)助回帖

★
余澤成(金幣+1): 謝謝參與應(yīng)助！ 2011-06-04 19:29:51

樓上我的很對(duì)，MATLAB直接就可以擬合的
要么你把數(shù)據(jù)貼上來(lái)，我給你擬合一下

贊一下(1人)

回復(fù)此樓

4樓2011-06-04 14:39:09

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

yangger

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1.1
紅花: 1
帖子: 20
在線(xiàn): 17.2小時(shí)
蟲(chóng)號(hào): 1226591
注冊(cè): 2011-03-09
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by smithjohn at 2011-06-03 23:28:10:
加上幾十上百組數(shù)不就over了？
為什么要一個(gè)一個(gè)的加啊，就算驗(yàn)證的話(huà)也不用這樣吧？

普通編程思路，對(duì)matlab不是很熟的后果。。。。。。。。。。。。。

贊一下

回復(fù)此樓

6樓2011-06-17 19:15:06

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

yangger

新蟲(chóng) (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 1.1
紅花: 1
帖子: 20
在線(xiàn): 17.2小時(shí)
蟲(chóng)號(hào): 1226591
注冊(cè): 2011-03-09
專(zhuān)業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by dubo at 2011-06-04 14:39:52:
或者你仿照這個(gè)代碼改改
function y = fitme_diheral( c,phi )
% A1=c(1);
% k1=c(2);
% x1=c(3);
% kBT=443*8.3145/1000;
aa=c(1);
bb=c(2);
cc=c(3);
dd=c(4);
ee=c(5);
ff=c(6);
psi=(phi-pi);
% ...

非常感謝，我試試

回復(fù)此樓

7樓2011-06-17 19:16:02

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

liuweihong1980

送鮮花一朵

支持

8樓2011-06-21 10:10:08

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

rns521

金蟲(chóng) (正式寫(xiě)手)

應(yīng)助: 0 (幼兒園)
金幣: 2049
散金: 10
紅花: 1
帖子: 377
在線(xiàn): 196.2小時(shí)
蟲(chóng)號(hào): 1115623
注冊(cè): 2010-10-07
性別: GG
專(zhuān)業(yè): 高分子材料結(jié)構(gòu)與性能

都是大牛啊

[ 發(fā)自手機(jī)版 http://www.gaoyang168.com/3g ]

回復(fù)此樓

9樓2011-06-21 11:32:29

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

liang201101

銅蟲(chóng) (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 169.3
帖子: 120
在線(xiàn): 25.2小時(shí)
蟲(chóng)號(hào): 1444020
注冊(cè): 2011-10-15
性別: GG
專(zhuān)業(yè): 制造系統(tǒng)與自動(dòng)化

引用回帖:

5樓: Originally posted by dubo at 2011-06-04 14:39:52
或者你仿照這個(gè)代碼改改
function y = fitme_diheral( c,phi )
% A1=c(1);
% k1=c(2);
% x1=c(3);
% kBT=443*8.3145/1000;
aa=c(1);
bb=c(2);
cc=c(3);
dd=c(4);
ee=c(5);
ff=c(6);
psi=(phi-pi);
% y ...

你好，我有個(gè)相似的問(wèn)題，希望大俠幫幫忙，非常感謝，一共有8組數(shù)據(jù)，abcd是自變量，F(xiàn)是因變量，需要利用這些數(shù)據(jù)用matlab編程擬合成一條曲線(xiàn)，得出一個(gè)含有abcd變量的F=f（abcd）函數(shù)，由于初學(xué)不知道該怎么入手，希望能有源代碼程序，這樣我比較好理解一點(diǎn)，非常感謝！拜托大俠！
         F             a                b             c             d
      4             30             98          330          16
      4             50             80          300          14
      6             70             85          300          14
      9             90             82          300          14
      9          130             61.9       300          14
      12          150             75          300          14
      16          182             90          300          14
      20          200          115          330          16

贊一下

回復(fù)此樓

希望能做出自己的東西，加油!

10樓2013-05-23 21:44:09

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

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

返回列表

下一頁(yè)

普通表情龍兔虎貓

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

[考研] 278求調(diào)劑 +5	煙火先于春 2026-03-17	5/250	2026-03-18 08:43 by 星空星月
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-16	8/400	2026-03-17 19:03 by Wangjingyue
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 本人考085602 化學(xué)工程專(zhuān)碩 +16	不知道叫什么！ 2026-03-15	18/900	2026-03-17 17:05 by ruiyingmiao
[考研] 有沒(méi)有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3	TqlXswl 2026-03-16	7/350	2026-03-17 15:23 by TqlXswl
[考研] 304求調(diào)劑 +4	ahbd 2026-03-14	4/200	2026-03-16 16:48 by 我的船我的海
[考研] 085600調(diào)劑 +5	漾漾123sun 2026-03-12	6/300	2026-03-16 15:58 by 漾漾123sun
[考研] 機(jī)械專(zhuān)碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 328求調(diào)劑 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 331求調(diào)劑（0703有機(jī)化學(xué) +5	ZY-05 2026-03-13	6/300	2026-03-14 10:51 by Jy?
[考研] 304求調(diào)劑 +6	Mochaaaa 2026-03-12	7/350	2026-03-13 22:18 by 星空星月
[考研] 四川大學(xué)085601材料工程專(zhuān)碩初試294求調(diào)劑 +4	祝我們好在冬天 2026-03-11	4/200	2026-03-13 21:39 by peike
[考研] （081700）化學(xué)工程與技術(shù)-298分求調(diào)劑 +12	11啦啦啦 2026-03-11	35/1750	2026-03-13 21:25 by JourneyLucky
[考研] 329求調(diào)劑 +3	miaodesi 2026-03-12	4/200	2026-03-13 20:53 by 18595523086
[考研] 材料專(zhuān)碩350 求調(diào)劑 +4	王金科 2026-03-12	4/200	2026-03-13 16:02 by ruiyingmiao
[考研] 工科278分求調(diào)劑 +5	周慢熱啊 2026-03-12	7/350	2026-03-13 15:49 by JourneyLucky
[考研] 289求調(diào)劑 +3	李政瑩 2026-03-12	3/150	2026-03-13 11:02 by 求調(diào)劑zz
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-12	4/200	2026-03-12 19:33 by 求調(diào)劑zz

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

[求助] 急切求助用matlab曲線(xiàn)擬合色散方程

» 本帖已獲得的紅花（最新10朵）

» 猜你喜歡

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

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

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