| 查看: 1105 | 回复: 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编程绘图 |
» 猜你喜欢
filecode，4个jtjc了
已经有11人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有117人回复
-
filecode=后面第一个是大写字母
已经有12人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有4人回复
-
今天放榜吗？
已经有13人回复
-
今天基金会出结果吗？20260819
已经有8人回复
-
2027广东省杰青
已经有10人回复
-
朋友圈看到的
已经有8人回复
-
时间戳变了，能看出什么问题？
已经有14人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有30人回复
zhaoqian59
捐助贵宾 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1514
- 散金: 8
- 红花: 5
- 帖子: 194
- 在线: 57.4小时
- 虫号: 2656296
- 注册: 2013-09-15
- 性别: MM
- 专业: 爆炸与冲击动力学
7楼2015-08-12 09:40:23
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68166.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.3小时
- 虫号: 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
- 金币: 68166.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.3小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
4楼2015-08-11 15:08:12
