24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2346  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

jjliu2012

新虫 (小有名气)

[求助] 新手matlab求助,谢谢各位

最优控制中的最速下降法,若是求不出函数几乎函数梯度的解析表达式,最速下降法的程序该怎么样处理呢?
以下是我们平时用的最速下降法程序

function [x,val,k]=grad(fun,gfun,x0)
%功能: 用最速下降法求解无约束问题: min f(x)
%输入: x0是初始点, fun, gfun分别是目标函数和梯度
%输出: x, val分别是近似最优点和最优值, k是迭代次数.
maxk=5000; %最大迭代次数
rho=0.5;sigma=0.4;
k=0; epsilon=1e-5;
while(k¡maxk)
g=feval(gfun,x0); %计算梯度
d=-g; %计算搜索方向
if(norm(d)¡epsilon), break; end
m=0; mk=0;
while(m¡20) %Armijo搜索
if(feval(fun,x0+rho^m*d)¡feval(fun,x0)+sigma*rho^m*g’*d)
mk=m; break;
end
m=m+1;
end
x0=x0+rho^mk*d;
k=k+1;
end
x=x0;
val=feval(fun,x0);

%这里的fun是函数,f是函数,g是f函数的梯度,这里的f和g都有解析的表达式
function f=fun(x)
f=100*(x(1)^2-x(2))^2+(x(1)-1)^2;
function g=gfun(x)
g=[400*x(1)*(x(1)^2-x(2))+2*(x(1)-1), -200*(x(1)^2-x(2))]’;
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
你这里fun是定积分,不是t的函数,换个例子来吧

PS:如果换成不定积分fun(t)=int[x(s)^2.ds,a,t], 那么gfun(t)=x(t)^2=int[f,a0,t]^2
7楼2013-09-26 20:50:51
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

buhely

银虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
如果你的fun是用明确的数学语言描述的,但不能直接算出gfun,可以考虑用matlab的符号变量编程,然后算出gfun的表达式;
再不行貌似就只能再添加各种启发式搜素算法了
2楼2013-09-26 10:00:19
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

jjliu2012

新虫 (小有名气)

引用回帖:
2楼: Originally posted by buhely at 2013-09-26 10:00:19
如果你的fun是用明确的数学语言描述的,但不能直接算出gfun,可以考虑用matlab的符号变量编程,然后算出gfun的表达式;
再不行貌似就只能再添加各种启发式搜素算法了

这里的x满足的是一个非线性的微分方程,x就没有解析解,fun是关于x的,也没有解析表达式啊
3楼2013-09-26 10:23:53
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)

引用回帖:
3楼: Originally posted by jjliu2012 at 2013-09-25 21:23:53
这里的x满足的是一个非线性的微分方程,x就没有解析解,fun是关于x的,也没有解析表达式啊...

你能不能把你得问题说得再明白点,贴的程序貌似和你的问题也不相关
4楼2013-09-26 11:24:25
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 385分 生物学(071000)求调剂 +4 qf626 2026-04-01 4/200 2026-04-02 06:44 by ilovexiaobin
[考研] 270求调剂 +7 小杰pp 2026-03-31 8/400 2026-04-01 21:36 by 七度不信任
[考研] 310分求调剂 +4 成功上岸wang 2026-04-01 4/200 2026-04-01 20:35 by liu823948201
[考研] 286求调剂 +5 Sa67890. 2026-04-01 7/350 2026-04-01 19:50 by 6781022
[考研] 一志愿085600中科院宁波所276分求调剂 +22 材料学257求调剂 2026-03-28 23/1150 2026-04-01 19:03 by 客尔美德
[考研] 314求调剂 +4 溪云珂 2026-03-26 4/200 2026-04-01 17:00 by oooqiao
[考研] 311求调剂 +11 蓝月亮亮 2026-03-30 11/550 2026-04-01 16:33 by 七度不信任
[考研] 348求调剂 +9 zzzzyk123 2026-04-01 9/450 2026-04-01 14:24 by w虫虫123
[考研] 一志愿中农0710生物学,微生物方向总分338求调剂 +3 柒xxxx. 2026-03-26 3/150 2026-04-01 12:30 by 冰乌龙
[考研] 土木304求调剂 +5 兔突突突, 2026-03-31 6/300 2026-04-01 09:37 by JourneyLucky
[考研] 【调剂】一志愿厦大生物与医药调剂 +3 Echo虾米 2026-03-31 3/150 2026-04-01 08:40 by JourneyLucky
[考研] 一志愿西交大080500材料学硕349 +6 jqx1258 2026-03-31 7/350 2026-03-31 21:08 by yuq
[考研] 334求调剂 +7 Trying] 2026-03-31 7/350 2026-03-31 12:33 by 无际的草原
[考研] 一志愿浙江大学工科动力工程370,数一121,专业课135,现在能去哪里 +3 080700调剂 2026-03-30 4/200 2026-03-31 12:00 by KLMY666
[考研] 285求调剂 +6 AZMK 2026-03-29 9/450 2026-03-30 21:02 by dophin1985
[考研] 322求调剂 +10 宋明欣 2026-03-27 10/500 2026-03-30 18:47 by 544594351
[考研] 356求调剂 +3 gysy?s?a 2026-03-28 3/150 2026-03-29 00:33 by 544594351
[考研] 材料与化工(0856)304求B区调剂 +8 邱gl 2026-03-27 8/400 2026-03-28 12:42 by 唐沐儿
[考研] 330一志愿中国海洋大学 化学工程 085602 有读博意愿 求调剂 +3 wywy.. 2026-03-27 4/200 2026-03-28 03:32 by fmesaito
[考研] 调剂 +3 李嘉图·S·路 2026-03-27 3/150 2026-03-27 11:19 by wangjy2002
信息提示
请填处理意见