| 1 | 1/1 | 返回列表 |
| 查看: 2570 | 回復: 0 | ||
colored銅蟲 (初入文壇)
|
[求助]
關于 fsolve 函數(shù)的使用
|
|
樓主現(xiàn)在正在做一個仿真,遇到一些困難,請各位路過的高手幫幫忙! 問題如下:主函數(shù)產(chǎn)生六個參數(shù)lamdad, miud, etad,lamdae,miue,etae,function 需要這六個參數(shù)中的一些,請問要怎么讀取這些值,而且在主函數(shù)中還要用fsolve函數(shù)來求解一個三元的非線性方程組(方程組為F)。 ------------------------------------------------------------------------------------------------------ 主函數(shù): for ii = 1 : cd +dd [lamdad, miud, etad] = fun_myunified(ii, ad, bd, cd, dd); for jj = 1 : ce + de [lamdae,miue,etae] = fun_myunified(jj, ae, be, ce, de); ------------------------------------------------------------------------------------------------------ function output = fun(*,*,*) F = [x(1)*gamma(x(2)+1/x(3))/(x(2)^(1/x(3))*gamma(x(2))) - 2* miue.*miud .* (betae./etad).^((miud-miue)./2) .* average_gg((miud+miue)/2 , miud-miue, xi); x(1)^2*gamma(x(2)+2/x(3))/(x(2)^(2/x(3))*gamma(x(2))) - 2* miue.*miud .* (betae./etad).^((miud-miue)./2) .* average_gg((miud+miue)/2+1 , miud-miue, xi); x(1)^3*gamma(x(2)+3/x(3))/(x(2)^(3/x(3))*gamma(x(2))) - 2* miue.*miud .* (betae./etad).^((miud-miue)./2) .* average_gg((miud+miue)/2+2 , miud-miue, xi)]; ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- 在主函數(shù)中又需要用fsolve解出來的三個參數(shù)。 ------------------------------------------------------------------------------------------------------- 主函數(shù): for ii = 1 : cd +dd [lamdad, miud, etad] = fun_myunified(ii, ad, bd, cd, dd); for jj = 1 : ce + de [lamdae,miue,etae] = fun_myunified(jj, ae, be, ce, de); x0 = [1; 1; 1]; % Make a starting guess at the solution options = optimoptions('fsolve','Display','iter'); % Option to display output [x1,fval1] = fsolve(@fun,x0,options); % Call solver --------------------------------------------------------------------------------------------------------- 謝謝! |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 302求調(diào)劑 +6 | 呼呼呼。。。。 2026-03-17 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +5 | 慕寒mio 2026-03-16 | 5/250 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 085600材料與化工求調(diào)劑 +5 | 緒幸與子 2026-03-17 | 5/250 |
|
|
[碩博家園] 深圳大學碩士招生(2026秋,傳感器方向,僅錄取第一志愿) +4 | xujiaoszu 2026-03-11 | 9/450 |
|
|
[基金申請]
今年的國基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
[考研] 080500,材料學碩302分求調(diào)劑學校 +4 | 初識可樂 2026-03-14 | 5/250 |
|
|
[考研] 265求調(diào)劑 +4 | 威化餅07 2026-03-12 | 4/200 |
|
|
[基金申請] 現(xiàn)在如何回避去年的某一個專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機器學習 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
[考研] 290求調(diào)劑 +9 | ADT 2026-03-11 | 9/450 |
|
|
[考研] 26調(diào)劑/材料科學與工程/總分295/求收留 +9 | 2026調(diào)劑俠 2026-03-12 | 9/450 |
|
|
[考研] 材料與化工085600調(diào)劑求老師收留 +9 | jiaanl 2026-03-11 | 9/450 |
|
|
[考研] 求調(diào)劑 +7 | 18880831720 2026-03-11 | 7/350 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|