应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » ode45
41/1返回列表
查看: 600  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

DIScreet123

新虫 (初入文坛)

[求助] ode45

用oed45如何设定非线性系统的随机输入，此随机输入要随时间变化
回复此楼

» 猜你喜欢
1楼 2015-11-15 14:19:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
% 定义X=[ …, …, …,…]，
% 写出微分方程
dX=inline(……);

% 初始条件，可以是变量
X0=[…, …, …,…];

% 微分方程求解
[t,X] = ode45(dX, [t0, t1], X0);
赞一下 回复此楼
2楼2015-11-15 16:24:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

DIScreet123

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by feixiaolin at 2015-11-15 16:24:33
% 定义X=，
% 写出微分方程
dX=inline(……);

% 初始条件，可以是变量
X0=;

% 微分方程求解
= ode45(dX, , X0);

谢谢您的答复 我想问下
function dx=suijishuru(t,x)
       dx(1)=x(2);
       dx(2)=normrnd(25,1,1,1)-25.*x(1)-x(1).^3-0.5.*x(2);
       dx=dx(;
end
x0=[0,0];
[t,x]=ode45('suijishuru',[0 20],x0);
plot(t,x(:,1))
为什么得到的时域图是衰减的？？ 按道理我给了他一个随机输入 他应该是受迫振动，到最终会稳定的啊？  请问我的程序有什么问题的啊？？万分感谢
赞一下 回复此楼
3楼2015-11-15 18:17:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

feixiaolin

荣誉版主 (文坛精英)

优秀版主
引用回帖:
3楼: Originally posted by DIScreet123 at 2015-11-15 18:17:49
谢谢您的答复 我想问下
function dx=suijishuru(t,x)
       dx(1)=x(2);
       dx(2)=normrnd(25,1,1,1)-25.*x(1)-x(1).^3-0.5.*x(2);
       dx=dx(;
end
x0=;
=ode45('suijishuru',,x0);
plot(t, ...

随机振动为噪声信号，用于系统辨识。

发自小木虫Android客户端
赞一下 回复此楼
4楼2015-11-15 19:03:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 DIScreet123 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