24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 2343  |  回复: 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材料调剂
[考研] 材料求调剂 +8 呢呢妮妮 2026-04-01 10/500 2026-04-01 22:24 by salamander`
[考研] 379求调剂 +3 ?苦瓜不苦 2026-04-01 3/150 2026-04-01 20:09 by 暮云清寒
[考研] 085600,材料与化工321分,求调剂 +11 大馋小子 2026-03-27 11/550 2026-04-01 16:10 by cymywx
[考研] 336求调剂 +3 kiyy 2026-04-01 3/150 2026-04-01 12:46 by wxiongid
[考研] 土木304求调剂 +5 兔突突突, 2026-03-31 6/300 2026-04-01 09:37 by JourneyLucky
[考研] 0856材料与化工调剂,339 +10 10213207 2026-03-31 10/500 2026-04-01 08:28 by Linda Hu
[考研] 求调剂:085600材料与化工,考材科基,总分319 +17 678lucky 2026-03-31 21/1050 2026-04-01 01:40 by 1018329917
[考研] 一志愿:西北大学,英一数一408-284分求调剂 +7 12.27 2026-03-27 7/350 2026-03-31 21:59 by lbsjt
[考研] 材料工程085601数二英一335求调剂 +5 双马尾痞老板2 2026-03-31 5/250 2026-03-31 19:07 by Wang200018
[考研] 合肥区域性重点一本招收调剂 +4 6266jl 2026-03-30 8/400 2026-03-31 18:43 by 6266jl
[考研] 085601英二数二求调剂 总分325 +4 余航航 2026-03-31 4/200 2026-03-31 17:38 by 唐沐儿
[考研] 080500-315分复试调剂 +9 上岸3821 2026-03-31 9/450 2026-03-31 17:29 by 唐沐儿
[考研] 生医工0831调剂求推荐 +4 小熊睿睿_s 2026-03-27 6/300 2026-03-31 17:12 by 记事本2026
[考研] 289求调剂 +3 Acesczlo 2026-03-29 4/200 2026-03-31 14:48 by 热情沙漠
[考研] 327求调剂 +5 小卡不卡. 2026-03-29 5/250 2026-03-30 19:30 by Wang200018
[考研] 0703化学求调剂 +6 丹青奶盖 2026-03-26 8/400 2026-03-30 18:33 by 探123
[考研] 328求调剂 +8 嗯滴的基本都 2026-03-27 8/400 2026-03-30 17:20 by Wang200018
[考研] 298调剂 +3 jiyingjie123 2026-03-27 3/150 2026-03-27 11:57 by wxiongid
[考研] 考研调剂 +9 小蜡新笔 2026-03-26 9/450 2026-03-27 11:10 by 不吃魚的貓
信息提示
请填处理意见