當(dāng)前只顯示滿(mǎn)足指定條件的回帖，點(diǎn)擊這里查看本話(huà)題的所有回帖

在水一方801

應(yīng)助: 1 (幼兒園)
金幣: 12.5
帖子: 18
在線(xiàn): 12.8小時(shí)
蟲(chóng)號(hào): 2796876
注冊(cè): 2013-11-12
專(zhuān)業(yè): 光學(xué)

[求助] 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ò)程寫(xiě)的詳細(xì)一些！謝謝！
為了防止上述方程描述錯(cuò)誤  我在下面附件用word把方程寫(xiě)出來(lái)了！

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶(hù)自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶(hù)本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : matlab解超越方程組求助.docx

2015-05-26 10:29:20, 12.03 K

» 猜你喜歡

2026年循環(huán)經(jīng)濟(jì)功能材料國(guó)際會(huì)議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國(guó)際會(huì)議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有112人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請(qǐng)問(wèn)還有沒(méi)有用Latex寫(xiě)文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有30人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運(yùn)，祝面上順利！已經(jīng)有61人回復(fù)

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

新人求助解MATLAB超越方程，方程中含有復(fù)數(shù) 已經(jīng)有9人回復(fù)
matlab對(duì)含貝塞爾函數(shù)的超越方程求根已經(jīng)有5人回復(fù)
求助：超越方程的解法已經(jīng)有13人回復(fù)
求解一個(gè)方程組，急！��！已經(jīng)有19人回復(fù)
matlab如何求解多元多次方程，包括繪圖已經(jīng)有7人回復(fù)
matlab解方程組求助攻已經(jīng)有4人回復(fù)
多元一次方程的matlab解法，要求解為正整數(shù) 已經(jīng)有10人回復(fù)
求助用matlab或者（1st0pt)編程解方程組已經(jīng)有3人回復(fù)
求助，用matlab解微分方程組，希望高手能給予指點(diǎn) 已經(jīng)有17人回復(fù)
matlab解偏微分方程求助已經(jīng)有12人回復(fù)
matlab求解高階方程組已經(jīng)有3人回復(fù)
有償求助：matlab怎么解聯(lián)立方程組？已經(jīng)有11人回復(fù)
matlab求解非線(xiàn)性方程組，錯(cuò)誤提示怎么解決已經(jīng)有5人回復(fù)
求助 MATLAB解方程組-fslove 已經(jīng)有7人回復(fù)
求助matlab---fsolve解非線(xiàn)性方程組已經(jīng)有6人回復(fù)
matlab求解色散方程，求助已經(jīng)有12人回復(fù)
求助Matlab解方程問(wèn)題已經(jīng)有6人回復(fù)
MATLAB解方程出現(xiàn)??? Subscripted assignment dimension mismatch. 已經(jīng)有3人回復(fù)
【求助】matlab求解非線(xiàn)性方程組，并畫(huà)圖處理。要求y,z是實(shí)數(shù)解！已經(jīng)有18人回復(fù)
【求助】matlab 解矩陣方程重金懸賞解決方案已經(jīng)有10人回復(fù)
【求助】復(fù)超越方程的解法已經(jīng)有19人回復(fù)
【求助】求教matlab解非線(xiàn)性方程組已經(jīng)有9人回復(fù)

1樓 2015-05-26 10:30:27

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

在水一方801

新蟲(chóng) (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 12.5
帖子: 18
在線(xiàn): 12.8小時(shí)
蟲(chóng)號(hào): 2796876
注冊(cè): 2013-11-12
專(zhuān)業(yè): 光學(xué)

引用回帖:

2樓: Originally posted by 月只藍(lán) at 2015-05-26 19:00:43
function solve_odes
clear all;clc

y0=;
tspan=0:1e-9:3e-8;
=ode45(@myfun,tspan,y0);
y1=y(:,1)
y2=y(:,2)
y3=y(:,3)
figure(1)
plot(t,y(:,1),'--r<')
figure(2)
plot(t,y(:,2),'b-o')
figure ...

我把原來(lái)的式子、參數(shù)、和應(yīng)該得到的圖形上傳上來(lái)了，并附我寫(xiě)的程序，麻煩請(qǐng)看一下哪里出了問(wèn)題，為什么得不到想要的結(jié)果，請(qǐng)大神寫(xiě)清具體步驟！求賜教！并附全部金幣，只有8個(gè)了。

參數(shù)1.png

參數(shù)2.jpg

參數(shù)3.jpg

待解式子.png

結(jié)果圖（應(yīng)該為I的分布圖）.png

贊一下

回復(fù)此樓

» 本帖附件資源列表

歡迎監(jiān)督和反饋：小木蟲(chóng)僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶(hù)自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶(hù)本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com
附件 1 : 參數(shù)和代碼.docx

2015-05-27 19:01:30, 13.03 K

4樓2015-05-27 19:04:19

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

查看全部 6 個(gè)回答

月只藍(lán)

主管區(qū)長(zhǎng) (職業(yè)作家)

專(zhuān)家經(jīng)驗(yàn): +1059
計(jì)算強(qiáng)帖: 8
應(yīng)助: 1712 (講師)
貴賓: 8.888
金幣: 68127.7
散金: 1938
紅花: 443
沙發(fā): 4
帖子: 4373
在線(xiàn): 3291.4小時(shí)
蟲(chóng)號(hào): 1122189
注冊(cè): 2010-10-14
專(zhuān)業(yè): 宇宙學(xué)
管轄: 計(jì)算模擬區(qū)

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
在水一方801: 金幣+7, ★★★★★最佳答案, 大蝦的回答很有對(duì)我啟示，雖然沒(méi)有得到我想要的結(jié)果，原因可能是我表達(dá)的不夠明確。把昨天全部的7個(gè)金幣給您。下面我把原來(lái)的式子參數(shù)發(fā)一下，并把今天的8個(gè)金附上，希望高手蟲(chóng)友幫忙解答一下，有可能是我對(duì)式子不到位而引起的錯(cuò)誤。 2015-05-27 16:45:29

CODE:

function solve_odes

clear all;clc

y0=[1 0 0];

tspan=0:1e-9:3e-8;

[t y]=ode45(@myfun,tspan,y0);

y1=y(:,1)

y2=y(:,2)

y3=y(:,3)

figure(1)

plot(t,y(:,1),'--r<')

figure(2)

plot(t,y(:,2),'b-o')

figure(3),

plot(t,y(:,3),'k<-')

function dydt=myfun(t,y)

dydt = [  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)];

