 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 394 | 回復(fù): 1 | |||
missqiushui
|
[交流]
matlab牛頓法如何調(diào)用問題，調(diào)用部分總是出錯 已有1人參與
|
|
牛頓法的程序： function [minx,minf]=example3g(f,x0,var,eps) %目標(biāo)函數(shù)：f %初始點：x0 %自變量向量var %精度：eps %目標(biāo)函數(shù)取最小值時的自變量值：minx; %目標(biāo)函數(shù)的最小值：minf format long; if nargin==4 eps=1.0e-6; end tol=1; x0=transpose(x0); while tol>eps gradutef=jacobian(f,var); %梯度方向 %其大概意思是：函數(shù)jacobian就是雅可比矩陣, %命令jacobian（f，v）是計算標(biāo)量函數(shù)或矢量函數(shù)f關(guān)于矢量v的雅可比矩陣。 %其結(jié)果就相當(dāng)于df/dv。需要注意的是，當(dāng)f是標(biāo)量時，就相當(dāng)于是求f的梯度。 jacf=Jacobian(gradutef,var); %雅克比矩陣 v=Funval(gradutef,var,x0); tol=norm(v); pv=Funval(jacf,var,x0); p=-inv(pv)*transpose(v); %搜索方向 x1=x0+p; x0=x1; end minx=x1; minf=Funval(f,var,minx); format short; 想優(yōu)化的函數(shù)： function f=op(x) f=561+7.92*x(1)+0.001562*x(1)^2+310+7.85*x(2)+0.00194*x(2)^2+78+7.97*x(3)+0.00482*x(3)^2+x(4)*(800-x(1)-x(2)-x(3)); 調(diào)用部分，不知道var部分如何初始化 syms x x0; x0=[300 200 300 10]; [minx,minf]=example3g(op,x0,[x(1) x(2) x(3) x(4)],1e-6) |
回復(fù)此樓» 猜你喜歡
材料考研調(diào)劑
已經(jīng)有3人回復(fù)
-
346求調(diào)劑[0856]
已經(jīng)有6人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有18人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有9人回復(fù)
-
288求調(diào)劑
已經(jīng)有4人回復(fù)
-
328求調(diào)劑，英語六級551，有科研經(jīng)歷
已經(jīng)有12人回復(fù)
-
一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑
已經(jīng)有5人回復(fù)
-
085601材料工程專碩求調(diào)劑
已經(jīng)有7人回復(fù)
-
321求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿中海洋材料工程專碩330分求調(diào)劑
已經(jīng)有7人回復(fù)
48038074
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 147
- 散金: 2
- 帖子: 37
- 在線: 38.4小時
- 蟲號: 4032798
- 注冊: 2015-08-21
- 性別: GG
- 專業(yè): 電力系統(tǒng)
2樓2015-11-23 00:40:35
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 missqiushui 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 294求調(diào)劑材料與化工專碩 +9 | 陌の森林 2026-03-18 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 304求調(diào)劑 +6 | 司空. 2026-03-18 | 6/300 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +4 | 材化逐夢人 2026-03-18 | 4/200 |
|
|
|
[考研] 材料專碩274一志愿陜西師范大學(xué)求調(diào)劑 +6 | 薛云鵬 2026-03-13 | 6/300 |
|
|
|
[考研] 312求調(diào)劑 +8 | 陌宸希 2026-03-16 | 9/450 |
|
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨； 2026-03-16 | 3/150 |
|
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請 +6 | zlingli 2026-03-13 | 6/300 |
|
|
|
[考研] 289求調(diào)劑 +4 | 這么名字咋樣 2026-03-14 | 6/300 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|
|
|
[考研] 266求調(diào)劑 +4 | 學(xué)員97LZgn 2026-03-13 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +8 | 困于星晨 2026-03-12 | 10/500 |
|
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|