应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 仿真模拟 » MATLAB » matlab求数值解，一个编程求助！
101/1返回列表
查看: 1104  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

匿名

用户注销 (正式写手)
本帖仅楼主可见

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2016-01-06 11:35:26
已阅   同方向广播   申请仿真EPI   回复此楼   编辑   查看我的主页

文刀星

木虫 (小有名气)
不就是一个微分方程的数值解吗？help一下ode45，很容易就编出来了！
赞一下 回复此楼
2楼2016-01-06 14:28:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)
本帖仅楼主可见
赞一下
3楼2016-01-06 14:35:08
已阅   申请仿真EPI   回复此楼   编辑   查看我的主页

lishouyin

铁杆木虫 (著名写手)

【答案】应助回帖

感谢参与，应助指数 +1
用c语言做for循环。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
4楼2016-01-06 14:47:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)
本帖仅楼主可见
赞一下
5楼2016-01-06 14:52:47
已阅   申请仿真EPI   回复此楼   编辑   查看我的主页

lishouyin

铁杆木虫 (著名写手)

【答案】应助回帖

应该是能，有相关语法。但matlab做循环效率低。

发自小木虫Android客户端
赞一下 回复此楼
6楼2016-01-06 14:59:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)
本帖仅楼主可见
7楼2016-01-06 15:15:54
已阅   申请仿真EPI   回复此楼   编辑   查看我的主页

lishouyin

铁杆木虫 (著名写手)

【答案】应助回帖

最好自己学下吧，很好学的。图片上的内容我也不太懂，很容易抄错。

发自小木虫Android客户端
赞一下 回复此楼
8楼2016-01-06 16:23:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)
本帖仅楼主可见
赞一下
9楼2016-01-06 16:28:32
已阅   申请仿真EPI   回复此楼   编辑   查看我的主页

文刀星

木虫 (小有名气)

【答案】应助回帖

感谢参与，应助指数 +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');
赞一下(1人) 回复此楼
10楼2016-01-06 16:53:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xingyunhu4 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