版塊導(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)版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

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

當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖

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

» 猜你喜歡

一志愿西電085401數(shù)一英一299求調(diào)劑六級(jí)521 已經(jīng)有3人回復(fù)
生物考研337分求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑：085600材料與化工，考材科基，總分319 已經(jīng)有14人回復(fù)
285求調(diào)劑已經(jīng)有10人回復(fù)
化學(xué)工程085602 305分求調(diào)劑已經(jīng)有29人回復(fù)
理學(xué)07化學(xué) 303求調(diào)劑已經(jīng)有5人回復(fù)
081200-11408-276學(xué)碩求調(diào)劑已經(jīng)有3人回復(fù)
土木304求調(diào)劑已經(jīng)有5人回復(fù)
一志愿中國(guó)科學(xué)院大學(xué)265求調(diào)劑已經(jīng)有7人回復(fù)
286求調(diào)劑已經(jīng)有8人回復(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的回帖

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的回帖

查看全部 11 個(gè)回答

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的回帖

查看全部 11 個(gè)回答

普通表情龍兔虎貓

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

[考研] 一志愿中國(guó)科學(xué)院大學(xué)265求調(diào)劑 +6	恬淡ye 2026-03-31	7/350	2026-03-31 11:28 by 恬淡ye
[考研] 313求調(diào)劑 +6	賣(mài)個(gè)關(guān)子吧 2026-03-31	6/300	2026-03-31 10:58 by Jaylen.
[考研] 282求調(diào)劑 +3	呼吸都是減肥 2026-03-31	3/150	2026-03-31 10:00 by wangjy2002
[考研] 一志愿西交288化工專(zhuān)碩求調(diào)劑 +4	好運(yùn)好運(yùn)接接 2026-03-24	4/200	2026-03-31 09:18 by 西京學(xué)院土木工?/a>
[考研] 085602化工求調(diào)劑（331分） +8	111@127 2026-03-30	8/400	2026-03-30 21:23 by 研究僧導(dǎo)導(dǎo)
[考研] 一志愿華東師范大學(xué)有機(jī)化學(xué)專(zhuān)業(yè)，初試351分，復(fù)試被刷求調(diào)劑! +7	真名有冰 2026-03-29	8/400	2026-03-30 21:08 by maaj02
[考研] 309求調(diào)劑 +15	誰(shuí)不是少年 2026-03-29	15/750	2026-03-30 19:28 by 2026材料調(diào)劑
[考研] 311求調(diào)劑 +7	藍(lán)月亮亮 2026-03-30	7/350	2026-03-30 19:24 by 無(wú)際的草原
[考研] 332求調(diào)劑 +6	@MZB382400 2026-03-28	6/300	2026-03-30 16:57 by 無(wú)際的草原
[考研] 化學(xué)308分調(diào)劑 +7	你好明天你好 2026-03-30	8/400	2026-03-30 13:02 by 你好明天你好
[考研] 375求調(diào)劑 +6	雨夏整夜 2026-03-29	6/300	2026-03-30 10:21 by herarysara
[考研] 295求調(diào)劑 +5	wei-5 2026-03-26	5/250	2026-03-30 08:34 by 探123
[考研] 一志愿南航 335分 \| 0856 \| GPA 4.07 \| 有科研經(jīng)歷 +8	cccchenso 2026-03-29	8/400	2026-03-29 23:53 by 我是小康
[考研] 304求調(diào)劑 +6	曼殊2266 2026-03-27	6/300	2026-03-28 14:10 by 唐沐兒
[考研] 調(diào)劑 +3	好好讀書(shū)。 2026-03-28	3/150	2026-03-28 12:04 by 王保杰33
[考研] 341求調(diào)劑 +7	青檸檬1 2026-03-26	7/350	2026-03-27 00:19 by wxiongid
[考研] 調(diào)劑 +4	柚柚yoyo 2026-03-26	4/200	2026-03-26 20:43 by fmesaito
[考研] 環(huán)境專(zhuān)碩324分求調(diào)劑推薦 +5	軒小寧—— 2026-03-26	5/250	2026-03-26 12:05 by i_cooler
[考研] 考研一志愿蘇州大學(xué)初始315（英一）求調(diào)劑 +3	sbdksD 2026-03-24	4/200	2026-03-25 18:16 by xcjcqu
[考研] 302求調(diào)劑 +4	錦衣衛(wèi)藤椒 2026-03-25	4/200	2026-03-25 16:29 by 功夫瘋狂