 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1445 | 回復(fù): 5 | ||
yangfan832銅蟲 (初入文壇)
|
[求助]
matlab下bvp4c求解非線方程問題
|
|
請問各位大俠，我自己編了個用bvp4c求解器求解非線性邊值條件方程組的程序。當(dāng)方程中系數(shù)n=1時，即系數(shù)簡化的情況下，程序能較好的運行并算出較好的數(shù)值結(jié)果。當(dāng)n不等于1，程序提示出現(xiàn)奇異性，不能求解方程。請問這是給定初值不合適的原因嗎還是有別的問題，能幫忙解決一下嗎？謝謝！ 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ù)此樓» 猜你喜歡
081700 調(diào)劑 267分
已經(jīng)有8人回復(fù)
-
請教下大家 2026年國家基金申請是雙盲審嗎？
已經(jīng)有5人回復(fù)
-
276求調(diào)劑。有半年電池和半年高分子實習(xí)經(jīng)歷
已經(jīng)有10人回復(fù)
-
一志愿南航材料專317分求調(diào)劑
已經(jīng)有4人回復(fù)
-
求調(diào)劑
已經(jīng)有4人回復(fù)
-
生物學(xué)學(xué)碩求調(diào)劑
已經(jīng)有5人回復(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-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑
已經(jīng)有12人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- matlab用多元線性最小二乘法求解多元線性方程的問題
已經(jīng)有11人回復(fù)
- matlab數(shù)值求解邊界條件微分方程組
已經(jīng)有7人回復(fù)
- 求matlab求解一個方程組的問題
已經(jīng)有6人回復(fù)
- 一個微分方程組，求解
已經(jīng)有18人回復(fù)
- 求高人指點用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求解微分方程組遇到的一個問題
已經(jīng)有21人回復(fù)
- 【求助】matlab求解非線性方程組，并畫圖處理。要求y,z是實數(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小時
- 蟲號: 291104
- 注冊: 2006-10-28
2樓2012-01-17 23:14:07
yangfan832
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 187.8
- 帖子: 26
- 在線: 103.7小時
- 蟲號: 667949
- 注冊: 2008-12-04
- 性別: MM
- 專業(yè): 力學(xué)
3樓2012-01-18 11:00:31
yangfan832
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 187.8
- 帖子: 26
- 在線: 103.7小時
- 蟲號: 667949
- 注冊: 2008-12-04
- 性別: MM
- 專業(yè): 力學(xué)
4樓2012-01-18 11:01:15
yangfan832
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 187.8
- 帖子: 26
- 在線: 103.7小時
- 蟲號: 667949
- 注冊: 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小時
- 蟲號: 291104
- 注冊: 2006-10-28
6樓2012-01-18 13:52:28
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 yangfan832 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +10 | 梁梁校校 2026-03-17 | 10/500 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)求調(diào)劑 +4 | 奶油草莓. 2026-03-22 | 5/250 |
|
|
|
[考研] 333求調(diào)劑 +3 | ALULU4408 2026-03-23 | 3/150 |
|
|
|
[考研] 352求調(diào)劑 +3 | 大米飯！ 2026-03-22 | 3/150 |
|
|
|
[考研] 311求調(diào)劑 +6 | 冬十三 2026-03-18 | 6/300 |
|
|
|
[考研] 307求調(diào)劑 +11 | 冷笙123 2026-03-17 | 11/550 |
|
|
|
[考研] 306求調(diào)劑 +5 | 來好運來來來 2026-03-22 | 5/250 |
|
|
|
[考研] 289求調(diào)劑 +7 | 懷瑾握瑜l 2026-03-20 | 7/350 |
|
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
|
[考研] 一志愿深大，0703化學(xué)，總分302，求調(diào)劑 +4 | 七月-七七 2026-03-21 | 4/200 |
|
|
|
[考研] 297求調(diào)劑 +11 | 戲精丹丹丹 2026-03-17 | 12/600 |
|
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[基金申請]
學(xué)校已經(jīng)提交到NSFC，還能修改嗎？
 40+4
|
babangida 2026-03-19 | 9/450 |
|
|
|
[考研] 268求調(diào)劑 +9 | 簡單點0 2026-03-17 | 9/450 |
|
|
|
[考研] 求調(diào)劑 +3 | Ma_xt 2026-03-17 | 3/150 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +5 | 材化逐夢人 2026-03-18 | 5/250 |
|
|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8 | 小材化本科 2026-03-18 | 8/400 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +4 | xbxudjdn 2026-03-18 | 4/200 |
|
|
|
[考研] 【同濟軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|