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

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

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

匿名

用戶注銷 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 193.9
帖子: 112
在線: 37.6小時(shí)
蟲號(hào): 0
注冊(cè): 2012-09-14
專業(yè): 防災(zāi)工程

本帖僅樓主可見

» 猜你喜歡

中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博士后入站、德國兩年全職已經(jīng)有0人回復(fù)
面向全球申請(qǐng)人的中德博士后2026項(xiàng)目5月30日截止- 國內(nèi)博后入站德國兩年全職已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費(fèi)? 已經(jīng)有214人回復(fù)
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額已經(jīng)有2人回復(fù)
重慶交大材料學(xué)院26年碩士招生復(fù)試實(shí)施細(xì)則已出。歡迎加入光子學(xué)微結(jié)構(gòu)與器件課題組。已經(jīng)有0人回復(fù)
北工大宋曉艷教授團(tuán)隊(duì)招收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ù)
計(jì)算與實(shí)驗(yàn)科學(xué)工程國際會(huì)議(8月7-11日,香港) 二維材料分會(huì)S5與聲學(xué)換能器分會(huì)S8 已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有2人回復(fù)

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

請(qǐng)教下可以用UDF編寫第二類邊界條件么？已經(jīng)有4人回復(fù)
第二類邊界條件的追趕法已經(jīng)有5人回復(fù)
熱傳導(dǎo)第二類邊界條件（傅里葉定理）已經(jīng)有6人回復(fù)
求助：偏微分方程求解已經(jīng)有18人回復(fù)
有限元第二類邊界條件 matlab仿真怎么編已經(jīng)有3人回復(fù)

1樓 2014-05-09 10:48:18

已閱同方向廣播申請(qǐng)計(jì)算強(qiáng)帖回復(fù)此樓編輯查看我的主頁

回帖支持 ( 顯示支持度最高的前 50 名 )

xcenxcen

銀蟲 (小有名氣)

應(yīng)助: 10 (幼兒園)
金幣: 442.3
紅花: 3
帖子: 53
在線: 107.1小時(shí)
蟲號(hào): 3038190
注冊(cè): 2014-03-11
性別: MM
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

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

★
fegg7502: 金幣+1, 應(yīng)助指數(shù)+1, 鼓勵(lì)交流 2014-05-13 09:17:07

function [y err] = soodesv(q1,q2,q3,alpha,beta,a,b,dx)
N = (b-a)/dx+1;
y = zeros(N,1);
D = zeros(N);
f = zeros(N,1);
for ir=1:N
x = a+(ir-1)*dx;
if ir==1
      D(ir,ir) = -q2(x)-2.0/(dx^2);
      D(ir,ir+1) = 2.0/(dx^2);
      f(ir) = 2.0*alpha/dx+q1(x)*alpha+q3(x);
end
if ir>1 && ir<N
      D(ir,ir) = -2.0/(dx^2)-q2(x);
      D(ir,ir+1) = 1.0/(dx^2)-q1(x)/(2.0*dx);
      D(ir,ir-1) = 1.0/(dx^2)+q1(x)/(2.0*dx);
      f(ir) = q3(x);
end
if ir==N
      D(ir,ir) = -q2(x)-2.0/(dx^2);
      D(ir,ir-1) = 2.0/(dx^2);
      f(ir) = q1(x)*beta+q3(x)-2.0*beta/dx;
end
end
y = D\f;
t = a:dx:b;
err = norm(D*y-f);
fprintf('2rd order ode solver error: ||D*y-f||_{2}=%e\n',norm(D*y-f));
plot(t,y);

例子：
y''-3.0*y'+2.0*y-x=0 y'(0)=1 y'(1)=2
>>soodesv(@(x)3.0,@(x)(-2.0),@(x)x,1,2,0,1,0.01);

解的圖像.jpg

結(jié)果.jpg

贊一下(2人)

回復(fù)此樓

待我長發(fā)及腰，遮住一身肥膘�？v然虎背熊腰，也要高冷傲嬌。

6樓2014-05-12 10:16:45

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

普通回帖

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗(yàn): +4
計(jì)算強(qiáng)帖: 12
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

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

