| 查看: 930 | 回复: 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编程绘图 |
» 猜你喜欢
材料与化工专硕调剂
已经有6人回复
-
086500 325 求调剂
已经有3人回复
-
材料080500调剂求收留
已经有7人回复
-
一志愿南昌大学，327分，材料与化工085600
已经有6人回复
-
环境工程调剂
已经有9人回复
-
工科材料085601 279求调剂
已经有9人回复
-
287求调剂
已经有5人回复
-
22408 344分 求调剂 一志愿 华电计算机技术
已经有3人回复
-
289求调剂
已经有3人回复
-
317求调剂
已经有10人回复
zhaoqian59
捐助贵宾 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1514
- 散金: 8
- 红花: 5
- 帖子: 194
- 在线: 57.4小时
- 虫号: 2656296
- 注册: 2013-09-15
- 性别: MM
- 专业: 爆炸与冲击动力学
7楼2015-08-12 09:40:23
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68127.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
- 金币: 68127.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3291.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
4楼2015-08-11 15:08:12
50