| 查看: 937 | 回复: 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编程绘图 |
» 猜你喜欢
324分 085600材料化工求调剂
已经有4人回复
-
一志愿华中科技大学，080502，354分求调剂
已经有5人回复
-
一志愿武理材料305分求调剂
已经有7人回复
-
271材料工程求调剂
已经有8人回复
-
【考研调剂】化学专业 281分，一志愿四川大学，诚心求调剂
已经有8人回复
-
一志愿华南师大 070300（化学）304分求调剂
已经有3人回复
-
一志愿 西北大学 ，070300化学学硕，总分287，双非一本，求调剂。
已经有3人回复
-
一志愿重庆大学085700资源与环境专硕，总分308求调剂
已经有3人回复
-
354求调剂
已经有8人回复
-
304求调剂
已经有6人回复
zhaoqian59
捐助贵宾 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1514
- 散金: 8
- 红花: 5
- 帖子: 194
- 在线: 57.4小时
- 虫号: 2656296
- 注册: 2013-09-15
- 性别: MM
- 专业: 爆炸与冲击动力学
3楼2015-08-11 15:03:26
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
2楼2015-08-11 10:37:11
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
4楼2015-08-11 15:08:12
zhaoqian59
捐助贵宾 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1514
- 散金: 8
- 红花: 5
- 帖子: 194
- 在线: 57.4小时
- 虫号: 2656296
- 注册: 2013-09-15
- 性别: MM
- 专业: 爆炸与冲击动力学
5楼2015-08-11 15:14:44
10