★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
fegg7502: 金幣+3, 應(yīng)助指數(shù)+1, 專家考核, 3ks 2014-05-13 09:16:37

這個(gè)可以考慮用bvp解算器做
BVP解算器
solinit = bvpinit(x, yinit, params)
sol = bvpsolver(odefun,bcfun,solinit,options)
由于邊值問題可能有多解，為了便于我們確定那個(gè)解是我們需要的，所以必須使用bvpinit函數(shù)對(duì)初值進(jìn)行估計(jì)
解算器(bvpsolver)：Matlab中提供了bvp4c和bvp5c，后者誤差控制更好些
輸入?yún)?shù)：
x：需要計(jì)算的網(wǎng)格點(diǎn)，相當(dāng)于ode**的tspan
yinit：猜測(cè)的值，可以是具體值，也可以是函數(shù)，類似與 ode**的 x0
params：其它未知參數(shù)，也是一個(gè)猜測(cè)值
odefun：描述邊值問題微分方程的函數(shù)句柄

bcfun：邊值函數(shù)，一般是雙邊值(x 的上下限即認(rèn)為兩個(gè)邊界)，但也支持多邊值(具體看幫助)
bc(y(a),y(b))=0
舉個(gè)例子，y(0)=0,y(4)=-2 ---bcfun可寫成 res =@(ya,yb) [ ya(1) ; yb(1) + 2];
y'(0)=y'(pi)=0,y(0)=1-----------bcfun可寫成res = @(ya,yb) [ ya(2); yb(2);ya(1)-1 ];
solinit：由bvpinit生成的初始化網(wǎng)格
solinit = bvpinit(linspace(pa,pb,10),[1 0]);
bvpinit的第一個(gè)參數(shù)表示求解域及中心節(jié)點(diǎn)，第二個(gè)參數(shù)表示猜測(cè)初始值，基本都寫成[1 0]，它的變化一般對(duì)解無影響
options：BVP解算器優(yōu)化參數(shù)，可以通過bvpset設(shè)置，具體參數(shù)查看幫助

降階法表達(dá)二階的常微分方程，ya(1)表示函數(shù)在求解域[a,b]的左邊界a上的值，ya(2)表示函數(shù)的一階導(dǎo)數(shù)在a上的值，“ya(1)-5”實(shí)際上是“ya(1)-5=0”的簡寫，表示y(x=a)=5，即函數(shù)在左邊界上的值為5。yb是右邊界，其他與ya相同。
用bvp4c求得的解sol是一個(gè)結(jié)構(gòu)體，其中的子項(xiàng)y是一個(gè)2行多列的矩陣，第一行表示函數(shù)值，第二行表示函數(shù)的導(dǎo)數(shù)。另一個(gè)子項(xiàng)yp也是一個(gè)2行多列的矩陣，第一行表示函數(shù)一階導(dǎo)數(shù)值，第二行表示函數(shù)二階導(dǎo)數(shù)值。用deval求出的y與sol.y相同

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

2樓2014-05-09 19:47:44

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

匿名

用戶注銷 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 193.9
帖子: 112
在線: 37.6小時(shí)
蟲號(hào): 0
注冊(cè): 2012-09-14
專業(yè): 防災(zāi)工程

本帖僅樓主可見

贊一下

3樓2014-05-09 21:24:53

已閱申請(qǐng)計(jì)算強(qiáng)帖回復(fù)此樓編輯查看我的主頁

dbb627

榮譽(yù)版主 (著名寫手)

專家經(jīng)驗(yàn): +4
計(jì)算強(qiáng)帖: 12
應(yīng)助: 289 (大學(xué)生)
貴賓: 0.589
金幣: 24640.4
散金: 551
紅花: 61
沙發(fā): 1
帖子: 1246
在線: 1794.8小時(shí)
蟲號(hào): 149791
注冊(cè): 2005-12-29
性別: GG
專業(yè): 污染控制化學(xué)
管轄: 計(jì)算模擬

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

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
478107371: 金幣+60, ★★★★★最佳答案, 大神，請(qǐng)接受我的膝蓋，分全是你一人兒的 2014-05-10 00:11:25
fegg7502: 金幣+2, 應(yīng)助指數(shù)+1, 專家考核, 3ks 2014-05-13 09:16:54

