 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 1 | 1/1 | 返回列表 |
| 查看: 539 | 回復: 0 | ||
engine9817銅蟲 (小有名氣)
|
[求助]
求設計一個被積分函數(shù)中需要解方程的積分程序
|
|
如果有一個被積方程為G(x,y),積分區(qū)間假設為（1，2） 其中y的值與x是由另一個函數(shù)決定如f(x,y)=0; 這個式子不能直接得到y(tǒng)的解析式，復雜點的話，可能一個x有多個y值，但只有一個合適，如正的或?qū)嵉摹?br /> 這樣在求積G(x,y)時，對應于每個x都要求解一個y值，還可能要對y值進行篩選。 我一開始在function文件里面設計了G(x,y)的程序.方便起見，我們假設 f(x,y)=x-y^2=0; 假設這里把y解析表示，要求用fsolve解出y function integrandfun=G(x) y=fsolve(@(t) x-t.^2, 1); integrandfun=x.^2+y; 主函數(shù)中調(diào)用 這個方程積分如 In=quadgk(@(x) integrandfun, 0,1); 結(jié)果是： 3.558078204809361與直接積分x.^2+sqrt(x)的結(jié)果3.552284749830793有差別 提示： 警告: Trust-region-dogleg algorithm of FSOLVE cannot handle non-square systems; using Levenberg-Marquardt algorithm instead. > In fsolve at 285 In fkk at 3 In @(x)fkk(x) In quadgk>evalFun at 330 In quadgk>f1 at 348 In quadgk>vadapt at 249 In quadgk at 188 In figures at 12 No solution found. fsolve stopped because the last step was ineffective. However, the vector of function values is not near zero, as measured by the default value of the function tolerance. <stopping criteria details> 請問有進一步的優(yōu)化嗎？事實上如果把f(x,y)變得復雜點，積分存在更多問題。 |
回復此樓» 猜你喜歡
新疆大學地質(zhì)與礦業(yè)工程學院招生
已經(jīng)有32人回復
-
0856材料與化工353分求調(diào)劑
已經(jīng)有4人回復
-
282分材料專業(yè)求調(diào)劑院校
已經(jīng)有13人回復
-
315食品工程求調(diào)劑
已經(jīng)有5人回復
-
環(huán)境調(diào)劑
已經(jīng)有3人回復
-
接受26屆調(diào)劑生
已經(jīng)有20人回復
-
考研材料與化工，求調(diào)劑
已經(jīng)有3人回復
-
歡迎211本科同學，過A區(qū)國家線，A區(qū)非偏遠一本，交叉學科課題組
已經(jīng)有35人回復
-
304求調(diào)劑（085602一志愿985）
已經(jīng)有4人回復
-
294 英二數(shù)二物化 求調(diào)劑
已經(jīng)有6人回復
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 engine9817 的主題更新
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料調(diào)劑，307分 +8 | 張泳銘1 2026-03-09 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 調(diào)劑 +4 | 調(diào)劑的考研學生 2026-03-09 | 4/200 |
|
|
[考研] 334求調(diào)劑 +5 | Trying] 2026-03-07 | 5/250 |
|
|
|
[考研] 【求調(diào)劑】293分環(huán)境工程求調(diào)劑材料/化工，服從調(diào)劑，抗壓能力強！ +13 | xiiiia 2026-03-04 | 14/700 |
|
|
|
[考研] 考研調(diào)劑，一志愿山東大學材料與化工，328分，政治51 +5 | 關(guān)你西紅柿929 2026-03-08 | 6/300 |
|
|
|
[考研] 313求調(diào)劑 +4 | Yyt楊1 2026-03-07 | 5/250 |
|
|
|
[考博] 求材料讀博院校 +8 | yanglei131 2026-03-08 | 8/400 |
|
|
|
[考研] 一志愿武理314求調(diào)劑 +4 | （￣～￣;） 2026-03-08 | 5/250 |
|
|
|
[考博] 26年博士申請 +4 | 科研狗111 2026-03-07 | 4/200 |
|
|
|
[考研] 286求調(diào)劑 +12 | Faune 2026-03-06 | 14/700 |
|
|
|
[考研] 303求調(diào)劑 +8 | forgman95 2026-03-05 | 10/500 |
|
|
|
[考研] 化學290求調(diào)劑 +3 | 兩顆 西柚 2026-03-07 | 4/200 |
|
|
|
[考研] 0703 東華大學 理學化學均過 a 線，總分 281 求調(diào)劑 +6 | 香香 Lu 2026-03-03 | 11/550 |
|
|
|
[考研] 材料085601一志愿哈工大317 +4 | 壓迫感行 2026-03-04 | 4/200 |
|
|
|
[考研] 材料調(diào)劑 +4 | L9370 2026-03-05 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-04 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +4 | xuyuan111 2026-03-04 | 5/250 |
|
|
|
[考研] 347分電子信息專碩求調(diào)劑 +3 | 咕嘟咕嘟1 2026-03-03 | 3/150 |
|
|
|
[考研] 322,求調(diào)劑 +3 | 菜菜愛玩 2026-03-04 | 3/150 |
|
|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|