 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 223 | 回復(fù): 3 | |||
我要飛鐵蟲 (正式寫手)
|
[求助]
為何出錯？
|
 請大家看看這個程序為何出錯？如何修改？謝謝了！function [T, Y]=FOChen(parameters, orders, TSim, Y0) % format long e % Numerical Solution of the Fractional-Order Chen's System % % D^q1 x(t) = a(y(t)-x(t)) % D^q2 y(t) = dx(t) - x(t)z(t) + cy(t) % D^q3 z(t) = x(t)y(t) - bz(t) % % function [T, Y] = FOChen(parameters, orders, TSim, Y0) % % Input: parameters - model parameters [a, b, c, d] % orders - derivatives orders [q1, q2, q3] % TSim - simulation time (0 - TSim) in sec % Y0 - initial conditions [Y0(1), Y0(2), Y0(3)] % % Output: T - simulation time (0 : Tstep : TSim) % Y - solution of the system (x=Y(1), y=Y(2), z=Y(3)) % % Author: (c) Ivo Petras (ivo.petras@tuke.sk), 2010. % % time step: h=0.005; % number of calculated mesh points: n=100/h; %orders of derivatives, respectively: q1=0.99; q2=0.96; q3=0.92; % constants of Chen's system: a=10; b=1; c=1; d=5; % binomial coefficients calculation: cp1=1; cp2=1; cp3=1; for j=1:n c1(j)=(1-(1+q1)/j)*cp1; c2(j)=(1-(1+q2)/j)*cp2; c3(j)=(1-(1+q3)/j)*cp3; cp1=c1(j); cp2=c2(j); cp3=c3(j); end % initial conditions setting: x(1)=0.5; y(1)=0.6; z(1)=0.1; % calculation of phase portraits /numerical solution/: for i=2:n x(i)=(a*(y(i-1)-x(i-1)))*h^q1 - memo(x, c1, i); y(i)=(-d*x(i)-x(i)*z(i-1)+c*y(i-1))*h^q2 - memo(y, c2, i); z(i)=(x(i)*y(i)-b*z(i-1))*h^q3 - memo(z, c3, i); end for j=1:n Y(j,1)=x(j); Y(j,2)=y(j); Y(j,3)=z(j); end T=h:h:TSim; % [ Last edited by 我要飛 on 2011-11-18 at 10:00 ] |
回復(fù)此樓» 猜你喜歡
調(diào)劑
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費?
已經(jīng)有285人回復(fù)
-
中德博士后2026項目5月30日截止- 國內(nèi)博士后入站、德國兩年全職
已經(jīng)有0人回復(fù)
-
面向全球申請人的中德博士后2026項目5月30日截止- 國內(nèi)博后入站德國兩年全職
已經(jīng)有0人回復(fù)
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有2人回復(fù)
-
重慶交大材料學(xué)院26年碩士招生復(fù)試實施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復(fù)
-
北工大宋曉艷教授團隊招收2026年9月入學(xué)博士【學(xué)博和專博名額均有】！
已經(jīng)有1人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收10人
已經(jīng)有0人回復(fù)
-
天津城建大學(xué)理學(xué)院化學(xué)、環(huán)境工程、材料物理等方向接收調(diào)劑，10人
已經(jīng)有0人回復(fù)
dbb627
榮譽版主 (著名寫手)
-
專家經(jīng)驗: +4 - 計算強帖: 12
- 應(yīng)助: 289 (大學(xué)生)
- 貴賓: 0.589
- 金幣: 24640.4
- 散金: 551
- 紅花: 61
- 沙發(fā): 1
- 帖子: 1246
- 在線: 1794.8小時
- 蟲號: 149791
- 注冊: 2005-12-29
- 性別: GG
- 專業(yè): 污染控制化學(xué)
- 管轄: 計算模擬
2樓2011-11-18 10:18:17
我要飛
鐵蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1515.3
- 散金: 144
- 紅花: 5
- 帖子: 458
- 在線: 591.2小時
- 蟲號: 671362
- 注冊: 2008-12-09
- 專業(yè): 計算數(shù)學(xué)與科學(xué)工程計算
3樓2011-11-18 10:25:27
dbb627
榮譽版主 (著名寫手)
-
專家經(jīng)驗: +4 - 計算強帖: 12
- 應(yīng)助: 289 (大學(xué)生)
- 貴賓: 0.589
- 金幣: 24640.4
- 散金: 551
- 紅花: 61
- 沙發(fā): 1
- 帖子: 1246
- 在線: 1794.8小時
- 蟲號: 149791
- 注冊: 2005-12-29
- 性別: GG
- 專業(yè): 污染控制化學(xué)
- 管轄: 計算模擬
4樓2011-11-18 10:30:05
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 我要飛 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +5 | 腳滑的守法公民 2026-03-27 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 266分求材料化工冶金礦業(yè)等專業(yè)的調(diào)劑 +4 | 哇呼哼呼哼 2026-03-26 | 4/200 |
|
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
|
[考研] 考研化學(xué)308分求調(diào)劑 +10 | 你好明天你好 2026-03-23 | 12/600 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +7 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
|
[考研] 284求調(diào)劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
|
[考研] 【雙一流院校新能源、環(huán)境材料，材料加工與模擬招收大量調(diào)劑】 +4 | Higraduate 2026-03-22 | 8/400 |
|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-26 | 6/300 |
|
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(333分) +9 | holy J 2026-03-21 | 9/450 |
|
|
|
[考研] 291 求調(diào)劑 +7 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 8/400 |
|
|
|
[考研] 一志愿哈工大，085400，320，求調(diào)劑 +4 | gdlf9999 2026-03-24 | 4/200 |
|
|
|
[考研] 299求調(diào)劑 +4 | 15188958825 2026-03-25 | 4/200 |
|
|
|
[考研] 296求調(diào)劑 +4 | 汪��？！ 2026-03-25 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
|
[考研] 一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑 +3 | 開開心心沒煩惱 2026-03-23 | 4/200 |
|
|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
|
[考研] 384求調(diào)劑 +3 | 子系博 2026-03-22 | 6/300 |
|
|
|
[考研]
 一志愿070300浙大化學(xué)358分，求調(diào)劑！
+4
|
酥酥魚.. 2026-03-21 | 4/200 |
|
|
|
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|