| 查看: 846 | 回复: 7 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
zhaoqian59捐助贵宾 (小有名气)
|
[求助]
这个程序是不是哪里有问题？初始值跟实验值是一样的，但是模拟结果不对已有1人参与
|
||
|
f(x) pl=1000; phyd=110000; P0=29.39E9; a0=0.0341; v=1.5; x1=x(1); x2=x(2); pg=P0*(a0/x1)^(3*v); y=[x2 (1/pl*(pg-phyd)-3/2*x2^2)/x1]'; Ts = 0.001; % 仿真步长 Tn = 1; % 仿真终止时间 t = 0:Ts:Tn; % 仿真时间范围 N = length(t); x=[0.034 0]'; % 系统初值：x=[a a的导数] %% Variable Declaration and Initialization a=0.034*ones(1,N); a_dt=zeros(1,N); %x = [a a_d]'; % intial states %% Simulation Iteration for n = 1:N a(n)=x(1); a_dt(n)=x(2); k1 = fx(x); k2 = fx(x+Ts/2*k1); k3 = fx(x+Ts/2*k2); k4 = fx( x+Ts *k3); x = x + Ts/6 * ( k1 + 2*k2 + 2*k3 + k4 ); end %% Plot Response plot(t,a) xlabel('t/(s)') ylabel('a') 是一个二阶非线性常微分方程 |
回复此楼» 收录本帖的淘帖专辑推荐
 matlab编程绘图 |
» 猜你喜欢
三甲基碘化亚砜的氧化反应
已经有4人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有5人回复
-
孩子确诊有中度注意力缺陷
已经有12人回复
-
2025冷门绝学什么时候出结果
已经有3人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有4人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
-
AI论文写作工具：是科研加速器还是学术作弊器？
已经有3人回复
-
论文投稿，期刊推荐
已经有4人回复
-
硕士和导师闹得不愉快
已经有13人回复
zhaoqian59
捐助贵宾 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1514
- 散金: 8
- 红花: 5
- 帖子: 194
- 在线: 57.4小时
- 虫号: 2656296
- 注册: 2013-09-15
- 性别: MM
- 专业: 爆炸与冲击动力学
7楼2015-08-12 09:40:23
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68120.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2015-08-11 10:37:11
zhaoqian59
捐助贵宾 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1514
- 散金: 8
- 红花: 5
- 帖子: 194
- 在线: 57.4小时
- 虫号: 2656296
- 注册: 2013-09-15
- 性别: MM
- 专业: 爆炸与冲击动力学
3楼2015-08-11 15:03:26
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68120.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
4楼2015-08-11 15:08:12