舉個(gè)例子吧
y''-3y'+2y=x， y‘(0)=1,y’(1)=2

CODE:

%解析解

 yy=dsolve('D2y-3*Dy+2*y=x','Dy(0)=1','Dy(1)=2','x');

x=0:0.01:1;hold on;plot(x,eval(yy),'bo')

%bvp 數(shù)值解

 f1=@(x,y)[y(2);x-2*y(1)+3*y(2)];

 solinit = bvpinit(linspace(0,1,20),[-1 1]);

 f2=@(ya,yb)[ya(2)-1;yb(2)-2];

 sol = bvp4c(f1,f2,solinit);

xint = linspace(0,1,25);

yint = deval(sol,xint);

hold on; plot(xint,yint(1,:),'rp-');

贊一下

回復(fù)此樓

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

4樓2014-05-09 21:51:36

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

匿名

用戶注銷 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 193.9
帖子: 112
在線: 37.6小時(shí)
蟲號(hào): 0
注冊(cè): 2012-09-14
專業(yè): 防災(zāi)工程

本帖僅樓主可見

贊一下

5樓2014-05-10 01:25:56

已閱申請(qǐng)計(jì)算強(qiáng)帖回復(fù)此樓編輯查看我的主頁

匿名

用戶注銷 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 193.9
帖子: 112
在線: 37.6小時(shí)
蟲號(hào): 0
注冊(cè): 2012-09-14
專業(yè): 防災(zāi)工程

本帖僅樓主可見

贊一下

7樓2014-05-12 10:39:14

已閱申請(qǐng)計(jì)算強(qiáng)帖回復(fù)此樓編輯查看我的主頁

xcenxcen

銀蟲 (小有名氣)

應(yīng)助: 10 (幼兒園)
金幣: 442.3
紅花: 3
帖子: 53
在線: 107.1小時(shí)
蟲號(hào): 3038190
注冊(cè): 2014-03-11
性別: MM
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

★
fegg7502: 金幣+1, 鼓勵(lì)交流 2014-05-13 09:17:18

引用回帖:

7樓: Originally posted by 478107371 at 2014-05-12 10:39:14
大神，你這個(gè)是用什么方法寫的代碼，能不能稍微解釋一下子...

內(nèi)部網(wǎng)格用的二階中心差分
(y(j+1)-2*y(j)+y(j-1))/(dx^2)-q1j*(y(j+1)-y(j-1))/(2*dx)-q2j*yj-q3=0
j=1,....,N-1
邊界條件用的二階中心差分，引入虛網(wǎng)格點(diǎn)
左邊邊界條件

y1-y(-1))/(2.0*dx)=alpha
y(-1)=y1-2*alpha*dx
利用方程
(y(1)-2*y(0)+y(-1))/(dx^2)-q10*alpha-q20*y0-q3=0
把y(-1)帶入到上式
每個(gè)網(wǎng)格點(diǎn)都有一個(gè)方程
然后形成一個(gè)(N+1)×(N+1)的方程組
Dy=f
解出來就完了

贊一下

回復(fù)此樓

待我長發(fā)及腰，遮住一身肥膘。縱然虎背熊腰，也要高冷傲嬌。

8樓2014-05-12 12:49:14

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

匿名

用戶注銷 (小有名氣)

應(yīng)助: 0 (幼兒園)
金幣: 193.9
帖子: 112
在線: 37.6小時(shí)
蟲號(hào): 0
注冊(cè): 2012-09-14
專業(yè): 防災(zāi)工程

本帖僅樓主可見

贊一下

9樓2014-05-12 19:05:12

已閱申請(qǐng)計(jì)算強(qiáng)帖回復(fù)此樓編輯查看我的主頁

xcenxcen

銀蟲 (小有名氣)

應(yīng)助: 10 (幼兒園)
金幣: 442.3
紅花: 3
帖子: 53
在線: 107.1小時(shí)
蟲號(hào): 3038190
注冊(cè): 2014-03-11
性別: MM
專業(yè): 計(jì)算數(shù)學(xué)與科學(xué)工程計(jì)算

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

