| 查看: 1070 | 回复: 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 但是，求解以后不对。希望得到高手帮助。 |
回复此楼» 猜你喜欢
欢迎发来filecode的Mz6后的代码验证其规律
已经有118人回复
-
今天放榜吗？
已经有15人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有3人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
filecode，4个jtjc了
已经有11人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有4人回复
-
今天基金会出结果吗？20260819
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- IOP的审稿期限给的多长？审稿人没有看到邀请信会多久自动切换状态？
已经有13人回复
- 两道有关固液体系 相变，要设计过程求G,H,S状态函数变化的题目
已经有23人回复
- CrystEngComm投稿后的状态变化请教
已经有7人回复
- 用MATLAB ode45求解2阶微分方程
已经有4人回复
- 关于 journal of power sources 投稿状态和时间的变化
已经有9人回复
- 【求助】求解释：荧光光谱法研究变性剂对蛋白质构象影响，峰强变而峰位不发生变化的可
已经有7人回复
- 2010年和2011年国家社科申报指南的一些变化（语言学）
已经有8人回复
- 【求助完毕】用ode逆向求解微分方程问题
已经有4人回复
- 解析under review状态中日期（status date）的变化
已经有116人回复
- 【交流】电离属于化学变化吗？
已经有7人回复
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
2楼2011-04-20 22:14:15
Hongtan
木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 5798.8
- 散金: 1044
- 红花: 3
- 帖子: 1468
- 在线: 843.3小时
- 虫号: 439487
- 注册: 2007-10-14
- 性别: GG
- 专业: 控制理论与方法
3楼2011-04-20 22:31:45
信彼南山
木虫 (著名写手)
- 应助: 33 (小学生)
- 金币: 4142.9
- 散金: 1221
- 红花: 16
- 帖子: 1178
- 在线: 233.5小时
- 虫号: 1133529
- 注册: 2010-10-27
- 专业: 导航、制导与传感技术
4楼2011-04-20 23:15:15