附圖1.png

附圖2.png

附圖3.png

贊一下(1人)

回復(fù)此樓

MATLAB、MS小問(wèn)題、普通問(wèn)題請(qǐng)發(fā)帖求助！時(shí)間精力有限，恕不接受無(wú)償私信求助。

2樓2015-05-26 19:00:43

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

在水一方801

新蟲(chóng) (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 12.5
帖子: 18
在線(xiàn): 12.8小時(shí)
蟲(chóng)號(hào): 2796876
注冊(cè): 2013-11-12
專(zhuān)業(yè): 光學(xué)

引用回帖:

非常感謝您的幫忙，這個(gè)規(guī)律的圖我也畫(huà)出來(lái)著，我糾結(jié)的是這些參數(shù)得不出文獻(xiàn)中所給的圖形，會(huì)不會(huì)是所取精度不夠。規(guī)律圖類(lèi)似梳妝函數(shù)，或者脈沖波形。我自己解決不了，所以過(guò)來(lái)請(qǐng)求大神！

趨勢(shì)圖1.png

趨勢(shì)圖2.png

贊一下

回復(fù)此樓

3樓2015-05-26 21:53:10

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

在水一方801

新蟲(chóng) (初入文壇)

應(yīng)助: 1 (幼兒園)
金幣: 12.5
帖子: 18
在線(xiàn): 12.8小時(shí)
蟲(chóng)號(hào): 2796876
注冊(cè): 2013-11-12
專(zhuān)業(yè): 光學(xué)

再貼一下文件中代碼，方面觀看。

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)蟲(chóng)友幫忙！

贊一下

回復(fù)此樓

5樓2015-05-27 19:08:59

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

最具人氣熱帖推薦 [查看全部]	作者	回/看	最后發(fā)表

[考研] 一志愿西安交通大學(xué)材料工程專(zhuān)業(yè) 282分求調(diào)劑 +5	楓橋ZL 2026-03-18	7/350	2026-03-19 14:52 by 功夫瘋狂
[考研] 材料080500調(diào)劑求收留 +4	一顆meteor 2026-03-13	4/200	2026-03-19 10:32 by 30660438
[考研] 材料專(zhuān)碩英一數(shù)二306 +5	z1z2z3879 2026-03-18	5/250	2026-03-19 07:43 by BruceLiu320
[考研] 0703化學(xué)調(diào)劑 +8	妮妮ninicgb 2026-03-15	12/600	2026-03-19 05:55 by anny19840123
[考研] 330求調(diào)劑 +3	小材化本科 2026-03-18	3/150	2026-03-18 21:55 by 無(wú)懈可擊111
[考研] 328求調(diào)劑，英語(yǔ)六級(jí)551，有科研經(jīng)歷 +3	生物工程調(diào)劑 2026-03-17	7/350	2026-03-18 20:41 by Wangjingyue
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 085600材料與化工 +5	安全上岸！ 2026-03-16	5/250	2026-03-18 15:33 by cmz0325
[考研] 070300化學(xué)319求調(diào)劑 +6	錦鯉0909 2026-03-17	6/300	2026-03-18 13:22 by Iveryant
[考研] 生物學(xué)071000 329分求調(diào)劑 +3	我愛(ài)生物生物愛(ài)?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[考研] 268求調(diào)劑 +8	一定有學(xué)上- 2026-03-14	9/450	2026-03-17 17:47 by laoshidan
[考研] 326求調(diào)劑 +5	上岸的小葡 2026-03-15	6/300	2026-03-17 17:26 by ruiyingmiao
[考研] 材料專(zhuān)碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 275求調(diào)劑 +4	太陽(yáng)花天天開(kāi)心 2026-03-16	4/200	2026-03-17 10:53 by 功夫瘋狂
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[碩博家園] 085600 260分求調(diào)劑 +3	天空還下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空還下雨么