| 6 | 1/1 | 返回列表 |
| 查看: 1058 | 回復(fù): 5 | |||
在水一方801新蟲 (初入文壇)
|
[求助]
matlab解超越方程組求助!
|
|
dy = [ y(1)*(2*7.6e-23*y(2)*0.006-2*2.3e-21*y(3)*0.00058-2*1e-20*(1.2e22-y(3))*0.00058-2.64e-20*y(1)-(log(1/0.85)+0.1))/1.56e-9; -4*7.6e-23*y(1)*y(2)/(1.88e-19); 4*y(1)*(-2.3e-21*y(3)+1e-20*(1.2e22-y(3)))/(1.88e-19)]; dy是對(duì)t求導(dǎo), t的范圍在0-1e-7或者0-1e-8; 限制:y(1) = max(y(1),1); 初始條件:y(1) = 1; y(2) = 0; y(3) = 0; 我整理方程得到了這種方程組,不會(huì)解 麻煩高手們把Matlab程序等解題過(guò)程寫的詳細(xì)一些!謝謝! 為了防止上述方程描述錯(cuò)誤 我在下面附件用word把方程寫出來(lái)了! |
主管區(qū)長(zhǎng) (職業(yè)作家)
![]() |
專家經(jīng)驗(yàn): +1059 |

新蟲 (初入文壇)
新蟲 (初入文壇)
新蟲 (初入文壇)
|
再貼一下文件中代碼,方面觀看。 lc clear close all y0 = [1e9;4e20;0]; tspan=[0 1e-7]; tic [t,y] = ode23('rate_eq',tspan,y0); toc y(:,1) = max(y(:,1),1); figure subplot(3,1,1); plot(t,y(:,1)); xlabel(' 時(shí)間(s)'); ylabel('光子數(shù)密度(m^{-3})'); subplot(3,1,2); plot(t,y(:,2)); xlabel('時(shí)間(s)'); ylabel('反轉(zhuǎn)粒子束密度(m^{-3})'); subplot(3,1,3); plot(t,y(:,3)); xlabel('時(shí)間(s)'); ylabel('基態(tài)粒子數(shù)密度{-3})'); function dy = rate_eq(t,y) sigma = 7.6e-23; sigma_o = 1e-20; sigma_s = 2.3e-21; n1 = 1.972; n2 = 3.48; L = 0.1; l = 0.006; d = 0.00058; R = 0.85; c = 2.997963e8; lc = n1*l+n2*d+(0.23-l-d); tr = 2*lc/c; n0 = 1.2e22; beta = 2.6e-10;h = 6.623e-34;lamda = 1.06e-6;niu = c/lamda; B = 6*beta*h*niu*c*(0.3)^2;%此項(xiàng)可以為0; y(1) = max(y(1),1); dy = [ y(1)*(2*sigma*y(2)*l-2*sigma_s*y(3)*d-2*sigma_o*(n0-y(3))*d-B*y(1)-(log(1/R)+L))/tr; -4*sigma*y(1)*y(2)/(h*niu); 4*y(1)*(sigma_o*(n0-y(3))-sigma_s*y(3))/(h*niu)]; 請(qǐng)蟲友幫忙! |
新蟲 (初入文壇)
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
[考博] 環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名 +3 | QGZDSYS 2026-03-13 | 5/250 |
|
|
[考研] 296求調(diào)劑 +5 | 大口吃飯 身體健 2026-03-13 | 5/250 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個(gè)二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
[考研] 0703 物理化學(xué)調(diào)劑 +3 | 我可以上岸的對(duì)?/a> 2026-03-13 | 5/250 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
[考研] 308 085701 四六級(jí)已過(guò)求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
[考研] 材料與化工(0856)304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|
|
[考研] 328化工專碩求調(diào)劑 +4 | 。,。,。,。i 2026-03-12 | 4/200 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級(jí)已過(guò) 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|