| 2 | 1/1 | 返回列表 |
| 查看: 393 | 回復: 1 | |||
missqiushui
|
[交流]
matlab牛頓法如何調用問題,調用部分總是出錯 已有1人參與
|
|
牛頓法的程序: function [minx,minf]=example3g(f,x0,var,eps) %目標函數(shù):f %初始點:x0 %自變量向量var %精度:eps %目標函數(shù)取最小值時的自變量值:minx; %目標函數(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)是計算標量函數(shù)或矢量函數(shù)f關于矢量v的雅可比矩陣。 %其結果就相當于df/dv。需要注意的是,當f是標量時,就相當于是求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)); 調用部分,不知道var部分如何初始化 syms x x0; x0=[300 200 300 10]; [minx,minf]=example3g(op,x0,[x(1) x(2) x(3) x(4)],1e-6) |
銅蟲 (初入文壇)

| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿985,本科211,0817化學工程與技術319求調劑 +8 | Liwangman 2026-03-15 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 274求調劑 +6 | S.H1 2026-03-18 | 6/300 |
|
|
[考研] 本科鄭州大學物理學院,一志愿華科070200學碩,346求調劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
[考研] 085410人工智能專碩317求調劑(0854都可以) +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
[考研] 311求調劑 +4 | 冬十三 2026-03-18 | 4/200 |
|
|
[考研] 一志愿西南交大,求調劑 +4 | 材化逐夢人 2026-03-18 | 4/200 |
|
|
[考研] 070300化學319求調劑 +6 | 錦鯉0909 2026-03-17 | 6/300 |
|
|
[考研] 0703化學調劑 +3 | 妮妮ninicgb 2026-03-17 | 3/150 |
|
|
[考研] 301求調劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考博] 26申博 +4 | 八6八68 2026-03-16 | 4/200 |
|
|
[考研] 一志愿南京大學,080500材料科學與工程,調劑 +4 | Jy? 2026-03-16 | 4/200 |
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
[考研] 283求調劑 +3 | 聽風就是雨; 2026-03-16 | 3/150 |
|
|
[考研] 中科院材料273求調劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
[考研] 0703化學調劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 0856專碩279求調劑 +5 | 加油加油!? 2026-03-15 | 5/250 |
|
|
[考研] 22408總分284求調劑 +3 | InAspic 2026-03-13 | 3/150 |
|
|
[考研] 304求調劑 +6 | Mochaaaa 2026-03-12 | 7/350 |
|
|
[考研] 289求調劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|