24小时热门版块排行榜    

查看: 3442  |  回复: 9

893473348

银虫 (初入文坛)

[求助] 有关L1自适应控制投影算子Proj的Matlab编程问题

求助:有关L1自适应控制投影算子Proj的Matlab编程问题。现在需要编程实现下图所示的算子,编程错误,不能够保证输出量是有界的而导致控制系统发散,急求牛人指点迷津!非常感谢!!!
有关L1自适应控制投影算子Proj的Matlab编程问题
Proj.jpg
回复此楼
Nothing is impossible!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunshine819

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★
1592203609: 金币+1, 应助指数+1, 谢谢回帖 2013-05-30 17:43:59
893473348: 金币+5, 有帮助, 多谢 2013-06-17 09:58:36
如果在一维情况下,这个投影算子就会退化为限制信号上下界的算子,如果没有其他更好的办法,用这种办法至少可以保证不发散
日光之下皆覆辙 月光之下皆旧梦
2楼2013-05-30 16:27:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

893473348

银虫 (初入文坛)

引用回帖:
2楼: Originally posted by sunshine819 at 2013-05-30 16:27:00
如果在一维情况下,这个投影算子就会退化为限制信号上下界的算子,如果没有其他更好的办法,用这种办法至少可以保证不发散

你好,这是我在Matlab中编的算子程序,其中gradient=detaf(x,b)失球函数f的梯度的,求指点,非常感谢!

function z=Proj(x,y,a,b)    % 投影算子    % a,b分别表示所求参数的下、上界

% a=-10;    % 下界
% b=10;     % 上界
% x=-12;y=-2;
epsilon_x=0.01;    % 算子的扰动误差边界(以使f满足所设凸集的定义)
% if a>x           % 限定变量的范围([-10,10]只是theta的一组值,这个值对不同参数是不同的,据已知定)
%     x=a;
% elseif b<x
%     x=b;
% end

xmax=b;       % 此处参数还有待进一步确定,先用参数边界代替   %%%%%%(可能出现问题的地方,xmax的定义未知!!!)
f=(((epsilon_x+1)*(x^2))-xmax^2)/(epsilon_x*(xmax^2));   % 定义凸函数(函数值)

gradient=detaf(x,b);   % 调用函数detaf求上述函数的梯度向量(函数形式)
FX=eval(gradient);   % 计算梯度值
nabla_f=norm(FX);    % 求梯度值的范数

if f < 0
    z=y;
elseif f >= 0 && FX'*y <= 0
    z=y;
elseif f >= 0 && FX'*y > 0
    z=y-(FX/nabla_f)*(FX'*y/nabla_f)*f;
end
Nothing is impossible!
3楼2013-05-31 08:48:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunshine819

金虫 (小有名气)

【答案】应助回帖

引用回帖:
3楼: Originally posted by 893473348 at 2013-05-31 08:48:27
你好,这是我在Matlab中编的算子程序,其中gradient=detaf(x,b)失球函数f的梯度的,求指点,非常感谢!

function z=Proj(x,y,a,b)    % 投影算子    % a,b分别表示所求参数的下、上界

% a=-10;    % 下界
...

您的这个程序计算的 f(theta)不满足书里面给出的取值范围吧
日光之下皆覆辙 月光之下皆旧梦
4楼2013-07-16 12:52:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lwk223442

新虫 (初入文坛)

引用回帖:
4楼: Originally posted by sunshine819 at 2013-07-16 12:52:13
您的这个程序计算的 f(theta)不满足书里面给出的取值范围吧...

求一个proj的程序
谢谢
5楼2017-11-13 21:49:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

难忘今宵233

新虫 (初入文坛)

引用回帖:
5楼: Originally posted by lwk223442 at 2017-11-13 21:49:03
求一个proj的程序
谢谢...

q请问您现在有proj的程序吗?
6楼2018-01-22 15:48:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
本帖内容被屏蔽

7楼2021-07-20 09:52:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
本帖内容被屏蔽

8楼2021-07-20 09:53:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
本帖内容被屏蔽

9楼2021-07-20 09:54:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
本帖内容被屏蔽

10楼2021-07-20 09:54:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 893473348 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 h4CP7TrQR8Lg 2026-08-22 3/150 2026-08-23 16:30 by lXgU0RvNTDC9
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 DpSrDtM079iu 2026-08-22 3/150 2026-08-23 15:33 by KM5EcsNQRBPn
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 DpSrDtM079iu 2026-08-22 4/200 2026-08-23 15:30 by KM5EcsNQRBPn
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-21 6/300 2026-08-23 14:45 by KM5EcsNQRBPn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-21 5/250 2026-08-23 14:33 by KM5EcsNQRBPn
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 2JOx3r2CYEgw 2026-08-21 8/400 2026-08-23 14:30 by KM5EcsNQRBPn
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 h4CP7TrQR8Lg 2026-08-22 5/250 2026-08-23 12:09 by LR9qGULyN2ew
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[论文投稿] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 2JOx3r2CYEgw 2026-08-22 5/250 2026-08-23 03:52 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[基金申请] 人气不行了 +8 fansofjerry 2026-08-21 8/400 2026-08-22 16:30 by zyqchem
[基金申请] 今天基金会出结果吗?20260819 +16 kkkl_v 2026-08-19 17/850 2026-08-22 16:12 by 阿布Abu
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 放榜前的不淡定 40+4 snowwithsea 2026-08-19 14/700 2026-08-21 23:51 by cratir
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 我面上完蛋了 +7 且听虎啸 2026-08-20 8/400 2026-08-21 12:31 by 酷酷墨镜
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
信息提示
请填处理意见