| 查看: 3620 | 回复: 5 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
mathsxdg金虫 (正式写手)
|
[求助]
如何用MATLAB求解脉冲微分方程并画图
|
|
有以下脉冲微分方程，请高手赐教如何用Matlab求解并画图，多谢了！！！ l:  其中，最后一个式子表示脉冲线l, b0 为一常数。 [ Last edited by mathsxdg on 2013-10-15 at 21:23 ] |
回复此楼» 猜你喜欢
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有15人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有14人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有16人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有15人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有15人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有12人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有12人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助matlab求解三元超越方程
已经有4人回复
- 中心差分求解二维热传导的matlab程序
已经有12人回复
- 求助matlab画图，急急急
已经有11人回复
- 用matlab如何解方程？
已经有6人回复
- MATLAB微分方程参数拟合问题，求大神
已经有7人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- 请教mathematic问题？关于解微分方程组的和画图的一些问题。
已经有19人回复
- 如何用matlab求解矩阵系数的二阶微分方程
已经有13人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- 求助：如何用MATLAB的偏微分方程工具箱画出不同区域内的温度分布图？
已经有5人回复
- matlab非线性方程组求解，急！！！
已经有5人回复
- matlab有限差分法求解复杂的偏微分方程组
已经有19人回复
- 求助——求解常系数二次微分方程
已经有6人回复
- matlab数值求解边界条件微分方程组
已经有7人回复
- 求matlab求解一个方程组的问题
已经有6人回复
- MATLAB非线性方程组求解~~~反应釜模拟
已经有9人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！
已经有18人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 求助一个差分方程的MATLAB实现
已经有3人回复
- 【求助】先解常微分方程 然后画图【已完结】
已经有11人回复
- 【求助】利用matlab求解微分方程得出解函数中含0的意思
已经有3人回复
mathsxdg
金虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 1985.4
- 散金: 620
- 红花: 3
- 帖子: 328
- 在线: 177.4小时
- 虫号: 1383339
- 注册: 2011-08-29
- 性别: GG
- 专业: 人类遗传学
5楼2013-10-26 08:53:18
【答案】应助回帖
|
你这里 l 是一点吧，可以这样求： 编一个循环 每个循环长度为 l ，在其中求解ode45的方程，然后将解得最后一点作为下一个循环的初值。 即可。 程序： function draw() l=10; x10=1;x20=2;b0=0.2; for j=1:10 [t,x]=ode45(@test,[(j-1)*l,j*l],[x10,x20]); [mm,nn]=size(x); plot(t,x);hold on; x10=x(nn,1)+(sqrt(3)/2-1)*x(nn,1)-1/2*b0*x(nn,2); y10=x(nn,2)+1/2*b0*x(nn,1)+(sqrt(3)/2-1)*x(nn,2); x=[]; t=[]; end plot(t,x); function dydt = test(t,x) %%以下为参数，可修改。 b0=1; %%以上为参数，可修改。 x1=x(1);x2=x(2); dydt =[2*x1-x2 x1+2*x2]; |
2楼2013-10-26 08:21:08
3楼2013-10-26 08:22:26
4楼2013-10-26 08:24:06
