 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 3267 | 回復(fù): 0 | ||
Figo1937銀蟲(chóng) (初入文壇)
|
[求助]
bvp4c解邊界值微分方程
|
|
用MATLAB的bvp4c遇到“奇異性雅可比行列式”無(wú)法求解！報(bào)錯(cuò)： 錯(cuò)誤使用 bvp4c (line 251) 無(wú)法求解排列方程式 - 遇到了奇異性雅可比行列式。 出錯(cuò) mybvp1>exmpbh_solve (line 28) sol=bvp4c(@ode,@bc,s0); 出錯(cuò) mybvp1 (line 13) S=exmpbh_solve(s0,ph0,phi); 代碼如下： function S=mybvp1 x0=linspace(0,1,10); %自變量初始網(wǎng)格點(diǎn)數(shù)組 根據(jù)實(shí)際修改自變量上下線0~10,及網(wǎng)格點(diǎn)個(gè)數(shù)5 v0=[0.7,0,0.68,0.0,0.0]; %各網(wǎng)格點(diǎn)處初值猜測(cè) 若自變量n等分，微分方程組有m個(gè)，則v0為m*（n+1）的二維數(shù)組 %第一行n的初值： %第二行y的初值： %第三行y\'的初值： %第四行y"的初值： %第五行y"\'的初值： phi=0.243; %已知參數(shù) d0=0.0; %參數(shù)1初值猜測(cè) ph0=d0; s0=bvpinit(x0,v0,ph0); %產(chǎn)生初始猜測(cè)構(gòu)架 S=exmpbh_solve(s0,ph0,phi); end function sol=exmpbh_solve(s0,ph,q) %q用來(lái)傳遞已知參量 ph為d，q為phi phi =q; %————已知量———— D=0.28*10^(-3);L=8*10^(-3);alpha=0.5; %表一 E=1.15*10^11;mu=0.7;I=7.7*10^(-9); %表二 kx=1.65*10^8;p0=1.05*10^4;p1=375;f=1.5*10^3; %表三 %————代入式———— beta=L/sqrt(2)*(kx/E/I)^(1/4); pi1=(p0+(p1-p0)*exp(-alpha*f*L*ph))/kx/L; pi2=p1/kx/L; sol=bvp4c(@ode,@bc,s0); function dydx=ode(x,y,ph) %以下微分方程公式中所有phi用q替代;d用ph替代;y用y(2)替代;y的二階導(dǎo)數(shù)用y(4)替代 dydx=zeros(5,1); %確保為列向量 %————條件1———— if x>=0 && x<=alpha*(1-ph) dydx(1)=4*beta^4*(pi1+mu*(D/L)*abs(y(2)));% n\' dydx(2)=y(3);%y\' dydx(3)=y(4);%y\'\' dydx(4)=y(5);%y\'\'\' dydx(5)=y(1)*y(4)-4*beta^4*y(2);%y\'\'\'\' %————條件2———— elseif x>alpha*(1-ph) && x<alpha dydx(1)=0; dydx(2)=y(3); dydx(3)=y(4); dydx(4)=y(5); dydx(5)=y(1)*y(4); %————條件3———— elseif x>=alpha && x<1 dydx(1)=4*beta^4*(pi2+mu*(D/L)*abs(y(2)-alpha*phi*ph/sqrt(1-phi*phi))); dydx(2)=y(3); dydx(3)=y(4); dydx(4)=y(5); dydx(5)=y(1)*y(4)-4*beta^4*(y(2)-alpha*phi*ph/sqrt(1-phi*phi)*L/D); end end function res=bc(ya,yb,ph) %ya表示自變量取初值時(shí)狀態(tài)量的值，yb表示自變量取末值時(shí)的狀態(tài)量的值 res=zeros(6,1); %確保為列向量 res(1)=ya(1); %n\'(0)=0 res(2)=ya(2); %y（0）=0 res(3)=yb(2)-alpha*phi*ph/sqrt(1-phi*phi)*(L/D);%y（1）=... res(4)=yb(3); %y\'(1)=0 res(5)=ya(4); %y"(0)=0 end end mx326D4.png |
回復(fù)此樓» 猜你喜歡
材料科學(xué)與工程求調(diào)劑
已經(jīng)有8人回復(fù)
-
一志愿大連理工大學(xué)材料求調(diào)劑
已經(jīng)有3人回復(fù)
-
一志愿北化085600材料專(zhuān)碩275|有文章專(zhuān)利｜求調(diào)劑
已經(jīng)有9人回復(fù)
-
材料求調(diào)劑 一志愿哈工大總分298分，前三科223分
已經(jīng)有8人回復(fù)
-
274求調(diào)劑
已經(jīng)有3人回復(fù)
-
合肥區(qū)域性重點(diǎn)一本招收調(diào)劑
已經(jīng)有4人回復(fù)
-
085600材料與化工調(diào)劑
已經(jīng)有10人回復(fù)
-
22408 359分調(diào)劑
已經(jīng)有7人回復(fù)
-
抱歉
已經(jīng)有4人回復(fù)
-
085602化工求調(diào)劑（331分）
已經(jīng)有8人回復(fù)
找到一些相關(guān)的精華帖子，希望有用哦~
邊界條件微分方程用bvp4c求解，猜測(cè)值不同有不同曲線解的原因以及邊界條件明顯不符
已經(jīng)有4人回復(fù)
- matlab求解邊界值問(wèn)題
已經(jīng)有8人回復(fù)
- MATLAB求解數(shù)值微分方程的問(wèn)題
已經(jīng)有3人回復(fù)
- 求高手幫忙編一個(gè)matlab的求解參數(shù)的語(yǔ)句
已經(jīng)有6人回復(fù)
- Matlab數(shù)值求解二階常微分方程
已經(jīng)有9人回復(fù)
- MATLAB 解偏微分方程時(shí)的邊界條件傳遞問(wèn)題
已經(jīng)有4人回復(fù)
- CFX求解時(shí)程序出錯(cuò)問(wèn)題（急�。。。�
已經(jīng)有5人回復(fù)
- 請(qǐng)教MATLAB中的LMI求解問(wèn)題，急�。�！
已經(jīng)有7人回復(fù)
- 求高手解答Matlab求解微分方程的初始值和邊界值問(wèn)題
已經(jīng)有5人回復(fù)
- Matlab: 用bvp解二階微分方程（邊界值）遇到困難，請(qǐng)高手指點(diǎn)
已經(jīng)有5人回復(fù)
- matlab數(shù)值求解邊界條件微分方程組
已經(jīng)有7人回復(fù)
- 解微分方程邊界條件問(wèn)題：什么是open boundary conditions
已經(jīng)有1人回復(fù)
點(diǎn)擊這里搜索更多相關(guān)資源科研從小木蟲(chóng)開(kāi)始，人人為我，我為人人
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 Figo1937 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 274求調(diào)劑 +3 | xiao愛(ài)同學(xué) 2026-03-30 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 08工科求調(diào)劑286 +4 | tgs_001 2026-03-28 | 4/200 |
|
|
|
[考研] 一志愿鄭大材料工程290求調(diào)劑 +11 | Youth_ 2026-03-30 | 11/550 |
|
|
|
[考研] 一志愿：西北大學(xué)，英一數(shù)一408-284分求調(diào)劑 +5 | 12.27 2026-03-27 | 5/250 |
|
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +14 | 北道巷 2026-03-24 | 14/700 |
|
|
|
[考研] 材料專(zhuān)碩調(diào)劑 +11 | 椰椰。 2026-03-29 | 11/550 |
|
|
|
[考研] 求調(diào)劑 +7 | 青春裁為三截 2026-03-29 | 7/350 |
|
|
|
[考研] 南京大學(xué)化學(xué)調(diào)劑 +10 | 景隨風(fēng) 2026-03-29 | 15/750 |
|
|
|
[考研] 332求調(diào)劑 +6 | Lyy930824@ 2026-03-29 | 6/300 |
|
|
|
[考研] 305求調(diào)劑 +8 | RuiFairyrui 2026-03-28 | 8/400 |
|
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 266分，求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
|
[考研] 286求調(diào)劑 +12 | PolarBear11 2026-03-26 | 12/600 |
|
|
|
[考研] 081200-314 +3 | LILIQQ 2026-03-27 | 4/200 |
|
|
|
[考研] 286求調(diào)劑 +4 | 丟掉懶惰 2026-03-27 | 7/350 |
|
|
|
[考研] 復(fù)試調(diào)劑，一志愿南農(nóng)083200食品科學(xué)與工程 +5 | XQTJZ 2026-03-26 | 5/250 |
|
|
|
[考研] 085600，材料與化工321分調(diào)劑 +4 | 大饞小子 2026-03-27 | 6/300 |
|
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
|
[考研] 一志愿吉林大學(xué)材料與化工303分求調(diào)劑 +4 | 為學(xué)666 2026-03-24 | 4/200 |
|