24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2345  |  回复: 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的回帖

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的回帖
查看全部 11 个回答

buhely

银虫 (小有名气)

【答案】应助回帖

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

jjliu2012

新虫 (小有名气)

比如 dx/dt=f(x(t)),这里的f是一个非线性的函数。而fun=int(x^2,a,b),这种情况x(t)没有解析表达式
5楼2013-09-26 15:19:50
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂 +13 zchqwer 2026-04-01 15/750 2026-04-01 22:30 by 2026材料调剂
[考研] 材料调剂 +11 一样YWY 2026-03-31 11/550 2026-04-01 22:25 by zhouyuwinner
[考研] 296求调剂 +4 汪!?! 2026-03-31 7/350 2026-04-01 22:04 by 客尔美德
[考研] 284求调剂 +12 小熊~~ 2026-03-31 12/600 2026-04-01 20:23 by 花??
[考研] 0710生物学考研调剂 +3 李多米lee. 2026-03-27 4/200 2026-04-01 16:21 by zzchen2000
[考研] 318一志愿吉林大学生物与医药 求调剂 +6 笃行致远. 2026-03-28 6/300 2026-04-01 09:28 by oooqiao
[考研] 335求调剂 +3 321* 2026-03-31 4/200 2026-04-01 00:00 by 321*
[考研] 合肥区域性重点一本招收调剂 +4 6266jl 2026-03-30 8/400 2026-03-31 18:43 by 6266jl
[考研] 299求调剂 +8 嗯嗯嗯嗯2 2026-03-27 8/400 2026-03-31 18:23 by lizhi8172
[考研] 322求调剂:一志愿湖南大学 材料与化工(085600),已过六级。 +10 XX小邓 2026-03-29 10/500 2026-03-31 16:46 by 不吃魚的貓
[考研] 英一数一总分334求调剂 +4 陈阳坤 2026-03-31 4/200 2026-03-31 14:22 by 记事本2026
[考研] 298求调剂 +3 种圣赐 2026-03-29 3/150 2026-03-29 12:06 by longlotian
[考研] 086000生物与医药调剂 +5 Feisty。 2026-03-28 9/450 2026-03-29 12:02 by longlotian
[考研] 调剂考研 +3 王杰一 2026-03-29 3/150 2026-03-29 08:09 by fmesaito
[硕博家园] 招收生物学/细胞生物学调剂 +4 IceGuo 2026-03-26 5/250 2026-03-29 01:25 by griffith2014
[考研] 求调剂 +7 争取九点睡 2026-03-28 8/400 2026-03-28 21:07 by 争取九点睡
[考研] 085405 考的11408求各位老师带走 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 09:19 by 乐呵呵的追梦人
[考研] 材料求调剂一志愿哈工大324 +7 闫旭东 2026-03-28 9/450 2026-03-28 08:51 by Xu de nuo
[考研] 330一志愿中国海洋大学 化学工程 085602 有读博意愿 求调剂 +3 wywy.. 2026-03-27 4/200 2026-03-28 03:32 by fmesaito
[考研] 298调剂 +3 jiyingjie123 2026-03-27 3/150 2026-03-27 11:57 by wxiongid
信息提示
请填处理意见