| 5 | 1/1 | 返回列表 |
| 查看: 930 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
[求助]
求助高手幫忙用matlab接下面這個(gè)算式,具體見(jiàn)附件
|
|||
|
本人在科研過(guò)程中找到了這個(gè)公式,并做試驗(yàn)驗(yàn)證,公式的計(jì)算規(guī)則如下,不知道要怎么用matlab計(jì)算,請(qǐng)高手幫忙算一下: AN=((0.02105*xN^0.9208)/(1- xN)^0.9208-aN)*0.0002778 1式 BN=(1.515*0.2842* xN /(1+0.2842* xN)-bN)*0.0002778 2式 M=11.8554*A+2.0937*B 3式 總式為0.0010652*(3600*0.8/300-3168*cN/298)-M=0 4式 ①N=1時(shí),a1 =0.007,b1=0.048,x1 =0.0903 將a1、b1、x1代入解出4式中的cN值為c1; ②N=2時(shí)a2=0.007+ A1 b2=0.048+ B1 x2= c1 代入取出值c2 ③ N=3時(shí)a3=0.007+ A1+ A2 b3=0.048+ B1+ B2 x3= c2 代入取出值c3 ④ N=4時(shí)a4=0.007+ A1+ A2+ A3 b4=0.048+ B1+ B2+ B3 x4= c3 代入取出值c4 ⑤ N=5時(shí)a5=0.007+ A1+ A2+ A3+ A4 b5=0.048+ B1+ B2+ B3+ B4 x5= c4 代入取出值c5 以此類推。 請(qǐng)高手給出m文件并將計(jì)算的結(jié)果繪制成曲線,cN-N 或者給出表格形式的cN,N 若不理解請(qǐng)看附件[ Last edited by likuihei on 2011-10-3 at 22:05 ] |

鐵桿木蟲(chóng) (著名寫手)
![]() |
鐵桿木蟲(chóng) (著名寫手)
方丈大師
|
N = 10; n = 1:N; a = zeros(1,N); b = zeros(1,N); c = zeros(1,N); A = zeros(1,N); B = zeros(1,N); x = zeros(1,N); a(1) = 0.007; b(1) = 0.048; x(1) =0.0903; for i=1:N A(i) = ((0.02105*x(i)^0.9208)/(1- x(i))^0.9208-a(i))*0.0002778; B(i) = (1.515*0.2842* x(i) /(1+0.2842* x(i))-b(i))*0.0002778; M = 11.8554*A(i)+2.0937*B(i); c(i) = (3600*0.8/300-M/0.0010652)*298/3168; a(i+1) = a(i)+A(i); b(i+1) = b(i)+B(i); x(i+1) = c(i); end plot(n,c,'bo-') |

鐵桿木蟲(chóng) (著名寫手)
方丈大師
|
N = 50; a = zeros(1,N); b = zeros(1,N); c = zeros(1,N); A = zeros(1,N); B = zeros(1,N); x = zeros(1,N); M = zeros(1,N); S = zeros(1,N); a(1) = 0.007; b(1) = 0.048; x(1) =0.0903; for i=1:N A(i) = ((0.02105*x(i)^0.9208)/(1- x(i))^0.9208-a(i))*0.0002778; B(i) = (1.515*0.2842* x(i) /(1+0.2842* x(i))-b(i))*0.0002778; M(i) = 11.8554*A(i)+2.0937*B(i); S(i) = sum(M(1:i)); c(i) = (3600*0.8/300-S(i)/0.0010652)*298/3168; a(i+1) = a(i)+A(i); b(i+1) = b(i)+B(i); x(i+1) = c(i); end plot(1:N,c,'ro-') |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +7 | 某某某某位 2026-03-21 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 284求調(diào)劑 +10 | Zhao anqi 2026-03-22 | 10/500 |
|
|
[考研] 一志愿國(guó)科過(guò)程所081700,274求調(diào)劑 +3 | 三水研0水立方 2026-03-23 | 3/150 |
|
|
[考研]
|
pk3725069 2026-03-19 | 16/800 |
|
|
[考研] 265求調(diào)劑 +10 | 梁梁校校 2026-03-17 | 10/500 |
|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +7 | 生物工程調(diào)劑 2026-03-17 | 12/600 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問(wèn)津 2026-03-23 | 5/250 |
|
|
[考研] 306求調(diào)劑 +9 | chuanzhu川燭 2026-03-18 | 9/450 |
|
|
[考研] 北科281學(xué)碩材料求調(diào)劑 +8 | tcxiaoxx 2026-03-20 | 9/450 |
|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
[考研] 材料工程專碩 348分求調(diào)劑 +3 | 冬辭. 2026-03-17 | 5/250 |
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
[考研] 0703化學(xué)297求調(diào)劑 +3 | Daisy☆ 2026-03-20 | 3/150 |
|
|
[考研] 299求調(diào)劑 +6 | △小透明* 2026-03-17 | 6/300 |
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 一志愿 南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩 +5 | @taotao 2026-03-20 | 5/250 |
|
|
[考研]
|
不想起名字112 2026-03-19 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|