 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 3484 | 回復(fù): 10 | ||
| 本帖產(chǎn)生 1 個 程序強(qiáng)帖 ，點擊這里進(jìn)行查看 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
[求助]
急切求助用matlab曲線擬合色散方程
|
||
|
現(xiàn)在有如下色散方程： y^2 = A + B/(x^2-C) - D x^2 ----------------Eq. 1 其中A，B, C, D為待定系數(shù)，現(xiàn)在有(x,y)的一組實驗值； 如何得到最佳的A、B、C、D 我的思路如下： step1 通過四對(x_i,y_i) 得到 A，B，C，D 值， step2 代入Eq.1 輸入x實驗值，得到y(tǒng)擬合值， step3 判斷y擬合值和y實驗值差值， step4 如果差值大于10^-4則調(diào)節(jié)A,B,C,D; step5 用新的ABCD，返回到step2 現(xiàn)在遇到的問題是 1.不收斂 2. 調(diào)節(jié)過的A，B，C，D有時候使得Eq.1方程右邊小于零 求助：要是用matlab如何處理， [ Last edited by yangger on 2011-6-3 at 12:53 ] |
回復(fù)此樓» 本帖已獲得的紅花（最新10朵）
liuweihong1980» 猜你喜歡
291 求調(diào)劑
已經(jīng)有5人回復(fù)
-
298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有11人回復(fù)
-
一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑
已經(jīng)有7人回復(fù)
-
316求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料與化工考研調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿重慶大學(xué)085700資源與環(huán)境，總分308求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿211 初試270分 求調(diào)劑
已經(jīng)有4人回復(fù)
-
08工學(xué)調(diào)劑
已經(jīng)有11人回復(fù)
-
350求調(diào)劑
已經(jīng)有6人回復(fù)
-
求調(diào)劑一志愿武漢理工大學(xué)材料工程（085601）
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 高金求助matlab解微分方程組
已經(jīng)有12人回復(fù)
- 【求助】matlab 求解微分方程中的未知參數(shù)
已經(jīng)有20人回復(fù)
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題
已經(jīng)有21人回復(fù)
- 【求助】matlab 解矩陣方程 重金懸賞解決方案
已經(jīng)有10人回復(fù)
- 【求助】積分微分方程matlab求解
已經(jīng)有6人回復(fù)
- 【原創(chuàng)】多項式和非多項式曲線擬合的MATLAB實現(xiàn)視頻教程
已經(jīng)有10人回復(fù)
- 求助曲線擬合
已經(jīng)有6人回復(fù)
- 【求助】求教matlab解非線性方程組
已經(jīng)有9人回復(fù)
- 求助：用Origin做微分方程的曲線擬合
已經(jīng)有7人回復(fù)
3樓2011-06-04 13:36:21
smithjohn
金蟲 (正式寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 1334.9
- 散金: 372
- 紅花: 4
- 帖子: 673
- 在線: 148.8小時
- 蟲號: 1030430
- 注冊: 2010-05-28
2樓2011-06-03 23:28:10
dubo
金蟲 (著名寫手)
- 程序強(qiáng)帖: 4
- 應(yīng)助: 23 (小學(xué)生)
- 貴賓: 0.779
- 金幣: 569.2
- 散金: 3220
- 紅花: 31
- 沙發(fā): 1
- 帖子: 1821
- 在線: 349.5小時
- 蟲號: 559371
- 注冊: 2008-05-17
- 性別: GG
- 專業(yè): 高分子組裝與超分子結(jié)構(gòu)
4樓2011-06-04 14:39:09
dubo
金蟲 (著名寫手)
- 程序強(qiáng)帖: 4
- 應(yīng)助: 23 (小學(xué)生)
- 貴賓: 0.779
- 金幣: 569.2
- 散金: 3220
- 紅花: 31
- 沙發(fā): 1
- 帖子: 1821
- 在線: 349.5小時
- 蟲號: 559371
- 注冊: 2008-05-17
- 性別: GG
- 專業(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
|
或者你仿照這個代碼改改 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ā)表 | |
|---|---|---|---|---|
|
[考研] 08工學(xué)調(diào)劑 +7 | 用戶573181 2026-03-20 | 11/550 |
|
|---|---|---|---|---|
|
|
[考研] 招08考數(shù)學(xué) +6 | laoshidan 2026-03-20 | 14/700 |
|
|
[考研] 263求調(diào)劑 +6 | yqdszhdap－ 2026-03-22 | 9/450 |
|
|
|
[考研] 317求調(diào)劑 +12 | 申子申申 2026-03-19 | 18/900 |
|
|
|
[考研] 材料與化工085600，總分304，本科有兩篇sci參與，求調(diào)劑 +4 | 幸運的醬醬 2026-03-22 | 5/250 |
|
|
|
[考研] 315分，誠求調(diào)劑，材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
|
[考研] 289材料與化工（085600）B區(qū)求調(diào)劑 +3 | 這么名字咋樣 2026-03-22 | 4/200 |
|
|
|
[考研] 08工科 320總分 求調(diào)劑 +11 | 梨花珞晚風(fēng) 2026-03-17 | 11/550 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +5 | 我愛生物生物愛?/a> 2026-03-17 | 5/250 |
|
|
|
[考研] 269專碩求調(diào)劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|
|
[考研] 一志愿華中科技大學(xué)071000，求調(diào)劑 +4 | 沿岸有貝殼6 2026-03-21 | 4/200 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-20 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +3 | 勇敢的小吳 2026-03-20 | 3/150 |
|
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[考研] 265求調(diào)劑 +3 | Jack?k?y 2026-03-17 | 3/150 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑，總分315（英一） +5 | sbdksD 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
|
[考研] 材料學(xué)碩318求調(diào)劑 +5 | February_Feb 2026-03-19 | 5/250 |
|
|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|