| 查看: 855 | 回复: 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编程绘图 |
» 猜你喜欢
寻求一种能扛住强氧化性腐蚀性的容器密封件
已经有7人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
-
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入
已经有5人回复
-
2025冷门绝学什么时候出结果
已经有7人回复
-
请问有评职称，把科研教学业绩算分排序的高校吗
已经有6人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有7人回复
-
请问下大家为什么这个铃木偶联几乎不反应呢
已经有5人回复
-
康复大学泰山学者周祺惠团队招收博士研究生
已经有6人回复
zhaoqian59
捐助贵宾 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1514
- 散金: 8
- 红花: 5
- 帖子: 194
- 在线: 57.4小时
- 虫号: 2656296
- 注册: 2013-09-15
- 性别: MM
- 专业: 爆炸与冲击动力学
3楼2015-08-11 15:03:26
zhaoqian59
捐助贵宾 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1514
- 散金: 8
- 红花: 5
- 帖子: 194
- 在线: 57.4小时
- 虫号: 2656296
- 注册: 2013-09-15
- 性别: MM
- 专业: 爆炸与冲击动力学
5楼2015-08-11 15:14:44
zhaoqian59
捐助贵宾 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 1514
- 散金: 8
- 红花: 5
- 帖子: 194
- 在线: 57.4小时
- 虫号: 2656296
- 注册: 2013-09-15
- 性别: MM
- 专业: 爆炸与冲击动力学
7楼2015-08-12 09:40:23