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

登錄注冊

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

24小時熱門版塊排行榜

返回列表

benling1025

捐助貴賓 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 21.5
帖子: 5
在線: 4.1小時
蟲號: 3385937
注冊: 2014-08-28
專業(yè): 數(shù)論

[求助] 求助修改該matlab擬合帶參數(shù)的時滯微分方程程序

clear all
clc
tspan = [2 4 6 8 10 12 15 20 25 30 35 40 50 60 80 100 120 140 160 180];
x0=[300; 410];
k0 = [1 1 1 1 1 1 1 1]*1e-2;
lb = [0 0 0 0 0 0 0 0 0 0];
ub = [1 1 1 1 1 1 1 1 1 1]*1e5;

data=...
[
2  238  410
4  224  408
6  211  315
8  202  294
10  198  243
12  190  322
15  177  282
20  160  247
25  151  253
30  135  241
40  123  209
50  111  173
60  99  121
80  85  118
100  74  47
120  88  44
140  80  60
160  77  63
180  83  58
];
yexp = data(:,2:3);
[b,resnorm,residual,exitflag,output,lambda,jacobian] = lsqnonlin(@ObjFunc,k0,lb,ub,[],tspan,yexp);

function f = ObjFunc(b,tspan,yexp)          % 目標(biāo)函數(shù)
sol=dde23(@ddefun,[20],@ddehist,tspan);       % 20是時滯
f = [(ysim(:,1)-yexp(:,1)) (ysim(:,2)-yexp(:,2))];

function dydt=ddefun(t,y,Z)          % dde函數(shù)
ylag1=Z(:,1);
dGdt=-b(1)*y(1)-(b(4)*b(8)*y(1)*y(2))/(y(1)+b(8))+b(7);
dIdt=-b(2)*y(2)+b(6)*ylag1(1)
dydt=[dGdt; dIdt];

function S=ddehist(t)       % 初始值
S=[300;410];

這段程序是我根據(jù)論壇上類似的ode45的帖子中的程序段改編的，但是不知道為什么函數(shù)換成了dde23就沒法跑起來了，求高手指點！

回復(fù)此樓

» 猜你喜歡

哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有7人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有236人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有31人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有72人回復(fù)
華東師范大學(xué)芯片設(shè)計徐瓏真實水平咋樣已經(jīng)有1人回復(fù)
調(diào)劑已經(jīng)有0人回復(fù)

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

呼叫版主，在線緊急求助，關(guān)于matlab中微分方程組參數(shù)擬合得問題！已經(jīng)有12人回復(fù)
怎樣求滿足這些條件的平面曲線方程？已經(jīng)有10人回復(fù)
求助一個最小二乘法matlab擬合程序。已經(jīng)有7人回復(fù)
高分求助，用MATLAB模擬一曲線方程，已知曲線上散點已經(jīng)有3人回復(fù)
求助，matlab自定義微分方程擬合實驗數(shù)據(jù)來求方程中的參數(shù) 已經(jīng)有4人回復(fù)
緊急求助，利用Matlab對實驗數(shù)據(jù)進行擬合求解參數(shù)。已經(jīng)有27人回復(fù)
微分方程組參數(shù)擬合的問題（多參數(shù)動力系統(tǒng)參數(shù)估計）已經(jīng)有11人回復(fù)
微分方程組參數(shù)擬合的問題已經(jīng)有11人回復(fù)
MATLAB微分方程參數(shù)擬合問題，求大神已經(jīng)有7人回復(fù)
（十萬火急）切換退化時滯系統(tǒng)（帶約束條件的微分方程）的matlab作圖程序已經(jīng)有4人回復(fù)
matlab數(shù)據(jù)擬合求助�。∠Ｍ呤謳兔懸幌耺atlab的程序謝謝了�。� 已經(jīng)有15人回復(fù)
matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進行數(shù)值比較已經(jīng)有4人回復(fù)
偏微分方程數(shù)值解有基本理論和一些例子（化工類）MATLAB源程序已經(jīng)有296人回復(fù)
lsqnonlin函數(shù)擬合微分方程組參數(shù)擬合問題已經(jīng)有10人回復(fù)
急切求助用matlab曲線擬合色散方程已經(jīng)有10人回復(fù)
求助matlab隱函數(shù)參數(shù)擬合已經(jīng)有16人回復(fù)
【求助】用matlab最優(yōu)化方法進行參數(shù)擬合已經(jīng)有17人回復(fù)
【求助】Zero-order，F(xiàn)irst-order和Higuchi擴散方程的中文意思？已經(jīng)有3人回復(fù)

1樓 2014-08-31 00:29:44

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

塵年墨客

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 43.5
帖子: 14
在線: 3.5小時
蟲號: 2989469
注冊: 2014-02-24
性別: GG
專業(yè): 戰(zhàn)略管理

朋友您好，我最近在做時滯微分方程參數(shù)估計時也遇到了跟你一樣的問題，把ode45改成dde23后就是運行不了，不知道您最終解決了這一問題沒有？如果解決了的話，不知道可否不吝賜教一下，謝謝您！

贊一下

回復(fù)此樓

2樓2015-06-12 16:27:59