 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1444 | 回復(fù): 5 | ||
yangfan832銅蟲 (初入文壇)
|
[求助]
matlab下bvp4c求解非線方程問題
|
|
請(qǐng)問各位大俠，我自己編了個(gè)用bvp4c求解器求解非線性邊值條件方程組的程序。當(dāng)方程中系數(shù)n=1時(shí)，即系數(shù)簡(jiǎn)化的情況下，程序能較好的運(yùn)行并算出較好的數(shù)值結(jié)果。當(dāng)n不等于1，程序提示出現(xiàn)奇異性，不能求解方程。請(qǐng)問這是給定初值不合適的原因嗎還是有別的問題，能幫忙解決一下嗎？謝謝！ global b t E EI eta H l r Es EIz Q EIA m n ha; b=100/10^9;l=1000/10^9;t=50/10^9;r=0.89; E=76*10^9;Es=1.22; EI=E*b*t^3/12; H=2*(t+b)*r; EIz=EI+0.5*Es*b*t^2+Es*t^3/6; EIA=E*b*t+2*Es*(b+t); m=H*l^2/EIz; %n=EIA*l^2/EIz; n=1; ha=H/EIA; eta=t/l; Q=-1; %以上是給定常數(shù)參數(shù) myinit=[0 1 0 1 0.001 -0.1];%[0 1 0 0 7 -5],[0 1 0 1 0.1 -0.1] solint=bvpinit([0 0.25 0.5 0.75 1],myinit); sol=bvp4c(@ss01,@ss02,solint); function odefun1=ss01(x,y); %給出控制方程 global eta m n; odefun1=[y(2);-eta^2*y(4)*y(5);y(4);y(5);y(6);-m*y(5)-n*(y(2)+0.5*eta^2*y(4)^2)*y(5)]; return function bcfun1=ss02(ya,yb);%給出邊界條件 global Q n eta ha; bcfun1=[ya(1);ya(4);ya(3);yb(5);yb(2)+0.5*eta^2*yb(4)^2+ha;eta*(yb(2)+0.5*eta^2*yb(4)^2)*yb(4)+eta*yb(6)/n+ha*eta*yb(4)-Q]; return |
回復(fù)此樓» 猜你喜歡
284求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿山東大學(xué)藥學(xué)學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
07化學(xué)280分求調(diào)劑
已經(jīng)有4人回復(fù)
-
298-一志愿中國(guó)農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有12人回復(fù)
-
求材料，環(huán)境專業(yè)調(diào)劑
已經(jīng)有3人回復(fù)
-
335求調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿吉大化學(xué)322求調(diào)劑
已經(jīng)有4人回復(fù)
-
環(huán)境學(xué)碩288求調(diào)劑
已經(jīng)有8人回復(fù)
-
341求調(diào)劑(一志愿湖南大學(xué)070300)
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- matlab用多元線性最小二乘法求解多元線性方程的問題
已經(jīng)有11人回復(fù)
- matlab數(shù)值求解邊界條件微分方程組
已經(jīng)有7人回復(fù)
- 求matlab求解一個(gè)方程組的問題
已經(jīng)有6人回復(fù)
- 一個(gè)微分方程組，求解
已經(jīng)有18人回復(fù)
- 求高人指點(diǎn)用matlab求解非線性方程組，解決了追加100金幣；
已經(jīng)有11人回復(fù)
- matlab的fsove 命令求解非線性方程組
已經(jīng)有6人回復(fù)
- MATLAB非線性方程組求解~~~反應(yīng)釜模擬
已經(jīng)有9人回復(fù)
- 【求助】matlab 求解微分方程中的未知參數(shù)
已經(jīng)有20人回復(fù)
- 【求助】向各位大俠求助matlab求解微分方程組遇到的一個(gè)問題
已經(jīng)有21人回復(fù)
- 【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實(shí)數(shù)解！
已經(jīng)有18人回復(fù)
- 【求助】MATLAB 有限差分法（FDM）求解偏微分方程
已經(jīng)有22人回復(fù)
- 【求助】積分微分方程matlab求解
已經(jīng)有6人回復(fù)
- 【求助】求教matlab解非線性方程組
已經(jīng)有9人回復(fù)
dingd
鐵桿木蟲 (職業(yè)作家)
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 2006-10-28
2樓2012-01-17 23:14:07
yangfan832
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 187.8
- 帖子: 26
- 在線: 103.7小時(shí)
- 蟲號(hào): 667949
- 注冊(cè): 2008-12-04
- 性別: MM
- 專業(yè): 力學(xué)
3樓2012-01-18 11:00:31
yangfan832
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 187.8
- 帖子: 26
- 在線: 103.7小時(shí)
- 蟲號(hào): 667949
- 注冊(cè): 2008-12-04
- 性別: MM
- 專業(yè): 力學(xué)
4樓2012-01-18 11:01:15
yangfan832
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 187.8
- 帖子: 26
- 在線: 103.7小時(shí)
- 蟲號(hào): 667949
- 注冊(cè): 2008-12-04
- 性別: MM
- 專業(yè): 力學(xué)
5樓2012-01-18 11:01:42
dingd
鐵桿木蟲 (職業(yè)作家)
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 2006-10-28
6樓2012-01-18 13:52:28
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 yangfan832 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料/農(nóng)業(yè)專業(yè)，07/08開頭均可，過線就行 +3 | 呵唔哦豁 2026-03-23 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 361求調(diào)劑 +3 | Glack 2026-03-22 | 3/150 |
|
|
|
[考研] 上海電力大學(xué)材料防護(hù)與新材料重點(diǎn)實(shí)驗(yàn)室招收調(diào)劑研究生（材料、化學(xué)、電化學(xué)，環(huán)境） +3 | 我愛學(xué)電池 2026-03-23 | 3/150 |
|
|
|
[考研]
 材料與化工考研調(diào)劑
+4
|
孅華 2026-03-22 | 4/200 |
|
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
|
[考研] 306求調(diào)劑 +9 | chuanzhu川燭 2026-03-18 | 9/450 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-17 | 10/500 |
|
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長(zhǎng)期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +11 | 冷笙123 2026-03-17 | 11/550 |
|
|
|
[考研] 求調(diào)劑一志愿海大，0703化學(xué)學(xué)碩304分，有大創(chuàng)項(xiàng)目，四級(jí)已過 | 幸運(yùn)哩哩 2026-03-22 | 10/500 |
|
|
|
[考研] 311求調(diào)劑 +3 | 26研0 2026-03-20 | 3/150 |
|
|
|
[考研] 化學(xué)調(diào)劑 +5 | yzysaa 2026-03-21 | 5/250 |
|
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
|
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗(yàn)） +3 | 1911236844 2026-03-17 | 3/150 |
|
|
|
[考研] 274求調(diào)劑 +10 | S.H1 2026-03-18 | 10/500 |
|
|
|
[考研] 22408 344分 求調(diào)劑 一志愿 華電計(jì)算機(jī)技術(shù) +4 | solanXXX 2026-03-20 | 4/200 |
|
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
|
[考研] 一志愿西南交通 專碩 材料355 本科雙非 求調(diào)劑 +5 | 西南交通專材355 2026-03-19 | 5/250 |
|
|
|
[考研] 一志愿西安交通大學(xué) 學(xué)碩 354求調(diào)劑211或者雙一流 +3 | 我想要讀研究生 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +11 | Ymlll 2026-03-18 | 15/750 |
|