★
fegg7502: 金幣+1, 應(yīng)助指數(shù)+1, 3ks 2014-05-13 09:17:28

引用回帖:

9樓: Originally posted by 478107371 at 2014-05-12 19:05:12
邊界條件也是二階差分嗎，是不是一階差分？另外迭代過程呢？思路是不是用有限差分將二階方程變成線性方程組，然后用的牛頓迭代法求解？...

邊界條件也可以用二階的差分啊，當(dāng)然你也可以用一階的。
這是個(gè)線性方程組吧，直接消元就可以解了。你這個(gè)方程規(guī)模也不大，還是一維的問題。
真要迭代可以考慮用Jacobi迭代什么的。
牛頓迭代法不是解非線性方程的嗎？

贊一下

回復(fù)此樓

待我長發(fā)及腰，遮住一身肥膘。縱然虎背熊腰，也要高冷傲嬌。

10樓2014-05-13 08:31:02

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

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 478107371 的主題更新

返回列表

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

[考研] 346求調(diào)劑一志愿070303有機(jī)化學(xué) +5	蘿卜燉青菜 2026-03-28	6/300	2026-03-30 14:32 by limeifeng
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +13	丹青奶蓋 2026-03-26	15/750	2026-03-30 12:35 by fangnagu
[考研] 317分一志愿南理工材料工程本科湖工大求調(diào)劑 +11	芋泥小鈴鐺 2026-03-28	11/550	2026-03-30 12:23 by Wang200018
[考研] 310求調(diào)劑 +9	爭取九點(diǎn)睡 2026-03-30	9/450	2026-03-30 11:44 by Wang200018
[考研] 282求調(diào)劑 +4	wcq131415 2026-03-24	4/200	2026-03-30 10:39 by Delta2012
[考研] 337求調(diào)劑 +6	《樹》 2026-03-29	6/300	2026-03-30 10:15 by herarysara
[考博] 26申博自薦 +6	whh869393 2026-03-24	6/300	2026-03-29 21:05 by nxgogo
[考研] 289求調(diào)劑 +5	BrightLL 2026-03-29	5/250	2026-03-29 17:24 by zhyzzh
[考研] 一志愿雙一流機(jī)械285分求調(diào)劑 +4	幸運(yùn)的三木 2026-03-29	5/250	2026-03-29 14:49 by Miko19
[考研] 求調(diào)劑，一志愿南京航空航天大學(xué) ，080500材料科學(xué)與工程學(xué)碩，總分289分 +7	@taotao 2026-03-29	7/350	2026-03-29 12:03 by longlotian
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5	大饞小子 2026-03-28	5/250	2026-03-29 08:56 by qingfeng258
[考研] 求調(diào)劑推薦材料 304 +15	荷包蛋hyj 2026-03-26	15/750	2026-03-28 04:13 by fmesaito
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿求調(diào)劑 +3	wywy.. 2026-03-27	4/200	2026-03-28 03:32 by fmesaito
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4	崔wj 2026-03-26	4/200	2026-03-27 08:04 by chemisry
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-26	3/150	2026-03-26 19:57 by nihaoar
[考研] 材料調(diào)劑 5+4	想要一壺桃花水 2026-03-25	10/500	2026-03-26 19:56 by 不吃魚的貓
[考研] 07化學(xué)303求調(diào)劑 +5	睿08 2026-03-25	5/250	2026-03-25 22:46 by 418490947
[考研] 302求調(diào)劑 +4	錦衣衛(wèi)藤椒 2026-03-25	4/200	2026-03-25 16:29 by 功夫瘋狂
[考研] 材料調(diào)劑 +3	iwinso 2026-03-23	3/150	2026-03-25 11:29 by greychen00
[考研] 一志愿南航材料專317分求調(diào)劑 +5	炸呀炸呀炸薯?xiàng)l 2026-03-23	5/250	2026-03-24 16:52 by 星空星月

24小時(shí)熱門版塊排行榜

» 猜你喜歡

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

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

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

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

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

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