24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 884  |  回复: 3

Hongtan

木虫 (著名写手)

[求助] 求助:控制器u随状态变化而切换时的ode45求解

考虑的系统很简单,即

\dot x_1=x_2
\dot x_2=u

u=-U sign(x_1-\beta x_1m)

其中x_1m 为 \dot x_1=0 时的值,也即 \dot x_1=0 时, x_1m 就变化一次, 且初始值为 x_1m=x_1(0)

我的程序是

主程序

clear all
global i;
global x10;
global M;

i=1;
tfinal=5;
step=0.001;
t0=[0:step:tfinal];
x10=2;
x20=1;
x0=[x10,x20];

[t,x]=ode45('slide',t0,x0);

figure(1)
subplot(2,2,1)
plot(t,x(:,1),'k-',t,x(:,2),'k:','Linewidth',1);
xlabel('t');
ylabel('x');



子程序

function y=slide(t,x)
global i
global x10;
global M;

U=6;
beta=0.5;

M(1)=x10;
if abs(x(2))<=0.1;
i=i+1;
M(i)=x(1);
end
x_1m=M(i);

u=-U*sign(x(1)-beta*x_1m);   

f1=x(2);
f2=3*sin(x(1))+(2+sin(t))*u;
y=[f1;f2];
t

但是,求解以后不对。希望得到高手帮助。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)

程序怎么这么乱啊?
M(1)=x10;
if abs(x(2))<=0.1;
i=i+1;
M(i)=x(1);
end
x_1m=M(i);

这一段想做什么?
2楼2011-04-20 22:14:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Hongtan

木虫 (著名写手)

引用回帖:
Originally posted by 信彼南山 at 2011-04-20 22:14:15:
程序怎么这么乱啊?
M(1)=x10;
if abs(x(2))<=0.1;
i=i+1;
M(i)=x(1);
end
x_1m=M(i);

这一段想做什么?

这一段是切换,呵呵。x_10是x_1的初始值,作为x_1m的初始值。因为无法实现正好在时刻 \dot x_1(t)=0 的时候切换,所以就取满足|\dotx_1|=|x_2|<0.1时,x_1m切换到下一个值。
3楼2011-04-20 22:31:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信彼南山

木虫 (著名写手)

你那段if一共就执行1次,切换什么了?
我咋没弄明白呢?
4楼2011-04-20 23:15:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Hongtan 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 071000生物学调剂 +7 拉提桃 2026-04-06 7/350 2026-04-06 18:55 by 52305043001
[考研] 305求调剂 +3 77Qi 2026-04-06 3/150 2026-04-06 15:36 by going home
[考研] 070300化学学硕311分求调剂 +11 梁富贵险中求 2026-04-04 11/550 2026-04-06 10:43 by 蓝云思雨
[考研] 一志愿南京航空航天大学 材料与化工329分求调剂 +8 Mr. Z 2026-04-05 8/400 2026-04-06 09:24 by dongzh2009
[考研] 315求调剂 +5 &123456789 2026-04-05 5/250 2026-04-05 19:55 by nepu_uu
[考研] 301求调剂 +3 XYPLR 2026-04-05 4/200 2026-04-05 19:07 by XYPLR
[考研] 085601,一志愿厦大334复试被刷求调剂 +13 曾仰之 2026-04-03 15/750 2026-04-04 20:13 by dongzh2009
[考研] 302求调剂一志愿华中师范大学 +8 小江小江江江 2026-04-02 8/400 2026-04-04 19:50 by 蓝云思雨
[考研] 怎么删帖子啊 +3 缝曦1000 2026-04-04 3/150 2026-04-04 14:20 by 土木硕士招生
[考研] 0710生物学336分求调剂 +6 kiyy 2026-04-01 8/400 2026-04-04 10:10 by kiyy
[考研] 268求调剂 +8 你好tg 2026-04-03 9/450 2026-04-04 05:08 by gswylq
[考研] 化工调剂303分,过四级 +28 栖梧待风 2026-04-02 28/1400 2026-04-03 21:40 by qlm5820
[考研] 295求调剂 +3 尚偌呀 2026-04-03 4/200 2026-04-03 21:23 by zhq0425
[考研] 335求调剂 +7 沈清璃 2026-04-03 7/350 2026-04-03 18:55 by lijunpoly
[考研] 071000生物学调剂 +8 知昭蔓 2026-04-02 8/400 2026-04-03 10:36 by macy2011
[考研] 22408 266求调剂 +3 masss11222 2026-04-02 3/150 2026-04-02 18:11 by 笔落锦州
[考研] 一志愿北交大材料工程,总分358 +4 cs0106 2026-04-01 4/200 2026-04-02 07:42 by 尚水阁主
[考研] 085600 一志愿9 总分351 求调剂学校 +7 czhcz 2026-03-31 9/450 2026-04-01 19:24 by 无际的草原
[考研] 复试调剂 +7 双马尾痞老板2 2026-03-31 7/350 2026-03-31 19:49 by Dyhoer
[考研] 370求调剂 +3 080700调剂 2026-03-30 3/150 2026-03-31 01:09 by A_Zhe
信息提示
请填处理意见