| 查看: 3519 | 回复: 9 | ||
[求助]
有关L1自适应控制投影算子Proj的Matlab编程问题
|
求助：有关L1自适应控制投影算子Proj的Matlab编程问题。现在需要编程实现下图所示的算子，编程错误，不能够保证输出量是有界的而导致控制系统发散，急求牛人指点迷津！非常感谢！！！ Proj.jpg |
回复此楼
sunshine819
金虫 (小有名气)
- 应助: 32 (小学生)
- 金币: 1170.9
- 散金: 7
- 红花: 2
- 帖子: 137
- 在线: 42.6小时
- 虫号: 2254567
- 注册: 2013-01-22
- 性别: GG
- 专业: 控制理论与方法
2楼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 |
3楼2013-05-31 08:48:27
sunshine819
金虫 (小有名气)
- 应助: 32 (小学生)
- 金币: 1170.9
- 散金: 7
- 红花: 2
- 帖子: 137
- 在线: 42.6小时
- 虫号: 2254567
- 注册: 2013-01-22
- 性别: GG
- 专业: 控制理论与方法
4楼2013-07-16 12:52:13
5楼2017-11-13 21:49:03
6楼2018-01-22 15:48:45
|
本帖内容被屏蔽 |
7楼2021-07-20 09:52:36
|
本帖内容被屏蔽 |
8楼2021-07-20 09:53:29
|
本帖内容被屏蔽 |
9楼2021-07-20 09:54:05
|
本帖内容被屏蔽 |
10楼2021-07-20 09:54:36