24小时热门版块排行榜

返回列表

mathsxdg

金虫 (正式写手)

应助: 5 (幼儿园)
金币: 1985.4
散金: 620
红花: 3
帖子: 328
在线: 177.4小时
虫号: 1383339
注册: 2011-08-29
性别: GG
专业: 人类遗传学

[求助] 如何用MATLAB求解脉冲微分方程并画图

有以下脉冲微分方程，请高手赐教如何用Matlab求解并画图，多谢了！！！

其中，最后一个式子表示脉冲线l, b0 为一常数。

[ Last edited by mathsxdg on 2013-10-15 at 21:23 ]

回复此楼

» 猜你喜欢

想换工作。大多数高校都是评职称时认可5年内在原单位取得的成果吗？已经有9人回复
求个博导看看已经有16人回复
上海工程技术大学张培磊教授团队招收博士生已经有4人回复
上海工程技术大学【激光智能制造】课题组招收硕士已经有5人回复
求助院士们，这个如何合成呀已经有4人回复
临港实验室与上科大联培博士招生1名已经有9人回复
需要合成515-64-0，50g，能接单的留言已经有4人回复
自荐读博已经有4人回复
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？已经有6人回复
带资进组求博导收留已经有10人回复

» 本主题相关价值贴推荐，对您同样有帮助:

求助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人回复

生命科学

1楼 2013-10-15 21:19:29

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

niubenhit

金虫 (正式写手)

应助: 69 (初中生)
金币: 970.7
散金: 1813
红花: 9
帖子: 517
在线: 193.3小时
虫号: 2413993
注册: 2013-04-12
专业: 力学中的基本问题和方法

【答案】应助回帖

你这里 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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

niubenhit

金虫 (正式写手)

应助: 69 (初中生)
金币: 970.7
散金: 1813
红花: 9
帖子: 517
在线: 193.3小时
虫号: 2413993
注册: 2013-04-12
专业: 力学中的基本问题和方法

引用回帖:

2楼: Originally posted by niubenhit at 2013-10-26 08:21:08
你这里 l 是一点吧，可以这样求：
编一个循环
每个循环长度为 l ，在其中求解ode45的方程，然后将解得最后一点作为下一个循环的初值。
即可。

程序：
function draw()
l=10;
x10=1;x20=2;b0=0.2;
for j=1 ...

没看到最后一行，不好意思，楼主，上面程序有误，删了吧。

赞一下

回复此楼

有眼大如天，山高月更阔

3楼2013-10-26 08:22:26

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