 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 3200 | 回復: 17 | ||||
musejianglin銀蟲 (小有名氣)
|
[求助]
求助，用matlab解微分方程組，希望高手能給予指點
|
|||
|
微分方程的形式如下所示： d2y/dx2 =a*sinh(y) 邊界條件為y(0)=b,dy/dx(x=0處)=c。 求y與x的曲線關(guān)系。 因為已經(jīng)知道y 不存在解析解，只存在數(shù)值解，而且文獻中提到是用matlab中的bvp4c路徑做的，所以請matlab高手和數(shù)學高手駐足，幫忙解決一下，不勝感激。 |
回復此樓» 猜你喜歡
計算與實驗科學工程國際會議(8月7-11日,香港) 二維材料分會S5與聲學換能器分會S8
已經(jīng)有0人回復
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有4人回復
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有211人回復
-
歡迎調(diào)劑到贛南師范大學智能制造與未來能源學院
已經(jīng)有29人回復
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有4人回復
-
歡迎調(diào)劑到贛南師范大學智能制造與未來能源學院
已經(jīng)有26人回復
-
重慶交大26年碩士生招生擬調(diào)劑通知已出！歡迎加入光子學微結(jié)構(gòu)與器件課題組。
已經(jīng)有0人回復
-
UJN物理學專業(yè)調(diào)劑
已經(jīng)有7人回復
-
0702一志愿吉大B區(qū)求調(diào)劑 本科期間發(fā)表一篇Sci
已經(jīng)有0人回復
-
北京納米能源與系統(tǒng)研究所王中林院士/曹南穎研究員課題組招收2026碩士研究生1名
已經(jīng)有0人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 請問這個4階的非線性偏微分方程組在PDE中怎么解？等待高手解答。
已經(jīng)有17人回復
- 求Matlab高手解決線性方程組的迭代求解問題
已經(jīng)有17人回復
- 如何用matlab求解矩陣系數(shù)的二階微分方程
已經(jīng)有13人回復
- matlab求解偏微分方程組
已經(jīng)有3人回復
- 求助 MATLAB解方程組-fslove
已經(jīng)有7人回復
- 求助matlab---fsolve解非線性方程組
已經(jīng)有6人回復
- 高手進來用matlab求解常微分方程，十分感謝！
已經(jīng)有17人回復
- matlab有限差分法求解復雜的偏微分方程組
已經(jīng)有19人回復
- matlab微分方程組參數(shù)擬合，以周為單位求解，匯總后以年為單位進行數(shù)值比較
已經(jīng)有4人回復
- matlab解微分方程
已經(jīng)有10人回復
- Matlab求解偏微分方程組
已經(jīng)有22人回復
- matlab數(shù)值求解邊界條件微分方程組
已經(jīng)有7人回復
- 【求助】matlab 求解微分方程中的未知參數(shù)
已經(jīng)有20人回復
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個問題
已經(jīng)有21人回復
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已經(jīng)有22人回復
- 【求助】積分微分方程matlab求解
已經(jīng)有6人回復
- 【求助】matlab怎么求解偏微分方程組啊，先謝謝了
已經(jīng)有13人回復
- 【求助】有限元法求解二階偏微分方程組
已經(jīng)有6人回復
- 【求助】求教matlab解非線性方程組
已經(jīng)有9人回復
xxz903
金蟲 (小有名氣)
- 應(yīng)助: 33 (小學生)
- 金幣: 1045.6
- 散金: 190
- 紅花: 4
- 帖子: 253
- 在線: 66.7小時
- 蟲號: 1360336
- 注冊: 2011-08-03
- 性別: GG
- 專業(yè): 化工系統(tǒng)工程
2樓2013-05-14 23:20:11
月只藍
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗: +1059 - 計算強帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68128.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學
- 管轄: 計算模擬區(qū)
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+2, 謝謝 2013-05-15 17:44:25
musejianglin: 金幣+15 2013-05-16 21:19:25
musejianglin: 金幣+3 2013-05-20 09:05:12
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+2, 謝謝 2013-05-15 17:44:25
musejianglin: 金幣+15 2013-05-16 21:19:25
musejianglin: 金幣+3 2013-05-20 09:05:12
|
你的方程少個邊界條件吧，是不是少個"dy/dx(x=某數(shù)值)= ”的邊界條件？ 不然按你目前的邊界條件，都是x=0處的初始值，這個微分方程已經(jīng)不是邊值問題了，而是初值問題了。在現(xiàn)有條件下，解該問題用ode可解，程序如下，復制進一個m文件，F(xiàn)5運行即可，根據(jù)需要輸入a，b，c的數(shù)值。a=b=c=1的結(jié)果見附圖1，該方程組有剛性，不知道你希望x什么多大的區(qū)間內(nèi)求解。補充完邊界條件后，可再發(fā)上來。 function solvebvp333 clear all;clc global a b c a=input('參數(shù)a='); b=input('參數(shù)b='); c=input('參數(shù)c='); y0=[1 1]; xspan=0:0.02:1; [x,y]=ode15s(@fun332,xspan,y0); figure plot(x,y(:,1),'bo-'); function dydx=fun332(x,y) global a dydx=[y(2);a*sinh(y(1))]; 附圖1.jpg |
6樓2013-05-15 09:24:36
月只藍
主管區(qū)長 (職業(yè)作家)
-
專家經(jīng)驗: +1059 - 計算強帖: 8
- 應(yīng)助: 1712 (講師)
- 貴賓: 8.888
- 金幣: 68128.7
- 散金: 1938
- 紅花: 443
- 沙發(fā): 4
- 帖子: 4373
- 在線: 3291.4小時
- 蟲號: 1122189
- 注冊: 2010-10-14
- 專業(yè): 宇宙學
- 管轄: 計算模擬區(qū)
7樓2013-05-15 09:32:00
arising2010
銅蟲 (正式寫手)
- 應(yīng)助: 32 (小學生)
- 金幣: 1433.4
- 紅花: 11
- 帖子: 996
- 在線: 412.6小時
- 蟲號: 2464689
- 注冊: 2013-05-15
- 專業(yè): 可再生與替代能源利用中的
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
musejianglin: 金幣+2, ★有幫助 2013-05-15 09:06:33
fegg7502: 金幣+2, 應(yīng)助指數(shù)+1, 3ks 2013-05-18 14:31:46
musejianglin: 金幣+3 2013-05-20 09:04:56
感謝參與，應(yīng)助指數(shù) +1
musejianglin: 金幣+2, ★有幫助 2013-05-15 09:06:33
fegg7502: 金幣+2, 應(yīng)助指數(shù)+1, 3ks 2013-05-18 14:31:46
musejianglin: 金幣+3 2013-05-20 09:04:56
|
http://wenku.baidu.com/view/953dce370912a21614792967.html http://wenku.baidu.com/view/690deff2770bf78a65295436.html 可以參考一下上面兩個連接 |
3樓2013-05-15 01:26:56
musejianglin
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1301
- 帖子: 91
- 在線: 16.6小時
- 蟲號: 1225342
- 注冊: 2011-03-08
- 性別: MM
- 專業(yè): 膠體與界面化學
4樓2013-05-15 09:03:36
musejianglin
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1301
- 帖子: 91
- 在線: 16.6小時
- 蟲號: 1225342
- 注冊: 2011-03-08
- 性別: MM
- 專業(yè): 膠體與界面化學
5樓2013-05-15 09:06:08
somomo91
專家顧問 (職業(yè)作家)
-
專家經(jīng)驗: +3 - 計算強帖: 1
- 應(yīng)助: 144 (高中生)
- 金幣: 19353
- 紅花: 115
- 帖子: 3854
- 在線: 1053.2小時
- 蟲號: 2399301
- 注冊: 2013-04-03
- 專業(yè): 動力學與控制
- 管轄: 計算模擬
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+1, 謝謝 2013-05-16 17:09:32
musejianglin: 金幣+15, ★★★很有幫助 2013-05-16 21:19:06
感謝參與，應(yīng)助指數(shù) +1
csgt0: 金幣+1, 謝謝 2013-05-16 17:09:32
musejianglin: 金幣+15, ★★★很有幫助 2013-05-16 21:19:06
|
不必用1stopt， MATLAB 完全可以做到 這里是程序 ： =================================== % main program close all; clear all; clc global a a = 0.1; t0 = 0; tend=5000; tspan = t0:0.5:tend; v_y_ini = [1 0]; [t,Y] = ode45(@ode_sys, tspan, v_y_ini); v = Y(:, 1); y = Y(:, 2); =================================== % ode system function [ dy ] = ode_sys( ~, Y ) global a v = Y(1); y = Y(2); dy = zeros(2,1); dy(1) = a * sinh( y ); dy(2) = v ; end =================================== |
8樓2013-05-15 18:28:25
somomo91
專家顧問 (職業(yè)作家)
-
專家經(jīng)驗: +3 - 計算強帖: 1
- 應(yīng)助: 144 (高中生)
- 金幣: 19353
- 紅花: 115
- 帖子: 3854
- 在線: 1053.2小時
- 蟲號: 2399301
- 注冊: 2013-04-03
- 專業(yè): 動力學與控制
- 管轄: 計算模擬
9樓2013-05-15 18:29:29
musejianglin
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1301
- 帖子: 91
- 在線: 16.6小時
- 蟲號: 1225342
- 注冊: 2011-03-08
- 性別: MM
- 專業(yè): 膠體與界面化學
10樓2013-05-16 21:46:54
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 musejianglin 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856初試324分求調(diào)劑 +6 | 想上學求調(diào) 2026-04-01 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 266分，一志愿電氣工程，本科材料，求材料專業(yè)調(diào)劑 +10 | 哇呼哼呼哼 2026-04-01 | 11/550 |
|
|
[考研] 0710生物學求調(diào)劑 +9 | manman511 2026-04-01 | 9/450 |
|
|
|
[考研] 材料調(diào)劑 +12 | 一樣YWY 2026-04-01 | 12/600 |
|
|
|
[考研] 一志愿南昌大學324求調(diào)劑 +12 | hanamiko 2026-03-27 | 12/600 |
|
|
|
[考研] 286求調(diào)劑 +16 | PolarBear11 2026-03-26 | 16/800 |
|
|
|
[考研] 一志愿西安交大材料學碩（英一數(shù)二）347，求調(diào)劑到高分子/材料相關(guān)專業(yè) +7 | zju51 2026-03-31 | 9/450 |
|
|
|
[考研] 349求調(diào)劑 +6 | 吃的不少 2026-04-01 | 6/300 |
|
|
|
[考研] 化學308分調(diào)劑 +13 | 你好明天你好 2026-03-30 | 14/700 |
|
|
|
[考研] 一志愿鄭大085600，310分求調(diào)劑 +6 | 李瀟可 2026-03-26 | 6/300 |
|
|
|
[考研] 348求調(diào)劑 +9 | zzzzyk123 2026-04-01 | 9/450 |
|
|
|
[考研] 0855機械初試280求調(diào)劑 +3 | kazenotori 2026-03-31 | 3/150 |
|
|
|
[考研] 070300化學354求調(diào)劑 +15 | 101次希望 2026-03-28 | 15/750 |
|
|
|
[考研] 254材料與化工求調(diào)劑 +3 | 翰冬林楠 2026-03-30 | 4/200 |
|
|
|
[考研] 266分，求材料相關(guān)專業(yè)調(diào)劑 +10 | 哇呼哼呼哼 2026-03-30 | 12/600 |
|
|
|
[考研] 哈爾濱工業(yè)大學材料與化工專碩378求調(diào)劑 +3 | 塔比烏斯 2026-03-30 | 3/150 |
|
|
|
[考研] 環(huán)境科學與工程334分求調(diào)劑 +6 | 王一一依依 2026-03-30 | 8/400 |
|
|
|
[考研] 藥學105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
|
[考研] 324求調(diào)劑 +5 | hanamiko 2026-03-26 | 5/250 |
|