24小时热门版块排行榜

>论坛更新日志 (1621)
>虫友互识 (237)
>基金申请 (58)
>论文投稿 (39)
>博后之家 (24)
>导师招生 (20)
>休闲灌水 (20)
>文献求助 (16)
>硕博家园 (14)
>教师之家 (9)
>考博 (9)
>考研 (8)
>药学 (7)
>公派出国 (7)
>论文道贺祈福 (3)
>找工作 (3)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

匿名

用户注销 (正式写手)

应助: 0 (幼儿园)
金币: 291.3
散金: 509
红花: 3
帖子: 430
在线: 161.1小时
虫号: 0
注册: 2011-08-09
性别: MM
专业: 心肌细胞/血管细胞损伤、修

本帖仅楼主可见

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

求助一个物理问题的 matlab 编程已经有0人回复

1楼 2016-01-06 11:35:26

已阅同方向广播申请仿真EPI 回复此楼编辑查看我的主页

文刀星

木虫 (小有名气)

应助: 13 (小学生)
金币: 3651.2
散金: 500
红花: 5
帖子: 60
在线: 100.9小时
虫号: 3693767
注册: 2015-02-22
性别: GG
专业: 导航、制导与传感技术

【答案】应助回帖

感谢参与，应助指数 +1

function dm=mymfun(z,m)

R=0.00008314;
T=298.15;
N=1.326963e6;
H=1.230375;
C=0.0156;
P=(1+0.00001*994*9.8*z);%表达式看清，自己填补

r=(3*m*R*T/(4*pi*P))^1.5;

if r<7e-4
Vb=4474*r^1.357;
else
Vb=0.23;
end

if r<6.67*-4
Kl=0.6*r;
else
Kl=4e-4;
end

dm=-Kl*(H*P-C)*4*pi*r^2*N/Vb;
%%%%%%%%%%%%%%%%%%
把上面的保存在一个m文件中然后再运行下面的代码
[z,m]=ode45(@mymfun,22:-0.5:0,0.22019232);

plot(z,m);
grid on
xlabel('z');
ylabel('m');
title('z-m');