 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3485 | 回復(fù): 10 | ||
| 本帖產(chǎn)生 1 個(gè) 程序強(qiáng)帖 ，點(diǎn)擊這里進(jìn)行查看 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
急切求助用matlab曲線擬合色散方程
|
||
|
現(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» 猜你喜歡
086003食品工程求調(diào)劑
已經(jīng)有3人回復(fù)
-
【考研調(diào)劑】化學(xué)專(zhuān)業(yè) 281分，一志愿四川大學(xué)，誠(chéng)心求調(diào)劑
已經(jīng)有12人回復(fù)
-
一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑
已經(jīng)有12人回復(fù)
-
材料調(diào)劑
已經(jīng)有5人回復(fù)
-
081700 調(diào)劑 267分
已經(jīng)有8人回復(fù)
-
請(qǐng)教下大家 2026年國(guó)家基金申請(qǐng)是雙盲審嗎？
已經(jīng)有5人回復(fù)
-
276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷
已經(jīng)有10人回復(fù)
-
一志愿南航材料專(zhuān)317分求調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)學(xué)碩求調(diào)劑
已經(jīng)有5人回復(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)式曲線擬合的MATLAB實(shí)現(xiàn)視頻教程
已經(jīng)有10人回復(fù)
- 求助曲線擬合
已經(jīng)有6人回復(fù)
- 【求助】求教matlab解非線性方程組
已經(jīng)有9人回復(fù)
- 求助：用Origin做微分方程的曲線擬合
已經(jīng)有7人回復(fù)
dubo
金蟲(chóng) (著名寫(xiě)手)
- 程序強(qiáng)帖: 4
- 應(yīng)助: 23 (小學(xué)生)
- 貴賓: 0.779
- 金幣: 569.2
- 散金: 3220
- 紅花: 31
- 沙發(fā): 1
- 帖子: 1821
- 在線: 349.5小時(shí)
- 蟲(chóng)號(hào): 559371
- 注冊(cè): 2008-05-17
- 性別: GG
- 專(zhuān)業(yè): 高分子組裝與超分子結(jié)構(gòu)
4樓2011-06-04 14:39:09
smithjohn
金蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 1334.9
- 散金: 372
- 紅花: 4
- 帖子: 673
- 在線: 148.8小時(shí)
- 蟲(chóng)號(hào): 1030430
- 注冊(cè): 2010-05-28
2樓2011-06-03 23:28:10
3樓2011-06-04 13:36:21
dubo
金蟲(chóng) (著名寫(xiě)手)
- 程序強(qiáng)帖: 4
- 應(yīng)助: 23 (小學(xué)生)
- 貴賓: 0.779
- 金幣: 569.2
- 散金: 3220
- 紅花: 31
- 沙發(fā): 1
- 帖子: 1821
- 在線: 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
余澤成(金幣+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 |
5樓2011-06-04 14:39:52
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 276求調(diào)劑。有半年電池和半年高分子實(shí)習(xí)經(jīng)歷 +9 | 材料學(xué)257求調(diào)劑 2026-03-23 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +4 | 奶油草莓. 2026-03-22 | 5/250 |
|
|
|
[考研] 生物學(xué)一志愿985，分?jǐn)?shù)349求調(diào)劑 +6 | zxts12 2026-03-21 | 9/450 |
|
|
|
[考研] 291求調(diào)劑 +5 | 孅華 2026-03-22 | 5/250 |
|
|
|
[考研] 352求調(diào)劑 +3 | 大米飯！ 2026-03-22 | 3/150 |
|
|
|
[考研] 0854電子信息求調(diào)劑 +3 | α____ 2026-03-22 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +5 | Zhangbod 2026-03-21 | 7/350 |
|
|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-20 | 7/350 |
|
|
|
[考研] 求助 +5 | 夢(mèng)里的無(wú)言 2026-03-21 | 6/300 |
|
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
|
[考研] 材料 271求調(diào)劑 +5 | 展信悅_ 2026-03-21 | 5/250 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂(lè)29 2026-03-18 | 5/250 |
|
|
|
[考研] 求調(diào)劑，一志愿:南京航空航天大學(xué)大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +4 | @taotao 2026-03-19 | 4/200 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒(méi)變，給編輯又寫(xiě)郵件了，沒(méi)回復(fù)
 10+3
|
wangf9518 2026-03-17 | 4/200 |
|